- Timestamp:
- 2020-11-05T15:18:53+01:00 (4 years ago)
- Location:
- NEMO/branches/2020/dev_12905_xios_restart
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2020/dev_12905_xios_restart
- 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_restart/src/TOP/PISCES/P4Z/p4zfechem.F90
r12377 r13727 16 16 USE p4zche ! chemical model 17 17 USE p4zbc ! Boundary conditions from sediments 18 USE prtctl _trc! print control for debugging18 USE prtctl ! print control for debugging 19 19 USE iom ! I/O manager 20 20 … … 33 33 !! * Substitutions 34 34 # include "do_loop_substitute.h90" 35 # include "domzgr_substitute.h90" 35 36 !!---------------------------------------------------------------------- 36 37 !! NEMO/TOP 4.0 , NEMO Consortium (2018) … … 91 92 ! Chemistry is supposed to be fast enough to be at equilibrium 92 93 ! ------------------------------------------------------------ 93 DO_3D _11_11(1, jpkm1 )94 DO_3D( 1, 1, 1, 1, 1, jpkm1 ) 94 95 zTL1(ji,jj,jk) = ztotlig(ji,jj,jk) 95 96 zkeq = fekeq(ji,jj,jk) … … 106 107 107 108 zdust = 0. ! if no dust available 108 DO_3D _11_11(1, jpkm1 )109 DO_3D( 1, 1, 1, 1, 1, jpkm1 ) 109 110 ! Scavenging rate of iron. This scavenging rate depends on the load of particles of sea water. 110 111 ! This parameterization assumes a simple second order kinetics (k[Particles][Fe]). … … 117 118 ! 118 119 zfeequi = zFe3(ji,jj,jk) * 1E-9 119 zhplus = max( rtrn, hi(ji,jj,jk) )120 fe3sol = fesol(ji,jj,jk,1) * ( zhplus**3 + fesol(ji,jj,jk,2) * zhplus**2 &121 & + fesol(ji,jj,jk,3) * zhplus + fesol(ji,jj,jk,4) &122 & + fesol(ji,jj,jk,5) / zhplus )123 120 zfecoll = 0.5 * zFeL1(ji,jj,jk) * 1E-9 124 121 ! precipitation of Fe3+, creation of nanoparticles … … 176 173 IF( ln_ligand ) THEN 177 174 ! 178 DO_3D _11_11(1, jpkm1 )175 DO_3D( 1, 1, 1, 1, 1, jpkm1 ) 179 176 zlam1a = ( 0.369 * 0.3 * tr(ji,jj,jk,jpdoc,Kbb) + 102.4 * tr(ji,jj,jk,jppoc,Kbb) ) * xdiss(ji,jj,jk) & 180 177 & + ( 114. * 0.3 * tr(ji,jj,jk,jpdoc,Kbb) ) … … 221 218 IF(sn_cfctl%l_prttrc) THEN ! print mean trends (used for debugging) 222 219 WRITE(charout, FMT="('fechem')") 223 CALL prt_ctl_ trc_info(charout)224 CALL prt_ctl _trc(tab4d=tr(:,:,:,:,Krhs), mask=tmask, clinfo=ctrcnm)220 CALL prt_ctl_info( charout, cdcomp = 'top' ) 221 CALL prt_ctl(tab4d_1=tr(:,:,:,:,Krhs), mask1=tmask, clinfo=ctrcnm) 225 222 ENDIF 226 223 !
Note: See TracChangeset
for help on using the changeset viewer.