Changeset 12495 for NEMO/branches/2020/dev_r12472_ASINTER-05_Masson_CurrentFeedback/src/TOP/TRP/trdmxl_trc.F90
- Timestamp:
- 2020-03-02T09:10:34+01:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2020/dev_r12472_ASINTER-05_Masson_CurrentFeedback/src/TOP/TRP/trdmxl_trc.F90
r12377 r12495 401 401 DO jn = 1, jptra 402 402 IF( ln_trdtrc(jn) ) THEN 403 !-- Compute total trends (use rdttrc instead of rdt ???)403 !-- Compute total trends 404 404 IF ( ln_trcadv_muscl .OR. ln_trcadv_muscl2 ) THEN ! EULER-FORWARD schemes 405 ztmltot(:,:,jn) = ( tml_trc(:,:,jn) - tmlbn_trc(:,:,jn) )/r dt405 ztmltot(:,:,jn) = ( tml_trc(:,:,jn) - tmlbn_trc(:,:,jn) )/rn_Dt 406 406 ELSE ! LEAP-FROG schemes 407 ztmltot(:,:,jn) = ( tml_trc(:,:,jn) - tmlbn_trc(:,:,jn) + tmlb_trc(:,:,jn) - tmlbb_trc(:,:,jn))/(2.*r dt)407 ztmltot(:,:,jn) = ( tml_trc(:,:,jn) - tmlbn_trc(:,:,jn) + tmlb_trc(:,:,jn) - tmlbb_trc(:,:,jn))/(2.*rn_Dt) 408 408 ENDIF 409 409 … … 439 439 IF( ln_trdtrc(jn) ) THEN 440 440 tml_sum_trc(:,:,jn) = tmlbn_trc(:,:,jn) + 2 * ( tml_sum_trc(:,:,jn) - tml_trc(:,:,jn) ) + tml_trc(:,:,jn) 441 ztmltot2 (:,:,jn) = ( tml_sum_trc(:,:,jn) - tml_sumb_trc(:,:,jn) ) / ( 2.*r dt ) ! now tracer unit is /sec441 ztmltot2 (:,:,jn) = ( tml_sum_trc(:,:,jn) - tml_sumb_trc(:,:,jn) ) / ( 2.*rn_Dt ) ! now tracer unit is /sec 442 442 ENDIF 443 443 END DO … … 852 852 CALL ctl_stop( 'STOP', 'trd_mxl_trc : this was never checked. Comment this line to proceed...' ) 853 853 ENDIF 854 zsto = nn_trd_trc * r dt854 zsto = nn_trd_trc * rn_Dt 855 855 clop = "inst("//TRIM(clop)//")" 856 856 # else 857 857 IF( ln_trdmxl_trc_instant ) THEN 858 zsto = r dt ! inst. diags : we use IOIPSL time averaging858 zsto = rn_Dt ! inst. diags : we use IOIPSL time averaging 859 859 ELSE 860 zsto = nn_trd_trc * r dt ! mean diags : we DO NOT use any IOIPSL time averaging860 zsto = nn_trd_trc * rn_Dt ! mean diags : we DO NOT use any IOIPSL time averaging 861 861 ENDIF 862 862 clop = "ave("//TRIM(clop)//")" 863 863 # endif 864 zout = nn_trd_trc * r dt864 zout = nn_trd_trc * rn_Dt 865 865 iiter = nittrc000 - 1 866 866 … … 869 869 ! II.2 Compute julian date from starting date of the run 870 870 ! ------------------------------------------------------ 871 CALL ymds2ju( nyear, nmonth, nday, r dt, zjulian )871 CALL ymds2ju( nyear, nmonth, nday, rn_Dt, zjulian ) 872 872 zjulian = zjulian - adatrj ! set calendar origin to the beginning of the experiment 873 873 IF(lwp) WRITE(numout,*)' ' … … 901 901 CALL dia_nam( clhstnam, nn_trd_trc, csuff ) 902 902 CALL histbeg( clhstnam, jpi, glamt, jpj, gphit, & 903 & 1, jpi, 1, jpj, iiter, zjulian, r dt, nh_t(jn), nidtrd(jn), domain_id=nidom, snc4chunks=snc4set )903 & 1, jpi, 1, jpj, iiter, zjulian, rn_Dt, nh_t(jn), nidtrd(jn), domain_id=nidom, snc4chunks=snc4set ) 904 904 905 905 !-- Define the ML depth variable … … 938 938 CALL histdef(nidtrd(jn), trim(clvar)//trim(ctrd_trc(jl,2)), clmxl//" "//clvar//ctrd_trc(jl,1), & 939 939 & cltrcu, jpi, jpj, nh_t(jn), 1 , 1, 1 , -99 , 32, clop, zsto, zout ) ! IOIPSL: time mean 940 END DO ! if zsto=r dt above940 END DO ! if zsto=rn_Dt above 941 941 942 942 CALL histdef(nidtrd(jn), trim(clvar)//trim(ctrd_trc(jpmxl_trc_radb,2)), clmxl//" "//clvar//ctrd_trc(jpmxl_trc_radb,1), &
Note: See TracChangeset
for help on using the changeset viewer.