Changeset 12880


Ignore:
Timestamp:
2020-05-06T16:07:44+02:00 (7 months ago)
Author:
petesykes
Message:

resolved issue of no MLD data for forecast runs

File:
1 edited

Legend:

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

    r9538 r12880  
    578578                     hmld_zint_25h(:,:,jn) = hmld_zint_25h(:,:,jn) + hmld_zint(:,:) 
    579579                  ENDIF 
     580                  IF (lwp) THEN 
     581                     IF ( jn .EQ. 1 ) WRITE(numout,*) 'PS :: i_cnt_25h,kt :: ',i_cnt_25h,kt 
     582                  ENDIF 
    580583                  IF( i_cnt_25h .EQ. 25 .AND.  MOD( kt, i_steps*24) == 0 .AND. kt .NE. nn_it000 ) THEN 
     584                     IF (lwp) THEN 
     585                        WRITE(numout,*) 'zdf_mxl_zint (25h) : Outputting 25h data at i_cnt_25h=',i_cnt_25h 
     586                     ENDIF 
    581587                     CALL iom_put( "mldzint25h_"//cmld , hmld_zint_25h(:,:,jn) / 25._wp   ) 
     588                     ! Reset array 
     589                     hmld_zint_25h(:,:,jn) = hmld_zint(:,:) 
     590                     ! Reset 25h counter on last mld_diag 
     591                     IF ( jn .EQ. nn_mld_diag ) i_cnt_25h = 1 
    582592                  ENDIF 
    583593               ENDIF 
     
    590600               IF (lwp) THEN 
    591601                  WRITE(numout,*) 'zdf_mxl_zint (25h) : Summed the following number of hourly values so far',i_cnt_25h 
    592           ENDIF 
     602               ENDIF 
    593603               i_cnt_25h = i_cnt_25h + 1 
    594604               IF( mld_25h_init ) mld_25h_init = .FALSE. 
    595             ENDIF 
    596             IF( i_cnt_25h .EQ. 25 .AND.  MOD( kt, i_steps*24) == 0 .AND. kt .NE. nn_it000 ) THEN 
    597                i_cnt_25h = 1  
    598                DO jn = 1, nn_mld_diag  
    599                      hmld_zint_25h(:,:,jn) = hmld_zint(:,:)  
    600                ENDDO  
    601605            ENDIF 
    602606         ENDIF 
Note: See TracChangeset for help on using the changeset viewer.