Changeset 12880

Timestamp:

2020-05-06T16:07:44+02:00 (3 years ago)

Author:

petesykes

Message:

resolved issue of no MLD data for forecast runs

File:

• branches/UKMO/AMM15_v3_6_STABLE_package_collate_PS44/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfmxl.F90 (modified) (2 diffs)

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

@@ -578 +578 @@
                      hmld_zint_25h(:,:,jn) = hmld_zint_25h(:,:,jn) + hmld_zint(:,:)
                   ENDIF
+                  IF (lwp) THEN
+                     IF ( jn .EQ. 1 ) WRITE(numout,*) 'PS :: i_cnt_25h,kt :: ',i_cnt_25h,kt
+                  ENDIF
                   IF( i_cnt_25h .EQ. 25 .AND.  MOD( kt, i_steps*24) == 0 .AND. kt .NE. nn_it000 ) THEN
+                     IF (lwp) THEN
+                        WRITE(numout,*) 'zdf_mxl_zint (25h) : Outputting 25h data at i_cnt_25h=',i_cnt_25h
+                     ENDIF
                      CALL iom_put( "mldzint25h_"//cmld , hmld_zint_25h(:,:,jn) / 25._wp   )
+                     ! Reset array
+                     hmld_zint_25h(:,:,jn) = hmld_zint(:,:)
+                     ! Reset 25h counter on last mld_diag
+                     IF ( jn .EQ. nn_mld_diag ) i_cnt_25h = 1
                   ENDIF
                ENDIF
@@ -590 +600 @@
                IF (lwp) THEN
                   WRITE(numout,*) 'zdf_mxl_zint (25h) : Summed the following number of hourly values so far',i_cnt_25h
-          ENDIF
+               ENDIF
                i_cnt_25h = i_cnt_25h + 1
                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