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