1.1 --- /dev/null Thu Jan 01 00:00:00 1970 +0000
1.2 +++ b/all/run_2-model_37cn_cont_42cn.csh Thu Mar 26 14:02:21 2009 -0400
1.3 @@ -0,0 +1,247 @@
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 = i01.37cn_cont_vs_i01.42cn
1.32 +
1.33 +#*******************************************************
1.34 +# user modification-(2)
1.35 +
1.36 +# model1
1.37 +set MODEL1 = i01.37cn_cont
1.38 +set DIR_M = /lustre/wolf-ddn/scratch/hof/cases/i01.37cn_cont_out/
1.39 +set FILE1 = i01.37cn_cont_out_ANN_climo.nc
1.40 +set FILE2 = i01.37cn_cont_out_MONS_climo.nc
1.41 +set FILE3 =
1.42 +set FILE4 = $FILE1
1.43 +set FILE5 = $FILE1
1.44 +set FILE6 = $FILE1
1.45 +set FILE7 = i01.37cn_cont_Fire_C_2976-3000_monthly.nc
1.46 +set FILE8 = i01.37cn_cont_ameriflux_2976-3000_monthly.nc
1.47 +set FILE9 = $FILE1
1.48 +set FILE10 = $FILE2
1.49 +set GRID = 1.9
1.50 +set BGC = cn
1.51 +set ENERGY = old
1.52 +
1.53 +# in the "CLAMP metric processing" section:
1.54 +# only 00.initial.ncl and 99.final.ncl are required,
1.55 +# user can comment out any one or more of the other ncl scripts,
1.56 +# e.g.
1.57 +# #ncl $INPUT_TEXT $DIR_SCRIPTS/10.fire.ncl
1.58 +
1.59 +# model surface data
1.60 +set DIR_S = /lustre/wolf-ddn/scratch/hof/clamp_data/surface_model/
1.61 +
1.62 +# observed data
1.63 +set DIR_O = /lustre/wolf-ddn/scratch/hof/clamp_data/observed/
1.64 +
1.65 +# directory for scripts, templates and ncl files
1.66 +set DIR_SCRIPTS = /lustre/wolf-ddn/scratch/hof/clamp/all/
1.67 +
1.68 +#********************************************************
1.69 +
1.70 +# create model1 and model1_vs_model2 directory by copying templates
1.71 +if ($FILE3 != "") then
1.72 + set TEMPLATE1 = template_1-model
1.73 + set TEMPLATE2 = template_2-model
1.74 +else
1.75 + set TEMPLATE1 = template_1-model_noCO2
1.76 + set TEMPLATE2 = template_2-model_noCO2
1.77 +endif
1.78 +#cp -r $DIR_SCRIPTS/$TEMPLATE1 $MODEL1
1.79 +#cp -r $DIR_SCRIPTS/$TEMPLATE2 $COMPARE
1.80 +rm -rf $MODEL1
1.81 +rm -rf $COMPARE
1.82 +tar xvpf $DIR_SCRIPTS/$TEMPLATE1.tar
1.83 +mv $TEMPLATE1 $MODEL1
1.84 +tar xvpf $DIR_SCRIPTS/$TEMPLATE2.tar
1.85 +mv $TEMPLATE2 $COMPARE
1.86 +
1.87 +# add quote, to be usesd in INPUT_TEXT
1.88 +set MODELQ = \"$MODEL1\"
1.89 +set DIRMQ = \"$DIR_M\"
1.90 +set F1 = \"$FILE1\"
1.91 +set F2 = \"$FILE2\"
1.92 +set F3 = \"$FILE3\"
1.93 +set F4 = \"$FILE4\"
1.94 +set F5 = \"$FILE5\"
1.95 +set F6 = \"$FILE6\"
1.96 +set F7 = \"$FILE7\"
1.97 +set F8 = \"$FILE8\"
1.98 +set F9 = \"$FILE9\"
1.99 +set F10 = \"$FILE10\"
1.100 +set GRIDQ = \"$GRID\"
1.101 +set BGCQ = \"$BGC\"
1.102 +set ENERGYQ = \"$ENERGY\"
1.103 +set DIRSQ = \"$DIR_S\"
1.104 +set DIROQ = \"$DIR_O\"
1.105 +set DIRCQ = \"$DIR_SCRIPTS\"
1.106 +
1.107 +set COMPAREQ = \"$COMPARE\"
1.108 +set MODELN = \"model1\"
1.109 +
1.110 +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"
1.111 +
1.112 +echo "`date` - Running diagnostics for model ${MODELQ}"
1.113 +# CLAMP metric processing of model1
1.114 +echo "`date` - Running 00.initial.ncl"
1.115 +ncl $INPUT_TEXT $DIR_SCRIPTS/00.initial.ncl
1.116 +echo "`date` - Running 01.npp.ncl"
1.117 +ncl $INPUT_TEXT $DIR_SCRIPTS/01.npp.ncl
1.118 +echo "`date` - Running 01.lai.ncl"
1.119 +ncl $INPUT_TEXT $DIR_SCRIPTS/02.lai.ncl
1.120 +
1.121 +if ($FILE3 != "") then
1.122 +echo "`date` - Running 03.co2.ncl"
1.123 +ncl $INPUT_TEXT $DIR_SCRIPTS/03.co2.ncl
1.124 +endif
1.125 +
1.126 +echo "`date` - Running 04.biomass.ncl"
1.127 +ncl $INPUT_TEXT $DIR_SCRIPTS/04.biomass.ncl
1.128 +echo "`date` - Running 06.fluxnet.ncl"
1.129 +ncl $INPUT_TEXT $DIR_SCRIPTS/06.fluxnet.ncl
1.130 +echo "`date` - Running 07.beta.ncl"
1.131 +ncl $INPUT_TEXT $DIR_SCRIPTS/07.beta.ncl
1.132 +echo "`date` - Running 08.turnover.ncl"
1.133 +ncl $INPUT_TEXT $DIR_SCRIPTS/08.turnover.ncl
1.134 +
1.135 +if ($BGC != "casa") then
1.136 +echo "`date` - Running 09.carbon_sink.ncl"
1.137 +ncl $INPUT_TEXT $DIR_SCRIPTS/09.carbon_sink.ncl
1.138 +else
1.139 +echo "`date` - Running 09x.carbon_sink.ncl"
1.140 +ncl $INPUT_TEXT $DIR_SCRIPTS/09x.carbon_sink.ncl
1.141 +endif
1.142 +
1.143 +if ($BGC != "casa") then
1.144 +echo "`date` - Running 10.fire.ncl"
1.145 +ncl $INPUT_TEXT $DIR_SCRIPTS/10.fire.ncl
1.146 +endif
1.147 +
1.148 +echo "`date` - Running 11.ameriflux.ncl"
1.149 +ncl $INPUT_TEXT $DIR_SCRIPTS/11.ameriflux.ncl
1.150 +echo "`date` - Running 99.final.ncl"
1.151 +ncl $INPUT_TEXT $DIR_SCRIPTS/99.final.ncl
1.152 +
1.153 +#*******************************************************
1.154 +# user modification-(3)
1.155 +
1.156 +# model2
1.157 +set MODEL2 = i01.42cn
1.158 +set DIR_M = /lustre/wolf-ddn/scratch/hof/cases/i01.42cn_out/
1.159 +set FILE1 = i01.42cn_out_ANN_climo.nc
1.160 +set FILE2 = i01.42cn_out_MONS_climo.nc
1.161 +set FILE3 =
1.162 +set FILE4 = $FILE1
1.163 +set FILE5 = $FILE1
1.164 +set FILE6 = $FILE1
1.165 +set FILE7 = i01.42cn_Fire_C_2976-3000_monthly.nc
1.166 +set FILE8 = i01.42cn_ameriflux_2976-3000_monthly.nc
1.167 +set FILE9 = $FILE1
1.168 +set FILE10 = $FILE2
1.169 +set GRID = 1.9
1.170 +set BGC = cn
1.171 +set ENERGY = old
1.172 +#*******************************************************
1.173 +# create model2 directory by copying templates
1.174 +if ($FILE3 != "") then
1.175 + set TEMPLATE1 = template_1-model
1.176 +else
1.177 + set TEMPLATE1 = template_1-model_noCO2
1.178 +endif
1.179 +#cp -r $DIR_SCRIPTS/$TEMPLATE1 $MODEL2
1.180 +rm -rf $MODEL2
1.181 +tar xvpf $DIR_SCRIPTS/$TEMPLATE1.tar
1.182 +mv $TEMPLATE1 $MODEL2
1.183 +
1.184 +# add quote, to be usesd in INPUT_TEXT
1.185 +set MODELQ = \"$MODEL2\"
1.186 +set DIRMQ = \"$DIR_M\"
1.187 +set F1 = \"$FILE1\"
1.188 +set F2 = \"$FILE2\"
1.189 +set F3 = \"$FILE3\"
1.190 +set F4 = \"$FILE4\"
1.191 +set F5 = \"$FILE5\"
1.192 +set F6 = \"$FILE6\"
1.193 +set F7 = \"$FILE7\"
1.194 +set F8 = \"$FILE8\"
1.195 +set F9 = \"$FILE9\"
1.196 +set F10 = \"$FILE10\"
1.197 +set GRIDQ = \"$GRID\"
1.198 +set BGCQ = \"$BGC\"
1.199 +set ENERGYQ = \"$ENERGY\"
1.200 +
1.201 +set COMPAREQ = \"$COMPARE\"
1.202 +set MODELN = \"model2\"
1.203 +
1.204 +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"
1.205 +
1.206 +echo "`date` - Running diagnostics for model ${MODELQ}"
1.207 +# CLAMP metric processing of model2
1.208 +echo "`date` - Running 00.inital.ncl"
1.209 +ncl $INPUT_TEXT $DIR_SCRIPTS/00.initial.ncl
1.210 +echo "`date` - Running 01.npp.ncl"
1.211 +ncl $INPUT_TEXT $DIR_SCRIPTS/01.npp.ncl
1.212 +echo "`date` - Running 02.lai.ncl"
1.213 +ncl $INPUT_TEXT $DIR_SCRIPTS/02.lai.ncl
1.214 +
1.215 +if ($FILE3 != "") then
1.216 +echo "`date` - Running 03.co2.ncl"
1.217 +ncl $INPUT_TEXT $DIR_SCRIPTS/03.co2.ncl
1.218 +endif
1.219 +
1.220 +echo "`date` - Running 04.biomass.ncl"
1.221 +ncl $INPUT_TEXT $DIR_SCRIPTS/04.biomass.ncl
1.222 +echo "`date` - Running 06.fluxnet.ncl"
1.223 +ncl $INPUT_TEXT $DIR_SCRIPTS/06.fluxnet.ncl
1.224 +echo "`date` - Running 07.beta.ncl"
1.225 +ncl $INPUT_TEXT $DIR_SCRIPTS/07.beta.ncl
1.226 +echo "`date` - Running 08.turnover.ncl"
1.227 +ncl $INPUT_TEXT $DIR_SCRIPTS/08.turnover.ncl
1.228 +
1.229 +if ($BGC != "casa") then
1.230 +echo "`date` - Running 09.carbon_sink.ncl"
1.231 +ncl $INPUT_TEXT $DIR_SCRIPTS/09.carbon_sink.ncl
1.232 +else
1.233 +echo "`date` - Running 09x.carbon_sink.ncl"
1.234 +ncl $INPUT_TEXT $DIR_SCRIPTS/09x.carbon_sink.ncl
1.235 +endif
1.236 +
1.237 +if ($BGC != "casa") then
1.238 +echo "`date` - Running 10.fire.ncl"
1.239 +ncl $INPUT_TEXT $DIR_SCRIPTS/10.fire.ncl
1.240 +endif
1.241 +
1.242 +echo "`date` - Running 11.ameriflux.ncl"
1.243 +ncl $INPUT_TEXT $DIR_SCRIPTS/11.ameriflux.ncl
1.244 +echo "`date` - Running 99.final.ncl"
1.245 +ncl $INPUT_TEXT $DIR_SCRIPTS/99.final.ncl
1.246 +
1.247 +# create a tar file from the final output
1.248 +echo "`date` - tarring up final output into all.tar"
1.249 +tar cf - $MODEL1 $MODEL2 $COMPARE > all.tar
1.250 +