- Timestamp:
- 2018-12-04T16:19:46+01:00 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2018/dev_r10164_HPC09_ESIWACE_PREP_MERGE/src/TOP/trcstp.F90
r10314 r10372 20 20 USE trdtrc_oce 21 21 USE trdmxl_trc 22 USE sms_pisces, ONLY : ln_check_mass 22 23 ! 23 24 USE prtctl_trc ! Print control for debbuging … … 73 74 cvol(:,:,jk) = e1e2t(:,:) * e3t_n(:,:,jk) * tmask(:,:,jk) 74 75 END DO 75 areatot = glob_sum( 'trcstp', cvol(:,:,:) ) 76 IF ( ln_ctl .OR. kt == nitrst .OR. ( ln_check_mass .AND. kt == nitend ) & 77 & .OR. iom_use( "pno3tot" ) .OR. iom_use( "ppo4tot" ) .OR. iom_use( "psiltot" ) & 78 & .OR. iom_use( "palktot" ) .OR. iom_use( "pfertot" ) ) & 79 & areatot = glob_sum( 'trcstp', cvol(:,:,:) ) 76 80 ENDIF 77 81 ! … … 105 109 ENDIF 106 110 ! 107 ztrai = 0._wp ! content of all tracers 108 DO jn = 1, jptra 109 ztrai = ztrai + glob_sum( 'trcstp', trn(:,:,:,jn) * cvol(:,:,:) ) 110 END DO 111 IF( lwp ) WRITE(numstr,9300) kt, ztrai / areatot 111 IF (ln_ctl ) THEN 112 ztrai = 0._wp ! content of all tracers 113 DO jn = 1, jptra 114 ztrai = ztrai + glob_sum( 'trcstp', trn(:,:,:,jn) * cvol(:,:,:) ) 115 END DO 116 IF( lwm ) WRITE(numstr,9300) kt, ztrai / areatot 117 ENDIF 112 118 9300 FORMAT(i10,D23.16) 113 119 !
Note: See TracChangeset
for help on using the changeset viewer.