Modifications to scoring and graphics production for the final version of code for the C-LAMP paper in GCB.
2 #-------------------------------------------------------------------
3 # note: user modifies ONLY the "user modification" section
6 # DIR_M : directory of model data
7 # DIR_O : directory of observed data
8 # DIR_S : directory of model surface data
9 # DIR_SCRIPTS : directory of run scripts
10 # FILE1 : time_mean climatology from CLM diagnostic package
11 # FILE2 : 12-monthly climatology from CLM diagnostic package
12 # FILE3 : 12-monthly climatology from ATM diagnostic package
13 # leave it blank, if no ATM file:
15 # FILE7 : timeseries file generated from ....
16 # FILE8 : timeseries file generated from ....
17 # GRID : T31, T42, or 1.9
19 # ENERGY: new or old (model data fields)
20 #********************************************************************
23 # model data (no ATM file)
25 set DIR_M = /fis/cgd/cseg/people/jeff/clamp_data/model/
26 set FILE1 = i01.10cn_1948-2004_ANN_climo.nc
27 set FILE2 = i01.10cn_1948-2004_MONS_climo.nc
29 set FILE4 = i01.06cn_1980-2004_ANN_climo.nc
30 set FILE5 = i01.10cn_1990-2004_ANN_climo.nc
31 set FILE6 = i01.07cn_1990-2004_ANN_climo.nc
32 set FILE7 = i01.10cn_Fire_C_1979-2004_monthly.nc
33 set FILE8 = i01.10cn_ameriflux_1990-2004_monthly.nc
38 # in the "CLAMP metric processing" section:
39 # only 00.initial.ncl and 99.final.ncl are required,
40 # user can comment out any one or more of the other ncl scripts,
42 # #ncl $INPUT_TEXT $DIR_SCRIPTS/10.fire.ncl
45 set DIR_S = /fis/cgd/cseg/people/jeff/clamp_data/surface_model/
48 set DIR_O = /fis/cgd/cseg/people/jeff/clamp_data/observed/
50 # directory for scripts, templates and ncl files
51 set DIR_SCRIPTS = /fis/cgd/cseg/people/jeff/clamp/all/
53 #********************************************************************
55 # create a directory for the model by copying from a template
56 if ($FILE3 != "") then
57 set TEMPLATE = template_1-model
59 set TEMPLATE = template_1-model_noCO2
61 #cp -r $DIR_SCRIPTS/$TEMPLATE $MODEL
63 tar xvpf $DIR_SCRIPTS/$TEMPLATE.tar
66 # add quote, to be usesd in INPUT_TEXT
67 set MODELQ = \"$MODEL\"
68 set DIRMQ = \"$DIR_M\"
79 set ENERGYQ = \"$ENERGY\"
80 set DIRSQ = \"$DIR_S\"
81 set DIROQ = \"$DIR_O\"
83 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"
85 # CLAMP metric processing
86 ncl $INPUT_TEXT $DIR_SCRIPTS/00.initial.ncl
87 ncl $INPUT_TEXT $DIR_SCRIPTS/01.npp.ncl
88 ncl $INPUT_TEXT $DIR_SCRIPTS/02.lai.ncl
90 if ($FILE3 != "") then
91 ncl $INPUT_TEXT $DIR_SCRIPTS/03.co2.ncl
94 ncl $INPUT_TEXT $DIR_SCRIPTS/04.biomass.ncl
95 ncl $INPUT_TEXT $DIR_SCRIPTS/06.fluxnet.ncl
96 ncl $INPUT_TEXT $DIR_SCRIPTS/07.beta.ncl
97 ncl $INPUT_TEXT $DIR_SCRIPTS/08.turnover.ncl
99 if ($BGC != "casa") then
100 ncl $INPUT_TEXT $DIR_SCRIPTS/09.carbon_sink.ncl
102 ncl $INPUT_TEXT $DIR_SCRIPTS/09x.carbon_sink.ncl
105 if ($BGC != "casa") then
106 ncl $INPUT_TEXT $DIR_SCRIPTS/10.fire.ncl
109 ncl $INPUT_TEXT $DIR_SCRIPTS/11.ameriflux.ncl
110 ncl $INPUT_TEXT $DIR_SCRIPTS/99.final.ncl
112 # create a tar file from the final output
113 tar cf $MODEL.tar $MODEL