Skip to content

MoorDyn - SeaState runtime improvements #3308

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
andrew-platt merged 7 commits into from
Apr 21, 2026
+181 −111
Merged

MoorDyn - SeaState runtime improvements #3308

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
7 commits
Select commit Hold shift + click to select a range
115fddd
Change MoorDyn warning string for output files to remove too specific…
RyanDavies19 Feb 22, 2026
d254da2
MD: Add ramp time for water kinematics
RyanDavies19 Nov 18, 2025
d60da3e
MD: Add line node acceleration output
RyanDavies19 Apr 8, 2026
d070a19
MD: runtime improvements for SeaState coupling
RyanDavies19 Apr 14, 2026
2f1252a
MD: remove forgotten print statements
RyanDavies19 Apr 14, 2026
3f50051
Merge remote-tracking branch 'OpenFAST/rc-5.0.1' into MD_improvements
andrew-platt Apr 21, 2026
e4e6e52
Update r-test: RM1 floating tank case
andrew-platt Apr 21, 2026
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
43 changes: 27 additions & 16 deletions modules/moordyn/src/MoorDyn.f90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -487,6 +487,8 @@ SUBROUTINE MD_Init(InitInp, u, p, x, xd, z, other, y, m, DTcoupling, InitOut, Er
read (OptValue,*) p%inertialF
else if ( OptString == 'INERTIALF_RAMPT') then
read (OptValue,*) p%inertialF_rampT
else if ( OptString == 'WAVEKIN_RAMPT') then
read (OptValue,*) p%waveKin_rampT
else if ( OptString == 'OUTSWITCH') then
read (OptValue,*) p%OutSwitch
else if ( OptString == 'DISABLEOUTTIME') then
Expand All @@ -501,6 +503,14 @@ SUBROUTINE MD_Init(InitInp, u, p, x, xd, z, other, y, m, DTcoupling, InitOut, Er
Line = NextLine(i)

END DO

! extract dtCoupling from input and warn if dtOut is invalid
p%dtCoupling = DTcoupling ! store coupling time step for use in updatestates
IF (p%dtOut > 0.0_DbKi .and. p%dtOut < DTcoupling) THEN
ErrStat2 = ErrID_Info
ErrMsg2 = 'MoorDyn dtOut is less than the coupling time step. Output will be written at the coupling time step.'
CALL CheckError( ErrStat2, ErrMsg2 )
END IF

if (p%writeLog > 1) then
write(p%UnLog, '(A)' ) " - Options List:"
Expand Down Expand Up @@ -1635,25 +1645,29 @@ SUBROUTINE MD_Init(InitInp, u, p, x, xd, z, other, y, m, DTcoupling, InitOut, Er
! process output flag characters (LineOutString) and set line output flag array (OutFlagList)
m%LineList(l)%OutFlagList = 0 ! first set array all to zero
! per node 3 component
IF ( scan( LineOutString, 'p') > 0 ) m%LineList(l)%OutFlagList(2) = 1 ! node position (p)
IF ( scan( LineOutString, 'v') > 0 ) m%LineList(l)%OutFlagList(3) = 1 ! node velocity (v)
IF ( scan( LineOutString, 'U') > 0 ) m%LineList(l)%OutFlagList(4) = 1 ! node displacement (U)
IF ( scan( LineOutString, 'D') > 0 ) m%LineList(l)%OutFlagList(5) = 1 ! node rotation (D)
IF ( scan( LineOutString, 'b') > 0 ) m%LineList(l)%OutFlagList(6) = 1 ! seabed contact forces (B)
IF ( scan( LineOutString, 'V') > 0 ) m%LineList(l)%OutFlagList(7) = 1 ! VIV forces
IF ( scan( LineOutString, 'p') > 0 ) m%LineList(l)%OutFlagList(2) = 1 ! node position (p)
IF ( scan( LineOutString, 'v') > 0 ) m%LineList(l)%OutFlagList(3) = 1 ! node velocity (v)
IF ( scan( LineOutString, 'a') > 0 ) m%LineList(l)%OutFlagList(4) = 1 ! node acceleration (a)
IF ( scan( LineOutString, 'U') > 0 ) m%LineList(l)%OutFlagList(5) = 1 ! wave velocity (U)
IF ( scan( LineOutString, 'D') > 0 ) m%LineList(l)%OutFlagList(6) = 1 ! hydrodynamic force (D)
IF ( scan( LineOutString, 'b') > 0 ) m%LineList(l)%OutFlagList(7) = 1 ! seabed contact forces (b)
IF ( scan( LineOutString, 'V') > 0 ) m%LineList(l)%OutFlagList(8) = 1 ! VIV forces
! per node 1 component
IF ( scan( LineOutString, 'W') > 0 ) m%LineList(l)%OutFlagList(8) = 1 ! node weight/buoyancy (positive up)
IF ( scan( LineOutString, 'K') > 0 ) m%LineList(l)%OutFlagList(9) = 1 ! curvature at node
IF ( scan( LineOutString, 'W') > 0 ) m%LineList(l)%OutFlagList(9) = 1 ! node weight/buoyancy (positive up)
IF ( scan( LineOutString, 'K') > 0 ) m%LineList(l)%OutFlagList(10) = 1 ! curvature at node
! per element 1 component
IF ( scan( LineOutString, 't') > 0 ) m%LineList(l)%OutFlagList(10) = 1 ! segment tension force (just EA)
IF ( scan( LineOutString, 'c') > 0 ) m%LineList(l)%OutFlagList(11) = 1 ! segment internal damping force
IF ( scan( LineOutString, 's') > 0 ) m%LineList(l)%OutFlagList(12) = 1 ! Segment strain
IF ( scan( LineOutString, 'd') > 0 ) m%LineList(l)%OutFlagList(13) = 1 ! Segment strain rate
IF ( scan( LineOutString, 'l') > 0 ) m%LineList(l)%OutFlagList(14) = 1 ! Segment stretched length
IF ( scan( LineOutString, 't') > 0 ) m%LineList(l)%OutFlagList(11) = 1 ! segment tension force (just EA)
IF ( scan( LineOutString, 'c') > 0 ) m%LineList(l)%OutFlagList(12) = 1 ! segment internal damping force
IF ( scan( LineOutString, 's') > 0 ) m%LineList(l)%OutFlagList(13) = 1 ! Segment strain
IF ( scan( LineOutString, 'd') > 0 ) m%LineList(l)%OutFlagList(14) = 1 ! Segment strain rate
IF ( scan( LineOutString, 'l') > 0 ) m%LineList(l)%OutFlagList(15) = 1 ! Segment stretched length

IF (SUM(m%LineList(l)%OutFlagList) > 0) m%LineList(l)%OutFlagList(1) = 1 ! this first entry signals whether to create any output file at all
! the above letter-index combinations define which OutFlagList entry corresponds to which output type

! set flag to store node accelerations if needed for VIV or acceleration output
m%LineList(l)%store_rdd = (m%LineTypeList(m%LineList(l)%PropsIdNum)%Cl > 0) .OR. (m%LineList(l)%OutFlagList(4) == 1)


if (p%writeLog > 1) then
write(p%UnLog, '(A)' ) " - Line"//trim(num2lstr(m%LineList(l)%IdNum))//":"
Expand Down Expand Up @@ -2930,9 +2944,6 @@ SUBROUTINE MD_Init(InitInp, u, p, x, xd, z, other, y, m, DTcoupling, InitOut, Er
END DO

end if ! InputFileDat%TMaxIC > 0


p%dtCoupling = DTcoupling ! store coupling time step for use in updatestates

other%dummy = 0
xd%dummy = 0
Expand Down
Loading
Loading