- Timestamp:
- 2017-02-21T16:46:07+01:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/trcstp.F90
r6234 r7709 55 55 !! Update the passive tracers 56 56 !!------------------------------------------------------------------- 57 58 USE dom_oce, ONLY: narea 59 57 60 INTEGER, INTENT( in ) :: kt ! ocean time-step index 58 61 INTEGER :: jk, jn ! dummy loop indices … … 129 132 ENDIF 130 133 ! 131 ztrai = 0._wp ! content of all tracers 132 DO jn = 1, jptra 133 ztrai = ztrai + glob_sum( trn(:,:,:,jn) * cvol(:,:,:) ) 134 END DO 135 IF( lwp ) WRITE(numstr,9300) kt, ztrai / areatot 136 9300 FORMAT(i10,e18.10) 134 IF (ln_ctl) THEN 135 ! The following code is very expensive since it involves multiple 136 ! reproducible global sums over all tracer fields and is potentially 137 ! called on every timestep. The results it produces are purely for 138 ! informational purposes and do not affect model evolution. 139 ! Hence we restrict its use by protecting it with the ln_ctl RTL 140 ! which should normally only be used under debugging conditions 141 ! and not in operational runs. We also need to restrict output 142 ! to the master PE since there's no point duplicating the same results 143 ! on all processors. 144 ztrai = 0._wp ! content of all tracers 145 DO jn = 1, jptra 146 ztrai = ztrai + glob_sum( trn(:,:,:,jn) * cvol(:,:,:) ) 147 END DO 148 IF( lwp ) WRITE(numstr,9300) kt, ztrai / areatot 149 9300 FORMAT(i10,e18.10) 150 ENDIF 137 151 ! 138 152 IF( nn_timing == 1 ) CALL timing_stop('trc_stp')
Note: See TracChangeset
for help on using the changeset viewer.