Mercurial > hgweb > burst-coupling

changeset 80:78e31cdb54fb 0.9.0

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Continue with timestepping during sediment spin-up (BUT it will be reset with next recoupling); for output filenames, don't rewind for diagdec, diagcen and diagmil.
author Marco van Hulten <marco@hulten.org>
date Fri, 31 May 2019 11:52:13 +0200
parents 3fb63810b3be
children 8383e92c1312
files mo_sedmnt_offline.F90 mod_dia.F
diffstat 2 files changed, 7 insertions(+), 2 deletions(-) [+]
line wrap: on
line diff
--- a/mo_sedmnt_offline.F90	Sun May 12 12:13:33 2019 +0200
+++ b/mo_sedmnt_offline.F90	Fri May 31 11:52:13 2019 +0200
@@ -347,8 +347,11 @@
 
 ! get new date
 nday=1
-if (nmonth == 1) then
+nmonth = nmonth + 1
+if (nmonth > 12) then
    nday_of_year=1
+   nmonth = 1
+   nyear = nyear + 1
    if (calendar(1:3) == 'sta') then
       if (mod(nyear,4)   == 0 .and.                                  &
        & (mod(nyear,100) /= 0 .or. mod(nyear,400) == 0)) then
--- a/mod_dia.F	Sun May 12 12:13:33 2019 +0200
+++ b/mod_dia.F	Fri May 31 11:52:13 2019 +0200
@@ -277,7 +277,9 @@
                write(fname,'(4a,i6.6,a)')
      .           runid(1:runid_len),'.micom.',ctag,'.',ny,'.nc'
             else
-               call pstdat(ny,nm,nd,1)
+               if (diagann) then
+                  call pstdat(ny,nm,nd,1)
+               endif ! don't rewind for diagdec, diagcen and diagmil
                write(fname,'(4a,i4.4,a)')
      .           runid(1:runid_len),'.micom.',ctag,'.',ny,'.nc'
             endif