Commit f69f65512063b2b378322de24d1f28c6799cddad

Authored by Dany Dumont
1 parent 7818ae3b
Exists in master and in 1 other branch snow

Ajout du script sverdrup.jnl dans le cas test ghent

Showing 1 changed file with 105 additions and 0 deletions   Show diff stats
test_cases/ghent/sverdrup.jnl 0 → 100644
... ... @@ -0,0 +1,105 @@
  1 +let zooc=zoo*84
  2 +let phytoC=phy*79.5
  3 +let detc=det*79.5
  4 +let bacc=bac*79.5
  5 +
  6 +!TOTF=total food= r1*Phytoc+r2*Bacc+r3*Detc
  7 +
  8 +
  9 +!Parameters definition.
  10 +
  11 +define symbol r1=0.55
  12 +define symbol r2=0.4
  13 +define symbol r3=0.05
  14 +
  15 +let TOTF=($r1)*phytoc+($r2)*bacc+($r3)*detc
  16 +
  17 +
  18 +!GRAZ=gmax*Zooc*r1*phytoc/(k3+TOTF)
  19 +
  20 +define symbol gmax=0.7
  21 +define symbol k3=1.0
  22 +
  23 +let GRAZ=($gmax)*zooc*($r1)*phy/(($k3)+TOTF)
  24 +
  25 +
  26 +!MORT=mu1*phytoc
  27 +
  28 +define symbol mu1=0.05
  29 +
  30 +let MORT=($mu1)*phytoc
  31 +
  32 +!EXUD=gamma*phytoc
  33 +
  34 +define symbol gamma=0.05
  35 +
  36 +let EXUD=($gamma)*phytoc
  37 +
  38 +let LOSS=(GRAZ+MORT+EXUD)
  39 +
  40 +let loss_day=loss[l=@sbx:24]
  41 +
  42 +!Define the primary production.
  43 +
  44 +!Define Ps.
  45 +
  46 +define symbol alpha=0.1
  47 +define symbol inib=0.001
  48 +define symbol vp=0.8
  49 +
  50 +define symbol ps=($vp)/((($alpha)/(($alpha)+($inib)))*(($inib)/(($alpha)+($inib))))^(($inib)/($alpha))
  51 +
  52 +
  53 +!Photosynthesis equation.(Platt 1980)
  54 +
  55 +let prod=($ps)*(1-exp(-1*($alpha)*par/($ps)))*exp(-1*($inib)*par/($ps))
  56 +
  57 +
  58 +!Estimating daily production (no night effect)
  59 +
  60 +let prod_day=prod[l=@sbx:24]
  61 +
  62 +
  63 +!Define the primary production, total, new, and regenerated.
  64 +
  65 +let phytoC=phy*79.5
  66 +
  67 +define symbol p0=0.0001
  68 +define symbol k1=1.0
  69 +define symbol k2=0.23
  70 +
  71 +let new_prod=prod_day*(phytoc+($p0))*(nit/(($k1)*(1+(nit/($k1))+(amm/($k2)))))
  72 +let reg_prod=prod_day*(phytoc+($p0))*(amm/(($k2)*(1+(nit/($k1))+(amm/($k2)))))
  73 +
  74 +let tot_prod = new_prod + reg_prod
  75 +
  76 +
  77 +
  78 +!Estimating the difference between production and losses, and integrating the values
  79 +
  80 +let prod_net=tot_prod-loss_day
  81 +
  82 +!let prod_intZ=tot_prod[k=@din] ! integrating in the z axis, to get mgC m-2 d-1
  83 +
  84 +!let loss_day_intZ=loss_day[k=@din]
  85 +
  86 +!let prod_net_intz=prod_intZ-loss_day_intZ
  87 +
  88 +! Z(c)
  89 +
  90 +let Zc = prod_net[k=@DIN] - prod_net[k=@IIN]
  91 +
  92 +let compens = prod_net[Z=@LOC:0] ! But integration is missing...?
  93 +
  94 +list/clobber/file=zc_pc1_0910_plus10.dat zc, prod_net, tot_prod, loss_day, zsbl! to write the values to a file
  95 +
  96 +
  97 +plot/nolabel/line=8 zsbl[l=1651:4128@sbx:24]
  98 +plot/nolabel/over/line=7 zc[l=@sbx:3]
  99 +!frame/file=Zc_Zt_0910_E1.gif
  100 +
  101 +
  102 +let taus = (x_taus^2 + y_taus^2)^0.5
  103 +plot/nolabel/axes=0,1,1,0/vlimit=0:1.8:0.2/line=10 taus[l=1650:4180]
  104 +!frame/file=WindStress_0910.gif
  105 +
... ...