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