There seems to be a bug in the DEF_DIAGS at the point where the output time step for the diagnostics file is determined after a restart of ROMS. In my code (which is from August this year) lines 721-734 read as:
Code: Select all
!
! Set unlimited time record dimension to the appropriate value.
!
IF (nRST(ng).eq.nDIA(ng)) THEN
IF (ndefDIA(ng).gt.0) THEN
DIA(ng)%Rindex=((ntstart(ng)-1)- &
& ndefDIA(ng)*((ntstart(ng)-1)/ndefDIA(ng)))/ &
& nDIA(ng)
ELSE
DIA(ng)%Rindex=(ntstart(ng)-1)/nDIA(ng)
END IF
ELSE
DIA(ng)%Rindex=rec_size
END IF
Best,
Fabian