Changeset 10416 for NEMO/trunk/src/TOP/PISCES/P4Z/p4zligand.F90
- Timestamp:
- 2018-12-19T12:45:43+01:00 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/trunk/src/TOP/PISCES/P4Z/p4zligand.F90
r10362 r10416 25 25 REAL(wp), PUBLIC :: rlig !: Remin ligand production 26 26 REAL(wp), PUBLIC :: prlgw !: Photochemical of weak ligand 27 REAL(wp), PUBLIC :: rfep !: Dissolution rate of FeP28 27 29 28 !!---------------------------------------------------------------------- … … 43 42 ! 44 43 INTEGER :: ji, jj, jk 45 REAL(wp) :: zlgwp, zlgwpr, zlgwr, zlablgw , zrfepa, zfepr44 REAL(wp) :: zlgwp, zlgwpr, zlgwr, zlablgw 46 45 REAL(wp), DIMENSION(jpi,jpj,jpk) :: zligrem, zligpr, zrligprod 47 46 REAL(wp), ALLOCATABLE, DIMENSION(:,:,:) :: zw3d … … 71 70 zligpr(ji,jj,jk) = zlgwpr 72 71 zrligprod(ji,jj,jk) = zlgwp 73 !74 ! ----------------------------------------------------------75 ! Dissolution of nanoparticle Fe76 ! ----------------------------------------------------------77 ! dissolution rate is maximal in the presence of light and78 ! lower in the aphotici zone79 ! ! 25 Wm-2 constant80 zrfepa = rfep * ( 1. - EXP( -1. * etot(ji,jj,jk) / 25. ) ) * (1.- fr_i(ji,jj))81 zrfepa = MAX( (zrfepa / 10.0), zrfepa ) ! min of 10 days lifetime82 zfepr = rfep * xstep * trb(ji,jj,jk,jpfep)83 tra(ji,jj,jk,jpfep) = tra(ji,jj,jk,jpfep) - zfepr84 tra(ji,jj,jk,jpfer) = tra(ji,jj,jk,jpfer) + zfepr85 72 ! 86 73 END DO … … 130 117 INTEGER :: ios ! Local integer 131 118 ! 132 NAMELIST/nampislig/ rlgw, prlgw, rlgs, r fep, rlig119 NAMELIST/nampislig/ rlgw, prlgw, rlgs, rlig 133 120 !!---------------------------------------------------------------------- 134 121 ! … … 148 135 IF(lwp) THEN ! control print 149 136 WRITE(numout,*) ' Namelist : nampislig' 150 WRITE(numout,*) ' Dissolution rate of FeP rfep =', rfep151 137 WRITE(numout,*) ' Lifetime (years) of weak ligands rlgw =', rlgw 152 138 WRITE(numout,*) ' Remin ligand production per unit C rlig =', rlig
Note: See TracChangeset
for help on using the changeset viewer.