Skip to content

Commit 1f74491

Browse files
andrew-plattandrew-platt
authored
Merge pull request #3007 from luwang00/f/GridInterp
Cubic interpolation for SeaState and HydroDyn
2 parents 3d19110 + 7ce0149 commit 1f74491

36 files changed

+1391
-766
lines changed

‎modules/hydrodyn/src/Morison.f90‎

Lines changed: 8 additions & 19 deletions
Original file line numberDiff line numberDiff line change
@@ -3482,7 +3482,7 @@ SUBROUTINE Morison_CalcOutput( Time, u, p, x, xd, z, OtherState, y, m, errStat,
34823482
REAL(ReKi) :: sinBeta, sinBeta1, sinBeta2
34833483
REAL(ReKi) :: cosBeta, cosBeta1, cosBeta2
34843484
REAL(ReKi) :: CMatrix(3,3), CMatrix1(3,3), CMatrix2(3,3), CTrans(3,3) ! Direction cosine matrix for element, and its transpose
3485-
REAL(ReKi) :: l, z1, z2, zMid, r1, r2, r1b, r2b, r1In, r2In, rMidIn, rn, rn1, rn2, z_hi, zFillGroup
3485+
REAL(ReKi) :: l, z1, z2, zMid, r1, r2, r1b, r2b, r1In, r2In, rMidIn, z_hi, zFillGroup
34863486
REAL(ReKi) :: Sa1, Sa2, Sa1b, Sa2b, SaMidb, Sa1In, Sa2In, SaMidIn
34873487
REAL(ReKi) :: Sb1, Sb2, Sb1b, Sb2b, SbMidb, Sb1In, Sb2In, SbMidIn
34883488
REAL(ReKi) :: dRdl_mg, dSadl_mg, dSbdl_mg ! shorthand for taper including marine growth of element i
@@ -3779,16 +3779,10 @@ SUBROUTINE Morison_CalcOutput( Time, u, p, x, xd, z, OtherState, y, m, errStat,
37793779
ELSE IF (mem%MHstLMod == 2) THEN ! Alternative hydrostatic load calculation
37803780
! Get free surface elevation and normal at the element midpoint (both assumed constant over the element)
37813781
posMid = 0.5 * (pos1+pos2)
3782-
! rn is only used to estimate free surface normal numerically
3783-
IF (mem%MSecGeom == MSecGeom_Cyl) THEN
3784-
rn = 0.5 * (r1b +r2b )
3785-
ELSE IF (mem%MSecGeom == MSecGeom_Rec) THEN
3786-
rn = MAX( 0.5*(Sa1b+Sa2b), 0.5*(Sb1b+Sb2b) )
3787-
END IF
37883782
IF (p%WaveField%WaveStMod > 0) THEN
37893783
CALL GetTotalWaveElev( Time, posMid, ZetaMid, ErrStat2, ErrMsg2 )
37903784
CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
3791-
CALL GetFreeSurfaceNormal( Time, posMid, rn, n_hat, ErrStat2, ErrMsg2 )
3785+
CALL GetFreeSurfaceNormal( Time, posMid, n_hat, ErrStat2, ErrMsg2 )
37923786
CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
37933787
FSPt = (/posMid(1),posMid(2),ZetaMid/) ! Reference point on the free surface
37943788
ELSE
@@ -4563,22 +4557,18 @@ SUBROUTINE Morison_CalcOutput( Time, u, p, x, xd, z, OtherState, y, m, errStat,
45634557
if (mem%MSecGeom==MSecGeom_Cyl) then
45644558
r1 = mem%RMGB( 1)
45654559
r2 = mem%RMGB(N+1)
4566-
rn1 = r1
4567-
rn2 = r2
45684560
else if (mem%MSecGeom==MSecGeom_Rec) then
45694561
Sa1 = mem%SaMGB( 1)
45704562
Sa2 = mem%SaMGB(N+1)
45714563
Sb1 = mem%SbMGB( 1)
45724564
Sb2 = mem%SbMGB(N+1)
4573-
rn1 = MAX(Sa1,Sb1)
4574-
rn2 = MAX(Sa2,Sb2)
45754565
end if
45764566
if (mem%i_floor == 0) then ! both ends above or at seabed
45774567
! Compute loads on the end plate of node 1
45784568
IF (p%WaveField%WaveStMod > 0) THEN
45794569
CALL GetTotalWaveElev( Time, pos1, Zeta1, ErrStat2, ErrMsg2 )
45804570
CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
4581-
CALL GetFreeSurfaceNormal( Time, pos1, rn1, n_hat, ErrStat2, ErrMsg2 )
4571+
CALL GetFreeSurfaceNormal( Time, pos1, n_hat, ErrStat2, ErrMsg2 )
45824572
CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
45834573
FSPt = (/pos1(1),pos1(2),Zeta1/) ! Reference point on the free surface
45844574
ELSE
@@ -4606,7 +4596,7 @@ SUBROUTINE Morison_CalcOutput( Time, u, p, x, xd, z, OtherState, y, m, errStat,
46064596
IF (p%WaveField%WaveStMod > 0) THEN
46074597
CALL GetTotalWaveElev( Time, pos2, Zeta2, ErrStat2, ErrMsg2 )
46084598
CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
4609-
CALL GetFreeSurfaceNormal( Time, pos2, rn2, n_hat, ErrStat2, ErrMsg2 )
4599+
CALL GetFreeSurfaceNormal( Time, pos2, n_hat, ErrStat2, ErrMsg2 )
46104600
CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
46114601
FSPt = (/pos2(1),pos2(2),Zeta2/) ! Reference point on the free surface
46124602
ELSE
@@ -4635,7 +4625,7 @@ SUBROUTINE Morison_CalcOutput( Time, u, p, x, xd, z, OtherState, y, m, errStat,
46354625
IF (p%WaveField%WaveStMod > 0) THEN
46364626
CALL GetTotalWaveElev( Time, pos2, Zeta2, ErrStat2, ErrMsg2 )
46374627
CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
4638-
CALL GetFreeSurfaceNormal( Time, pos2, rn2, n_hat, ErrStat2, ErrMsg2 )
4628+
CALL GetFreeSurfaceNormal( Time, pos2, n_hat, ErrStat2, ErrMsg2 )
46394629
CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
46404630
FSPt = (/pos2(1),pos2(2),Zeta2/) ! Reference point on the free surface
46414631
ELSE
@@ -4769,10 +4759,9 @@ SUBROUTINE GetTotalWaveElev( Time, pos, Zeta, ErrStat, ErrMsg )
47694759

47704760
END SUBROUTINE GetTotalWaveElev
47714761

4772-
SUBROUTINE GetFreeSurfaceNormal( Time, pos, r, n, ErrStat, ErrMsg)
4762+
SUBROUTINE GetFreeSurfaceNormal( Time, pos, n, ErrStat, ErrMsg)
47734763
REAL(DbKi), INTENT( In ) :: Time
47744764
REAL(ReKi), INTENT( In ) :: pos(*) ! Position at which free-surface normal is to be calculated. Third entry ignored if present.
4775-
REAL(ReKi), INTENT( In ) :: r ! Distance for central differencing
47764765
REAL(ReKi), INTENT( OUT ) :: n(3) ! Free-surface normal vector
47774766
INTEGER(IntKi), INTENT( OUT ) :: ErrStat ! Error status of the operation
47784767
CHARACTER(*), INTENT( OUT ) :: ErrMsg ! Error message if errStat /= ErrID_None
@@ -4782,7 +4771,7 @@ SUBROUTINE GetFreeSurfaceNormal( Time, pos, r, n, ErrStat, ErrMsg)
47824771
ErrStat = ErrID_None
47834772
ErrMsg = ""
47844773

4785-
CALL WaveField_GetNodeWaveNormal( p%WaveField, m%WaveField_m, Time, pos, r, n, ErrStat2, ErrMsg2 )
4774+
CALL WaveField_GetNodeWaveNormal( p%WaveField, m%WaveField_m, Time, pos, n, ErrStat2, ErrMsg2 )
47864775
CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
47874776

47884777
END SUBROUTINE GetFreeSurfaceNormal
@@ -5716,7 +5705,7 @@ SUBROUTINE getElementHstLds_Mod1( Time, pos1, pos2, Zeta1, Zeta2, k_hat, r1, r2,
57165705
rh = r1 + h0*dRdl
57175706
! Estimate the free-surface normal at the free-surface intersection, n_hat
57185707
IF ( p%WaveField%WaveStMod > 0_IntKi ) THEN ! If wave stretching is enabled, compute free surface normal
5719-
CALL GetFreeSurfaceNormal( Time, FSInt, rh, n_hat, ErrStat2, ErrMsg2 )
5708+
CALL GetFreeSurfaceNormal( Time, FSInt, n_hat, ErrStat2, ErrMsg2 )
57205709
CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
57215710
ELSE ! Without wave stretching, use the normal of the SWL
57225711
n_hat = (/0.0_ReKi,0.0_ReKi,1.0_ReKi/)

‎modules/hydrodyn/src/Morison.txt‎

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -462,7 +462,7 @@ typedef ^ ^ ReKi
462462
typedef ^ ^ ReKi V_rel_n_HiPass {:} - - "High-pass filtered normal relative flow velocity at joints" m/s
463463
typedef ^ ^ ReKi zFillGroup {:} - - "Instantaneous highest point of each filled group" m
464464
typedef ^ ^ MeshMapType VisMeshMap - - - "Mesh mapping for visualization mesh" -
465-
typedef ^ ^ SeaSt_WaveField_MiscVarType WaveField_m - - - "misc var information from the SeaState Interpolation module" -
465+
typedef ^ ^ GridInterp_MiscVarType WaveField_m - - - "misc var information from the Grid Interpolation module" -
466466

467467
# ..... Parameters ................................................................................................................
468468
# Define parameters here:

‎modules/hydrodyn/src/Morison_Types.f90‎

Lines changed: 5 additions & 5 deletions
Original file line numberDiff line numberDiff line change
@@ -530,7 +530,7 @@ MODULE Morison_Types
530530
REAL(ReKi) , DIMENSION(:), ALLOCATABLE :: V_rel_n_HiPass !< High-pass filtered normal relative flow velocity at joints [m/s]
531531
REAL(ReKi) , DIMENSION(:), ALLOCATABLE :: zFillGroup !< Instantaneous highest point of each filled group [m]
532532
TYPE(MeshMapType) :: VisMeshMap !< Mesh mapping for visualization mesh [-]
533-
TYPE(SeaSt_WaveField_MiscVarType) :: WaveField_m !< misc var information from the SeaState Interpolation module [-]
533+
TYPE(GridInterp_MiscVarType) :: WaveField_m !< misc var information from the Grid Interpolation module [-]
534534
END TYPE Morison_MiscVarType
535535
! =======================
536536
! ========= Morison_ParameterType =======
@@ -4693,7 +4693,7 @@ subroutine Morison_CopyMisc(SrcMiscData, DstMiscData, CtrlCode, ErrStat, ErrMsg)
46934693
call NWTC_Library_CopyMeshMapType(SrcMiscData%VisMeshMap, DstMiscData%VisMeshMap, CtrlCode, ErrStat2, ErrMsg2)
46944694
call SetErrStat(ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName)
46954695
if (ErrStat >= AbortErrLev) return
4696-
call SeaSt_WaveField_CopyMisc(SrcMiscData%WaveField_m, DstMiscData%WaveField_m, CtrlCode, ErrStat2, ErrMsg2)
4696+
call GridInterp_CopyMisc(SrcMiscData%WaveField_m, DstMiscData%WaveField_m, CtrlCode, ErrStat2, ErrMsg2)
46974697
call SetErrStat(ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName)
46984698
if (ErrStat >= AbortErrLev) return
46994699
end subroutine
@@ -4780,7 +4780,7 @@ subroutine Morison_DestroyMisc(MiscData, ErrStat, ErrMsg)
47804780
end if
47814781
call NWTC_Library_DestroyMeshMapType(MiscData%VisMeshMap, ErrStat2, ErrMsg2)
47824782
call SetErrStat(ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName)
4783-
call SeaSt_WaveField_DestroyMisc(MiscData%WaveField_m, ErrStat2, ErrMsg2)
4783+
call GridInterp_DestroyMisc(MiscData%WaveField_m, ErrStat2, ErrMsg2)
47844784
call SetErrStat(ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName)
47854785
end subroutine
47864786

@@ -4821,7 +4821,7 @@ subroutine Morison_PackMisc(RF, Indata)
48214821
call RegPackAlloc(RF, InData%V_rel_n_HiPass)
48224822
call RegPackAlloc(RF, InData%zFillGroup)
48234823
call NWTC_Library_PackMeshMapType(RF, InData%VisMeshMap)
4824-
call SeaSt_WaveField_PackMisc(RF, InData%WaveField_m)
4824+
call GridInterp_PackMisc(RF, InData%WaveField_m)
48254825
if (RegCheckErr(RF, RoutineName)) return
48264826
end subroutine
48274827

@@ -4868,7 +4868,7 @@ subroutine Morison_UnPackMisc(RF, OutData)
48684868
call RegUnpackAlloc(RF, OutData%V_rel_n_HiPass); if (RegCheckErr(RF, RoutineName)) return
48694869
call RegUnpackAlloc(RF, OutData%zFillGroup); if (RegCheckErr(RF, RoutineName)) return
48704870
call NWTC_Library_UnpackMeshMapType(RF, OutData%VisMeshMap) ! VisMeshMap
4871-
call SeaSt_WaveField_UnpackMisc(RF, OutData%WaveField_m) ! WaveField_m
4871+
call GridInterp_UnpackMisc(RF, OutData%WaveField_m) ! WaveField_m
48724872
end subroutine
48734873

48744874
subroutine Morison_CopyParam(SrcParamData, DstParamData, CtrlCode, ErrStat, ErrMsg)

‎modules/hydrodyn/src/SS_Excitation.txt‎

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -43,7 +43,7 @@ typedef ^ ^ SS_Exc_ContinuousStateType
4343
# Define any data that are used only for efficiency purposes (these variables are not associated with time):
4444
# e.g. indices for searching in an array, large arrays that are local variables in any routine called multiple times, etc.
4545
typedef ^ MiscVarType INTEGER LastIndWave - 1 - "last used index in the WaveTime array" -
46-
typedef ^ ^ SeaSt_WaveField_MiscVarType WaveField_m - - - "misc var information from the SeaState Interpolation module" -
46+
typedef ^ ^ GridInterp_MiscVarType WaveField_m - - - "misc var information from the Grid Interpolation module" -
4747

4848

4949
# ..... Parameters .........................

‎modules/hydrodyn/src/SS_Excitation_Types.f90‎

Lines changed: 5 additions & 5 deletions
Original file line numberDiff line numberDiff line change
@@ -73,7 +73,7 @@ MODULE SS_Excitation_Types
7373
! ========= SS_Exc_MiscVarType =======
7474
TYPE, PUBLIC :: SS_Exc_MiscVarType
7575
INTEGER(IntKi) :: LastIndWave = 1 !< last used index in the WaveTime array [-]
76-
TYPE(SeaSt_WaveField_MiscVarType) :: WaveField_m !< misc var information from the SeaState Interpolation module [-]
76+
TYPE(GridInterp_MiscVarType) :: WaveField_m !< misc var information from the Grid Interpolation module [-]
7777
END TYPE SS_Exc_MiscVarType
7878
! =======================
7979
! ========= SS_Exc_ParameterType =======
@@ -494,7 +494,7 @@ subroutine SS_Exc_CopyMisc(SrcMiscData, DstMiscData, CtrlCode, ErrStat, ErrMsg)
494494
ErrStat = ErrID_None
495495
ErrMsg = ''
496496
DstMiscData%LastIndWave = SrcMiscData%LastIndWave
497-
call SeaSt_WaveField_CopyMisc(SrcMiscData%WaveField_m, DstMiscData%WaveField_m, CtrlCode, ErrStat2, ErrMsg2)
497+
call GridInterp_CopyMisc(SrcMiscData%WaveField_m, DstMiscData%WaveField_m, CtrlCode, ErrStat2, ErrMsg2)
498498
call SetErrStat(ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName)
499499
if (ErrStat >= AbortErrLev) return
500500
end subroutine
@@ -508,7 +508,7 @@ subroutine SS_Exc_DestroyMisc(MiscData, ErrStat, ErrMsg)
508508
character(*), parameter :: RoutineName = 'SS_Exc_DestroyMisc'
509509
ErrStat = ErrID_None
510510
ErrMsg = ''
511-
call SeaSt_WaveField_DestroyMisc(MiscData%WaveField_m, ErrStat2, ErrMsg2)
511+
call GridInterp_DestroyMisc(MiscData%WaveField_m, ErrStat2, ErrMsg2)
512512
call SetErrStat(ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName)
513513
end subroutine
514514

@@ -518,7 +518,7 @@ subroutine SS_Exc_PackMisc(RF, Indata)
518518
character(*), parameter :: RoutineName = 'SS_Exc_PackMisc'
519519
if (RF%ErrStat >= AbortErrLev) return
520520
call RegPack(RF, InData%LastIndWave)
521-
call SeaSt_WaveField_PackMisc(RF, InData%WaveField_m)
521+
call GridInterp_PackMisc(RF, InData%WaveField_m)
522522
if (RegCheckErr(RF, RoutineName)) return
523523
end subroutine
524524

@@ -528,7 +528,7 @@ subroutine SS_Exc_UnPackMisc(RF, OutData)
528528
character(*), parameter :: RoutineName = 'SS_Exc_UnPackMisc'
529529
if (RF%ErrStat /= ErrID_None) return
530530
call RegUnpack(RF, OutData%LastIndWave); if (RegCheckErr(RF, RoutineName)) return
531-
call SeaSt_WaveField_UnpackMisc(RF, OutData%WaveField_m) ! WaveField_m
531+
call GridInterp_UnpackMisc(RF, OutData%WaveField_m) ! WaveField_m
532532
end subroutine
533533

534534
subroutine SS_Exc_CopyParam(SrcParamData, DstParamData, CtrlCode, ErrStat, ErrMsg)

‎modules/hydrodyn/src/WAMIT.f90‎

Lines changed: 1 addition & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -241,7 +241,7 @@ SUBROUTINE WAMIT_Init( InitInp, u, p, x, xd, z, OtherState, y, m, Interval, ErrS
241241

242242
! Set up wave excitation grid - Can no longer use the WaveField parameters due to different headings
243243
! Copy WaveField grid parameters
244-
call SeaSt_WaveField_CopyParam(p%WaveField%GridParams, p%ExctnGridParams, 0, ErrStat2, ErrMsg2); CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName)
244+
call GridInterp_CopyParam(p%WaveField%VolGridParams, p%ExctnGridParams, 0, ErrStat2, ErrMsg2); CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName)
245245
if ( p%ExctnDisp == 0 ) then
246246
p%ExctnGridParams%n(2:3) = 1_IntKi
247247
p%ExctnGridParams%delta(2:3) = 0.0_SiKi
@@ -257,7 +257,6 @@ SUBROUTINE WAMIT_Init( InitInp, u, p, x, xd, z, OtherState, y, m, Interval, ErrS
257257
p%ExctnGridParams%pZero(4) = -Pi
258258
end if
259259
p%ExctnGridParams%n(4) = p%NExctnHdg+1
260-
p%ExctnGridParams%Z_depth = -1.0 ! Set to Z_depth to a negative value to indicate uniform "z" grid for platform heading
261260

262261
! This module's implementation requires that if NBodyMod = 2 or 3, then there is one instance of a WAMIT module for each body, therefore, HydroDyn may have NBody > 1, but this WAMIT module will have NBody = 1
263262
if ( (p%NBodyMod > 1) .and. (p%NBody > 1) ) then

‎modules/hydrodyn/src/WAMIT.txt‎

Lines changed: 3 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -17,6 +17,7 @@ usefrom Conv_Radiation.txt
1717
usefrom SS_Radiation.txt
1818
usefrom SS_Excitation.txt
1919
usefrom SeaSt_WaveField.txt
20+
usefrom GridInterp.txt
2021

2122
typedef WAMIT/WAMIT InitInputType INTEGER NBody - - - "[>=1; only used when PotMod=1. If NBodyMod=1, the WAMIT data contains a vector of size 6*NBody x 1 and matrices of size 6*NBody x 6*NBody; if NBodyMod>1, there are NBody sets of WAMIT data each with a vector of size 6 x 1 and matrices of size 6 x 6]" -
2223
typedef ^ ^ INTEGER NBodyMod - - - "Body coupling model {1: include coupling terms between each body and NBody in HydroDyn equals NBODY in WAMIT, 2: neglect coupling terms between each body and NBODY=1 with XBODY=0 in WAMIT, 3: Neglect coupling terms between each body and NBODY=1 with XBODY=/0 in WAMIT} (switch) [only used when PotMod=1]" -
@@ -95,7 +96,7 @@ typedef ^ ^ SS_Exc_Outp
9596
typedef ^ ^ Conv_Rdtn_MiscVarType Conv_Rdtn - - - "" -
9697
typedef ^ ^ Conv_Rdtn_InputType Conv_Rdtn_u - - - "" -
9798
typedef ^ ^ Conv_Rdtn_OutputType Conv_Rdtn_y - - - "" -
98-
typedef ^ ^ SeaSt_WaveField_MiscVarType WaveField_m - - - "misc var information from the SeaState Interpolation module" -
99+
typedef ^ ^ GridInterp_MiscVarType WaveField_m - - - "misc var information from the Grid Interpolation module" -
99100

100101
# ..... Parameters ................................................................................................................
101102
# Define parameters here:
@@ -120,7 +121,7 @@ typedef ^ ^ SS_Exc_Para
120121
typedef ^ ^ DbKi DT - - - "" -
121122
typedef ^ ^ SeaSt_WaveFieldType *WaveField - - - "Pointer to wave field"
122123
typedef ^ ^ INTEGER PtfmYMod - - - "Large yaw model" -
123-
typedef ^ ^ SeaSt_WaveField_ParameterType ExctnGridParams - - - "Parameters of WaveExctnGrid" -
124+
typedef ^ ^ GridInterp_ParameterType ExctnGridParams - - - "Parameters of WaveExctnGrid" -
124125
#
125126
#
126127
# ..... Inputs ....................................................................................................................

‎modules/hydrodyn/src/WAMIT2.f90‎

Lines changed: 1 addition & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -2646,7 +2646,7 @@ SUBROUTINE CheckInitInput( InitInp, p, MnDriftData, NewmanAppData, DiffQTFData,
26462646

26472647
! Set up 2nd-order wave excitation grid
26482648
! Copy WaveField grid parameters
2649-
call SeaSt_WaveField_CopyParam(InitInp%WaveField%GridParams, p%Exctn2GridParams, 0, ErrStatTmp, ErrMsgTmp); CALL SetErrStat( ErrStatTmp, ErrMsgTmp, ErrStat, ErrMsg, RoutineName)
2649+
call GridInterp_CopyParam(InitInp%WaveField%VolGridParams, p%Exctn2GridParams, 0, ErrStatTmp, ErrMsgTmp); CALL SetErrStat( ErrStatTmp, ErrMsgTmp, ErrStat, ErrMsg, RoutineName)
26502650
! x and y grids are currently not used for second-order wave excitation
26512651
p%Exctn2GridParams%n(2:3) = 1_IntKi
26522652
p%Exctn2GridParams%delta(2:3) = 0.0_SiKi
@@ -2662,7 +2662,6 @@ SUBROUTINE CheckInitInput( InitInp, p, MnDriftData, NewmanAppData, DiffQTFData,
26622662
p%Exctn2GridParams%pZero(4) = -Pi
26632663
end if
26642664
p%Exctn2GridParams%n(4) = p%NExctnHdg+1
2665-
p%Exctn2GridParams%Z_depth = -1.0 ! Set to Z_depth to a negative value to indicate uniform "z" grid for platform heading
26662665

26672666

26682667
!> 1. Check that we only specified one of MnDrift, NewmanApp, or DiffQTF

0 commit comments

Comments
 (0)