Modifications to scoring and graphics production for the final version of code for the C-LAMP paper in GCB.
2 #-------------------------------------------------------------------
3 # note: user modifies ONLY the "user modification" section
5 # COMPARE: model1 vs model2
7 # DIR_M : directory of model data
8 # DIR_O : directory of observed data
9 # DIR_S : directory of model surface data
10 # DIR_SCRIPTS : directory of run scripts
11 # FILE1 : time_mean climatology from CLM diagnostic package
12 # FILE2 : 12-monthly climatology from CLM diagnostic package
13 # FILE3 : 12-monthly climatology from ATM diagnostic package
14 # leave it blank, if no ATM file:
16 # FILE7 : timeseries file generated from ....
17 # FILE8 : timeseries file generated from ....
18 # GRID : T31, T42, or 1.9
20 # ENERGY : new or old (fields in model data)
21 #-------------------------------------------------------------------
23 #*******************************************************
24 # user modification-(1)
26 # directory name of model comparison
27 #et COMPARE = b30.061n_vs_b30.061m
28 set COMPARE = CN_vs_CASA
30 #*******************************************************
31 # user modification-(2)
35 set DIR_M = /lustre/wolf-ddn/scratch/hof/clamp_data/model/
36 #set FILE1 = i01.10casa_1948-2004_ANN_climo.nc
37 set FILE1 = i01.10casa_1975-2000_ANN_climo.nc
38 set FILE2 = i01.10casa_1948-2004_MONS_climo.nc
39 set FILE3 = f02.03casa_1876-1900_MONS_climo_atm.nc
40 set FILE4 = i01.06casa_1980-2004_ANN_climo.nc
41 set FILE5 = i01.10casa_1990-2004_ANN_climo.nc
42 set FILE6 = i01.07casa_1990-2004_ANN_climo.nc
43 set FILE7 = i01.10casa_Fire_C_1979-2004_monthly.nc
44 set FILE8 = i01.10casa_ameriflux_1990-2004_monthly.nc
45 set FILE9 = i01.10casa_2000-2004_ANN_climo.nc
46 set FILE10 = i01.10casa_2000-2004_MONS_climo.nc
47 set FILE11 = i01.10casa_2000_ANN_climo.nc
52 # in the "CLAMP metric processing" section:
53 # only 00.initial.ncl and 99.final.ncl are required,
54 # user can comment out any one or more of the other ncl scripts,
56 # #ncl $INPUT_TEXT $DIR_SCRIPTS/10.fire.ncl
59 set DIR_S = /lustre/wolf-ddn/scratch/hof/clamp_data/surface_model/
62 set DIR_O = /lustre/wolf-ddn/scratch/hof/clamp_data/observed/
64 # directory for scripts, templates and ncl files
65 set DIR_SCRIPTS = /lustre/wolf-ddn/scratch/hof/clamp/all/
67 #********************************************************
69 # create model1 and model1_vs_model2 directory by copying templates
70 if ($FILE3 != "") then
71 set TEMPLATE1 = template_1-model
72 set TEMPLATE2 = template_2-model
74 set TEMPLATE1 = template_1-model_noCO2
75 set TEMPLATE2 = template_2-model_noCO2
77 #cp -r $DIR_SCRIPTS/$TEMPLATE1 $MODEL1
78 #cp -r $DIR_SCRIPTS/$TEMPLATE2 $COMPARE
81 tar xvpf $DIR_SCRIPTS/$TEMPLATE1.tar
83 tar xvpf $DIR_SCRIPTS/$TEMPLATE2.tar
84 mv $TEMPLATE2 $COMPARE
86 # add quote, to be usesd in INPUT_TEXT
87 set MODELQ = \"$MODEL1\"
88 set DIRMQ = \"$DIR_M\"
100 set GRIDQ = \"$GRID\"
102 set ENERGYQ = \"$ENERGY\"
103 set DIRSQ = \"$DIR_S\"
104 set DIROQ = \"$DIR_O\"
105 set DIRCQ = \"$DIR_SCRIPTS\"
107 set COMPAREQ = \"$COMPARE\"
108 set MODELN = \"model1\"
110 set INPUT_TEXT = "model_name=$MODELQ model_grid=$GRIDQ dirm=$DIRMQ film1=$F1 film2=$F2 film3=$F3 film4=$F4 film5=$F5 film6=$F6 film7=$F7 film8=$F8 film9=$F9 film10=$F10 film11=$F11 BGC=$BGCQ ENERGY=$ENERGYQ dirs=$DIRSQ diro=$DIROQ dirscript=$DIRCQ modeln=$MODELN compare=$COMPAREQ"
112 echo "`date` - Running diagnostics for model ${MODELQ}"
113 # CLAMP metric processing of model1
114 echo "`date` - Running 00.initial.ncl"
115 ncl $INPUT_TEXT $DIR_SCRIPTS/00.initial.ncl
116 echo "`date` - Running 01.npp.ncl"
117 ncl $INPUT_TEXT $DIR_SCRIPTS/01.npp.ncl
118 echo "`date` - Running 01.lai.ncl"
119 ncl $INPUT_TEXT $DIR_SCRIPTS/02.lai.ncl
121 if ($FILE3 != "") then
122 echo "`date` - Running 03.co2.ncl"
123 ncl $INPUT_TEXT $DIR_SCRIPTS/03.co2.ncl
126 echo "`date` - Running 04.biomass.ncl"
127 ncl $INPUT_TEXT $DIR_SCRIPTS/04.biomass.ncl
128 echo "`date` - Running 06.fluxnet.ncl"
129 ncl $INPUT_TEXT $DIR_SCRIPTS/06.fluxnet.ncl
130 echo "`date` - Running 07.beta.ncl"
131 ncl $INPUT_TEXT $DIR_SCRIPTS/07.beta.ncl
132 echo "`date` - Running 08.turnover.ncl"
133 ncl $INPUT_TEXT $DIR_SCRIPTS/08.turnover.ncl
135 if ($BGC != "casa") then
136 echo "`date` - Running 09.carbon_sink.ncl"
137 ncl $INPUT_TEXT $DIR_SCRIPTS/09.carbon_sink.ncl
139 echo "`date` - Running 09x.carbon_sink.ncl"
140 ncl $INPUT_TEXT $DIR_SCRIPTS/09x.carbon_sink.ncl
143 if ($BGC != "casa") then
144 echo "`date` - Running 10.fire.ncl"
145 ncl $INPUT_TEXT $DIR_SCRIPTS/10.fire.ncl
148 echo "`date` - Running 11.ameriflux.ncl"
149 ncl $INPUT_TEXT $DIR_SCRIPTS/11.ameriflux.ncl
150 echo "`date` - Running 99.final.ncl"
151 ncl $INPUT_TEXT $DIR_SCRIPTS/99.final.ncl
153 #*******************************************************
154 # user modification-(3)
158 set DIR_M = /lustre/wolf-ddn/scratch/hof/clamp_data/model/
159 #set FILE1 = i01.10cn_1948-2004_ANN_climo.nc
160 set FILE1 = i01.10cn_1975-2000_ANN_climo.nc
161 set FILE2 = i01.10cn_1948-2004_MONS_climo.nc
162 set FILE3 = f02.03cn_1901-1925_MONS_climo_atm.nc
163 set FILE4 = i01.06cn_1980-2004_ANN_climo.nc
164 set FILE5 = i01.10cn_1990-2004_ANN_climo.nc
165 set FILE6 = i01.07cn_1990-2004_ANN_climo.nc
166 set FILE7 = i01.10cn_Fire_C_1979-2004_monthly.nc
167 set FILE8 = i01.10cn_ameriflux_1990-2004_monthly.nc
168 set FILE9 = i01.10cn_2000-2004_ANN_climo.nc
169 set FILE10 = i01.10cn_2000-2004_MONS_climo.nc
170 set FILE11 = i01.10cn_2000_ANN_climo.nc
175 #*******************************************************
176 # create model2 directory by copying templates
177 if ($FILE3 != "") then
178 set TEMPLATE1 = template_1-model
180 set TEMPLATE1 = template_1-model_noCO2
182 #cp -r $DIR_SCRIPTS/$TEMPLATE1 $MODEL2
184 tar xvpf $DIR_SCRIPTS/$TEMPLATE1.tar
185 mv $TEMPLATE1 $MODEL2
187 # add quote, to be usesd in INPUT_TEXT
188 set MODELQ = \"$MODEL2\"
189 set DIRMQ = \"$DIR_M\"
199 set F10 = \"$FILE10\"
200 set F11 = \"$FILE11\"
201 set GRIDQ = \"$GRID\"
203 set ENERGYQ = \"$ENERGY\"
205 set COMPAREQ = \"$COMPARE\"
206 set MODELN = \"model2\"
208 set INPUT_TEXT = "model_name=$MODELQ model_grid=$GRIDQ dirm=$DIRMQ film1=$F1 film2=$F2 film3=$F3 film4=$F4 film5=$F5 film6=$F6 film7=$F7 film8=$F8 film9=$F9 film10=$F10 film11=$F11 BGC=$BGCQ ENERGY=$ENERGYQ dirs=$DIRSQ diro=$DIROQ dirscript=$DIRCQ modeln=$MODELN compare=$COMPAREQ"
210 echo "`date` - Running diagnostics for model ${MODELQ}"
211 # CLAMP metric processing of model2
212 echo "`date` - Running 00.inital.ncl"
213 ncl $INPUT_TEXT $DIR_SCRIPTS/00.initial.ncl
214 echo "`date` - Running 01.npp.ncl"
215 ncl $INPUT_TEXT $DIR_SCRIPTS/01.npp.ncl
216 echo "`date` - Running 02.lai.ncl"
217 ncl $INPUT_TEXT $DIR_SCRIPTS/02.lai.ncl
219 if ($FILE3 != "") then
220 echo "`date` - Running 03.co2.ncl"
221 ncl $INPUT_TEXT $DIR_SCRIPTS/03.co2.ncl
224 echo "`date` - Running 04.biomass.ncl"
225 ncl $INPUT_TEXT $DIR_SCRIPTS/04.biomass.ncl
226 echo "`date` - Running 06.fluxnet.ncl"
227 ncl $INPUT_TEXT $DIR_SCRIPTS/06.fluxnet.ncl
228 echo "`date` - Running 07.beta.ncl"
229 ncl $INPUT_TEXT $DIR_SCRIPTS/07.beta.ncl
230 echo "`date` - Running 08.turnover.ncl"
231 ncl $INPUT_TEXT $DIR_SCRIPTS/08.turnover.ncl
233 if ($BGC != "casa") then
234 echo "`date` - Running 09.carbon_sink.ncl"
235 ncl $INPUT_TEXT $DIR_SCRIPTS/09.carbon_sink.ncl
237 echo "`date` - Running 09x.carbon_sink.ncl"
238 ncl $INPUT_TEXT $DIR_SCRIPTS/09x.carbon_sink.ncl
241 if ($BGC != "casa") then
242 echo "`date` - Running 10.fire.ncl"
243 ncl $INPUT_TEXT $DIR_SCRIPTS/10.fire.ncl
246 echo "`date` - Running 11.ameriflux.ncl"
247 ncl $INPUT_TEXT $DIR_SCRIPTS/11.ameriflux.ncl
248 echo "`date` - Running 99.final.ncl"
249 ncl $INPUT_TEXT $DIR_SCRIPTS/99.final.ncl
251 # create a tar file from the final output
252 echo "`date` - tarring up final output into all.tar"
253 tar cf - $MODEL1 $MODEL2 $COMPARE > all.tar