Changeset 3594 for branches/2012/dev_r3452_UKMO9_RESTART/NEMOGCM/NEMO/OPA_SRC/DOM/domain.F90

Timestamp:

2012-11-19T13:28:55+01:00 (11 years ago)

Author:

rfurner

Message:

code not tested through SETTEE, builds and runs, but has not been thoroughly tested, so will not be included in 2012 merge, however submitted back to keep record of work done for 2013 developments

File:

• branches/2012/dev_r3452_UKMO9_RESTART/NEMOGCM/NEMO/OPA_SRC/DOM/domain.F90 (modified) (3 diffs)

branches/2012/dev_r3452_UKMO9_RESTART/NEMOGCM/NEMO/OPA_SRC/DOM/domain.F90

@@ -123 +123 @@
       !!----------------------------------------------------------------------
       USE ioipsl
+
+      INTEGER  ::   js               ! dummy loop indice
+      CHARACTER(LEN=10)   ::   cltimes  ! restart dump times needed
+
       NAMELIST/namrun/ nn_no   , cn_exp    , cn_ocerst_in, cn_ocerst_out, ln_rstart , nn_rstctl,   &
          &             nn_it000, nn_itend  , nn_date0    , nn_leapy     , nn_istate , nn_stock ,   &
@@ -152 +156 @@
          WRITE(numout,*) '      leap year calendar (0/1)        nn_leapy   = ', nn_leapy
          WRITE(numout,*) '      initial state output            nn_istate  = ', nn_istate
-         WRITE(numout,*) '      frequency of restart file       nn_stock   = ', nn_stock
+         IF ( nn_stock(1) < 0 ) THEN
+            WRITE(numout,*) '      frequency of restart file       nn_stock   = ', ABS(nn_stock(1))
+         ELSE
+            WRITE(numout,*) '      list of restart times           nn_stock   = ', nn_stock
+         ENDIF
          WRITE(numout,*) '      frequency of output file        nn_write   = ', nn_write
          WRITE(numout,*) '      multi file dimgout              ln_dimgnnn = ', ln_dimgnnn
@@ -168 +176 @@
       nleapy = nn_leapy
       ninist = nn_istate
-      nstock = nn_stock
       nwrite = nn_write
 
 
       !                             ! control of output frequency
-      IF ( nstock == 0 .OR. nstock > nitend ) THEN
-         WRITE(ctmp1,*) 'nstock = ', nstock, ' it is forced to ', nitend
+      IF ( ALL(nn_stock == 0) ) THEN
+         WRITE(ctmp1,*) 'No restart output times specified, restart forced to output at end of run, kt = ', nitend
          CALL ctl_warn( ctmp1 )
-         nstock = nitend
-      ENDIF
+         nn_stock(1) = nitend
+      ELSE
+         IF ( nn_stock(1) > 0 ) THEN
+            DO js=1,size(nn_stock)-1
+               IF ( nn_stock(js+1) <= nn_stock(js) .AND. nn_stock(js+1) .NE. 0.0 ) THEN
+                  WRITE(ctmp1,*) 'Restart times in nn_stock not monotonically increasing, some values ignored'
+                  CALL ctl_warn( ctmp1 )
+               ENDIF
+            ENDDO
+            IF ( ANY(nn_stock > nitend) ) THEN
+               WRITE(ctmp1,*) 'some values of nn_stock exceed run length they are forced to end of run, kt = ', nitend
+               CALL ctl_warn( ctmp1 )
+               WHERE (nn_stock > nitend ) nn_stock=nitend
+            ENDIF
+         ELSEIF ( nn_stock(1) < 0 ) THEN
+            IF ( ABS(nn_stock(1)) > nitend ) THEN
+               WRITE(ctmp1,*) 'Attempting to output restarts at frequency greater than run legth, restart forced to output at end of run, kt = ', nitend
+               CALL ctl_warn( ctmp1 )
+               nn_stock(1) = nitend
+            ELSE
+               DO js=1,SIZE(nn_stock) 
+                  nn_stock(js) = MIN( (ABS(nn_stock(1)) * js  +  nit000 - 1), nitend)
+               ENDDO
+               IF ( MAXVAL(nn_stock) < nitend ) THEN
+                  WRITE(cltimes,FMT='(i10)') (nitend - nit000 + 1)/nn_stock(1) + 1   
+                  CALL ctl_stop( 'rst_opn:', &   
+                  'Too many restart dump times to store in the array', &   
+                  'Increase jpstocks to ' // cltimes  ) 
+               ENDIF
+            ENDIF
+         ENDIF
+      ENDIF
+        
       IF ( nwrite == 0 ) THEN
          WRITE(ctmp1,*) 'nwrite = ', nwrite, ' it is forced to ', nitend