Changeset 10000


Ignore:
Timestamp:
2018-07-25T13:43:27+02:00 (2 years ago)
Author:
kingr
Message:

Pete's correction to 25h diagnostic.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/AMM15_v3_6_STABLE_package_withNOOS/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfmxl.F90

    r9474 r10000  
    581581                  IF( i_cnt_25h .EQ. 25 .AND.  MOD( kt, i_steps*24) == 0 .AND. kt .NE. nn_it000 ) THEN 
    582582                     CALL iom_put( "mldzint25h_"//cmld , hmld_zint_25h(:,:,jn) / 25._wp   ) 
     583                     ! PS - now we need to reset the 25h counter and re-initialise the array 
     584                     IF( jn .EQ. nn_mld_diag ) THEN 
     585                        ! only reset counter once all the MLDs have been written 
     586                        i_cnt_25h = 1 
     587                     ENDIF 
     588                     hmld_zint_25h(:,:,jn) = hmld_zint(:,:) 
    583589                  ENDIF 
    584590               ENDIF 
     
    586592            ENDIF 
    587593         END DO 
    588                    
     594          
    589595         IF(  mld_25h_write  ) THEN 
    590596            IF( ( MOD( kt, i_steps ) == 0 ) .OR.  mld_25h_init ) THEN 
     
    595601               IF( mld_25h_init ) mld_25h_init = .FALSE. 
    596602            ENDIF 
    597             IF( i_cnt_25h .EQ. 25 .AND.  MOD( kt, i_steps*24) == 0 .AND. kt .NE. nn_it000 ) THEN 
    598                i_cnt_25h = 1  
    599                DO jn = 1, nn_mld_diag 
    600                      hmld_zint_25h(:,:,jn) = hmld_zint(:,:) 
    601                ENDDO 
    602             ENDIF 
    603603         ENDIF 
    604604                   
Note: See TracChangeset for help on using the changeset viewer.