Commit 7818ae3b authored by Dany Dumont's avatar Dany Dumont

Nettoyage du cas test ghent et deplacement dans test_cases

parent ed3ed5c0
! 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=30]
let photic = light_fraction[z=@loc:0.01]
fill/nolabel/key=continuous/palette=rnb3/level=(10,100,5)(inf) 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) 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=-30: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/over/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=-30:0/axes=0,1,1,1 phy
frame/file=ztzc.gif
cancel viewport
! Nitrate et nitracline
let d_nit=nit[k=@ddc]
let nitracline=d_nit[k=@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
exit
plot/nolabel/thick=2 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=@din]
plot/over/thick=2/dash=(0.1,0.1,0.1,0.1)/color=black bac[k=@din]
plot/over/thick=2/dash=(0.3,0.1,0.1,0.1)/color=black zoo[k=@din]
!plot/over/thick=2 det[k=@din]
frame/file=pon_cycle.gif
! DN cycle
plot/thick=2/nolabel totaln[k=@din]
plot/over/thick=2/dash=(0.1,0.1,0.1,0.1)/color=black nit[k=@din]
plot/over/thick=2/dash=(0.3,0.1,0.1,0.1)/color=black amm[k=@din]
plot/over/thick=2/dash=(0.3,0.1,0.3,0.1)/color=black ldn[k=@din]
frame/file=dn_cycle.gif
! Autres compartiments
fill/nolabel/pal=orange_red/key=continuous zoo; frame/file=zoo.gif
fill/nolabel/key=continuous bac; frame/file=bac.gif
fill/nolabel/key=continuous/palette=no_green_centered amm; frame/file=amm.gif
fill/nolabel/key=continuous 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=-20: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
! NOAA/PMEL TMAP
! FERRET v6
! Linux(g77) 2.6.9-22. - 08/25/06
! 5-Mar-08 11:10
use out.nc
set region/l=1:24
set memory/size=200
set window/aspect=0.5/size=2000
ppl dfltfnt sr
ppl axlsze 0.18 0.18
! Physique
define viewport/ylim=0,1/xlim=0.00,0.17 v1
define viewport/ylim=0,1/xlim=0.17,0.33 v2
define viewport/ylim=0,1/xlim=0.33,0.50 v3
define viewport/ylim=0,1/xlim=0.50,0.66 v4
define viewport/ylim=0,1/xlim=0.66,0.83 v5
define viewport/ylim=0,1/xlim=0.83,1.00 v6
set viewport v1
go margins 0.5 0.5 0.2 0.2
!ppl axlint 1 2
!ppl axatic 0 11
fill/nolabel/nokey/axes=1,1,0,0/title=""/pal=rnb3 par
set viewport v2
go margins 0.5 0.5 0.2 0.2
!ppl axlint 1 8
!ppl color 2 70 70 70
fill/nolabel/nokey/axes=1,1,0,0/title=""/pal=no_green_centered nit
label/nouser -1, 0.5, 0, 90, 0.24, W m^-^2
set viewport v3
go margins 0.5 0.5 0.2 0.2
!ppl axlint 1 2
!ppl axatic 0 11
let phy = dia + fla
fill/nolabel/nokey/axes=1,1,0,0/title=""/level=(0.2,1,0.1)(inf)/pal=no_red phy
set viewport v4
go margins 0.5 0.5 0.2 0.2
!ppl axlint 1 4
let zoo = mcz + msz
fill/nolabel/nokey/axes=1,1,0,0/pal=orange_red zoo
!label/nouser -1, 0.5, 0, 90, 0.24, m^2 s^-^2
set viewport v5
go margins 0.5 0.5 0.2 0.2
!ppl axlint 0 4
fill/nolabel/nokey/axes=1,1,0,0/pal=orange_brown det
frame/file=snapshot.gif
cancel viewport
exit
#!/bin/bash
if [ $# -ne 1 ]
then
echo "Usage: $0 <ncfile>"
exit
fi
NCFILE=$1.nc
JNLFILE=anim_tmp.jnl
GIFFILE=$1.gif
if [ -f $JNLFILE ]
then
rm -f $JNLFILE
fi
cat >>$JNLFILE<< !
use $NCFILE
!set region/l=1:24
let nframe = `temp,return=lsize`
set memory/size=200
set window/aspect=0.5/size=2000
!ppl dfltfnt sr
!ppl axlsze 0.18 0.18
define viewport/ylim=0,1/xlim=0.0,0.2 v1
define viewport/ylim=0,1/xlim=0.2,0.4 v2
define viewport/ylim=0,1/xlim=0.4,0.6 v3
define viewport/ylim=0,1/xlim=0.6,0.8 v4
define viewport/ylim=0,1/xlim=0.8,1.0 v5
let phy = dia + fla
let zoo = mcz + msz
REPEAT/RANGE=1:100:5/NAME=m ( \\
set viewport v1; \\
go margins 0.5 0.5 0.2 0.2; \\
fill/nolabel/nokey/axes=1,1,0,0/title=""/pal=rnb3/lev=(0,400,10)(inf) par; \\
set viewport v2; \\
go margins 0.5 0.5 0.2 0.2; \\
fill/nolabel/nokey/axes=1,1,0,0/title=""/pal=no_green_centered/lev=(0,25,1)(inf) nit; \\
set viewport v3; \\
go margins 0.5 0.5 0.2 0.2; \\
fill/nolabel/nokey/axes=1,1,0,0/title=""/level=(0.2,1,0.1)(inf)/pal=no_red phy; \\
set viewport v4; \\
go margins 0.5 0.5 0.2 0.2; \\
fill/nolabel/nokey/axes=1,1,0,0/pal=orange_red zoo; \\
set viewport v5; \\
go margins 0.5 0.5 0.2 0.2; \\
fill/nolabel/nokey/axes=1,1,0,0/pal=orange_brown det; \\
frame/file=`m`.gif; \\
cancel viewport )
exit
!
echo 'YOUHOU'
ferret -gif -script $JNLFILE
set gifsicle_options = "-d 30 -l"
gifsicle $gifsicle_options `ls cn.gif.~?~| sort -n` > movie.gif
rm -f $JNLFILE
exit
/home/dumoda01/git/gotm_ismer/src/gotm_prod_IFORT
\ No newline at end of file
This diff is collapsed.
This diff is collapsed.
1.00000
0.80000
0.04000
0.04000
0.00600
0.00600
1.00000
0.80000
1.00000
0.80000
0.05000
0.05000
0.20000
0.05000
-0.05000
-0.20000
1.00000
1.00000
1.00000
0.62500
0.30000
0.30000
0.20000
0.10000
0.70000
0.55000
0.40000
0.05000
0.10000
0.50000
0.30000
0.05000
0.15000
0.24000
0.50000
0.03000
0.00000
0.02000
0.00000
-5.00000
0.03000
3.61000
14.58000
3.26500
24.92300
********
********
1.00000
1.10000
0.40000
0.04000
0.50000
This diff is collapsed.
......@@ -77,14 +77,14 @@
! only for output
!-------------------------------------------------------------------------------
&airsea
calc_fluxes= .true.
calc_fluxes= .false.
meteo_file= 'meteo.dat'
wet_mode= 1
heat_method= 0
const_swr= 200.0
const_heat= 0.0
heatflux_file= 'heatflux.dat'
momentum_method= 1
momentum_method= 0
const_tx= 0.1
const_ty= 0.0
momentumflux_file='momentumflux.dat'
......
......@@ -46,7 +46,7 @@
!-------------------------------------------------------------------------------
&bio_nml
bio_calc= .true.
bio_model= 7
bio_model= 4
bio_eulerian= .true.
cnpar= 1.0
w_adv_discr= 4
......
......@@ -88,7 +88,7 @@
out_fmt= 2
out_dir= "."
out_fn= out
nsave= 10
nsave= 30
diagnostics= .false.
mld_method= 2
diff_k= 1.e-5
......
......@@ -55,7 +55,7 @@
z_s2= 35.
s_2= 32.
s_obs_NN= 2.56e-4
s_prof_file= 'franklin_sprof_ctd.dat'
s_prof_file= 'sprof.dat'
SRelaxTauM= 1209600.
SRelaxTauB= 1209600.
SRelaxTauS= 1209600.
......@@ -110,7 +110,7 @@
n_1= 3.0
z_n2= 70.
n_2= 15.0
n_prof_file= 'nprof_ctd.dat'
n_prof_file= 'nprof.dat'
NRelaxTauM= 1.e15
NRelaxTauB= 1.e15
NRelaxTauS= 1.e15
......@@ -224,7 +224,7 @@
z_t2= 35.
t_2= 0.0
t_obs_NN= 2.56e-4
t_prof_file= 'franklin_tprof_ctd.dat'
t_prof_file= 'tprof.dat'
TRelaxTauM= 1209600.
TRelaxTauB= 1209600.
TRelaxTauS= 1209600.
......@@ -339,18 +339,18 @@
!
!-------------------------------------------------------------------------------
&ext_pressure
ext_press_method=0
ext_press_method=1
ext_press_mode= 0
ext_press_file= 'pressure.dat'
PressConstU= 0.0
PressConstU= 1.e-6
PressConstV= 0.0
PressHeight= 0.0
PeriodM= 0.0
AmpMu= 0.0
PeriodM= 44712.0
AmpMu= 1.e-5
AmpMv= 0.0
PhaseMu= 0.0
PhaseMv= 0.0
PeriodS= 0.0
PeriodS= 43200.0
AmpSu= 0.0
AmpSv= 0.0
PhaseSu= 0.0
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment