Changeset 7177 for branches/2016/dev_r7012_ROBUST5_CNRS/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zligand.F90
- Timestamp:
- 2016-11-03T16:29:51+01:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2016/dev_r7012_ROBUST5_CNRS/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zligand.F90
r7174 r7177 63 63 ! This is based on the idea that as LGW is lower 64 64 ! there is a larger fraction of refractory OM 65 zlgwr = max( rlgs , rlgw * exp( -2 * (tr n(ji,jj,jk,jplgw)*1e9) ) ) ! years66 zlgwr = 1. / zlgwr * tgfunc(ji,jj,jk) * ( xstep / nyear_len(1) ) * tr n(ji,jj,jk,jplgw)65 zlgwr = max( rlgs , rlgw * exp( -2 * (trb(ji,jj,jk,jplgw)*1e9) ) ) ! years 66 zlgwr = 1. / zlgwr * tgfunc(ji,jj,jk) * ( xstep / nyear_len(1) ) * trb(ji,jj,jk,jplgw) 67 67 ! photochem loss of weak ligand 68 zlgwpr = prlgw * xstep * etot(ji,jj,jk) * tr n(ji,jj,jk,jplgw) * (1. - fr_i(ji,jj))68 zlgwpr = prlgw * xstep * etot(ji,jj,jk) * trb(ji,jj,jk,jplgw) * (1. - fr_i(ji,jj)) 69 69 tra(ji,jj,jk,jplgw) = tra(ji,jj,jk,jplgw) + zlgwp - zlgwr - zlgwpr 70 70 END DO … … 84 84 zrfepa = rfep * ( 1. - EXP( -1. * etot(ji,jj,jk) / 25. ) ) * (1.- fr_i(ji,jj)) 85 85 zrfepa = MAX( (zrfepa / 10.0), zrfepa ) ! min of 10 days lifetime 86 zfepr = rfep * xstep * tr n(ji,jj,jk,jpfep)86 zfepr = rfep * xstep * trb(ji,jj,jk,jpfep) 87 87 tra(ji,jj,jk,jpfep) = tra(ji,jj,jk,jpfep) - zfepr 88 88 tra(ji,jj,jk,jpfer) = tra(ji,jj,jk,jpfer) + zfepr
Note: See TracChangeset
for help on using the changeset viewer.