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.csh
author Forrest Hoffman <forrest@climatemodeling.org>
Mon, 26 Jan 2009 22:08:20 -0500
changeset 0 0c6405ab2ff4
child 1 4be95183fbcd
permissions -rwxr-xr-x
Initial commit of C-LAMP Diagnostics from Jeff Lee 
     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 = CN

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

    36 set FILE1  = i01.10cn_1948-2004_ANN_climo.nc

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

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

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

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

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

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

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

    44 set GRID   = T42

    45 set BGC    = cn

    46 set ENERGY = new

    47 

    48 # in the "CLAMP metric processing" section:

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

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

    51 #  e.g. 

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

    53 

    54 # model surface data

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

    56 

    57 # observed data

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

    59 

    60 # directory for scripts, templates and ncl files

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

    62 

    63 #********************************************************

    64 

    65 # create model1 and model1_vs_model2 directory by copying templates

    66 if ($FILE3 != "") then

    67    set TEMPLATE1 = template_1-model

    68    set TEMPLATE2 = template_2-model

    69 else

    70    set TEMPLATE1 = template_1-model_noCO2

    71    set TEMPLATE2 = template_2-model_noCO2

    72 endif

    73 cp -r $DIR_SCRIPTS/$TEMPLATE1 $MODEL1

    74 cp -r $DIR_SCRIPTS/$TEMPLATE2 $COMPARE

    75 

    76 # add quote, to be usesd in INPUT_TEXT

    77 set MODELQ = \"$MODEL1\"

    78 set DIRMQ  = \"$DIR_M\"

    79 set F1  = \"$FILE1\"

    80 set F2  = \"$FILE2\"

    81 set F3  = \"$FILE3\"

    82 set F4  = \"$FILE4\"

    83 set F5  = \"$FILE5\"

    84 set F6  = \"$FILE6\"

    85 set F7  = \"$FILE7\"

    86 set F8  = \"$FILE8\"

    87 set GRIDQ   = \"$GRID\"

    88 set BGCQ    = \"$BGC\"

    89 set ENERGYQ = \"$ENERGY\"

    90 set DIRSQ   = \"$DIR_S\"

    91 set DIROQ   = \"$DIR_O\"

    92 set DIRCQ   = \"$DIR_SCRIPTS\"

    93 

    94 set COMPAREQ = \"$COMPARE\"

    95 set MODELN   = \"model1\"

    96 

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

    98  

    99 # CLAMP metric processing of model1

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

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

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

   103 

   104 if ($FILE3 != "") then

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

   106 endif

   107 

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

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

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

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

   112 

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

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

   115 else

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

   117 endif

   118 

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

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

   121 endif

   122 

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

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

   125 

   126 #*******************************************************

   127 # user modification-(3)

   128 

   129 # model2

   130 set MODEL2 = CASA

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

   132 set FILE1  = i01.10casa_1948-2004_ANN_climo.nc

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

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

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

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

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

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

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

   140 set GRID   = T42

   141 set BGC    = casa

   142 set ENERGY = new

   143 #*******************************************************

   144 # create model2 directory by copying templates

   145 if ($FILE3 != "") then

   146    set TEMPLATE1 = template_1-model

   147 else

   148    set TEMPLATE1 = template_1-model_noCO2

   149 endif

   150 cp -r $DIR_SCRIPTS/$TEMPLATE1 $MODEL2

   151 

   152 # add quote, to be usesd in INPUT_TEXT

   153 set MODELQ  = \"$MODEL2\"

   154 set DIRMQ   = \"$DIR_M\"

   155 set F1  = \"$FILE1\"

   156 set F2  = \"$FILE2\"

   157 set F3  = \"$FILE3\"

   158 set F4  = \"$FILE4\"

   159 set F5  = \"$FILE5\"

   160 set F6  = \"$FILE6\"

   161 set F7  = \"$FILE7\"

   162 set F8  = \"$FILE8\"

   163 set GRIDQ   = \"$GRID\"

   164 set BGCQ    = \"$BGC\"

   165 set ENERGYQ = \"$ENERGY\"

   166 

   167 set COMPAREQ = \"$COMPARE\"

   168 set MODELN   = \"model2\"

   169 

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

   171  

   172 # CLAMP metric processing of model2

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

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

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

   176 

   177 if ($FILE3 != "") then

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

   179 endif

   180 

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

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

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

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

   185 

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

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

   188 else

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

   190 endif

   191 

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

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

   194 endif

   195 

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

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

   198 

   199 # create a tar file from the final output

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

   201
clamp