|
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 = CN |
|
35 set DIR_M = /lustre/wolf-ddn/scratch/hof/clamp_data/model/ |
|
36 set FILE1 = i01.10cn_1948-2004_ANN_climo.nc |
|
37 set FILE2 = i01.10cn_1948-2004_MONS_climo.nc |
|
38 set FILE3 = f02.03cn_1901-1925_MONS_climo_atm.nc |
|
39 set FILE4 = i01.06cn_1980-2004_ANN_climo.nc |
|
40 set FILE5 = i01.10cn_1990-2004_ANN_climo.nc |
|
41 set FILE6 = i01.07cn_1990-2004_ANN_climo.nc |
|
42 set FILE7 = i01.10cn_Fire_C_1979-2004_monthly.nc |
|
43 set FILE8 = i01.10cn_ameriflux_1990-2004_monthly.nc |
|
44 set GRID = T42 |
|
45 set BGC = cn |
|
46 set ENERGY = new |
|
47 |
|
48 # in the "CLAMP metric processing" section: |
|
49 # only 00.initial.ncl and 99.final.ncl are required, |
|
50 # user can comment out any one or more of the other ncl scripts, |
|
51 # e.g. |
|
52 # #ncl $INPUT_TEXT $DIR_SCRIPTS/10.fire.ncl |
|
53 |
|
54 # model surface data |
|
55 set DIR_S = /lustre/wolf-ddn/scratch/hof/clamp_data/surface_model/ |
|
56 |
|
57 # observed data |
|
58 set DIR_O = /lustre/wolf-ddn/scratch/hof/clamp_data/observed/ |
|
59 |
|
60 # directory for scripts, templates and ncl files |
|
61 set DIR_SCRIPTS = /lustre/wolf-ddn/scratch/hof/clamp/all/ |
|
62 |
|
63 #******************************************************** |
|
64 |
|
65 # create model1 and model1_vs_model2 directory by copying templates |
|
66 if ($FILE3 != "") then |
|
67 set TEMPLATE1 = template_1-model |
|
68 set TEMPLATE2 = template_2-model |
|
69 else |
|
70 set TEMPLATE1 = template_1-model_noCO2 |
|
71 set TEMPLATE2 = template_2-model_noCO2 |
|
72 endif |
|
73 cp -r $DIR_SCRIPTS/$TEMPLATE1 $MODEL1 |
|
74 cp -r $DIR_SCRIPTS/$TEMPLATE2 $COMPARE |
|
75 |
|
76 # add quote, to be usesd in INPUT_TEXT |
|
77 set MODELQ = \"$MODEL1\" |
|
78 set DIRMQ = \"$DIR_M\" |
|
79 set F1 = \"$FILE1\" |
|
80 set F2 = \"$FILE2\" |
|
81 set F3 = \"$FILE3\" |
|
82 set F4 = \"$FILE4\" |
|
83 set F5 = \"$FILE5\" |
|
84 set F6 = \"$FILE6\" |
|
85 set F7 = \"$FILE7\" |
|
86 set F8 = \"$FILE8\" |
|
87 set GRIDQ = \"$GRID\" |
|
88 set BGCQ = \"$BGC\" |
|
89 set ENERGYQ = \"$ENERGY\" |
|
90 set DIRSQ = \"$DIR_S\" |
|
91 set DIROQ = \"$DIR_O\" |
|
92 set DIRCQ = \"$DIR_SCRIPTS\" |
|
93 |
|
94 set COMPAREQ = \"$COMPARE\" |
|
95 set MODELN = \"model1\" |
|
96 |
|
97 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 BGC=$BGCQ ENERGY=$ENERGYQ dirs=$DIRSQ diro=$DIROQ dirscript=$DIRCQ modeln=$MODELN compare=$COMPAREQ" |
|
98 |
|
99 # CLAMP metric processing of model1 |
|
100 ncl $INPUT_TEXT $DIR_SCRIPTS/00.initial.ncl |
|
101 ncl $INPUT_TEXT $DIR_SCRIPTS/01.npp.ncl |
|
102 ncl $INPUT_TEXT $DIR_SCRIPTS/02.lai.ncl |
|
103 |
|
104 if ($FILE3 != "") then |
|
105 ncl $INPUT_TEXT $DIR_SCRIPTS/03.co2.ncl |
|
106 endif |
|
107 |
|
108 ncl $INPUT_TEXT $DIR_SCRIPTS/04.biomass.ncl |
|
109 ncl $INPUT_TEXT $DIR_SCRIPTS/06.fluxnet.ncl |
|
110 ncl $INPUT_TEXT $DIR_SCRIPTS/07.beta.ncl |
|
111 ncl $INPUT_TEXT $DIR_SCRIPTS/08.turnover.ncl |
|
112 |
|
113 if ($BGC != "casa") then |
|
114 ncl $INPUT_TEXT $DIR_SCRIPTS/09.carbon_sink.ncl |
|
115 else |
|
116 ncl $INPUT_TEXT $DIR_SCRIPTS/09x.carbon_sink.ncl |
|
117 endif |
|
118 |
|
119 if ($BGC != "casa") then |
|
120 ncl $INPUT_TEXT $DIR_SCRIPTS/10.fire.ncl |
|
121 endif |
|
122 |
|
123 ncl $INPUT_TEXT $DIR_SCRIPTS/11.ameriflux.ncl |
|
124 ncl $INPUT_TEXT $DIR_SCRIPTS/99.final.ncl |
|
125 |
|
126 #******************************************************* |
|
127 # user modification-(3) |
|
128 |
|
129 # model2 |
|
130 set MODEL2 = CASA |
|
131 set DIR_M = /lustre/wolf-ddn/scratch/hof/clamp_data/model/ |
|
132 set FILE1 = i01.10casa_1948-2004_ANN_climo.nc |
|
133 set FILE2 = i01.10casa_1948-2004_MONS_climo.nc |
|
134 set FILE3 = f02.03casa_1876-1900_MONS_climo_atm.nc |
|
135 set FILE4 = i01.06casa_1980-2004_ANN_climo.nc |
|
136 set FILE5 = i01.10casa_1990-2004_ANN_climo.nc |
|
137 set FILE6 = i01.07casa_1990-2004_ANN_climo.nc |
|
138 set FILE7 = i01.10casa_Fire_C_1979-2004_monthly.nc |
|
139 set FILE8 = i01.10casa_ameriflux_1990-2004_monthly.nc |
|
140 set GRID = T42 |
|
141 set BGC = casa |
|
142 set ENERGY = new |
|
143 #******************************************************* |
|
144 # create model2 directory by copying templates |
|
145 if ($FILE3 != "") then |
|
146 set TEMPLATE1 = template_1-model |
|
147 else |
|
148 set TEMPLATE1 = template_1-model_noCO2 |
|
149 endif |
|
150 cp -r $DIR_SCRIPTS/$TEMPLATE1 $MODEL2 |
|
151 |
|
152 # add quote, to be usesd in INPUT_TEXT |
|
153 set MODELQ = \"$MODEL2\" |
|
154 set DIRMQ = \"$DIR_M\" |
|
155 set F1 = \"$FILE1\" |
|
156 set F2 = \"$FILE2\" |
|
157 set F3 = \"$FILE3\" |
|
158 set F4 = \"$FILE4\" |
|
159 set F5 = \"$FILE5\" |
|
160 set F6 = \"$FILE6\" |
|
161 set F7 = \"$FILE7\" |
|
162 set F8 = \"$FILE8\" |
|
163 set GRIDQ = \"$GRID\" |
|
164 set BGCQ = \"$BGC\" |
|
165 set ENERGYQ = \"$ENERGY\" |
|
166 |
|
167 set COMPAREQ = \"$COMPARE\" |
|
168 set MODELN = \"model2\" |
|
169 |
|
170 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 BGC=$BGCQ ENERGY=$ENERGYQ dirs=$DIRSQ diro=$DIROQ dirscript=$DIRCQ modeln=$MODELN compare=$COMPAREQ" |
|
171 |
|
172 # CLAMP metric processing of model2 |
|
173 ncl $INPUT_TEXT $DIR_SCRIPTS/00.initial.ncl |
|
174 ncl $INPUT_TEXT $DIR_SCRIPTS/01.npp.ncl |
|
175 ncl $INPUT_TEXT $DIR_SCRIPTS/02.lai.ncl |
|
176 |
|
177 if ($FILE3 != "") then |
|
178 ncl $INPUT_TEXT $DIR_SCRIPTS/03.co2.ncl |
|
179 endif |
|
180 |
|
181 ncl $INPUT_TEXT $DIR_SCRIPTS/04.biomass.ncl |
|
182 ncl $INPUT_TEXT $DIR_SCRIPTS/06.fluxnet.ncl |
|
183 ncl $INPUT_TEXT $DIR_SCRIPTS/07.beta.ncl |
|
184 ncl $INPUT_TEXT $DIR_SCRIPTS/08.turnover.ncl |
|
185 |
|
186 if ($BGC != "casa") then |
|
187 ncl $INPUT_TEXT $DIR_SCRIPTS/09.carbon_sink.ncl |
|
188 else |
|
189 ncl $INPUT_TEXT $DIR_SCRIPTS/09x.carbon_sink.ncl |
|
190 endif |
|
191 |
|
192 if ($BGC != "casa") then |
|
193 ncl $INPUT_TEXT $DIR_SCRIPTS/10.fire.ncl |
|
194 endif |
|
195 |
|
196 ncl $INPUT_TEXT $DIR_SCRIPTS/11.ameriflux.ncl |
|
197 ncl $INPUT_TEXT $DIR_SCRIPTS/99.final.ncl |
|
198 |
|
199 # create a tar file from the final output |
|
200 tar cf - $MODEL1 $MODEL2 $COMPARE > all.tar |
|
201 |