Modes Change Log
================
v1.00, Unknown Date Oregon State University
This was the version obtained from Alan Wright.
v2.00, 30-Mar-2000, M. Buhl
I moved the eignevalue and eigenvector routines into Eigen.for.
I converted the remaining routines into free-form Fortran. I moved
the data that was stored in common blocks to modules in Modules.f90.
I cleaned up the code, eliminated unnecessary code and added some
performance improvements.
I added the program name, version, and the current date and time to
the output file.
I simplified the logic in Interp.
I redid the output file.
I made it so it is easy to convert the program to double precision. I
did not test it with double conversion.
v2.01, 11-Apr-2000, M. Buhl
Unknown bug fix.
v2.10, 21-Apr-2000, M. Buhl
I modified the code so that, for blades, we read in both the out-of-
plane and in-plane stiffnesses in addition to a pitch angle. We then
use the pitch angle to pitch calculate the in-plane and out-of-plane
stiffnesses from the flatwise and edgewise stiffnesses. We then
calculate two sets of mode shapes. In order to accomodate the above, I
enhanced the input file. I also reordered some of the input parameters.
I also changed the blade/tower flag to be a logical variable. I also
adjusted the output.
I dynamically allocated the BM and BS arrays.
v2.20, 24-Apr-2000, M. Buhl
I added twist distribution to the input file for blades. Using it, we
now transform the stiffness from the local principal axes (flatwise and
edgewise) to the out-of-plane/in-plane axis system before doing the
eigen analysis.
v2.21, 05-Apr-2001, M. Buhl
I changed the equation for TKK in Modes.f90\Modes() so the the gravity
term for the tower case is subtracted from this stiffness matrix instead
of being added. This lowers the tower frequencies slightly.
v2.22, 29-Apr-2002
At Craig Hansen's request, I modified the output statements so they use
F12.4 format instead of F12.3. This will reduce the chances of FAST
rejecting a mode shape because it's terms do not add up to 1.0 within a
plus or minus 0.001 error. Hopefully, we won't overflow the field with
coefficients greater than a million.