Ajout du cas test ghent

ghent/100m/fig_bio_ismer.jnl

+ ! NOAA/PMEL TMAP
+ ! FERRET v6  
+ ! Linux(g77) 2.6.9-22. - 08/25/06
+ !  5-Mar-08 11:10     
+
+use out.nc
+!set region/z=-120:0
+set memory/size=200
+set window/aspect=0.5/size=200
+ppl dfltfnt sr
+ppl axlsze 0.18 0.18
+
+! Symbols
+define symbol ns91=7.81  ! W m-2 (Nelson and Smith 1991)
+
+
+! ======================================================================
+! Zone photique
+let light_fraction = par/par[k=50]
+let photic = light_fraction[z=@loc:0.01]
+fill/nolabel/key=continuous/palette=rnb3/level=(10,100,5)(inf)/k=21:50 par; frame/file=par.gif
+
+define axis/t="1-jan-2009":"28-jun-2009":24/units=hours thours
+fill/nolabel/key=continuous/palette=rnb3/level=(10,100,5)(inf)/k=21:50 par[l=@sbx:72]; frame/file=par_sbx.gif
+
+
+! Sverdrup Critical Depth
+let light_int = par[l=@sbx:72]
+let zc  = light_int[k=@loc:($ns91)]
+let zt  = zsbl[l=@sbx:72]
+let ztzc_ratio = zt/zc
+
+plot zc
+frame/file=zc.gif
+plot zt
+frame/file=zt.gif
+
+
+let phy    = dia + fla
+let zoo    = msz + mcz
+let totaln = phy + bac + det + zoo + nit + amm + ldn
+
+! Schloss et al. (2014)
+set viewport upper
+go margins 1.0 0.0 1.5 1.5
+ppl axlint 0 5
+define axis/t="1-jan-2009":"28-jun-2009":24/units=hours thours
+fill/nolabel/key=continuous/palette=rnb2 num[gt=thours]
+plot/over/nolabel/thick=2/axes=0,0,1,1 zt
+
+set viewport lower
+go margins 0 1.0 1.5 1.5
+ppl axlint 1 5
+fill/nolabel/key=continuous/palette=no_red/level=(0.1,2,0.1)(inf)/vlimits=-40:0/axes=0,1,1,1 phy
+frame/file=zt_phy.gif
+cancel viewport
+
+set viewport upper
+go margins 1.0 0.0 1.5 1.5
+ppl axlint 0 5
+plot/nolabel/thick=2/axes=0,0,1,1/vlimits=0:3 ztzc_ratio
+
+set viewport lower
+go margins 0 1.0 1.5 1.5
+ppl axlint 1 5
+fill/nolabel/key=continuous/palette=no_red/level=(0.1,2,0.1)(inf)/vlimits=-40:0/axes=0,1,1,1 phy
+frame/file=ztzc.gif
+cancel viewport
+
+! Nitrate et nitracline
+let d_nit=nit[k=1:80@ddc]
+let nitracline=d_nit[k=1:50@loc:-0.02]
+fill/nolabel/key=continuous/palette=no_green_centered nit; frame/file=nit.gif
+plot nitracline; frame/file=nitracline.gif
+
+! Maximum de chlorophylle
+let phymax = phy[k=@max]
+let zero_at_phymax = phy - phymax
+let depth_of_phymax = zero_at_phymax[z=@loc:0]
+let depth_of_phymax = if (phymax gt 0.1) then zero_at_phymax[z=@loc:0] else 999
+set variable/bad=999 depth_of_phymax
+
+
+! Phytoplancton
+define axis/t="1-jan-2009":"28-jun-2009":24/units=hours thours
+fill/nolabel/key=continuous/palette=no_red/level=(0.1,2,0.1)(inf) phy[gt=thours]
+plot/over nitracline
+plot/over photic
+frame/file=phy.gif
+
+plot/nolabel/thick=2/vlimits=-120:0 depth_of_phymax
+plot/over/dash/color=black/thick=2 photic
+frame/file=phymax.gif
+
+! PON cycle
+!ppl axlint 1 5
+!plot/thick=2/vlimits=0:30/nolabel phy[k=1:80@din]
+!plot/over/thick=2/dash=(0.1,0.1,0.1,0.1)/color=black bac[k=1:80@din]
+!plot/over/thick=2/dash=(0.3,0.1,0.1,0.1)/color=black zoo[k=1:80@din]
+!plot/over/thick=2 det[k=1:80@din]
+!frame/file=pon_cycle.gif
+
+! DN cycle
+!plot/thick=2/nolabel totaln[k=1:80@din]
+!plot/over/thick=2/dash=(0.1,0.1,0.1,0.1)/color=black nit[k=1:80@din]
+!plot/over/thick=2/dash=(0.3,0.1,0.1,0.1)/color=black amm[k=1:80@din]
+!plot/over/thick=2/dash=(0.3,0.1,0.3,0.1)/color=black ldn[k=1:80@din]
+!frame/file=dn_cycle.gif
+
+
+! Autres compartiments
+fill/nolabel/pal=orange_red/key=continuous zoo;  frame/file=zoo.gif
+fill/nolabel/pal=red_purple/key=continuous bac;  frame/file=bac.gif
+fill/nolabel/key=continuous/pal=no_green_centered amm;  frame/file=amm.gif
+fill/nolabel/key=continuous/pal=reds_cmyk ldn;  frame/file=ldn.gif
+fill/nolabel/key=continuous/pal=orange_brown det;  frame/file=det.gif
+
+! ======================================================================
+! Hydrographie
+
+fill/nolabel/key=continuous/palette=rnb3/level=(-inf)(30,34,0.2)(inf) salt; contour/over/size=0.12/level=(30,34,0.2) salt
+frame/file=salt.gif
+
+fill/nolabel/key=continuous/palette=ocean_temp temp; contour/over/size=0.12 temp
+frame/file=temp.gif
+
+!plot sst;  frame/file=sst.gif
+!plot sss;  frame/file=sss.gif
+
+! Radiation
+plot/nolabel/thick=2 swr;  frame/file=swr.gif
+!plot int_swr; frame/file=int_swr.gif
+
+! Chaleur
+plot/nolabel/thick=2/title="Surface heat flux (W/m^2)" heat; frame/file=heat.gif
+!plot int_heat; frame/file=int_heat.gif
+
+! Vent
+let taus = (x_taus^2 + y_taus^2)^0.5
+set variable/unit="Pa"/title="surface wind stress" taus
+
+set viewport upper
+go margins 1.0 0.0 1.5 1.5
+ppl axlint 0 5
+plot/nolabel/thick=2/axes=0,0,1,1 taus
+
+set viewport lower
+go margins 0 1.0 1.5 1.5
+ppl axlint 1 5
+fill/nolabel/key=continuous/palette=rnb3/vlimits=-40:0/axes=0,1,1,1 num
+frame/file=taus_num.gif
+cancel viewport
+
+set viewport upper
+go margins 1.0 0.0 1.5 1.5
+ppl axlint 0 5
+plot/nolabel/thick=2/axes=0,0,1,1 taus
+
+set viewport lower
+go margins 0 1.0 1.5 1.5
+ppl axlint 1 5
+plot/nolabel/thick=2/axes=0,1,1,1/vlimits=-40:0 zsbl
+frame/file=taus_zsbl.gif
+cancel viewport
+
+
+set viewport upper
+go margins 1.0 0.0 1.5 1.5
+ppl axlint 0 5
+fill/nolabel/key=continuous/pal=warm_cmyk (u^2+v^2)^0.5
+
+set viewport lower
+go margins 0 1.0 1.5 1.5
+ppl axlint 1 5
+fill/nolabel/key=continuous/pal=rnb2/axes=0,1,1,1/vlimits=-100:-60/lev=(-5,-2,0.05)(inf) log(num)
+frame/file=u_num_bot.gif
+cancel viewport
+
+
+! Turbulence
+!fill/nolabel u;   frame/file=u.gif
+!fill/nolabel v;   frame/file=v.gif
+!fill/nolabel uu;  frame/file=uu.gif
+!fill/nolabel vv;  frame/file=vv.gif
+!fill/nolabel ww;  frame/file=ww.gif
+!fill/nolabel p;   frame/file=p.gif
+!fill/nolabel tke; frame/file=tke.gif
+!fill/nolabel eps; frame/file=eps.gif
+plot/nolabel/thick=2/vlimits=-60:0 zsbl; frame/file=zsbl.gif
+!plot/nolabel zbbl; frame/file=zbbl.gif
+
+! Non-local fluxes
+!fill/nolabel/key=continuous/palette=redscale gamh; frame/file=gamh.gif
+!fill/nolabel/key=continuous/palette=redscale gams; frame/file=gams.gif
+!fill/nolabel/key=continuous/palette=redscale gamu; frame/file=gamu.gif
+!fill/nolabel/key=continuous/palette=redscale gamv; frame/file=gamv.gif
+
+! ======================================================================
+! Résolution verticale
+set window/aspect=1.8
+plot/nolabel/l=1/symbol=27/thick=2/size=0.15 h; frame/file=h.gif
+
+exit