- Timestamp:
- 2020-11-27T17:26:33+01:00 (4 years ago)
- Location:
- NEMO/branches/2020/tickets_icb_1900
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2020/tickets_icb_1900
- Property svn:externals
-
NEMO/branches/2020/tickets_icb_1900/src/TOP/PISCES/P2Z/p2zexp.F90
r13237 r13899 17 17 USE p2zsed 18 18 USE lbclnk 19 USE prtctl _trc! Print control for debbuging19 USE prtctl ! Print control for debbuging 20 20 USE trd_oce 21 21 USE trdtrc … … 82 82 ! LAYERS IS DETERMINED BY DMIN3 DEFINED IN sms_p2z.F90 83 83 ! ---------------------------------------------------------------------- 84 DO_3D _00_00(1, jpkm1 )84 DO_3D( 0, 0, 0, 0, 1, jpkm1 ) 85 85 ze3t = 1. / e3t(ji,jj,jk,Kmm) 86 86 tr(ji,jj,jk,jpno3,Krhs) = tr(ji,jj,jk,jpno3,Krhs) + ze3t * dmin3(ji,jj,jk) * xksi(ji,jj) … … 93 93 zgeolpoc = 0.e0 ! Initialization 94 94 ! Release of nutrients from the "simple" sediment 95 DO_2D _00_0095 DO_2D( 0, 0, 0, 0 ) 96 96 ikt = mbkt(ji,jj) 97 97 tr(ji,jj,ikt,jpno3,Krhs) = tr(ji,jj,ikt,jpno3,Krhs) + sedlam * sedpocn(ji,jj) / e3t(ji,jj,ikt,Kmm) … … 103 103 END_2D 104 104 105 DO_2D _00_00105 DO_2D( 0, 0, 0, 0 ) 106 106 tr(ji,jj,1,jpno3,Krhs) = tr(ji,jj,1,jpno3,Krhs) + zgeolpoc * cmask(ji,jj) / areacot / e3t(ji,jj,1,Kmm) 107 107 END_2D … … 121 121 ELSE 122 122 ! 123 DO_2D _11_11123 DO_2D( 1, 1, 1, 1 ) 124 124 zsedpocd = zsedpoca(ji,jj) - 2. * sedpocn(ji,jj) + sedpocb(ji,jj) ! time laplacian on tracers 125 125 sedpocb(ji,jj) = sedpocn(ji,jj) + rn_atfp * zsedpocd ! sedpocb <-- filtered sedpocn … … 140 140 IF(sn_cfctl%l_prttrc) THEN ! print mean trends (used for debugging) 141 141 WRITE(charout, FMT="('exp')") 142 CALL prt_ctl_ trc_info(charout)143 CALL prt_ctl _trc(tab4d=tr(:,:,:,:,Krhs), mask=tmask, clinfo=ctrcnm)142 CALL prt_ctl_info( charout, cdcomp = 'top' ) 143 CALL prt_ctl(tab4d_1=tr(:,:,:,:,Krhs), mask1=tmask, clinfo=ctrcnm) 144 144 ENDIF 145 145 ! … … 174 174 zdm0 = 0._wp 175 175 zrro = 1._wp 176 DO_3D _11_11(jpkb, jpkm1 )176 DO_3D( 1, 1, 1, 1, jpkb, jpkm1 ) 177 177 zfluo = ( gdepw(ji,jj,jk ,Kmm) / gdepw(ji,jj,jpkb,Kmm) )**xhr 178 178 zfluu = ( gdepw(ji,jj,jk+1,Kmm) / gdepw(ji,jj,jpkb,Kmm) )**xhr … … 191 191 dminl(:,:) = 0._wp 192 192 dmin3(:,:,:) = zdm0 193 DO_3D _11_11(1, jpk )193 DO_3D( 1, 1, 1, 1, 1, jpk ) 194 194 IF( tmask(ji,jj,jk) == 0._wp ) THEN 195 195 dminl(ji,jj) = dminl(ji,jj) + dmin3(ji,jj,jk) … … 198 198 END_3D 199 199 200 DO_2D _11_11200 DO_2D( 1, 1, 1, 1 ) 201 201 IF( tmask(ji,jj,1) == 0 ) dmin3(ji,jj,1) = 0._wp 202 202 END_2D … … 204 204 ! Coastal mask 205 205 cmask(:,:) = 0._wp 206 DO_2D _00_00206 DO_2D( 0, 0, 0, 0 ) 207 207 IF( tmask(ji,jj,1) /= 0. ) THEN 208 208 zmaskt = tmask(ji+1,jj,1) * tmask(ji-1,jj,1) * tmask(ji,jj+1,1) * tmask(ji,jj-1,1) … … 214 214 ! 215 215 IF( ln_rsttr ) THEN 216 CALL iom_get( numrtr, jpdom_auto glo, 'SEDB'//ctrcnm(jpdet), sedpocb(:,:) )217 CALL iom_get( numrtr, jpdom_auto glo, 'SEDN'//ctrcnm(jpdet), sedpocn(:,:) )216 CALL iom_get( numrtr, jpdom_auto, 'SEDB'//ctrcnm(jpdet), sedpocb(:,:) ) 217 CALL iom_get( numrtr, jpdom_auto, 'SEDN'//ctrcnm(jpdet), sedpocn(:,:) ) 218 218 ELSE 219 219 sedpocb(:,:) = 0._wp
Note: See TracChangeset
for help on using the changeset viewer.