1.1 --- /dev/null Thu Jan 01 00:00:00 1970 +0000
1.2 +++ b/all/run_2-model_FMH.csh Thu Mar 26 14:02:21 2009 -0400
1.3 @@ -0,0 +1,254 @@
1.4 +#!/bin/csh -f
1.5 +#-------------------------------------------------------------------
1.6 +# note: user modifies ONLY the "user modification" section
1.7 +#
1.8 +# COMPARE: model1 vs model2
1.9 +# MODELn : model name
1.10 +# DIR_M : directory of model data
1.11 +# DIR_O : directory of observed data
1.12 +# DIR_S : directory of model surface data
1.13 +# DIR_SCRIPTS : directory of run scripts
1.14 +# FILE1 : time_mean climatology from CLM diagnostic package
1.15 +# FILE2 : 12-monthly climatology from CLM diagnostic package
1.16 +# FILE3 : 12-monthly climatology from ATM diagnostic package
1.17 +# leave it blank, if no ATM file:
1.18 +# set FILE3 =
1.19 +# FILE7 : timeseries file generated from ....
1.20 +# FILE8 : timeseries file generated from ....
1.21 +# GRID : T31, T42, or 1.9
1.22 +# BGC : cn or casa
1.23 +# ENERGY : new or old (fields in model data)
1.24 +#-------------------------------------------------------------------
1.25 +
1.26 +#*******************************************************
1.27 +# user modification-(1)
1.28 +
1.29 +# directory name of model comparison
1.30 +#et COMPARE = b30.061n_vs_b30.061m
1.31 +set COMPARE = CN_vs_CASA
1.32 +
1.33 +#*******************************************************
1.34 +# user modification-(2)
1.35 +
1.36 +# model1
1.37 +set MODEL1 = CASA
1.38 +set DIR_M = /lustre/wolf-ddn/scratch/hof/clamp_data/model/
1.39 +#set FILE1 = i01.10casa_1948-2004_ANN_climo.nc
1.40 +set FILE1 = i01.10casa_1975-2000_ANN_climo.nc
1.41 +set FILE2 = i01.10casa_1948-2004_MONS_climo.nc
1.42 +set FILE3 = f02.03casa_1876-1900_MONS_climo_atm.nc
1.43 +set FILE4 = i01.06casa_1980-2004_ANN_climo.nc
1.44 +set FILE5 = i01.10casa_1990-2004_ANN_climo.nc
1.45 +set FILE6 = i01.07casa_1990-2004_ANN_climo.nc
1.46 +set FILE7 = i01.10casa_Fire_C_1979-2004_monthly.nc
1.47 +set FILE8 = i01.10casa_ameriflux_1990-2004_monthly.nc
1.48 +set FILE9 = i01.10casa_2000-2004_ANN_climo.nc
1.49 +set FILE10 = i01.10casa_2000-2004_MONS_climo.nc
1.50 +set FILE11 = i01.10casa_2000_ANN_climo.nc
1.51 +set GRID = T42
1.52 +set BGC = casa
1.53 +set ENERGY = new
1.54 +
1.55 +# in the "CLAMP metric processing" section:
1.56 +# only 00.initial.ncl and 99.final.ncl are required,
1.57 +# user can comment out any one or more of the other ncl scripts,
1.58 +# e.g.
1.59 +# #ncl $INPUT_TEXT $DIR_SCRIPTS/10.fire.ncl
1.60 +
1.61 +# model surface data
1.62 +set DIR_S = /lustre/wolf-ddn/scratch/hof/clamp_data/surface_model/
1.63 +
1.64 +# observed data
1.65 +set DIR_O = /lustre/wolf-ddn/scratch/hof/clamp_data/observed/
1.66 +
1.67 +# directory for scripts, templates and ncl files
1.68 +set DIR_SCRIPTS = /lustre/wolf-ddn/scratch/hof/clamp/all/
1.69 +
1.70 +#********************************************************
1.71 +
1.72 +# create model1 and model1_vs_model2 directory by copying templates
1.73 +if ($FILE3 != "") then
1.74 + set TEMPLATE1 = template_1-model
1.75 + set TEMPLATE2 = template_2-model
1.76 +else
1.77 + set TEMPLATE1 = template_1-model_noCO2
1.78 + set TEMPLATE2 = template_2-model_noCO2
1.79 +endif
1.80 +#cp -r $DIR_SCRIPTS/$TEMPLATE1 $MODEL1
1.81 +#cp -r $DIR_SCRIPTS/$TEMPLATE2 $COMPARE
1.82 +rm -rf $MODEL1
1.83 +rm -rf $COMPARE
1.84 +tar xvpf $DIR_SCRIPTS/$TEMPLATE1.tar
1.85 +mv $TEMPLATE1 $MODEL1
1.86 +tar xvpf $DIR_SCRIPTS/$TEMPLATE2.tar
1.87 +mv $TEMPLATE2 $COMPARE
1.88 +
1.89 +# add quote, to be usesd in INPUT_TEXT
1.90 +set MODELQ = \"$MODEL1\"
1.91 +set DIRMQ = \"$DIR_M\"
1.92 +set F1 = \"$FILE1\"
1.93 +set F2 = \"$FILE2\"
1.94 +set F3 = \"$FILE3\"
1.95 +set F4 = \"$FILE4\"
1.96 +set F5 = \"$FILE5\"
1.97 +set F6 = \"$FILE6\"
1.98 +set F7 = \"$FILE7\"
1.99 +set F8 = \"$FILE8\"
1.100 +set F9 = \"$FILE9\"
1.101 +set F10 = \"$FILE10\"
1.102 +set F11 = \"$FILE11\"
1.103 +set GRIDQ = \"$GRID\"
1.104 +set BGCQ = \"$BGC\"
1.105 +set ENERGYQ = \"$ENERGY\"
1.106 +set DIRSQ = \"$DIR_S\"
1.107 +set DIROQ = \"$DIR_O\"
1.108 +set DIRCQ = \"$DIR_SCRIPTS\"
1.109 +
1.110 +set COMPAREQ = \"$COMPARE\"
1.111 +set MODELN = \"model1\"
1.112 +
1.113 +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"
1.114 +
1.115 +echo "`date` - Running diagnostics for model ${MODELQ}"
1.116 +# CLAMP metric processing of model1
1.117 +echo "`date` - Running 00.initial.ncl"
1.118 +ncl $INPUT_TEXT $DIR_SCRIPTS/00.initial.ncl
1.119 +echo "`date` - Running 01.npp.ncl"
1.120 +ncl $INPUT_TEXT $DIR_SCRIPTS/01.npp.ncl
1.121 +echo "`date` - Running 01.lai.ncl"
1.122 +ncl $INPUT_TEXT $DIR_SCRIPTS/02.lai.ncl
1.123 +
1.124 +if ($FILE3 != "") then
1.125 +echo "`date` - Running 03.co2.ncl"
1.126 +ncl $INPUT_TEXT $DIR_SCRIPTS/03.co2.ncl
1.127 +endif
1.128 +
1.129 +echo "`date` - Running 04.biomass.ncl"
1.130 +ncl $INPUT_TEXT $DIR_SCRIPTS/04.biomass.ncl
1.131 +echo "`date` - Running 06.fluxnet.ncl"
1.132 +ncl $INPUT_TEXT $DIR_SCRIPTS/06.fluxnet.ncl
1.133 +echo "`date` - Running 07.beta.ncl"
1.134 +ncl $INPUT_TEXT $DIR_SCRIPTS/07.beta.ncl
1.135 +echo "`date` - Running 08.turnover.ncl"
1.136 +ncl $INPUT_TEXT $DIR_SCRIPTS/08.turnover.ncl
1.137 +
1.138 +if ($BGC != "casa") then
1.139 +echo "`date` - Running 09.carbon_sink.ncl"
1.140 +ncl $INPUT_TEXT $DIR_SCRIPTS/09.carbon_sink.ncl
1.141 +else
1.142 +echo "`date` - Running 09x.carbon_sink.ncl"
1.143 +ncl $INPUT_TEXT $DIR_SCRIPTS/09x.carbon_sink.ncl
1.144 +endif
1.145 +
1.146 +if ($BGC != "casa") then
1.147 +echo "`date` - Running 10.fire.ncl"
1.148 +ncl $INPUT_TEXT $DIR_SCRIPTS/10.fire.ncl
1.149 +endif
1.150 +
1.151 +echo "`date` - Running 11.ameriflux.ncl"
1.152 +ncl $INPUT_TEXT $DIR_SCRIPTS/11.ameriflux.ncl
1.153 +echo "`date` - Running 99.final.ncl"
1.154 +ncl $INPUT_TEXT $DIR_SCRIPTS/99.final.ncl
1.155 +
1.156 +#*******************************************************
1.157 +# user modification-(3)
1.158 +
1.159 +# model2
1.160 +set MODEL2 = CN
1.161 +set DIR_M = /lustre/wolf-ddn/scratch/hof/clamp_data/model/
1.162 +#set FILE1 = i01.10cn_1948-2004_ANN_climo.nc
1.163 +set FILE1 = i01.10cn_1975-2000_ANN_climo.nc
1.164 +set FILE2 = i01.10cn_1948-2004_MONS_climo.nc
1.165 +set FILE3 = f02.03cn_1901-1925_MONS_climo_atm.nc
1.166 +set FILE4 = i01.06cn_1980-2004_ANN_climo.nc
1.167 +set FILE5 = i01.10cn_1990-2004_ANN_climo.nc
1.168 +set FILE6 = i01.07cn_1990-2004_ANN_climo.nc
1.169 +set FILE7 = i01.10cn_Fire_C_1979-2004_monthly.nc
1.170 +set FILE8 = i01.10cn_ameriflux_1990-2004_monthly.nc
1.171 +set FILE9 = i01.10cn_2000-2004_ANN_climo.nc
1.172 +set FILE10 = i01.10cn_2000-2004_MONS_climo.nc
1.173 +set FILE11 = i01.10cn_2000_ANN_climo.nc
1.174 +set GRID = T42
1.175 +set BGC = cn
1.176 +set ENERGY = new
1.177 +
1.178 +#*******************************************************
1.179 +# create model2 directory by copying templates
1.180 +if ($FILE3 != "") then
1.181 + set TEMPLATE1 = template_1-model
1.182 +else
1.183 + set TEMPLATE1 = template_1-model_noCO2
1.184 +endif
1.185 +#cp -r $DIR_SCRIPTS/$TEMPLATE1 $MODEL2
1.186 +rm -rf $MODEL2
1.187 +tar xvpf $DIR_SCRIPTS/$TEMPLATE1.tar
1.188 +mv $TEMPLATE1 $MODEL2
1.189 +
1.190 +# add quote, to be usesd in INPUT_TEXT
1.191 +set MODELQ = \"$MODEL2\"
1.192 +set DIRMQ = \"$DIR_M\"
1.193 +set F1 = \"$FILE1\"
1.194 +set F2 = \"$FILE2\"
1.195 +set F3 = \"$FILE3\"
1.196 +set F4 = \"$FILE4\"
1.197 +set F5 = \"$FILE5\"
1.198 +set F6 = \"$FILE6\"
1.199 +set F7 = \"$FILE7\"
1.200 +set F8 = \"$FILE8\"
1.201 +set F9 = \"$FILE9\"
1.202 +set F10 = \"$FILE10\"
1.203 +set F11 = \"$FILE11\"
1.204 +set GRIDQ = \"$GRID\"
1.205 +set BGCQ = \"$BGC\"
1.206 +set ENERGYQ = \"$ENERGY\"
1.207 +
1.208 +set COMPAREQ = \"$COMPARE\"
1.209 +set MODELN = \"model2\"
1.210 +
1.211 +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"
1.212 +
1.213 +echo "`date` - Running diagnostics for model ${MODELQ}"
1.214 +# CLAMP metric processing of model2
1.215 +echo "`date` - Running 00.inital.ncl"
1.216 +ncl $INPUT_TEXT $DIR_SCRIPTS/00.initial.ncl
1.217 +echo "`date` - Running 01.npp.ncl"
1.218 +ncl $INPUT_TEXT $DIR_SCRIPTS/01.npp.ncl
1.219 +echo "`date` - Running 02.lai.ncl"
1.220 +ncl $INPUT_TEXT $DIR_SCRIPTS/02.lai.ncl
1.221 +
1.222 +if ($FILE3 != "") then
1.223 +echo "`date` - Running 03.co2.ncl"
1.224 +ncl $INPUT_TEXT $DIR_SCRIPTS/03.co2.ncl
1.225 +endif
1.226 +
1.227 +echo "`date` - Running 04.biomass.ncl"
1.228 +ncl $INPUT_TEXT $DIR_SCRIPTS/04.biomass.ncl
1.229 +echo "`date` - Running 06.fluxnet.ncl"
1.230 +ncl $INPUT_TEXT $DIR_SCRIPTS/06.fluxnet.ncl
1.231 +echo "`date` - Running 07.beta.ncl"
1.232 +ncl $INPUT_TEXT $DIR_SCRIPTS/07.beta.ncl
1.233 +echo "`date` - Running 08.turnover.ncl"
1.234 +ncl $INPUT_TEXT $DIR_SCRIPTS/08.turnover.ncl
1.235 +
1.236 +if ($BGC != "casa") then
1.237 +echo "`date` - Running 09.carbon_sink.ncl"
1.238 +ncl $INPUT_TEXT $DIR_SCRIPTS/09.carbon_sink.ncl
1.239 +else
1.240 +echo "`date` - Running 09x.carbon_sink.ncl"
1.241 +ncl $INPUT_TEXT $DIR_SCRIPTS/09x.carbon_sink.ncl
1.242 +endif
1.243 +
1.244 +if ($BGC != "casa") then
1.245 +echo "`date` - Running 10.fire.ncl"
1.246 +ncl $INPUT_TEXT $DIR_SCRIPTS/10.fire.ncl
1.247 +endif
1.248 +
1.249 +echo "`date` - Running 11.ameriflux.ncl"
1.250 +ncl $INPUT_TEXT $DIR_SCRIPTS/11.ameriflux.ncl
1.251 +echo "`date` - Running 99.final.ncl"
1.252 +ncl $INPUT_TEXT $DIR_SCRIPTS/99.final.ncl
1.253 +
1.254 +# create a tar file from the final output
1.255 +echo "`date` - tarring up final output into all.tar"
1.256 +tar cf - $MODEL1 $MODEL2 $COMPARE > all.tar
1.257 +