Ocean Modeling Discussion

ROMS/TOMS

Search for:
It is currently Thu Oct 19, 2017 8:00 pm




Post new topic Reply to topic  [ 1 post ] 

All times are UTC

Author Message
PostPosted: Wed Oct 04, 2017 4:36 pm 
Offline

Joined: Tue Jun 06, 2017 1:54 pm
Posts: 2
Location: Dalhousie University
Hello,

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:
!
!  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

In case nRST does not equal nDIA and the last written restart time stamp is before the last one in the diagnostics file, the old data in the diagnostics file are not overwritten at the correct position but the newly calculated data are appended to the file. I suppose changing the logic in DEF_DIAGS analogous to that of DEF_HIS should solve this issue.

Best,
Fabian


Top
 Profile  
Reply with quote  
Display posts from previous:  Sort by  
Post new topic Reply to topic  [ 1 post ] 

All times are UTC


Who is online

Users browsing this forum: No registered users and 2 guests


You cannot post new topics in this forum
You cannot reply to topics in this forum
You cannot edit your posts in this forum
You cannot delete your posts in this forum
You cannot post attachments in this forum

Search for:
Jump to:  
Powered by phpBB® Forum Software © phpBB Group