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_38cn_41cn.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 = i01.38cn_vs_i01.41cn

    29 

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

    31 # user modification-(2)

    32 

    33 # model1 

    34 set MODEL1 = i01.38cn

    35 set DIR_M  = /lustre/wolf-ddn/scratch/hof/cases1/i01.38cn_out/

    36 set FILE1  = i01.38cn_out_ANN_climo.nc

    37 set FILE2  = i01.38cn_out_MONS_climo.nc

    38 set FILE3  = 

    39 set FILE4  = $FILE1

    40 set FILE5  = $FILE1

    41 set FILE6  = $FILE1

    42 set FILE7  = i01.38cn_Fire_C_2151-2175_monthly.nc

    43 set FILE8  = i01.38cn_ameriflux_2151-2175_monthly.nc

    44 set FILE9  = $FILE1

    45 set FILE10 = $FILE2

    46 set GRID   = 1.9

    47 set BGC    = cn

    48 set ENERGY = old

    49 

    50 # in the "CLAMP metric processing" section:

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

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

    53 #  e.g. 

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

    55 

    56 # model surface data

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

    58 

    59 # observed data

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

    61 

    62 # directory for scripts, templates and ncl files

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

    64 

    65 #********************************************************

    66 

    67 # create model1 and model1_vs_model2 directory by copying templates

    68 if ($FILE3 != "") then

    69    set TEMPLATE1 = template_1-model

    70    set TEMPLATE2 = template_2-model

    71 else

    72    set TEMPLATE1 = template_1-model_noCO2

    73    set TEMPLATE2 = template_2-model_noCO2

    74 endif

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

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

    77 rm -rf $MODEL1

    78 rm -rf $COMPARE

    79 tar xvpf $DIR_SCRIPTS/$TEMPLATE1.tar

    80 mv $TEMPLATE1 $MODEL1

    81 tar xvpf $DIR_SCRIPTS/$TEMPLATE2.tar

    82 mv $TEMPLATE2 $COMPARE

    83 

    84 # add quote, to be usesd in INPUT_TEXT

    85 set MODELQ = \"$MODEL1\"

    86 set DIRMQ  = \"$DIR_M\"

    87 set F1  = \"$FILE1\"

    88 set F2  = \"$FILE2\"

    89 set F3  = \"$FILE3\"

    90 set F4  = \"$FILE4\"

    91 set F5  = \"$FILE5\"

    92 set F6  = \"$FILE6\"

    93 set F7  = \"$FILE7\"

    94 set F8  = \"$FILE8\"

    95 set F9  = \"$FILE9\"

    96 set F10 = \"$FILE10\"

    97 set GRIDQ   = \"$GRID\"

    98 set BGCQ    = \"$BGC\"

    99 set ENERGYQ = \"$ENERGY\"

   100 set DIRSQ   = \"$DIR_S\"

   101 set DIROQ   = \"$DIR_O\"

   102 set DIRCQ   = \"$DIR_SCRIPTS\"

   103 

   104 set COMPAREQ = \"$COMPARE\"

   105 set MODELN   = \"model1\"

   106 

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

   108  

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

   110 # CLAMP metric processing of model1

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

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

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

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

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

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

   117 

   118 if ($FILE3 != "") then

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

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

   121 endif

   122 

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

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

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

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

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

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

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

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

   131 

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

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

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

   135 else

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

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

   138 endif

   139 

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

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

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

   143 endif

   144 

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

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

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

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

   149 

   150 #*******************************************************

   151 # user modification-(3)

   152 

   153 # model2

   154 set MODEL2 = i01.41cn

   155 set DIR_M  = /lustre/wolf-ddn/scratch/hof/cases1/i01.41cn_out/

   156 set FILE1  = i01.41cn_out_ANN_climo.nc

   157 set FILE2  = i01.41cn_out_MONS_climo.nc

   158 set FILE3  = 

   159 set FILE4  = $FILE1

   160 set FILE5  = $FILE1

   161 set FILE6  = $FILE1

   162 set FILE7  = i01.41cn_Fire_C_2076-2100_monthly.nc

   163 set FILE8  = i01.41cn_ameriflux_2076-2100_monthly.nc

   164 set FILE9  = $FILE1

   165 set FILE10 = $FILE2

   166 set GRID   = 1.9

   167 set BGC    = cn

   168 set ENERGY = old

   169 #*******************************************************

   170 # create model2 directory by copying templates

   171 if ($FILE3 != "") then

   172    set TEMPLATE1 = template_1-model

   173 else

   174    set TEMPLATE1 = template_1-model_noCO2

   175 endif

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

   177 rm -rf $MODEL2

   178 tar xvpf $DIR_SCRIPTS/$TEMPLATE1.tar

   179 mv $TEMPLATE1 $MODEL2

   180 

   181 # add quote, to be usesd in INPUT_TEXT

   182 set MODELQ  = \"$MODEL2\"

   183 set DIRMQ   = \"$DIR_M\"

   184 set F1  = \"$FILE1\"

   185 set F2  = \"$FILE2\"

   186 set F3  = \"$FILE3\"

   187 set F4  = \"$FILE4\"

   188 set F5  = \"$FILE5\"

   189 set F6  = \"$FILE6\"

   190 set F7  = \"$FILE7\"

   191 set F8  = \"$FILE8\"

   192 set F9  = \"$FILE9\"

   193 set F10 = \"$FILE10\"

   194 set GRIDQ   = \"$GRID\"

   195 set BGCQ    = \"$BGC\"

   196 set ENERGYQ = \"$ENERGY\"

   197 

   198 set COMPAREQ = \"$COMPARE\"

   199 set MODELN   = \"model2\"

   200 

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

   202  

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

   204 # CLAMP metric processing of model2

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

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

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

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

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

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

   211 

   212 if ($FILE3 != "") then

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

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

   215 endif

   216 

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

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

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

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

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

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

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

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

   225 

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

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

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

   229 else

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

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

   232 endif

   233 

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

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

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

   237 endif

   238 

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

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

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

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

   243 

   244 # create a tar file from the final output

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

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

   247
clamp