Changeset 1334 for trunk/NEMO/TOP_SRC/TRP/trdmld_trc.F90
- Timestamp:
- 2009-03-03T15:07:48+01:00 (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/TOP_SRC/TRP/trdmld_trc.F90
r1317 r1334 406 406 !!---------------------------------------------------------------------- 407 407 INTEGER, INTENT( in ) :: kt ! ocean time-step index 408 INTEGER :: ji, jj, jk, jl, ik, it, jn 408 INTEGER :: ji, jj, jk, jl, ik, it, jn, itmod 409 409 REAL(wp) :: zavt, zfn, zfn2 410 410 !! … … 575 575 tmltrd_trc(:,:,:,:) = tmltrd_trc(:,:,:,:) * ucf_trc 576 576 577 it = kt - nit000 + 1 578 579 MODULO_NTRD : IF( MOD( it, ntrd_trc ) == 0 ) THEN ! nitend MUST be multiple of ntrd_trc 577 ! define time axis 578 it = kt 579 itmod = kt - nit000 + 1 580 581 MODULO_NTRD : IF( MOD( itmod, ntrd_trc ) == 0 ) THEN ! nitend MUST be multiple of ntrd_trc 580 582 ! 581 583 ztmltot (:,:,:) = 0.e0 ! reset arrays to zero … … 819 821 ! ---------------------------------- 820 822 821 IF( lwp .AND. MOD( it , ntrd_trc ) == 0 ) THEN823 IF( lwp .AND. MOD( itmod , ntrd_trc ) == 0 ) THEN 822 824 WRITE(numout,*) ' ' 823 825 WRITE(numout,*) 'trd_mld_trc : write passive tracer trends in the NetCDF file :' … … 930 932 # endif /* key_dimgout */ 931 933 932 IF( MOD( it , ntrd_trc ) == 0 ) THEN934 IF( MOD( itmod, ntrd_trc ) == 0 ) THEN 933 935 ! 934 936 ! Reset cumulative arrays to zero … … 999 1001 INTEGER, INTENT( in ) :: kt ! ocean time-step index 1000 1002 #if defined key_lobster 1001 INTEGER :: jl, it 1003 INTEGER :: jl, it, itmod 1002 1004 LOGICAL :: llwarn = .TRUE., lldebug = .TRUE. 1003 1005 REAL(wp), DIMENSION(jpi,jpj,jpdiabio) :: ztmltrdbio2 ! only needed for mean diagnostics … … 1059 1061 tmltrd_bio(:,:,:) = tmltrd_bio(:,:,:) * ucf_trc 1060 1062 1061 MODULO_NTRD : IF( MOD( kt, ntrd_trc ) == 0 ) THEN ! nitend MUST be multiple of ntrd 1063 ! define time axis 1064 it = kt 1065 itmod = kt - nit000 + 1 1066 1067 MODULO_NTRD : IF( MOD( itmod, ntrd_trc ) == 0 ) THEN ! nitend MUST be multiple of ntrd 1062 1068 ! 1063 1069 zfn = float(nmoymltrdbio) ; zfn2 = zfn * zfn … … 1132 1138 1133 1139 ! define time axis 1134 it = kt - nit000 + 1 1135 1136 IF( lwp .AND. MOD( it , ntrd_trc ) == 0 ) THEN 1140 it = kt 1141 itmod = kt - nit000 + 1 1142 1143 IF( lwp .AND. MOD( itmod , ntrd_trc ) == 0 ) THEN 1137 1144 WRITE(numout,*) ' ' 1138 1145 WRITE(numout,*) 'trd_mld_bio : write ML bio trends in the NetCDF file :' … … 1174 1181 # endif /* key_dimgout */ 1175 1182 1176 IF( MOD( it , ntrd_trc ) == 0 ) THEN1183 IF( MOD( itmod, ntrd_trc ) == 0 ) THEN 1177 1184 ! 1178 1185 ! III.5 Reset cumulative arrays to zero … … 1455 1462 CALL dia_nam( clhstnam, ntrd_trc, csuff ) 1456 1463 CALL histbeg( clhstnam, jpi, glamt, jpj, gphit, & 1457 & 1, jpi, 1, jpj, 0, zjulian, rdt, nh_t(jn), nidtrd(jn), domain_id=nidom )1464 & 1, jpi, 1, jpj, nit000-1, zjulian, rdt, nh_t(jn), nidtrd(jn), domain_id=nidom ) 1458 1465 1459 1466 !-- Define the ML depth variable … … 1468 1475 CALL dia_nam( clhstnam, ntrd_trc, 'trdbio' ) 1469 1476 CALL histbeg( clhstnam, jpi, glamt, jpj, gphit, & 1470 & 1, jpi, 1, jpj, 0, zjulian, rdt, nh_tb, nidtrdbio, domain_id=nidom )1477 & 1, jpi, 1, jpj, nit000-1, zjulian, rdt, nh_tb, nidtrdbio, domain_id=nidom ) 1471 1478 #endif 1472 1479
Note: See TracChangeset
for help on using the changeset viewer.