Changeset 14714
 Timestamp:
 20210414T18:55:43+02:00 (18 months ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

branches/UKMO/AMM15_v3_6_STABLE_package_reanalysis4/NEMOGCM/NEMO/OPA_SRC/DIA/diadct.F90
r11639 r14714 1 1 MODULE diadct 2 !! 3 !! JT/RJR: this version writes accumulated means 1hourly and 24hourly 4 !! NB 24hour values are means, not sums, of the 24 hourly values 5 !! 2 6 !!===================================================================== 3 7 !! *** MODULE diadct *** … … 177 181 178 182 !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 values183 !JT nn_dct=3600./rdt ! hard coded for NOOS transects, to give 25 hour means from hourly instantaneous values 180 184 181 185 !write out daily, 25hourly mean every day … … 394 398 395 399 ! 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 398 403 399 404 … … 419 424 ! Not 24 values, but 25! divide by ((nn_dctwri/nn_dct) +1) 420 425 !! 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)) 423 430 424 431 ! Sum over each class … … 453 460 454 461 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) 456 463 457 464
Note: See TracChangeset
for help on using the changeset viewer.