Changeset 783 for trunk/NEMO/OPA_SRC/TRD

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/TRD/trdmld_rst.F90 (modified) (4 diffs)

trunk/NEMO/OPA_SRC/TRD/trdmld_rst.F90

@@ -46 +46 @@
       !!--------------------------------------------------------------------------------
 
-      IF( kt == nitrst-1 ) THEN
-         IF( nitrst > 1.0e9 ) THEN   
-            WRITE(clkt,*) nitrst
-         ELSE
-            WRITE(clkt,'(i8.8)') nitrst
+      ! to get better performances with NetCDF format:
+      ! we open and define the ocean restart_mld file one time step before writing the data (-> at nitrst - 1)
+      ! except if we write ocean restart_mld files every time step or if an ocean restart_mld file was writen at nitend - 1
+      IF( kt == nitrst - 1 .OR. nstock == 1 .OR. ( kt == nitend .AND. MOD( nitend - 1, nstock ) == 0 ) ) THEN
+         ! beware of the format used to write kt (default is i8.8, that should be large enough...)
+         IF( nitrst > 999999999 ) THEN   ;   WRITE(clkt, *       ) nitrst
+         ELSE                            ;   WRITE(clkt, '(i8.8)') nitrst
          ENDIF
+         ! create the file
          clname = TRIM(cexper)//"_"//TRIM(ADJUSTL(clkt))//"_restart_mld"
-         IF(lwp) WRITE(numout,*) '             open ocean restart_mld NetCDF file: '//clname
+         IF(lwp) THEN
+            WRITE(numout,*)
+            SELECT CASE ( jprstlib )
+            CASE ( jprstdimg )   ;   WRITE(numout,*) '             open ocean restart_mld binary file: '//clname
+            CASE DEFAULT         ;   WRITE(numout,*) '             open ocean restart_mld NetCDF file: '//clname
+            END SELECT
+            IF( kt == nitrst - 1 ) THEN   ;   WRITE(numout,*) '             kt = nitrst - 1 = ', kt,' date= ', ndastp
+            ELSE                          ;   WRITE(numout,*) '             kt = '             , kt,' date= ', ndastp
+            ENDIF
+         ENDIF
+
          CALL iom_open( clname, nummldw, ldwrt = .TRUE., kiolib = jprstlib )
       ENDIF
@@ -59 +72 @@
       IF( kt == nitrst .AND. lwp ) THEN
          WRITE(numout,*)
-         WRITE(numout,*) 'trdmld_rst: output for ML diags. restart, with trd_mld_rst_write routine'
+         WRITE(numout,*) 'trdmld_rst: output for ML diags. restart, with trd_mld_rst_write routine kt =', kt
          WRITE(numout,*) '~~~~~~~~~~'
          WRITE(numout,*)
@@ -81 +94 @@
          CALL iom_rstput( kt, nitrst, nummldw, 'tml_sumb'        , tml_sumb        )
          DO jk = 1, jpltrd
-            IF( jk < 10 )   THEN
-               WRITE(charout,FMT="('tmltrd_csum_ub_', I1)") jk
-            ELSE
-               WRITE(charout,FMT="('tmltrd_csum_ub_', I2)") jk
+            IF( jk < 10 ) THEN   ;   WRITE(charout,FMT="('tmltrd_csum_ub_', I1)") jk
+            ELSE                 ;   WRITE(charout,FMT="('tmltrd_csum_ub_', I2)") jk
             ENDIF
             CALL iom_rstput( kt, nitrst, nummldw, charout,  tmltrd_csum_ub(:,:,jk) )
@@ -94 +105 @@
          CALL iom_rstput( kt, nitrst, nummldw, 'sml_sumb'        , sml_sumb        )
          DO jk = 1, jpltrd
-            IF( jk < 10 )   THEN
-               WRITE(charout,FMT="('smltrd_csum_ub_', I1)") jk
-            ELSE
-               WRITE(charout,FMT="('smltrd_csum_ub_', I2)") jk
+            IF( jk < 10 ) THEN   ;   WRITE(charout,FMT="('smltrd_csum_ub_', I1)") jk
+            ELSE                 ;   WRITE(charout,FMT="('smltrd_csum_ub_', I2)") jk
             ENDIF
             CALL iom_rstput( kt, nitrst, nummldw, charout , smltrd_csum_ub(:,:,jk) )