Commit 08ee118e705755a33f54094ba8b741a96a543ec4

Authored by dumoda01
1 parent 098bc3aa
Exists in master and in 1 other branch snow

Prise de nitrate modifiée pour celle du modèle de Fasham

Showing 1 changed file with 44 additions and 6 deletions   Show diff stats
src/extras/bio/bio_npzd.F90
... ... @@ -99,6 +99,8 @@
99 99 REALTYPE :: gmax=5.787037e-06
100 100 REALTYPE :: iv=1.1
101 101 REALTYPE :: alpha=0.3
  102 + REALTYPE :: inib=0.03 !DD
  103 + REALTYPE :: k1=0.5 !DD
102 104 REALTYPE :: rpn=1.157407e-07
103 105 REALTYPE :: rzn=1.157407e-07
104 106 REALTYPE :: rdn=3.472222e-08
... ... @@ -143,7 +145,7 @@
143 145 namelist /bio_npzd_nml/ numc, &
144 146 n_initial,p_initial,z_initial,d_initial, &
145 147 p0,z0,w_p,w_d,kc,i_min,rmax,gmax,iv,alpha,rpn, &
146   - rzn,rdn,rpdu,rpdl,rzd,aa,g2
  148 + rzn,rdn,rpdu,rpdl,rzd,aa,g2,inib,k1
147 149 !EOP
148 150 !-----------------------------------------------------------------------
149 151 !BOC
... ... @@ -273,19 +275,19 @@
273 275 !BOC
274 276 var_names(1) = 'nut'
275 277 var_units(1) = 'mmol/m**2'
276   - var_long(1) = 'nutrients'
  278 + var_long(1) = 'nutrients'
277 279  
278 280 var_names(2) = 'phy'
279 281 var_units(2) = 'mmol/m**2'
280   - var_long(2) = 'phytoplankton'
  282 + var_long(2) = 'phytoplankton'
281 283  
282 284 var_names(3) = 'zoo'
283 285 var_units(3) = 'mmol/m**2'
284   - var_long(3) = 'zooplankton'
  286 + var_long(3) = 'zooplankton'
285 287  
286 288 var_names(4) = 'det'
287 289 var_units(4) = 'mmol/m**2'
288   - var_long(4) = 'detritus'
  290 + var_long(4) = 'detritus'
289 291  
290 292 return
291 293 end subroutine var_info_npzd
... ... @@ -323,6 +325,41 @@
323 325 !-----------------------------------------------------------------------
324 326 !BOP
325 327 !
  328 +! !IROUTINE: Platt et al. (1980) and quota formulation for nutrient uptake
  329 +!
  330 +! !INTERFACE:
  331 + REALTYPE function fnp2(n,p,par)
  332 +!
  333 +! !DESCRIPTION:
  334 +! Here, the photoresponse of Platt et al. (1980) and quota for nutrient uptake
  335 +! is formulated.
  336 +!
  337 +! !USES:
  338 + IMPLICIT NONE
  339 +!
  340 +! !INPUT PARAMETERS:
  341 + REALTYPE, intent(in) :: n,p,par,iopt
  342 + REALTYPE :: Ps,ff,q1
  343 +!
  344 +! !REVISION HISTORY:
  345 +! Original author(s): Dany Dumont
  346 +!
  347 +!EOP
  348 +!-----------------------------------------------------------------------
  349 +!BOC
  350 +
  351 + Ps=rmax/((alpha/(alpha+inib))*(alpha/(alpha+inib))**(inib/alpha))
  352 + ff=Ps*(1.-exp(-1.*alpha*par/Ps))*exp(-1.*inib*par/Ps)
  353 + q1=n/k1/(1.+n/k1)
  354 +
  355 + fnp2=ff*q1*(p+p0)
  356 + return
  357 + end function fnp2
  358 +!EOC
  359 +
  360 +!-----------------------------------------------------------------------
  361 +!BOP
  362 +!
326 363 ! !IROUTINE: Ivlev formulation for zooplankton grazing on phytoplankton
327 364 !
328 365 ! !INTERFACE:
... ... @@ -491,7 +528,8 @@
491 528 rpd=rpdl
492 529 end if
493 530  
494   - dd(n,p,ci)=fnp(cc(n,ci),cc(p,ci),par(ci),iopt) ! snp
  531 + !dd(n,p,ci)=fnp(cc(n,ci),cc(p,ci),par(ci),iopt) ! snp
  532 + dd(n,p,ci)=fnp2(cc(n,ci),cc(p,ci),par(ci)) ! snp
495 533 dd(p,z,ci)=fpz(cc(p,ci),cc(z,ci)) ! spz
496 534 dd(p,n,ci)=rpn*cc(p,ci) ! spn
497 535 dd(z,n,ci)=rzn*cc(z,ci) ! szn
... ...