all/run_1-model.csh
changeset 0 0c6405ab2ff4
child 1 4be95183fbcd
     1.1 --- /dev/null	Thu Jan 01 00:00:00 1970 +0000
     1.2 +++ b/all/run_1-model.csh	Mon Jan 26 22:08:20 2009 -0500
     1.3 @@ -0,0 +1,112 @@
     1.4 +#!/bin/csh -f
     1.5 +#-------------------------------------------------------------------
     1.6 +# note: user modifies ONLY the "user modification" section 
     1.7 +#
     1.8 +#       MODEL : model name
     1.9 +#       DIR_M  : directory of model data
    1.10 +#       DIR_O  : directory of observed data
    1.11 +#       DIR_S  : directory of model surface data
    1.12 +#       DIR_SCRIPTS  : directory of run scripts
    1.13 +#       FILE1 : time_mean  climatology from CLM diagnostic package
    1.14 +#       FILE2 : 12-monthly climatology from CLM diagnostic package
    1.15 +#       FILE3 : 12-monthly climatology from ATM diagnostic package
    1.16 +#               leave it blank, if no ATM file:
    1.17 +#               set FILE3 =
    1.18 +#       FILE7 : timeseries file generated from .... 
    1.19 +#       FILE8 : timeseries file generated from .... 
    1.20 +#       GRID  : T31, T42, or 1.9
    1.21 +#       BGC   : cn or casa 
    1.22 +#       ENERGY: new or old (model data fields)
    1.23 +#********************************************************************
    1.24 +# user modification:
    1.25 +
    1.26 +# model data (no ATM file)
    1.27 +set MODEL  = CN
    1.28 +set DIR_M  = /fis/cgd/cseg/people/jeff/clamp_data/model/
    1.29 +set FILE1  = i01.10cn_1948-2004_ANN_climo.nc
    1.30 +set FILE2  = i01.10cn_1948-2004_MONS_climo.nc
    1.31 +set FILE3  =
    1.32 +set FILE4  = i01.06cn_1980-2004_ANN_climo.nc
    1.33 +set FILE5  = i01.10cn_1990-2004_ANN_climo.nc
    1.34 +set FILE6  = i01.07cn_1990-2004_ANN_climo.nc
    1.35 +set FILE7  = i01.10cn_Fire_C_1979-2004_monthly.nc
    1.36 +set FILE8  = i01.10cn_ameriflux_1990-2004_monthly.nc
    1.37 +set GRID   = T42
    1.38 +set BGC    = cn
    1.39 +set ENERGY = new
    1.40 +
    1.41 +# in the "CLAMP metric processing" section:
    1.42 +#  only 00.initial.ncl and 99.final.ncl are required,
    1.43 +#  user can comment out any one or more of the other ncl scripts,
    1.44 +#  e.g. 
    1.45 +#  #ncl $INPUT_TEXT $DIR_SCRIPTS/10.fire.ncl
    1.46 +
    1.47 +# model surface data
    1.48 +set DIR_S  = /fis/cgd/cseg/people/jeff/clamp_data/surface_model/
    1.49 +
    1.50 +# observed data
    1.51 +set DIR_O  = /fis/cgd/cseg/people/jeff/clamp_data/observed/
    1.52 +
    1.53 +# directory for scripts, templates and ncl files
    1.54 +set DIR_SCRIPTS = /fis/cgd/cseg/people/jeff/clamp/all/
    1.55 +
    1.56 +#********************************************************************
    1.57 +
    1.58 +# create a directory for the model by copying from a template
    1.59 +if ($FILE3 != "") then
    1.60 +   set TEMPLATE = template_1-model
    1.61 +else
    1.62 +   set TEMPLATE = template_1-model_noCO2
    1.63 +endif
    1.64 +cp -r $DIR_SCRIPTS/$TEMPLATE $MODEL
    1.65 +
    1.66 +# add quote, to be usesd in INPUT_TEXT
    1.67 +set MODELQ = \"$MODEL\"
    1.68 +set DIRMQ  = \"$DIR_M\"
    1.69 +set F1  = \"$FILE1\"
    1.70 +set F2  = \"$FILE2\"
    1.71 +set F3  = \"$FILE3\"
    1.72 +set F4  = \"$FILE4\"
    1.73 +set F5  = \"$FILE5\"
    1.74 +set F6  = \"$FILE6\"
    1.75 +set F7  = \"$FILE7\"
    1.76 +set F8  = \"$FILE8\"
    1.77 +set GRIDQ   = \"$GRID\"
    1.78 +set BGCQ    = \"$BGC\"
    1.79 +set ENERGYQ = \"$ENERGY\"
    1.80 +set DIRSQ   = \"$DIR_S\"
    1.81 +set DIROQ   = \"$DIR_O\"
    1.82 +
    1.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"
    1.84 + 
    1.85 +# CLAMP metric processing
    1.86 +ncl $INPUT_TEXT $DIR_SCRIPTS/00.initial.ncl
    1.87 +ncl $INPUT_TEXT $DIR_SCRIPTS/01.npp.ncl
    1.88 +ncl $INPUT_TEXT $DIR_SCRIPTS/02.lai.ncl
    1.89 +
    1.90 +if ($FILE3 != "") then
    1.91 +ncl $INPUT_TEXT $DIR_SCRIPTS/03.co2.ncl
    1.92 +endif
    1.93 +
    1.94 +ncl $INPUT_TEXT $DIR_SCRIPTS/04.biomass.ncl
    1.95 +ncl $INPUT_TEXT $DIR_SCRIPTS/06.fluxnet.ncl
    1.96 +ncl $INPUT_TEXT $DIR_SCRIPTS/07.beta.ncl
    1.97 +ncl $INPUT_TEXT $DIR_SCRIPTS/08.turnover.ncl
    1.98 +
    1.99 +if ($BGC != "casa") then
   1.100 +ncl $INPUT_TEXT $DIR_SCRIPTS/09.carbon_sink.ncl
   1.101 +else
   1.102 +ncl $INPUT_TEXT $DIR_SCRIPTS/09x.carbon_sink.ncl
   1.103 +endif
   1.104 +
   1.105 +if ($BGC != "casa") then
   1.106 +ncl $INPUT_TEXT $DIR_SCRIPTS/10.fire.ncl
   1.107 +endif
   1.108 +
   1.109 +ncl $INPUT_TEXT $DIR_SCRIPTS/11.ameriflux.ncl
   1.110 +ncl $INPUT_TEXT $DIR_SCRIPTS/99.final.ncl
   1.111 +
   1.112 +# create a tar file from the final output 
   1.113 +tar cf $MODEL.tar $MODEL
   1.114 +
   1.115 +