Changeset 276
- Timestamp:
- 2005-09-14T11:42:17+02:00 (19 years ago)
- Location:
- trunk/NEMO
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/OPA_SRC/step.F90
r268 r276 32 32 33 33 USE trcstp ! passive tracer time-stepping (trc_stp routine) 34 USE trcdia ! passive tracer diagnostics (trc_dia routine)35 USE trcrst ! restart for passive tracers36 34 37 35 USE dynhpg ! hydrostatic pressure grad. (dyn_hpg routine) … … 281 279 ! N.B. ua, va, ta, sa arrays are used as workspace in this section 282 280 !----------------------------------------------------------------------- 283 IF( kstp == nit000 ) CALL trc_dia( kstp, indic ) ! diagnostics init.284 281 285 282 CALL trc_stp( kstp ) ! time-stepping 286 283 287 CALL trc_wri( kstp ) ! outputs288 289 CALL trc_dia( kstp, indic ) ! diagnostics290 284 #endif 291 285 -
trunk/NEMO/TOP_SRC/TRP/trcstp.F90
r274 r276 13 13 USE trctrp ! passive tracers transport 14 14 USE trcsms ! passive tracers sources and sinks 15 USE trcdia ! passive tracer diagnostics (trc_dia routine) 16 USE trcrst ! restart for passive tracers 15 17 16 18 IMPLICIT NONE … … 45 47 REAL(wp) :: ztra 46 48 49 !! * local declarations 50 INTEGER :: indic ! error indicator if < 0 51 !! --------------------------------------------------------------------- 47 52 48 ! this ROUTINE is called only every ndttrc time step 49 IF( MOD( kt , ndttrc ) /= 0 ) RETURN 50 51 ! tracers: sink and source 53 indic = 1 ! reset to no error condition 52 54 53 55 54 CALL trc_sms( kt )55 56 56 DO jn = 1, jptra 57 IF(ln_ctl) THEN ! print mean trends (used for debugging) 58 ztra = SUM( tra(2:nictle,2:njctle,1:jpkm1,jn) * tmask(2:nictle,2:njctle,1:jpkm1) ) 59 WRITE(numout,*) ' trc/sms - ',ctrcnm(jn),' : ', ztra-tra_ctl(jn) 60 tra_ctl(jn) = ztra 61 ENDIF 62 END DO 57 ! diagnostics init 58 59 IF( kt == nit000 ) CALL trc_dia( kt, indic ) 60 61 ! these ROUTINES are called only every ndttrc time step 62 IF( MOD( kt , ndttrc ) == 0 ) THEN 63 64 ! tracers: sink and source 65 66 CALL trc_sms( kt ) 67 68 DO jn = 1, jptra 69 IF(ln_ctl) THEN ! print mean trends (used for debugging) 70 ztra = SUM( tra(2:nictle,2:njctle,1:jpkm1,jn) * tmask(2:nictle,2:njctle,1:jpkm1) ) 71 WRITE(numout,*) ' trc/sms - ',ctrcnm(jn),' : ', ztra-tra_ctl(jn) 72 tra_ctl(jn) = ztra 73 ENDIF 74 END DO 63 75 64 76 65 77 66 ! transport of passive tracers78 ! transport of passive tracers 67 79 68 CALL trc_trp( kt )80 CALL trc_trp( kt ) 69 81 70 DO jn = 1, jptra 71 IF(ln_ctl) THEN ! print mean trends (used for debugging) 72 ztra = SUM( tra(2:nictle,2:njctle,1:jpkm1,jn) * tmask(2:nictle,2:njctle,1:jpkm1) ) 73 WRITE(numout,*) ' trc/trp - ',ctrcnm(jn),' : ', ztra-tra_ctl(jn) 74 tra_ctl(jn) = ztra 75 ENDIF 76 END DO 82 DO jn = 1, jptra 83 IF(ln_ctl) THEN ! print mean trends (used for debugging) 84 ztra = SUM( tra(2:nictle,2:njctle,1:jpkm1,jn) * tmask(2:nictle,2:njctle,1:jpkm1) ) 85 WRITE(numout,*) ' trc/trp - ',ctrcnm(jn),' : ', ztra-tra_ctl(jn) 86 tra_ctl(jn) = ztra 87 ENDIF 88 END DO 89 90 ENDIF 91 92 93 94 ! outputs 95 96 CALL trc_wri( kt ) 97 98 99 100 ! diagnostics 101 102 CALL trc_dia( kt, indic ) 77 103 78 104
Note: See TracChangeset
for help on using the changeset viewer.