Changeset 14957 for NEMO/branches/2021
- Timestamp:
- 2021-06-07T14:26:34+02:00 (3 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2021/dev_r14383_PISCES_NEWDEV_PISCO/src/TOP/PISCES/P4Z/p4zfechem.F90
r14786 r14957 59 59 REAL(wp) :: zscave, zaggdfea, zaggdfeb, ztrc, zdust, zklight 60 60 REAL(wp) :: ztfe, zhplus, zxlam, zaggliga, zaggligb 61 REAL(wp) :: zprecip, zprecipno3, zvar61 REAL(wp) :: zprecip, zprecipno3, zconsfe, za1 62 62 REAL(wp) :: zrfact2 63 63 CHARACTER (len=25) :: charout … … 98 98 zTL1(ji,jj,jk) = ztotlig(ji,jj,jk) 99 99 zkeq = fekeq(ji,jj,jk) 100 zklight = 4.77E-7 * etot(ji,jj,jk) * 0.5 / 10**-6.3 100 zklight = 4.77E-7 * etot(ji,jj,jk) * 0.5 / ( 10**(-6.3) ) 101 zconsfe = consfe3(ji,jj,jk) / ( 10**(-6.3) ) 101 102 zfesatur = zTL1(ji,jj,jk) * 1E-9 102 103 ztfe = (1.0 + zklight) * tr(ji,jj,jk,jpfer,Kbb) 103 104 ! Fe' is the root of a 2nd order polynom 104 z var = 1. + zfesatur * zkeq + zklight + consfe3(ji,jj,jk)/10**-6.3- zkeq * tr(ji,jj,jk,jpfer,Kbb)105 zFe3 (ji,jj,jk) = ( -1 * z var + SQRT( zvar**2 + 4. * ztfe * zkeq) ) / ( 2. * zkeq + rtrn )105 za1 = 1. + zfesatur * zkeq + zklight + zconsfe - zkeq * tr(ji,jj,jk,jpfer,Kbb) 106 zFe3 (ji,jj,jk) = ( -1 * za1 + SQRT( za1**2 + 4. * ztfe * zkeq) ) / ( 2. * zkeq + rtrn ) 106 107 zFeL1(ji,jj,jk) = MAX( 0., tr(ji,jj,jk,jpfer,Kbb) - zFe3(ji,jj,jk) ) 107 108 END_3D … … 264 265 WRITE(numout,*) ' Scavenged iron that is added to POFe scaveff =', scaveff 265 266 ENDIF 266 ! 267 ! 267 268 END SUBROUTINE p4z_fechem_init 268 269
Note: See TracChangeset
for help on using the changeset viewer.