1.1 --- /dev/null Thu Jan 01 00:00:00 1970 +0000
1.2 +++ b/all/run_2-model.csh Mon Jan 26 22:08:20 2009 -0500
1.3 @@ -0,0 +1,201 @@
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 = CN
1.38 +set DIR_M = /lustre/wolf-ddn/scratch/hof/clamp_data/model/
1.39 +set FILE1 = i01.10cn_1948-2004_ANN_climo.nc
1.40 +set FILE2 = i01.10cn_1948-2004_MONS_climo.nc
1.41 +set FILE3 = f02.03cn_1901-1925_MONS_climo_atm.nc
1.42 +set FILE4 = i01.06cn_1980-2004_ANN_climo.nc
1.43 +set FILE5 = i01.10cn_1990-2004_ANN_climo.nc
1.44 +set FILE6 = i01.07cn_1990-2004_ANN_climo.nc
1.45 +set FILE7 = i01.10cn_Fire_C_1979-2004_monthly.nc
1.46 +set FILE8 = i01.10cn_ameriflux_1990-2004_monthly.nc
1.47 +set GRID = T42
1.48 +set BGC = cn
1.49 +set ENERGY = new
1.50 +
1.51 +# in the "CLAMP metric processing" section:
1.52 +# only 00.initial.ncl and 99.final.ncl are required,
1.53 +# user can comment out any one or more of the other ncl scripts,
1.54 +# e.g.
1.55 +# #ncl $INPUT_TEXT $DIR_SCRIPTS/10.fire.ncl
1.56 +
1.57 +# model surface data
1.58 +set DIR_S = /lustre/wolf-ddn/scratch/hof/clamp_data/surface_model/
1.59 +
1.60 +# observed data
1.61 +set DIR_O = /lustre/wolf-ddn/scratch/hof/clamp_data/observed/
1.62 +
1.63 +# directory for scripts, templates and ncl files
1.64 +set DIR_SCRIPTS = /lustre/wolf-ddn/scratch/hof/clamp/all/
1.65 +
1.66 +#********************************************************
1.67 +
1.68 +# create model1 and model1_vs_model2 directory by copying templates
1.69 +if ($FILE3 != "") then
1.70 + set TEMPLATE1 = template_1-model
1.71 + set TEMPLATE2 = template_2-model
1.72 +else
1.73 + set TEMPLATE1 = template_1-model_noCO2
1.74 + set TEMPLATE2 = template_2-model_noCO2
1.75 +endif
1.76 +cp -r $DIR_SCRIPTS/$TEMPLATE1 $MODEL1
1.77 +cp -r $DIR_SCRIPTS/$TEMPLATE2 $COMPARE
1.78 +
1.79 +# add quote, to be usesd in INPUT_TEXT
1.80 +set MODELQ = \"$MODEL1\"
1.81 +set DIRMQ = \"$DIR_M\"
1.82 +set F1 = \"$FILE1\"
1.83 +set F2 = \"$FILE2\"
1.84 +set F3 = \"$FILE3\"
1.85 +set F4 = \"$FILE4\"
1.86 +set F5 = \"$FILE5\"
1.87 +set F6 = \"$FILE6\"
1.88 +set F7 = \"$FILE7\"
1.89 +set F8 = \"$FILE8\"
1.90 +set GRIDQ = \"$GRID\"
1.91 +set BGCQ = \"$BGC\"
1.92 +set ENERGYQ = \"$ENERGY\"
1.93 +set DIRSQ = \"$DIR_S\"
1.94 +set DIROQ = \"$DIR_O\"
1.95 +set DIRCQ = \"$DIR_SCRIPTS\"
1.96 +
1.97 +set COMPAREQ = \"$COMPARE\"
1.98 +set MODELN = \"model1\"
1.99 +
1.100 +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"
1.101 +
1.102 +# CLAMP metric processing of model1
1.103 +ncl $INPUT_TEXT $DIR_SCRIPTS/00.initial.ncl
1.104 +ncl $INPUT_TEXT $DIR_SCRIPTS/01.npp.ncl
1.105 +ncl $INPUT_TEXT $DIR_SCRIPTS/02.lai.ncl
1.106 +
1.107 +if ($FILE3 != "") then
1.108 +ncl $INPUT_TEXT $DIR_SCRIPTS/03.co2.ncl
1.109 +endif
1.110 +
1.111 +ncl $INPUT_TEXT $DIR_SCRIPTS/04.biomass.ncl
1.112 +ncl $INPUT_TEXT $DIR_SCRIPTS/06.fluxnet.ncl
1.113 +ncl $INPUT_TEXT $DIR_SCRIPTS/07.beta.ncl
1.114 +ncl $INPUT_TEXT $DIR_SCRIPTS/08.turnover.ncl
1.115 +
1.116 +if ($BGC != "casa") then
1.117 +ncl $INPUT_TEXT $DIR_SCRIPTS/09.carbon_sink.ncl
1.118 +else
1.119 +ncl $INPUT_TEXT $DIR_SCRIPTS/09x.carbon_sink.ncl
1.120 +endif
1.121 +
1.122 +if ($BGC != "casa") then
1.123 +ncl $INPUT_TEXT $DIR_SCRIPTS/10.fire.ncl
1.124 +endif
1.125 +
1.126 +ncl $INPUT_TEXT $DIR_SCRIPTS/11.ameriflux.ncl
1.127 +ncl $INPUT_TEXT $DIR_SCRIPTS/99.final.ncl
1.128 +
1.129 +#*******************************************************
1.130 +# user modification-(3)
1.131 +
1.132 +# model2
1.133 +set MODEL2 = CASA
1.134 +set DIR_M = /lustre/wolf-ddn/scratch/hof/clamp_data/model/
1.135 +set FILE1 = i01.10casa_1948-2004_ANN_climo.nc
1.136 +set FILE2 = i01.10casa_1948-2004_MONS_climo.nc
1.137 +set FILE3 = f02.03casa_1876-1900_MONS_climo_atm.nc
1.138 +set FILE4 = i01.06casa_1980-2004_ANN_climo.nc
1.139 +set FILE5 = i01.10casa_1990-2004_ANN_climo.nc
1.140 +set FILE6 = i01.07casa_1990-2004_ANN_climo.nc
1.141 +set FILE7 = i01.10casa_Fire_C_1979-2004_monthly.nc
1.142 +set FILE8 = i01.10casa_ameriflux_1990-2004_monthly.nc
1.143 +set GRID = T42
1.144 +set BGC = casa
1.145 +set ENERGY = new
1.146 +#*******************************************************
1.147 +# create model2 directory by copying templates
1.148 +if ($FILE3 != "") then
1.149 + set TEMPLATE1 = template_1-model
1.150 +else
1.151 + set TEMPLATE1 = template_1-model_noCO2
1.152 +endif
1.153 +cp -r $DIR_SCRIPTS/$TEMPLATE1 $MODEL2
1.154 +
1.155 +# add quote, to be usesd in INPUT_TEXT
1.156 +set MODELQ = \"$MODEL2\"
1.157 +set DIRMQ = \"$DIR_M\"
1.158 +set F1 = \"$FILE1\"
1.159 +set F2 = \"$FILE2\"
1.160 +set F3 = \"$FILE3\"
1.161 +set F4 = \"$FILE4\"
1.162 +set F5 = \"$FILE5\"
1.163 +set F6 = \"$FILE6\"
1.164 +set F7 = \"$FILE7\"
1.165 +set F8 = \"$FILE8\"
1.166 +set GRIDQ = \"$GRID\"
1.167 +set BGCQ = \"$BGC\"
1.168 +set ENERGYQ = \"$ENERGY\"
1.169 +
1.170 +set COMPAREQ = \"$COMPARE\"
1.171 +set MODELN = \"model2\"
1.172 +
1.173 +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"
1.174 +
1.175 +# CLAMP metric processing of model2
1.176 +ncl $INPUT_TEXT $DIR_SCRIPTS/00.initial.ncl
1.177 +ncl $INPUT_TEXT $DIR_SCRIPTS/01.npp.ncl
1.178 +ncl $INPUT_TEXT $DIR_SCRIPTS/02.lai.ncl
1.179 +
1.180 +if ($FILE3 != "") then
1.181 +ncl $INPUT_TEXT $DIR_SCRIPTS/03.co2.ncl
1.182 +endif
1.183 +
1.184 +ncl $INPUT_TEXT $DIR_SCRIPTS/04.biomass.ncl
1.185 +ncl $INPUT_TEXT $DIR_SCRIPTS/06.fluxnet.ncl
1.186 +ncl $INPUT_TEXT $DIR_SCRIPTS/07.beta.ncl
1.187 +ncl $INPUT_TEXT $DIR_SCRIPTS/08.turnover.ncl
1.188 +
1.189 +if ($BGC != "casa") then
1.190 +ncl $INPUT_TEXT $DIR_SCRIPTS/09.carbon_sink.ncl
1.191 +else
1.192 +ncl $INPUT_TEXT $DIR_SCRIPTS/09x.carbon_sink.ncl
1.193 +endif
1.194 +
1.195 +if ($BGC != "casa") then
1.196 +ncl $INPUT_TEXT $DIR_SCRIPTS/10.fire.ncl
1.197 +endif
1.198 +
1.199 +ncl $INPUT_TEXT $DIR_SCRIPTS/11.ameriflux.ncl
1.200 +ncl $INPUT_TEXT $DIR_SCRIPTS/99.final.ncl
1.201 +
1.202 +# create a tar file from the final output
1.203 +tar cf - $MODEL1 $MODEL2 $COMPARE > all.tar
1.204 +