1 | #!/bin/ksh |
---|
2 | #set -vx |
---|
3 | |
---|
4 | #************************************************************** |
---|
5 | # Author: Nicolas Vuichard |
---|
6 | # Contact: Nicolas.Vuichard@lsce.jussieu.fr |
---|
7 | # $Date: |
---|
8 | # $Name$ |
---|
9 | # $Revision: 1 $ |
---|
10 | # History: |
---|
11 | # Modification: Nicolas Vuichard |
---|
12 | # |
---|
13 | #************************************************************** |
---|
14 | |
---|
15 | input=$@ |
---|
16 | paper=-p |
---|
17 | quality=cx1800 |
---|
18 | #cx800 cx1200 |
---|
19 | #frequency=1:12:6 |
---|
20 | echo $input |
---|
21 | #reference_file=RRRRRRR |
---|
22 | #validation directory=VVVVVVV |
---|
23 | debug= |
---|
24 | #-d |
---|
25 | |
---|
26 | # Ratio one day over time step |
---|
27 | one_day=86400 |
---|
28 | ratio=$(( one_day/1800 )) |
---|
29 | |
---|
30 | pwd |
---|
31 | |
---|
32 | cp ../ORCHIDEE_3dif.jnl . |
---|
33 | #---------------------- |
---|
34 | |
---|
35 | echo $input | ferretloop - ${debug} --ps2gif_option "-g $quality" -p -i "PFT distribution" -o pft ORCHIDEE_pft.jnl \ |
---|
36 | "\"PFT distribution\" RRRRRRR inputfile MAXVEGETFRAC" |
---|
37 | |
---|
38 | |
---|
39 | for freq in Day Month Year |
---|
40 | do |
---|
41 | if [ $freq = Day ] |
---|
42 | then |
---|
43 | ssbx=10 |
---|
44 | else |
---|
45 | ssbx=1 |
---|
46 | fi |
---|
47 | |
---|
48 | echo $input | ferretloop - ${debug} --ps2gif_option "-g $quality" -p -i "Sensible heat flux : $freq" -o fluxsens_$freq ORCHIDEE_3dif.jnl \ |
---|
49 | "\"Sensible heat flux : $freq\" VVVVVVV RRRRRRR inputfile H_f Qh $freq \"W/m2\" $ssbx" |
---|
50 | |
---|
51 | echo $input | ferretloop - ${debug} --ps2gif_option "-g $quality" -p -i "Latent heat flux : $freq" -o fluxlat_$freq ORCHIDEE_3dif.jnl \ |
---|
52 | "\"Latent heat flux : $freq\" VVVVVVV RRRRRRR inputfile LE_f Qle $freq \"W/m2\" $ssbx" |
---|
53 | |
---|
54 | echo $input | ferretloop - ${debug} --ps2gif_option "-g $quality" -p -i "Net radiation : $freq" -o netrad_$freq ORCHIDEE_3dif.jnl \ |
---|
55 | "\"Net radiation : $freq\" VVVVVVV RRRRRRR inputfile Rn_f swnet+lwnet $freq \"W/m2\" $ssbx" |
---|
56 | |
---|
57 | echo $input | ferretloop - ${debug} --ps2gif_option "-g $quality" -p -i "Ground heat flux : $freq" -o groundheat_$freq ORCHIDEE_3dif.jnl \ |
---|
58 | "\"Ground heat flux : $freq\" VVVVVVV RRRRRRR inputfile g_f Qg $freq \"W/m2\" $ssbx" |
---|
59 | |
---|
60 | |
---|
61 | echo $input | ferretloop - ${debug} --ps2gif_option "-g $quality" -p -i "Surface Temperature : $freq" -o tsoil_$freq ORCHIDEE_3dif.jnl \ |
---|
62 | "\"Surface Temperature : $freq\" VVVVVVV RRRRRRR inputfile Ts1_f (AvgSurfT-273.15) $freq \"degC\" $ssbx" |
---|
63 | |
---|
64 | echo $input | ferretloop - ${debug} --ps2gif_option "-g $quality" -p -i "Potential evapotranspiration : $freq" -o epot_$freq ORCHIDEE_3dif.jnl \ |
---|
65 | "\"Potential evapotranspiration : $freq\" VVVVVVV RRRRRRR inputfile (Epot_f*24) (PotEvap*${one_day}) $freq \"mm/d\" $ssbx" |
---|
66 | |
---|
67 | echo $input | ferretloop - ${debug} --ps2gif_option "-g $quality" -p -i "Soil moisture content 1 : $freq" -o humr1_$freq ORCHIDEE_3dif.jnl \ |
---|
68 | "\"Soil moisture content 1 : $freq\" VVVVVVV RRRRRRR inputfile SWC1_f SoilWet*100 $freq \"%\" $ssbx" |
---|
69 | |
---|
70 | ext_coef=0.5 |
---|
71 | echo $input | ferretloop - ${debug} --ps2gif_option "-g $quality" -p -i "LAI : $freq" -o lai_$freq ORCHIDEE_3dif.jnl \ |
---|
72 | "\"lai : $freq\" VVVVVVV RRRRRRR inputfile (Ts1_f/4) (0.-maxvegetfrac*ln(1.-vegetfrac/maxvegetfrac)/$ext_coef) $freq \"1\" $ssbx" |
---|
73 | |
---|
74 | echo $input | ferretloop - ${debug} --ps2gif_option "-g $quality" -p -i "LAI vs FAPAR : $freq" -o lai_fapar_$freq ORCHIDEE_3dif.jnl \ |
---|
75 | "\"lai vs FAPAR : $freq\" VVVVVVV RRRRRRR inputfile (-1.*ln(1.-FAPAR)/($ext_coef)) (0.-maxvegetfrac*ln(1.-vegetfrac/maxvegetfrac)/$ext_coef) $freq \"1\" $ssbx" |
---|
76 | |
---|
77 | |
---|
78 | echo $input | ferretloop - ${debug} --ps2gif_option "-g $quality" -p -i "Net Ecosystem Exchange : $freq" -o nee_$freq ORCHIDEE_3dif.jnl \ |
---|
79 | "\"Net Ecosystem Exchange : $freq\" VVVVVVV RRRRRRR inputfile (NEE*${ratio}) (NEE*${one_day}) $freq \"gC/m2/day\" $ssbx" |
---|
80 | |
---|
81 | echo $input | ferretloop - ${debug} --ps2gif_option "-g $quality" -p -i "Gross Primary Productivity : $freq" -o gpp_$freq ORCHIDEE_3dif.jnl \ |
---|
82 | "\"Gross Primary Productivity : $freq\" VVVVVVV RRRRRRR inputfile (GPP*${ratio}) (GPP*${one_day}) $freq \"gC/m2/day\" $ssbx" |
---|
83 | |
---|
84 | echo $input | ferretloop - ${debug} --ps2gif_option "-g $quality" -p -i "Total Ecosystem Respiration : $freq" -o reco_$freq ORCHIDEE_3dif.jnl \ |
---|
85 | "\"Total Ecosystem Respiration : $freq\" VVVVVVV RRRRRRR inputfile (Reco*${ratio}) ((maint_resp+hetero_resp+growth_resp)*${one_day}) $freq \"gC/m2/day\" $ssbx" |
---|
86 | |
---|
87 | done |
---|
88 | |
---|
89 | |
---|
90 | |
---|
91 | |
---|
92 | #Â Diurnal cycle |
---|
93 | |
---|
94 | let ANNEE=0 |
---|
95 | while [ $ANNEE -lt $DUREE ]; do |
---|
96 | let ANNEE2=$ANNEE+1 |
---|
97 | echo $input | ferretloop - ${debug} --ps2gif_option "-g $quality" -s ${ANNEE}:${ANNEE}:1 -p -i "Sensible heat flux : climatologie_${ANNEE2}" -o fluxsens_day ORCHIDEE_12day.jnl \ |
---|
98 | "\"Sensible heat flux\" VVVVVVV RRRRRRR inputfile H_f Qh sloop" |
---|
99 | |
---|
100 | echo $input | ferretloop - ${debug} --ps2gif_option "-g $quality" -s ${ANNEE}:${ANNEE}:1 -p -i "Latent heat flux : climatologie_${ANNEE2}" -o fluxlat_day ORCHIDEE_12day.jnl \ |
---|
101 | "\"Latent heat flux\" VVVVVVV RRRRRRR inputfile LE_f Qle sloop" |
---|
102 | |
---|
103 | echo $input | ferretloop - ${debug} --ps2gif_option "-g $quality" -s ${ANNEE}:${ANNEE}:1 -p -i "Net radiation : climatologie_${ANNEE2}" -o netrad_day ORCHIDEE_12day.jnl \ |
---|
104 | "\"Net radiation\" VVVVVVV RRRRRRR inputfile Rn_f swnet+lwnet sloop" |
---|
105 | |
---|
106 | echo $input | ferretloop - ${debug} --ps2gif_option "-g $quality" -s ${ANNEE}:${ANNEE}:1 -p -i "Ground heat flux : climatologie_${ANNEE2}" -o groundheat_day ORCHIDEE_12day.jnl \ |
---|
107 | "\"Ground heat flux\" VVVVVVV RRRRRRR inputfile g_f Qg sloop" |
---|
108 | |
---|
109 | |
---|
110 | echo $input | ferretloop - ${debug} --ps2gif_option "-g $quality" -s ${ANNEE}:${ANNEE}:1 -p -i "Surface Temperature : climatologie_${ANNEE2}" -o tsoil_day ORCHIDEE_12day.jnl \ |
---|
111 | "\"Surface Temperature\" VVVVVVV RRRRRRR inputfile Ts1_f (AvgSurfT-273.15) sloop" |
---|
112 | |
---|
113 | echo $input | ferretloop - ${debug} --ps2gif_option "-g $quality" -s ${ANNEE}:${ANNEE}:1 -p -i "gpp : climatologie_${ANNEE2}" -o gpp_day ORCHIDEE_12day.jnl \ |
---|
114 | "\"GPP\" VVVVVVV RRRRRRR inputfile (GPP*${ratio}) (GPP*${one_day}) sloop" |
---|
115 | |
---|
116 | echo $input | ferretloop - ${debug} --ps2gif_option "-g $quality" -s ${ANNEE}:${ANNEE}:1 -p -i "nee : climatologie_${ANNEE2}" -o nee_day ORCHIDEE_12day.jnl \ |
---|
117 | "\"nee\" VVVVVVV RRRRRRR inputfile (NEE*${ratio}) (NEE*${one_day}) sloop" |
---|
118 | |
---|
119 | echo $input | ferretloop - ${debug} --ps2gif_option "-g $quality" -s ${ANNEE}:${ANNEE}:1 -p -i "Ecosystem respiration : climatologie_${ANNEE2}" -o reco_day ORCHIDEE_12day.jnl \ |
---|
120 | "\"Ecosystem respiration\" VVVVVVV RRRRRRR inputfile (Reco*${ratio}) ((maint_resp+hetero_resp+growth_resp)*${one_day}) sloop" |
---|
121 | |
---|
122 | let ANNEE=$ANNEE+1 |
---|
123 | done |
---|