Changeset 14714 for branches


Ignore:
Timestamp:
2021-04-14T18:55:43+02:00 (4 months ago)
Author:
rrenshaw
Message:

section transports, hardwire calculation of 1-hour and 24-hour means

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/AMM15_v3_6_STABLE_package_reanalysis4/NEMOGCM/NEMO/OPA_SRC/DIA/diadct.F90

    r11639 r14714  
    11MODULE diadct 
     2  !! 
     3  !! JT/RJR: this version writes accumulated means 1-hourly and 24-hourly 
     4  !!         NB 24-hour values are means, not sums, of the 24 hourly values 
     5  !! 
    26  !!===================================================================== 
    37  !!                       ***  MODULE  diadct  *** 
     
    177181 
    178182        !Do calculation for daily, 25hourly mean every hour 
    179         nn_dct=3600./rdt         ! hard coded for NOOS transects, to give 25 hour means from hourly instantaneous values 
     183        !JT nn_dct=3600./rdt         ! hard coded for NOOS transects, to give 25 hour means from hourly instantaneous values 
    180184 
    181185        !write out daily, 25hourly mean every day 
     
    394398  
    395399        ! Compute transport and write only at nn_dctwri 
    396         IF ( MOD(kt,nn_dct)==0 .or. &               ! compute transport every nn_dct time steps 
    397             (ln_NOOS .and. kt==nn_it000 ) )  THEN   ! also include first time step when calculating NOOS 25 hour averages 
     400        !JT IF ( MOD(kt,nn_dct)==0 .or. &               ! compute transport every nn_dct time steps 
     401        !JT     (ln_NOOS .and. kt==nn_it000 ) )  THEN   ! also include first time step when calculating NOOS 25 hour averages 
     402        IF ( MOD(kt,nn_dct_h)==0 ) THEN            ! compute transport every nn_dct_h time steps 
    398403             
    399404             
     
    419424              ! Not 24 values, but 25! divide by ((nn_dctwri/nn_dct) +1) 
    420425              !! divide arrays by nn_dctwri/nn_dct  to obtain average 
    421               transports_3d(:,:,:,:)= transports_3d(:,:,:,:)/((nn_dctwri/nn_dct)+1.) 
    422               transports_2d(:,:,:)  = transports_2d(:,:,:)  /((nn_dctwri/nn_dct)+1.) 
     426              !JT transports_3d(:,:,:,:)= transports_3d(:,:,:,:)/((nn_dctwri/nn_dct)+1.) 
     427              !JT transports_2d(:,:,:)  = transports_2d(:,:,:)  /((nn_dctwri/nn_dct)+1.) 
     428              transports_3d(:,:,:,:)= transports_3d(:,:,:,:)/((nn_dctwri/nn_dct)) 
     429              transports_2d(:,:,:)  = transports_2d(:,:,:)  /((nn_dctwri/nn_dct)) 
    423430 
    424431              ! Sum over each class 
     
    453460                   
    454461                   
    455                   IF ( ln_NOOS ) CALL transport(secs(jsec),lldebug,jsec)  ! reinitialise for next 25 hour instantaneous average (overlapping values) 
     462                  !JT IF ( ln_NOOS ) CALL transport(secs(jsec),lldebug,jsec)  ! reinitialise for next 25 hour instantaneous average (overlapping values) 
    456463 
    457464 
Note: See TracChangeset for help on using the changeset viewer.