- required adding the SeaState misc vars to some HD types
- Added a warning if the writeOutputs request a position outside of the defined domain
- fix comments regarding where time is in the 4-D arrays (some comments had it in position 4, but it is in position 1)

- there was an extra call to SD_CalcOutput after the transfer and then a copy of some of the transfer codes, so I put the call to SD_CalcOutput inside `Transfer_Structure_to_Opt1Inputs` 
- also fixed some comments in SeaState

Use pointers in `Linear_HD_InputSolve_du`, 'Linear_HD_InputSolve_dy', and 'Linear_MAP_InputSolve_dy' to avoid duplicating code for a differene of 1 mesh

Also moved index routines for SD%y to make their placement consistent with the other modules

- Changed name of optional parameter to `ModName_GetOP` routine from `NeedLogMap` to `NeedPackedOrient` in an attempt to clarify what that does (returns the orientation matrix as a 3-valued vector instead of 9-valued one)
- Made the use of this value consistent in ElastoDyn (all meshes with orientation fields that are packed in the GetOP routine need to have the same value for their  call to `PackMotionMesh`)
- Added the optional `NeedPackedOrient` input to `SD_GetOP`
- Made ED, BD, and SD each pack their arrays the same way in `FAST_DiffInterpOutputs`; otherwise, the linear trim solution does not difference all of the outputs correctly to determine if it has converged

The field is getting set in `WAMIT2_CalcOutput`, but it is not allocated in the `MeshCreate` statement. Seems like that would have caused seg faults...

This could probably be improved by just plotting the SeaState grid at each time, but this works with the way it was done before.

- moved check for "has ice' from HD to SeaSt (SeaState must have WaveMod=0 and CurrMod=0) in that case
- don't check for WvHiCOff > WvLoCOff if WaveMod = 0,1 or 6
- warn if WvHiCOff is being modified internally
- fix module abbreviation for SeaState in FAST_Subs.f90
- removed some unused variables

SeaState/HydroDyn cleanup

…enFAST#1172)

- possibly segmentation faults due to arrays not being allocated before being used in SS_Exc_InputExtrapInterp
- also updated error handling

previous channel names did not have a leading 0 on the member #; this change allows the leading zero to be omitted

- added a few more outputs for MCF accelerations (used in validation vs other codes)
- removed unused variables
- combined some writeOutput variables and routines
- use SetOutParam routine generated by Matlab (instead of splitting this code into 3 other routines)
- added checks for invalid channel outputs

Do not print names of modules that are not used. Also, adds SeaState module to the list when/if it is used.

Reduces complexity of IF statements that are used each time step

- add `-module` argument to `manualRegressionTest.py` so it will run the module tests, too
- HD regression test now also copies the data it uses from `glue-codes\openfast\5MW_Baselin\HydroData`
- modules now return the error code if they exit on error after running an executable
- added `const` flag to optional logical arguments in `manualRegressionTest.py` and `executeHydrodynRegressionCase.py`

WaveMod = 7 was missing from WAMIT_Init, so it won't work with ExctnMod > 0.

… relative velocity filter with misc. variables.

Now HD and SeaSt data are output in the same file

Also note the .log files now have a different base name (haven't checked this for other drivers, just HD)

…hydro

Moved parts of driver into subroutines; also involved moving data into data types

- move linearization routines from main code to module
- move setting of HD inputs into a separate subroutine and attempt to simplify (also fix a bug on HD orientation inputs when prpInputsMod< 0
- destroy allocated variables
- add error checks on PRPin data (assumptions on time and length of arrays could have caused incorrect results or seg faults)
- linearization when PrpInputsMod == 2 or < 0 now uses correct HD inputs (otherwise inputs weren't set)
- updated some error checking

e.g., driver has only 1, and OpenFAST may have arrays of 2 or 3 elements.

return NaN when the number of output channels in the two files are not the same (there is no comparison for which channels are which)

The original equation for the white-noise spectrum is correct.

These outputs happen after moving on to the next step, so I wanted to clarify that in the code (and not confuse myself when I look at this again)

Was missing a transfer to the two meshes.
Also set time = 0 before calling SimStatus the first time

- renamed some Waves/Waves2 variables to avoid confusion with SeaSt variables of same name
- removed some unused variables (especially unused copies of data in different data structures)
- removed commented-out wave stretching code (old, from FAST7?)
- moved all of the SeaState summary file output into a subroutine
- fixed some allocation size errors
- fixed uninitialized variables: NWaveElev and ErrStat

The SeaSt module had the same name for the WaveElev and WaveKin arrays as its submodules, though at the SeaSt level they indicate values at user-defined points and at the submodule level, they indicate values on a regularly spaced grid. I renamed the submodule variables for less confusion (at least on my part)

Since the first nx*ny points of the WaveKinGrid*i arrays are identical to the WaveElevGrid*i arrays, I removed the WaveElevGrid*i arrays.

Code cleanup + HD driver linearization

- Bokeh must be at least v2.4
- executeSeaStateRegressionCase.py has same changes as other module scripts that were changed in dev
- manual test uses same keywords/arguments as execute*RegressionCase.py scripts

…t channels

this still doesn't check that the outputs are in the same order, so we could be comparing different channels against each other....