- Timestamp:
- 2018-12-04T16:19:46+01:00 (6 years ago)
- Location:
- NEMO/branches/2018/dev_r10164_HPC09_ESIWACE_PREP_MERGE/src/TOP
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2018/dev_r10164_HPC09_ESIWACE_PREP_MERGE/src/TOP/TRP/trcsbc.F90
r10170 r10372 113 113 114 114 ! 0. initialization 115 DO jn = 1, jptra 116 ! 117 IF( l_trdtrc ) ztrtrd(:,:,:) = tra(:,:,:,jn) ! save trends 118 ! 119 IF( nn_ice_tr == -1 ) THEN ! No tracers in sea ice (null concentration in sea ice) 115 IF( nn_ice_tr == -1 ) THEN ! No tracers in sea ice (null concentration in sea ice) 116 ! 117 DO jn = 1, jptra 120 118 DO jj = 2, jpj 121 119 DO ji = fs_2, fs_jpim1 ! vector opt. … … 123 121 END DO 124 122 END DO 125 ELSE 123 END DO 124 ! 125 ELSE 126 ! 127 DO jn = 1, jptra 126 128 DO jj = 2, jpj 127 129 DO ji = fs_2, fs_jpim1 ! vector opt. … … 143 145 END DO 144 146 END DO 145 ENDIF 146 ! 147 CALL lbc_lnk( 'trcsbc', sbc_trc(:,:,jn), 'T', 1. ) 148 ! Concentration dilution effect on tracers due to evaporation & precipitation 147 END DO 148 ENDIF 149 ! 150 CALL lbc_lnk( 'trcsbc', sbc_trc(:,:,:), 'T', 1. ) 151 ! Concentration dilution effect on tracers due to evaporation & precipitation 152 DO jn = 1, jptra 153 ! 154 IF( l_trdtrc ) ztrtrd(:,:,:) = tra(:,:,:,jn) ! save trends 155 ! 149 156 DO jj = 2, jpj 150 157 DO ji = fs_2, fs_jpim1 ! vector opt. -
NEMO/branches/2018/dev_r10164_HPC09_ESIWACE_PREP_MERGE/src/TOP/trcini.F90
r10345 r10372 72 72 CALL trc_ice_ini ! Tracers in sea ice 73 73 ! 74 IF(lw p) CALL ctl_opn( numstr, 'tracer.stat', 'REPLACE', 'FORMATTED', 'SEQUENTIAL', -1, numout, lwp , narea )74 IF(lwm) CALL ctl_opn( numstr, 'tracer.stat', 'REPLACE', 'FORMATTED', 'SEQUENTIAL', -1, numout, lwp , narea ) 75 75 ! 76 76 CALL trc_ini_state ! passive tracers initialisation : from a restart or from clim -
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.