|
1 #!/bin/csh -f |
|
2 #------------------------------------------------------------------- |
|
3 # note: user modifies ONLY the "user modification" section |
|
4 # |
|
5 # COMPARE: model1 vs model2 |
|
6 # MODELn : model name |
|
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: |
|
15 # set FILE3 = |
|
16 # FILE7 : timeseries file generated from .... |
|
17 # FILE8 : timeseries file generated from .... |
|
18 # GRID : T31, T42, or 1.9 |
|
19 # BGC : cn or casa |
|
20 # ENERGY : new or old (fields in model data) |
|
21 #------------------------------------------------------------------- |
|
22 |
|
23 #******************************************************* |
|
24 # user modification-(1) |
|
25 |
|
26 # directory name of model comparison |
|
27 #et COMPARE = b30.061n_vs_b30.061m |
|
28 set COMPARE = CN_vs_CASA |
|
29 |
|
30 #******************************************************* |
|
31 # user modification-(2) |
|
32 |
|
33 # model1 |
|
34 set MODEL1 = CASA |
|
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 |
|
48 set GRID = T42 |
|
49 set BGC = casa |
|
50 set ENERGY = new |
|
51 |
|
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, |
|
55 # e.g. |
|
56 # #ncl $INPUT_TEXT $DIR_SCRIPTS/10.fire.ncl |
|
57 |
|
58 # model surface data |
|
59 set DIR_S = /lustre/wolf-ddn/scratch/hof/clamp_data/surface_model/ |
|
60 |
|
61 # observed data |
|
62 set DIR_O = /lustre/wolf-ddn/scratch/hof/clamp_data/observed/ |
|
63 |
|
64 # directory for scripts, templates and ncl files |
|
65 set DIR_SCRIPTS = /lustre/wolf-ddn/scratch/hof/clamp/all/ |
|
66 |
|
67 #******************************************************** |
|
68 |
|
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 |
|
73 else |
|
74 set TEMPLATE1 = template_1-model_noCO2 |
|
75 set TEMPLATE2 = template_2-model_noCO2 |
|
76 endif |
|
77 #cp -r $DIR_SCRIPTS/$TEMPLATE1 $MODEL1 |
|
78 #cp -r $DIR_SCRIPTS/$TEMPLATE2 $COMPARE |
|
79 rm -rf $MODEL1 |
|
80 rm -rf $COMPARE |
|
81 tar xvpf $DIR_SCRIPTS/$TEMPLATE1.tar |
|
82 mv $TEMPLATE1 $MODEL1 |
|
83 tar xvpf $DIR_SCRIPTS/$TEMPLATE2.tar |
|
84 mv $TEMPLATE2 $COMPARE |
|
85 |
|
86 # add quote, to be usesd in INPUT_TEXT |
|
87 set MODELQ = \"$MODEL1\" |
|
88 set DIRMQ = \"$DIR_M\" |
|
89 set F1 = \"$FILE1\" |
|
90 set F2 = \"$FILE2\" |
|
91 set F3 = \"$FILE3\" |
|
92 set F4 = \"$FILE4\" |
|
93 set F5 = \"$FILE5\" |
|
94 set F6 = \"$FILE6\" |
|
95 set F7 = \"$FILE7\" |
|
96 set F8 = \"$FILE8\" |
|
97 set F9 = \"$FILE9\" |
|
98 set F10 = \"$FILE10\" |
|
99 set F11 = \"$FILE11\" |
|
100 set GRIDQ = \"$GRID\" |
|
101 set BGCQ = \"$BGC\" |
|
102 set ENERGYQ = \"$ENERGY\" |
|
103 set DIRSQ = \"$DIR_S\" |
|
104 set DIROQ = \"$DIR_O\" |
|
105 set DIRCQ = \"$DIR_SCRIPTS\" |
|
106 |
|
107 set COMPAREQ = \"$COMPARE\" |
|
108 set MODELN = \"model1\" |
|
109 |
|
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" |
|
111 |
|
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 |
|
120 |
|
121 if ($FILE3 != "") then |
|
122 echo "`date` - Running 03.co2.ncl" |
|
123 ncl $INPUT_TEXT $DIR_SCRIPTS/03.co2.ncl |
|
124 endif |
|
125 |
|
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 |
|
134 |
|
135 if ($BGC != "casa") then |
|
136 echo "`date` - Running 09.carbon_sink.ncl" |
|
137 ncl $INPUT_TEXT $DIR_SCRIPTS/09.carbon_sink.ncl |
|
138 else |
|
139 echo "`date` - Running 09x.carbon_sink.ncl" |
|
140 ncl $INPUT_TEXT $DIR_SCRIPTS/09x.carbon_sink.ncl |
|
141 endif |
|
142 |
|
143 if ($BGC != "casa") then |
|
144 echo "`date` - Running 10.fire.ncl" |
|
145 ncl $INPUT_TEXT $DIR_SCRIPTS/10.fire.ncl |
|
146 endif |
|
147 |
|
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 |
|
152 |
|
153 #******************************************************* |
|
154 # user modification-(3) |
|
155 |
|
156 # model2 |
|
157 set MODEL2 = CN |
|
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 |
|
171 set GRID = T42 |
|
172 set BGC = cn |
|
173 set ENERGY = new |
|
174 |
|
175 #******************************************************* |
|
176 # create model2 directory by copying templates |
|
177 if ($FILE3 != "") then |
|
178 set TEMPLATE1 = template_1-model |
|
179 else |
|
180 set TEMPLATE1 = template_1-model_noCO2 |
|
181 endif |
|
182 #cp -r $DIR_SCRIPTS/$TEMPLATE1 $MODEL2 |
|
183 rm -rf $MODEL2 |
|
184 tar xvpf $DIR_SCRIPTS/$TEMPLATE1.tar |
|
185 mv $TEMPLATE1 $MODEL2 |
|
186 |
|
187 # add quote, to be usesd in INPUT_TEXT |
|
188 set MODELQ = \"$MODEL2\" |
|
189 set DIRMQ = \"$DIR_M\" |
|
190 set F1 = \"$FILE1\" |
|
191 set F2 = \"$FILE2\" |
|
192 set F3 = \"$FILE3\" |
|
193 set F4 = \"$FILE4\" |
|
194 set F5 = \"$FILE5\" |
|
195 set F6 = \"$FILE6\" |
|
196 set F7 = \"$FILE7\" |
|
197 set F8 = \"$FILE8\" |
|
198 set F9 = \"$FILE9\" |
|
199 set F10 = \"$FILE10\" |
|
200 set F11 = \"$FILE11\" |
|
201 set GRIDQ = \"$GRID\" |
|
202 set BGCQ = \"$BGC\" |
|
203 set ENERGYQ = \"$ENERGY\" |
|
204 |
|
205 set COMPAREQ = \"$COMPARE\" |
|
206 set MODELN = \"model2\" |
|
207 |
|
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" |
|
209 |
|
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 |
|
218 |
|
219 if ($FILE3 != "") then |
|
220 echo "`date` - Running 03.co2.ncl" |
|
221 ncl $INPUT_TEXT $DIR_SCRIPTS/03.co2.ncl |
|
222 endif |
|
223 |
|
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 |
|
232 |
|
233 if ($BGC != "casa") then |
|
234 echo "`date` - Running 09.carbon_sink.ncl" |
|
235 ncl $INPUT_TEXT $DIR_SCRIPTS/09.carbon_sink.ncl |
|
236 else |
|
237 echo "`date` - Running 09x.carbon_sink.ncl" |
|
238 ncl $INPUT_TEXT $DIR_SCRIPTS/09x.carbon_sink.ncl |
|
239 endif |
|
240 |
|
241 if ($BGC != "casa") then |
|
242 echo "`date` - Running 10.fire.ncl" |
|
243 ncl $INPUT_TEXT $DIR_SCRIPTS/10.fire.ncl |
|
244 endif |
|
245 |
|
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 |
|
250 |
|
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 |
|
254 |