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.