Changeset 10000

Timestamp:

2018-07-25T13:43:27+02:00 (6 years ago)

Author:

kingr

Message:

Pete's correction to 25h diagnostic.

File:

• branches/UKMO/AMM15_v3_6_STABLE_package_withNOOS/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfmxl.F90 (modified) (3 diffs)

branches/UKMO/AMM15_v3_6_STABLE_package_withNOOS/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfmxl.F90

@@ -581 +581 @@
                   IF( i_cnt_25h .EQ. 25 .AND.  MOD( kt, i_steps*24) == 0 .AND. kt .NE. nn_it000 ) THEN
                      CALL iom_put( "mldzint25h_"//cmld , hmld_zint_25h(:,:,jn) / 25._wp   )
+                     ! PS - now we need to reset the 25h counter and re-initialise the array
+                     IF( jn .EQ. nn_mld_diag ) THEN
+                        ! only reset counter once all the MLDs have been written
+                        i_cnt_25h = 1
+                     ENDIF
+                     hmld_zint_25h(:,:,jn) = hmld_zint(:,:)
                   ENDIF
                ENDIF
@@ -586 +592 @@
             ENDIF
          END DO
-                  
+         
          IF(  mld_25h_write  ) THEN
             IF( ( MOD( kt, i_steps ) == 0 ) .OR.  mld_25h_init ) THEN
@@ -595 +601 @@
                IF( mld_25h_init ) mld_25h_init = .FALSE.
             ENDIF
-            IF( i_cnt_25h .EQ. 25 .AND.  MOD( kt, i_steps*24) == 0 .AND. kt .NE. nn_it000 ) THEN
-               i_cnt_25h = 1 
-               DO jn = 1, nn_mld_diag
-                     hmld_zint_25h(:,:,jn) = hmld_zint(:,:)
-               ENDDO
-            ENDIF
          ENDIF