|
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 = i01.39cn_vs_i01.40cn |
|
29 |
|
30 #******************************************************* |
|
31 # user modification-(2) |
|
32 |
|
33 # model1 |
|
34 set MODEL1 = i01.39cn |
|
35 set DIR_M = /lustre/wolf-ddn/scratch/hof/cases1/i01.39cn_out/ |
|
36 set FILE1 = i01.39cn_out_ANN_climo.nc |
|
37 set FILE2 = i01.39cn_out_MONS_climo.nc |
|
38 set FILE3 = |
|
39 set FILE4 = $FILE1 |
|
40 set FILE5 = $FILE1 |
|
41 set FILE6 = $FILE1 |
|
42 set FILE7 = i01.39cn_Fire_C_2101-2125_monthly.nc |
|
43 set FILE8 = i01.39cn_ameriflux_2101-2125_monthly.nc |
|
44 set FILE9 = $FILE1 |
|
45 set FILE10 = $FILE2 |
|
46 set GRID = 1.9 |
|
47 set BGC = cn |
|
48 set ENERGY = old |
|
49 |
|
50 # in the "CLAMP metric processing" section: |
|
51 # only 00.initial.ncl and 99.final.ncl are required, |
|
52 # user can comment out any one or more of the other ncl scripts, |
|
53 # e.g. |
|
54 # #ncl $INPUT_TEXT $DIR_SCRIPTS/10.fire.ncl |
|
55 |
|
56 # model surface data |
|
57 set DIR_S = /lustre/wolf-ddn/scratch/hof/clamp_data/surface_model/ |
|
58 |
|
59 # observed data |
|
60 set DIR_O = /lustre/wolf-ddn/scratch/hof/clamp_data/observed/ |
|
61 |
|
62 # directory for scripts, templates and ncl files |
|
63 set DIR_SCRIPTS = /lustre/wolf-ddn/scratch/hof/clamp/all/ |
|
64 |
|
65 #******************************************************** |
|
66 |
|
67 # create model1 and model1_vs_model2 directory by copying templates |
|
68 if ($FILE3 != "") then |
|
69 set TEMPLATE1 = template_1-model |
|
70 set TEMPLATE2 = template_2-model |
|
71 else |
|
72 set TEMPLATE1 = template_1-model_noCO2 |
|
73 set TEMPLATE2 = template_2-model_noCO2 |
|
74 endif |
|
75 #cp -r $DIR_SCRIPTS/$TEMPLATE1 $MODEL1 |
|
76 #cp -r $DIR_SCRIPTS/$TEMPLATE2 $COMPARE |
|
77 rm -rf $MODEL1 |
|
78 rm -rf $COMPARE |
|
79 tar xvpf $DIR_SCRIPTS/$TEMPLATE1.tar |
|
80 mv $TEMPLATE1 $MODEL1 |
|
81 tar xvpf $DIR_SCRIPTS/$TEMPLATE2.tar |
|
82 mv $TEMPLATE2 $COMPARE |
|
83 |
|
84 # add quote, to be usesd in INPUT_TEXT |
|
85 set MODELQ = \"$MODEL1\" |
|
86 set DIRMQ = \"$DIR_M\" |
|
87 set F1 = \"$FILE1\" |
|
88 set F2 = \"$FILE2\" |
|
89 set F3 = \"$FILE3\" |
|
90 set F4 = \"$FILE4\" |
|
91 set F5 = \"$FILE5\" |
|
92 set F6 = \"$FILE6\" |
|
93 set F7 = \"$FILE7\" |
|
94 set F8 = \"$FILE8\" |
|
95 set F9 = \"$FILE9\" |
|
96 set F10 = \"$FILE10\" |
|
97 set GRIDQ = \"$GRID\" |
|
98 set BGCQ = \"$BGC\" |
|
99 set ENERGYQ = \"$ENERGY\" |
|
100 set DIRSQ = \"$DIR_S\" |
|
101 set DIROQ = \"$DIR_O\" |
|
102 set DIRCQ = \"$DIR_SCRIPTS\" |
|
103 |
|
104 set COMPAREQ = \"$COMPARE\" |
|
105 set MODELN = \"model1\" |
|
106 |
|
107 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 BGC=$BGCQ ENERGY=$ENERGYQ dirs=$DIRSQ diro=$DIROQ dirscript=$DIRCQ modeln=$MODELN compare=$COMPAREQ" |
|
108 |
|
109 echo "`date` - Running diagnostics for model ${MODELQ}" |
|
110 # CLAMP metric processing of model1 |
|
111 echo "`date` - Running 00.initial.ncl" |
|
112 ncl $INPUT_TEXT $DIR_SCRIPTS/00.initial.ncl |
|
113 echo "`date` - Running 01.npp.ncl" |
|
114 ncl $INPUT_TEXT $DIR_SCRIPTS/01.npp.ncl |
|
115 echo "`date` - Running 01.lai.ncl" |
|
116 ncl $INPUT_TEXT $DIR_SCRIPTS/02.lai.ncl |
|
117 |
|
118 if ($FILE3 != "") then |
|
119 echo "`date` - Running 03.co2.ncl" |
|
120 ncl $INPUT_TEXT $DIR_SCRIPTS/03.co2.ncl |
|
121 endif |
|
122 |
|
123 echo "`date` - Running 04.biomass.ncl" |
|
124 ncl $INPUT_TEXT $DIR_SCRIPTS/04.biomass.ncl |
|
125 echo "`date` - Running 06.fluxnet.ncl" |
|
126 ncl $INPUT_TEXT $DIR_SCRIPTS/06.fluxnet.ncl |
|
127 echo "`date` - Running 07.beta.ncl" |
|
128 ncl $INPUT_TEXT $DIR_SCRIPTS/07.beta.ncl |
|
129 echo "`date` - Running 08.turnover.ncl" |
|
130 ncl $INPUT_TEXT $DIR_SCRIPTS/08.turnover.ncl |
|
131 |
|
132 if ($BGC != "casa") then |
|
133 echo "`date` - Running 09.carbon_sink.ncl" |
|
134 ncl $INPUT_TEXT $DIR_SCRIPTS/09.carbon_sink.ncl |
|
135 else |
|
136 echo "`date` - Running 09x.carbon_sink.ncl" |
|
137 ncl $INPUT_TEXT $DIR_SCRIPTS/09x.carbon_sink.ncl |
|
138 endif |
|
139 |
|
140 if ($BGC != "casa") then |
|
141 echo "`date` - Running 10.fire.ncl" |
|
142 ncl $INPUT_TEXT $DIR_SCRIPTS/10.fire.ncl |
|
143 endif |
|
144 |
|
145 echo "`date` - Running 11.ameriflux.ncl" |
|
146 ncl $INPUT_TEXT $DIR_SCRIPTS/11.ameriflux.ncl |
|
147 echo "`date` - Running 99.final.ncl" |
|
148 ncl $INPUT_TEXT $DIR_SCRIPTS/99.final.ncl |
|
149 |
|
150 #******************************************************* |
|
151 # user modification-(3) |
|
152 |
|
153 # model2 |
|
154 set MODEL2 = i01.40cn |
|
155 set DIR_M = /lustre/wolf-ddn/scratch/hof/cases1/i01.40cn_out/ |
|
156 set FILE1 = i01.40cn_out_ANN_climo.nc |
|
157 set FILE2 = i01.40cn_out_MONS_climo.nc |
|
158 set FILE3 = |
|
159 set FILE4 = $FILE1 |
|
160 set FILE5 = $FILE1 |
|
161 set FILE6 = $FILE1 |
|
162 set FILE7 = i01.40cn_Fire_C_2076-2100_monthly.nc |
|
163 set FILE8 = i01.40cn_ameriflux_2076-2100_monthly.nc |
|
164 set FILE9 = $FILE1 |
|
165 set FILE10 = $FILE2 |
|
166 set GRID = 1.9 |
|
167 set BGC = cn |
|
168 set ENERGY = old |
|
169 #******************************************************* |
|
170 # create model2 directory by copying templates |
|
171 if ($FILE3 != "") then |
|
172 set TEMPLATE1 = template_1-model |
|
173 else |
|
174 set TEMPLATE1 = template_1-model_noCO2 |
|
175 endif |
|
176 #cp -r $DIR_SCRIPTS/$TEMPLATE1 $MODEL2 |
|
177 rm -rf $MODEL2 |
|
178 tar xvpf $DIR_SCRIPTS/$TEMPLATE1.tar |
|
179 mv $TEMPLATE1 $MODEL2 |
|
180 |
|
181 # add quote, to be usesd in INPUT_TEXT |
|
182 set MODELQ = \"$MODEL2\" |
|
183 set DIRMQ = \"$DIR_M\" |
|
184 set F1 = \"$FILE1\" |
|
185 set F2 = \"$FILE2\" |
|
186 set F3 = \"$FILE3\" |
|
187 set F4 = \"$FILE4\" |
|
188 set F5 = \"$FILE5\" |
|
189 set F6 = \"$FILE6\" |
|
190 set F7 = \"$FILE7\" |
|
191 set F8 = \"$FILE8\" |
|
192 set F9 = \"$FILE9\" |
|
193 set F10 = \"$FILE10\" |
|
194 set GRIDQ = \"$GRID\" |
|
195 set BGCQ = \"$BGC\" |
|
196 set ENERGYQ = \"$ENERGY\" |
|
197 |
|
198 set COMPAREQ = \"$COMPARE\" |
|
199 set MODELN = \"model2\" |
|
200 |
|
201 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 BGC=$BGCQ ENERGY=$ENERGYQ dirs=$DIRSQ diro=$DIROQ dirscript=$DIRCQ modeln=$MODELN compare=$COMPAREQ" |
|
202 |
|
203 echo "`date` - Running diagnostics for model ${MODELQ}" |
|
204 # CLAMP metric processing of model2 |
|
205 echo "`date` - Running 00.inital.ncl" |
|
206 ncl $INPUT_TEXT $DIR_SCRIPTS/00.initial.ncl |
|
207 echo "`date` - Running 01.npp.ncl" |
|
208 ncl $INPUT_TEXT $DIR_SCRIPTS/01.npp.ncl |
|
209 echo "`date` - Running 02.lai.ncl" |
|
210 ncl $INPUT_TEXT $DIR_SCRIPTS/02.lai.ncl |
|
211 |
|
212 if ($FILE3 != "") then |
|
213 echo "`date` - Running 03.co2.ncl" |
|
214 ncl $INPUT_TEXT $DIR_SCRIPTS/03.co2.ncl |
|
215 endif |
|
216 |
|
217 echo "`date` - Running 04.biomass.ncl" |
|
218 ncl $INPUT_TEXT $DIR_SCRIPTS/04.biomass.ncl |
|
219 echo "`date` - Running 06.fluxnet.ncl" |
|
220 ncl $INPUT_TEXT $DIR_SCRIPTS/06.fluxnet.ncl |
|
221 echo "`date` - Running 07.beta.ncl" |
|
222 ncl $INPUT_TEXT $DIR_SCRIPTS/07.beta.ncl |
|
223 echo "`date` - Running 08.turnover.ncl" |
|
224 ncl $INPUT_TEXT $DIR_SCRIPTS/08.turnover.ncl |
|
225 |
|
226 if ($BGC != "casa") then |
|
227 echo "`date` - Running 09.carbon_sink.ncl" |
|
228 ncl $INPUT_TEXT $DIR_SCRIPTS/09.carbon_sink.ncl |
|
229 else |
|
230 echo "`date` - Running 09x.carbon_sink.ncl" |
|
231 ncl $INPUT_TEXT $DIR_SCRIPTS/09x.carbon_sink.ncl |
|
232 endif |
|
233 |
|
234 if ($BGC != "casa") then |
|
235 echo "`date` - Running 10.fire.ncl" |
|
236 ncl $INPUT_TEXT $DIR_SCRIPTS/10.fire.ncl |
|
237 endif |
|
238 |
|
239 echo "`date` - Running 11.ameriflux.ncl" |
|
240 ncl $INPUT_TEXT $DIR_SCRIPTS/11.ameriflux.ncl |
|
241 echo "`date` - Running 99.final.ncl" |
|
242 ncl $INPUT_TEXT $DIR_SCRIPTS/99.final.ncl |
|
243 |
|
244 # create a tar file from the final output |
|
245 echo "`date` - tarring up final output into all.tar" |
|
246 tar cf - $MODEL1 $MODEL2 $COMPARE > all.tar |
|
247 |