- Timestamp:
- 2016-09-23T11:57:08+02:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2015/nemo_v3_6_STABLE/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zlim.F90
r6204 r6943 44 44 REAL(wp), PUBLIC :: xkdoc !: 2nd half-sat. of DOC remineralization 45 45 REAL(wp), PUBLIC :: concbfe !: Fe half saturation for bacteria 46 REAL(wp), PUBLIC :: oxymin !: half saturation constant for anoxia 46 47 REAL(wp), PUBLIC :: qnfelim !: optimal Fe quota for nanophyto 47 48 REAL(wp), PUBLIC :: qdfelim !: optimal Fe quota for diatoms … … 121 122 zlim1 = xnanono3(ji,jj,jk) + xnanonh4(ji,jj,jk) 122 123 zlim2 = trb(ji,jj,jk,jppo4) / ( trb(ji,jj,jk,jppo4) + concbnh4 ) 123 zlim3 = trb(ji,jj,jk,jpfer) / ( concbfe + trb(ji,jj,jk,jpfer) )124 zlim3 = biron(ji,jj,jk) / ( concbfe + biron(ji,jj,jk) ) 124 125 zlim4 = trb(ji,jj,jk,jpdoc) / ( xkdoc + trb(ji,jj,jk,jpdoc) ) 125 126 xlimbacl(ji,jj,jk) = MIN( zlim1, zlim2, zlim3 ) … … 187 188 END DO 188 189 ! 190 DO jk = 1, jpkm1 191 DO jj = 1, jpj 192 DO ji = 1, jpi 193 ! denitrification factor computed from O2 levels 194 nitrfac(ji,jj,jk) = MAX( 0.e0, 0.4 * ( 6.e-6 - trb(ji,jj,jk,jpoxy) ) & 195 & / ( oxymin + trb(ji,jj,jk,jpoxy) ) ) 196 nitrfac(ji,jj,jk) = MIN( 1., nitrfac(ji,jj,jk) ) 197 END DO 198 END DO 199 END DO 189 200 ! 190 201 IF( lk_iomput .AND. knt == nrdttrc ) THEN ! save output diagnostics … … 216 227 NAMELIST/nampislim/ concnno3, concdno3, concnnh4, concdnh4, concnfer, concdfer, concbfe, & 217 228 & concbno3, concbnh4, xsizedia, xsizephy, xsizern, xsizerd, & 218 & xksi1, xksi2, xkdoc, qnfelim, qdfelim, caco3r 229 & xksi1, xksi2, xkdoc, qnfelim, qdfelim, caco3r, oxymin 219 230 INTEGER :: ios ! Local integer output status for namelist read 220 231 … … 249 260 WRITE(numout,*) ' Minimum size criteria for nanophyto xsizephy = ', xsizephy 250 261 WRITE(numout,*) ' Fe half saturation for bacteria concbfe = ', concbfe 262 WRITE(numout,*) ' halk saturation constant for anoxia oxymin =' , oxymin 251 263 WRITE(numout,*) ' optimal Fe quota for nano. qnfelim = ', qnfelim 252 264 WRITE(numout,*) ' Optimal Fe quota for diatoms qdfelim = ', qdfelim 253 265 ENDIF 254 266 ! 267 nitrfac (:,:,:) = 0._wp 268 ! 255 269 END SUBROUTINE p4z_lim_init 256 270
Note: See TracChangeset
for help on using the changeset viewer.