Changeset 783 for trunk/NEMO/OPA_SRC/DOM/domain.F90

Timestamp:

2008-01-07T15:29:41+01:00 (16 years ago)

Author:

smasson

Message:

write multiple restarts, see ticket:44

File:

• trunk/NEMO/OPA_SRC/DOM/domain.F90 (modified) (2 diffs)

trunk/NEMO/OPA_SRC/DOM/domain.F90

@@ -181 +181 @@
       ndastp = ndate0                ! Assign initial date to current date
 
-! ... Control of output frequency
-      IF ( nstock == 0 ) THEN
-          IF(lwp)WRITE(numout,cform_war)
-          IF(lwp)WRITE(numout,*) '           nstock = ', nstock, ' it is forced to ', nitend
-          nstock = nitend
-          nwarn = nwarn + 1
+      ! ... Control of output frequency
+      IF ( nstock == 0 .OR. nstock > nitend - nit000 + 1 ) THEN
+         WRITE(ctmp1,*) '           nstock = ', nstock, ' it is forced to ', nitend - nit000 + 1
+         CALL ctl_warn( ctmp1 )
+         nstock = nitend - nit000 + 1
       ENDIF
       IF ( nwrite == 0 ) THEN
-          IF(lwp)WRITE(numout,cform_war)
-          IF(lwp)WRITE(numout,*) '           nwrite = ', nwrite, ' it is forced to ', nitend
-          nwrite = nitend
-          nwarn = nwarn + 1
+         WRITE(ctmp1,*) '           nwrite = ', nwrite, ' it is forced to ', nitend
+         CALL ctl_warn( ctmp1 )
+         nwrite = nitend
       ENDIF
 
@@ -261 +259 @@
       IF( lk_ice_lim ) THEN
          IF(lwp) WRITE(numout,*) '           ice model coupling frequency      nfice  = ', nfice
+         IF( MOD( nitend - nit000 + 1, nfice) /= 0 ) THEN 
+            WRITE(ctmp1,*) 'experiment length (', nitend - nit000 + 1, ') is NOT a multiple of nfice (', nfice, ')'
+            CALL ctl_stop( ctmp1, 'Impossible to do proper restart files' )
+         ENDIF
+         IF( MOD( nstock             , nfice) /= 0 ) THEN 
+            WRITE(ctmp1,*) 'nstock ('           , nstock             , ') is NOT a multiple of nfice (', nfice, ')'
+            CALL ctl_stop( ctmp1, 'Impossible to do proper restart files' )
+         ENDIF
          nfbulk = nfice
-         IF( MOD( rday, nfice*rdt ) /= 0 ) THEN 
-            IF(lwp) WRITE(numout,*) ' '
-            IF(lwp) WRITE(numout,*) 'W A R N I N G :  nfice is NOT a multiple of the number of time steps in a day'
-            IF(lwp) WRITE(numout,*) ' '
-         ENDIF
+         IF( MOD( rday, nfice*rdt ) /= 0 )   CALL ctl_warn( 'nfice is NOT a multiple of the number of time steps in a day' )
          IF(lwp) WRITE(numout,*) '           bulk computation frequency       nfbulk  = ', nfbulk, ' = nfice if ice model used'
          IF(lwp) WRITE(numout,*) '           flag closed sea or not           nclosea = ', nclosea