Changeset 14714 for branches/UKMO/AMM15_v3_6_STABLE_package_reanalysis4
- Timestamp:
- 2021-04-14T18:55:43+02:00 (3 years 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 1-hourly and 24-hourly 4 !! NB 24-hour 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.