Mercurial Mercurial > repos > public_software > clamp / file revision
summary | shortlog | changelog | graph | tags | bookmarks | branches | files | changeset | file | latest | revisions | annotate | diff | comparison | raw | help
all/run_2-model_FMH.csh
author Forrest Hoffman <forrest@climatemodeling.org>
Thu, 26 Mar 2009 14:02:21 -0400
changeset 1 4be95183fbcd
permissions -rwxr-xr-x
Modifications to scoring and graphics production for the final version of code for the C-LAMP paper in GCB. 
     1 #!/bin/csh -f

     2 #-------------------------------------------------------------------

     3 # note: user modifies ONLY the "user modification" section 

     4 #

     5 #       COMPARE: model1 vs model2

     6 #       MODELn : model name

     7 #       DIR_M  : directory of model data

     8 #       DIR_O  : directory of observed data

     9 #       DIR_S  : directory of model surface data

    10 #       DIR_SCRIPTS  : directory of run scripts

    11 #       FILE1  : time_mean  climatology from CLM diagnostic package

    12 #       FILE2  : 12-monthly climatology from CLM diagnostic package

    13 #       FILE3  : 12-monthly climatology from ATM diagnostic package

    14 #                leave it blank, if no ATM file:

    15 #                set FILE3 =

    16 #       FILE7 : timeseries file generated from .... 

    17 #       FILE8 : timeseries file generated from .... 

    18 #       GRID   : T31, T42, or 1.9

    19 #       BGC    : cn or casa 

    20 #       ENERGY : new or old (fields in model data)

    21 #-------------------------------------------------------------------

    22 

    23 #*******************************************************

    24 # user modification-(1)

    25 

    26 # directory name of model comparison

    27 #et COMPARE = b30.061n_vs_b30.061m

    28 set COMPARE = CN_vs_CASA 

    29 

    30 #*******************************************************

    31 # user modification-(2)

    32 

    33 # model1 

    34 set MODEL1 = CASA

    35 set DIR_M  = /lustre/wolf-ddn/scratch/hof/clamp_data/model/

    36 #set FILE1  = i01.10casa_1948-2004_ANN_climo.nc

    37 set FILE1  = i01.10casa_1975-2000_ANN_climo.nc

    38 set FILE2  = i01.10casa_1948-2004_MONS_climo.nc

    39 set FILE3  = f02.03casa_1876-1900_MONS_climo_atm.nc

    40 set FILE4  = i01.06casa_1980-2004_ANN_climo.nc

    41 set FILE5  = i01.10casa_1990-2004_ANN_climo.nc

    42 set FILE6  = i01.07casa_1990-2004_ANN_climo.nc

    43 set FILE7  = i01.10casa_Fire_C_1979-2004_monthly.nc

    44 set FILE8  = i01.10casa_ameriflux_1990-2004_monthly.nc

    45 set FILE9  = i01.10casa_2000-2004_ANN_climo.nc

    46 set FILE10 = i01.10casa_2000-2004_MONS_climo.nc

    47 set FILE11 = i01.10casa_2000_ANN_climo.nc

    48 set GRID   = T42

    49 set BGC    = casa

    50 set ENERGY = new

    51 

    52 # in the "CLAMP metric processing" section:

    53 #  only 00.initial.ncl and 99.final.ncl are required,

    54 #  user can comment out any one or more of the other ncl scripts,

    55 #  e.g. 

    56 #  #ncl $INPUT_TEXT $DIR_SCRIPTS/10.fire.ncl

    57 

    58 # model surface data

    59 set DIR_S  = /lustre/wolf-ddn/scratch/hof/clamp_data/surface_model/

    60 

    61 # observed data

    62 set DIR_O  = /lustre/wolf-ddn/scratch/hof/clamp_data/observed/

    63 

    64 # directory for scripts, templates and ncl files

    65 set DIR_SCRIPTS = /lustre/wolf-ddn/scratch/hof/clamp/all/

    66 

    67 #********************************************************

    68 

    69 # create model1 and model1_vs_model2 directory by copying templates

    70 if ($FILE3 != "") then

    71    set TEMPLATE1 = template_1-model

    72    set TEMPLATE2 = template_2-model

    73 else

    74    set TEMPLATE1 = template_1-model_noCO2

    75    set TEMPLATE2 = template_2-model_noCO2

    76 endif

    77 #cp -r $DIR_SCRIPTS/$TEMPLATE1 $MODEL1

    78 #cp -r $DIR_SCRIPTS/$TEMPLATE2 $COMPARE

    79 rm -rf $MODEL1

    80 rm -rf $COMPARE

    81 tar xvpf $DIR_SCRIPTS/$TEMPLATE1.tar

    82 mv $TEMPLATE1 $MODEL1

    83 tar xvpf $DIR_SCRIPTS/$TEMPLATE2.tar

    84 mv $TEMPLATE2 $COMPARE

    85 

    86 # add quote, to be usesd in INPUT_TEXT

    87 set MODELQ = \"$MODEL1\"

    88 set DIRMQ  = \"$DIR_M\"

    89 set F1  = \"$FILE1\"

    90 set F2  = \"$FILE2\"

    91 set F3  = \"$FILE3\"

    92 set F4  = \"$FILE4\"

    93 set F5  = \"$FILE5\"

    94 set F6  = \"$FILE6\"

    95 set F7  = \"$FILE7\"

    96 set F8  = \"$FILE8\"

    97 set F9  = \"$FILE9\"

    98 set F10 = \"$FILE10\"

    99 set F11 = \"$FILE11\"

   100 set GRIDQ   = \"$GRID\"

   101 set BGCQ    = \"$BGC\"

   102 set ENERGYQ = \"$ENERGY\"

   103 set DIRSQ   = \"$DIR_S\"

   104 set DIROQ   = \"$DIR_O\"

   105 set DIRCQ   = \"$DIR_SCRIPTS\"

   106 

   107 set COMPAREQ = \"$COMPARE\"

   108 set MODELN   = \"model1\"

   109 

   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 film11=$F11 BGC=$BGCQ ENERGY=$ENERGYQ dirs=$DIRSQ diro=$DIROQ dirscript=$DIRCQ modeln=$MODELN compare=$COMPAREQ"

   111  

   112 echo "`date` - Running diagnostics for model ${MODELQ}"

   113 # CLAMP metric processing of model1

   114 echo "`date` - Running 00.initial.ncl"

   115 ncl $INPUT_TEXT $DIR_SCRIPTS/00.initial.ncl

   116 echo "`date` - Running 01.npp.ncl"

   117 ncl $INPUT_TEXT $DIR_SCRIPTS/01.npp.ncl

   118 echo "`date` - Running 01.lai.ncl"

   119 ncl $INPUT_TEXT $DIR_SCRIPTS/02.lai.ncl

   120 

   121 if ($FILE3 != "") then

   122 echo "`date` - Running 03.co2.ncl"

   123 ncl $INPUT_TEXT $DIR_SCRIPTS/03.co2.ncl

   124 endif

   125 

   126 echo "`date` - Running 04.biomass.ncl"

   127 ncl $INPUT_TEXT $DIR_SCRIPTS/04.biomass.ncl

   128 echo "`date` - Running 06.fluxnet.ncl"

   129 ncl $INPUT_TEXT $DIR_SCRIPTS/06.fluxnet.ncl

   130 echo "`date` - Running 07.beta.ncl"

   131 ncl $INPUT_TEXT $DIR_SCRIPTS/07.beta.ncl

   132 echo "`date` - Running 08.turnover.ncl"

   133 ncl $INPUT_TEXT $DIR_SCRIPTS/08.turnover.ncl

   134 

   135 if ($BGC != "casa") then

   136 echo "`date` - Running 09.carbon_sink.ncl"

   137 ncl $INPUT_TEXT $DIR_SCRIPTS/09.carbon_sink.ncl

   138 else

   139 echo "`date` - Running 09x.carbon_sink.ncl"

   140 ncl $INPUT_TEXT $DIR_SCRIPTS/09x.carbon_sink.ncl

   141 endif

   142 

   143 if ($BGC != "casa") then

   144 echo "`date` - Running 10.fire.ncl"

   145 ncl $INPUT_TEXT $DIR_SCRIPTS/10.fire.ncl

   146 endif

   147 

   148 echo "`date` - Running 11.ameriflux.ncl"

   149 ncl $INPUT_TEXT $DIR_SCRIPTS/11.ameriflux.ncl

   150 echo "`date` - Running 99.final.ncl"

   151 ncl $INPUT_TEXT $DIR_SCRIPTS/99.final.ncl

   152 

   153 #*******************************************************

   154 # user modification-(3)

   155 

   156 # model2

   157 set MODEL2 = CN

   158 set DIR_M  = /lustre/wolf-ddn/scratch/hof/clamp_data/model/

   159 #set FILE1  = i01.10cn_1948-2004_ANN_climo.nc

   160 set FILE1  = i01.10cn_1975-2000_ANN_climo.nc

   161 set FILE2  = i01.10cn_1948-2004_MONS_climo.nc

   162 set FILE3  = f02.03cn_1901-1925_MONS_climo_atm.nc

   163 set FILE4  = i01.06cn_1980-2004_ANN_climo.nc

   164 set FILE5  = i01.10cn_1990-2004_ANN_climo.nc

   165 set FILE6  = i01.07cn_1990-2004_ANN_climo.nc

   166 set FILE7  = i01.10cn_Fire_C_1979-2004_monthly.nc

   167 set FILE8  = i01.10cn_ameriflux_1990-2004_monthly.nc

   168 set FILE9  = i01.10cn_2000-2004_ANN_climo.nc

   169 set FILE10 = i01.10cn_2000-2004_MONS_climo.nc

   170 set FILE11 = i01.10cn_2000_ANN_climo.nc

   171 set GRID   = T42

   172 set BGC    = cn

   173 set ENERGY = new

   174 

   175 #*******************************************************

   176 # create model2 directory by copying templates

   177 if ($FILE3 != "") then

   178    set TEMPLATE1 = template_1-model

   179 else

   180    set TEMPLATE1 = template_1-model_noCO2

   181 endif

   182 #cp -r $DIR_SCRIPTS/$TEMPLATE1 $MODEL2

   183 rm -rf $MODEL2

   184 tar xvpf $DIR_SCRIPTS/$TEMPLATE1.tar

   185 mv $TEMPLATE1 $MODEL2

   186 

   187 # add quote, to be usesd in INPUT_TEXT

   188 set MODELQ  = \"$MODEL2\"

   189 set DIRMQ   = \"$DIR_M\"

   190 set F1  = \"$FILE1\"

   191 set F2  = \"$FILE2\"

   192 set F3  = \"$FILE3\"

   193 set F4  = \"$FILE4\"

   194 set F5  = \"$FILE5\"

   195 set F6  = \"$FILE6\"

   196 set F7  = \"$FILE7\"

   197 set F8  = \"$FILE8\"

   198 set F9  = \"$FILE9\"

   199 set F10 = \"$FILE10\"

   200 set F11 = \"$FILE11\"

   201 set GRIDQ   = \"$GRID\"

   202 set BGCQ    = \"$BGC\"

   203 set ENERGYQ = \"$ENERGY\"

   204 

   205 set COMPAREQ = \"$COMPARE\"

   206 set MODELN   = \"model2\"

   207 

   208 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"

   209  

   210 echo "`date` - Running diagnostics for model ${MODELQ}"

   211 # CLAMP metric processing of model2

   212 echo "`date` - Running 00.inital.ncl"

   213 ncl $INPUT_TEXT $DIR_SCRIPTS/00.initial.ncl

   214 echo "`date` - Running 01.npp.ncl"

   215 ncl $INPUT_TEXT $DIR_SCRIPTS/01.npp.ncl

   216 echo "`date` - Running 02.lai.ncl"

   217 ncl $INPUT_TEXT $DIR_SCRIPTS/02.lai.ncl

   218 

   219 if ($FILE3 != "") then

   220 echo "`date` - Running 03.co2.ncl"

   221 ncl $INPUT_TEXT $DIR_SCRIPTS/03.co2.ncl

   222 endif

   223 

   224 echo "`date` - Running 04.biomass.ncl"

   225 ncl $INPUT_TEXT $DIR_SCRIPTS/04.biomass.ncl

   226 echo "`date` - Running 06.fluxnet.ncl"

   227 ncl $INPUT_TEXT $DIR_SCRIPTS/06.fluxnet.ncl

   228 echo "`date` - Running 07.beta.ncl"

   229 ncl $INPUT_TEXT $DIR_SCRIPTS/07.beta.ncl

   230 echo "`date` - Running 08.turnover.ncl"

   231 ncl $INPUT_TEXT $DIR_SCRIPTS/08.turnover.ncl

   232 

   233 if ($BGC != "casa") then

   234 echo "`date` - Running 09.carbon_sink.ncl"

   235 ncl $INPUT_TEXT $DIR_SCRIPTS/09.carbon_sink.ncl

   236 else

   237 echo "`date` - Running 09x.carbon_sink.ncl"

   238 ncl $INPUT_TEXT $DIR_SCRIPTS/09x.carbon_sink.ncl

   239 endif

   240 

   241 if ($BGC != "casa") then

   242 echo "`date` - Running 10.fire.ncl"

   243 ncl $INPUT_TEXT $DIR_SCRIPTS/10.fire.ncl

   244 endif

   245 

   246 echo "`date` - Running 11.ameriflux.ncl"

   247 ncl $INPUT_TEXT $DIR_SCRIPTS/11.ameriflux.ncl

   248 echo "`date` - Running 99.final.ncl"

   249 ncl $INPUT_TEXT $DIR_SCRIPTS/99.final.ncl

   250 

   251 # create a tar file from the final output

   252 echo "`date` - tarring up final output into all.tar"

   253 tar cf - $MODEL1 $MODEL2 $COMPARE > all.tar

   254
clamp