- Timestamp:
- 2020-11-10T12:57:08+01:00 (4 years ago)
- Location:
- NEMO/branches/2020/dev_12905_xios_ancil
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2020/dev_12905_xios_ancil
- Property svn:externals
-
old new 3 3 ^/utils/build/mk@HEAD mk 4 4 ^/utils/tools@HEAD tools 5 ^/vendors/AGRIF/dev @HEADext/AGRIF5 ^/vendors/AGRIF/dev_r12970_AGRIF_CMEMS ext/AGRIF 6 6 ^/vendors/FCM@HEAD ext/FCM 7 7 ^/vendors/IOIPSL@HEAD ext/IOIPSL 8 8 9 9 # SETTE 10 ^/utils/CI/sette@ HEADsette10 ^/utils/CI/sette@13559 sette
-
- Property svn:externals
-
NEMO/branches/2020/dev_12905_xios_ancil/src/TOP/PISCES/P2Z/p2zexp.F90
r12489 r13766 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 … … 39 39 !! * Substitutions 40 40 # include "do_loop_substitute.h90" 41 # include "domzgr_substitute.h90" 41 42 !!---------------------------------------------------------------------- 42 43 !! NEMO/TOP 4.0 , NEMO Consortium (2018) … … 81 82 ! LAYERS IS DETERMINED BY DMIN3 DEFINED IN sms_p2z.F90 82 83 ! ---------------------------------------------------------------------- 83 DO_3D _00_00(1, jpkm1 )84 DO_3D( 0, 0, 0, 0, 1, jpkm1 ) 84 85 ze3t = 1. / e3t(ji,jj,jk,Kmm) 85 86 tr(ji,jj,jk,jpno3,Krhs) = tr(ji,jj,jk,jpno3,Krhs) + ze3t * dmin3(ji,jj,jk) * xksi(ji,jj) … … 92 93 zgeolpoc = 0.e0 ! Initialization 93 94 ! Release of nutrients from the "simple" sediment 94 DO_2D _00_0095 DO_2D( 0, 0, 0, 0 ) 95 96 ikt = mbkt(ji,jj) 96 97 tr(ji,jj,ikt,jpno3,Krhs) = tr(ji,jj,ikt,jpno3,Krhs) + sedlam * sedpocn(ji,jj) / e3t(ji,jj,ikt,Kmm) … … 102 103 END_2D 103 104 104 DO_2D _00_00105 DO_2D( 0, 0, 0, 0 ) 105 106 tr(ji,jj,1,jpno3,Krhs) = tr(ji,jj,1,jpno3,Krhs) + zgeolpoc * cmask(ji,jj) / areacot / e3t(ji,jj,1,Kmm) 106 107 END_2D 107 108 108 CALL lbc_lnk( 'p2zexp', sedpocn, 'T', 1. )109 CALL lbc_lnk( 'p2zexp', sedpocn, 'T', 1.0_wp ) 109 110 110 111 ! Oa & Ek: diagnostics depending on jpdia2d ! left as example … … 120 121 ELSE 121 122 ! 122 DO_2D _11_11123 DO_2D( 1, 1, 1, 1 ) 123 124 zsedpocd = zsedpoca(ji,jj) - 2. * sedpocn(ji,jj) + sedpocb(ji,jj) ! time laplacian on tracers 124 125 sedpocb(ji,jj) = sedpocn(ji,jj) + rn_atfp * zsedpocd ! sedpocb <-- filtered sedpocn … … 139 140 IF(sn_cfctl%l_prttrc) THEN ! print mean trends (used for debugging) 140 141 WRITE(charout, FMT="('exp')") 141 CALL prt_ctl_ trc_info(charout)142 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) 143 144 ENDIF 144 145 ! … … 173 174 zdm0 = 0._wp 174 175 zrro = 1._wp 175 DO_3D _11_11(jpkb, jpkm1 )176 DO_3D( 1, 1, 1, 1, jpkb, jpkm1 ) 176 177 zfluo = ( gdepw(ji,jj,jk ,Kmm) / gdepw(ji,jj,jpkb,Kmm) )**xhr 177 178 zfluu = ( gdepw(ji,jj,jk+1,Kmm) / gdepw(ji,jj,jpkb,Kmm) )**xhr … … 190 191 dminl(:,:) = 0._wp 191 192 dmin3(:,:,:) = zdm0 192 DO_3D _11_11(1, jpk )193 DO_3D( 1, 1, 1, 1, 1, jpk ) 193 194 IF( tmask(ji,jj,jk) == 0._wp ) THEN 194 195 dminl(ji,jj) = dminl(ji,jj) + dmin3(ji,jj,jk) … … 197 198 END_3D 198 199 199 DO_2D _11_11200 DO_2D( 1, 1, 1, 1 ) 200 201 IF( tmask(ji,jj,1) == 0 ) dmin3(ji,jj,1) = 0._wp 201 202 END_2D … … 203 204 ! Coastal mask 204 205 cmask(:,:) = 0._wp 205 DO_2D _00_00206 DO_2D( 0, 0, 0, 0 ) 206 207 IF( tmask(ji,jj,1) /= 0. ) THEN 207 208 zmaskt = tmask(ji+1,jj,1) * tmask(ji-1,jj,1) * tmask(ji,jj+1,1) * tmask(ji,jj-1,1) … … 209 210 END IF 210 211 END_2D 211 CALL lbc_lnk( 'p2zexp', cmask , 'T', 1. ) ! lateral boundary conditions on cmask (sign unchanged)212 CALL lbc_lnk( 'p2zexp', cmask , 'T', 1.0_wp ) ! lateral boundary conditions on cmask (sign unchanged) 212 213 areacot = glob_sum( 'p2zexp', e1e2t(:,:) * cmask(:,:) ) 213 214 ! 214 215 IF( ln_rsttr ) THEN 215 CALL iom_get( numrtr, jpdom_auto glo, 'SEDB'//ctrcnm(jpdet), sedpocb(:,:) )216 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(:,:) ) 217 218 ELSE 218 219 sedpocb(:,:) = 0._wp
Note: See TracChangeset
for help on using the changeset viewer.