forrest@1: #!/bin/csh -f
forrest@1: #-------------------------------------------------------------------
forrest@1: # note: user modifies ONLY the "user modification" section 
forrest@1: #
forrest@1: #       COMPARE: model1 vs model2
forrest@1: #       MODELn : model name
forrest@1: #       DIR_M  : directory of model data
forrest@1: #       DIR_O  : directory of observed data
forrest@1: #       DIR_S  : directory of model surface data
forrest@1: #       DIR_SCRIPTS  : directory of run scripts
forrest@1: #       FILE1  : time_mean  climatology from CLM diagnostic package
forrest@1: #       FILE2  : 12-monthly climatology from CLM diagnostic package
forrest@1: #       FILE3  : 12-monthly climatology from ATM diagnostic package
forrest@1: #                leave it blank, if no ATM file:
forrest@1: #                set FILE3 =
forrest@1: #       FILE7 : timeseries file generated from .... 
forrest@1: #       FILE8 : timeseries file generated from .... 
forrest@1: #       GRID   : T31, T42, or 1.9
forrest@1: #       BGC    : cn or casa 
forrest@1: #       ENERGY : new or old (fields in model data)
forrest@1: #-------------------------------------------------------------------
forrest@1: 
forrest@1: #*******************************************************
forrest@1: # user modification-(1)
forrest@1: 
forrest@1: # directory name of model comparison
forrest@1: #et COMPARE = b30.061n_vs_b30.061m
forrest@1: set COMPARE = i01.37cn_cont_vs_i01.43cn
forrest@1: 
forrest@1: #*******************************************************
forrest@1: # user modification-(2)
forrest@1: 
forrest@1: # model1 
forrest@1: set MODEL1 = i01.37cn_cont
forrest@1: set DIR_M  = /lustre/wolf-ddn/scratch/hof/cases/i01.37cn_cont_out/
forrest@1: set FILE1  = i01.37cn_cont_out_ANN_climo.nc
forrest@1: set FILE2  = i01.37cn_cont_out_MONS_climo.nc
forrest@1: set FILE3  = 
forrest@1: set FILE4  = $FILE1
forrest@1: set FILE5  = $FILE1
forrest@1: set FILE6  = $FILE1
forrest@1: set FILE7  = i01.37cn_cont_Fire_C_2976-3000_monthly.nc
forrest@1: set FILE8  = i01.37cn_cont_ameriflux_2976-3000_monthly.nc
forrest@1: set FILE9  = $FILE1
forrest@1: set FILE10 = $FILE2
forrest@1: set GRID   = 1.9
forrest@1: set BGC    = cn
forrest@1: set ENERGY = old
forrest@1: 
forrest@1: # in the "CLAMP metric processing" section:
forrest@1: #  only 00.initial.ncl and 99.final.ncl are required,
forrest@1: #  user can comment out any one or more of the other ncl scripts,
forrest@1: #  e.g. 
forrest@1: #  #ncl $INPUT_TEXT $DIR_SCRIPTS/10.fire.ncl
forrest@1: 
forrest@1: # model surface data
forrest@1: set DIR_S  = /lustre/wolf-ddn/scratch/hof/clamp_data/surface_model/
forrest@1: 
forrest@1: # observed data
forrest@1: set DIR_O  = /lustre/wolf-ddn/scratch/hof/clamp_data/observed/
forrest@1: 
forrest@1: # directory for scripts, templates and ncl files
forrest@1: set DIR_SCRIPTS = /lustre/wolf-ddn/scratch/hof/clamp/all/
forrest@1: 
forrest@1: #********************************************************
forrest@1: 
forrest@1: # create model1 and model1_vs_model2 directory by copying templates
forrest@1: if ($FILE3 != "") then
forrest@1:    set TEMPLATE1 = template_1-model
forrest@1:    set TEMPLATE2 = template_2-model
forrest@1: else
forrest@1:    set TEMPLATE1 = template_1-model_noCO2
forrest@1:    set TEMPLATE2 = template_2-model_noCO2
forrest@1: endif
forrest@1: #cp -r $DIR_SCRIPTS/$TEMPLATE1 $MODEL1
forrest@1: #cp -r $DIR_SCRIPTS/$TEMPLATE2 $COMPARE
forrest@1: rm -rf $MODEL1
forrest@1: rm -rf $COMPARE
forrest@1: tar xvpf $DIR_SCRIPTS/$TEMPLATE1.tar
forrest@1: mv $TEMPLATE1 $MODEL1
forrest@1: tar xvpf $DIR_SCRIPTS/$TEMPLATE2.tar
forrest@1: mv $TEMPLATE2 $COMPARE
forrest@1: 
forrest@1: # add quote, to be usesd in INPUT_TEXT
forrest@1: set MODELQ = \"$MODEL1\"
forrest@1: set DIRMQ  = \"$DIR_M\"
forrest@1: set F1  = \"$FILE1\"
forrest@1: set F2  = \"$FILE2\"
forrest@1: set F3  = \"$FILE3\"
forrest@1: set F4  = \"$FILE4\"
forrest@1: set F5  = \"$FILE5\"
forrest@1: set F6  = \"$FILE6\"
forrest@1: set F7  = \"$FILE7\"
forrest@1: set F8  = \"$FILE8\"
forrest@1: set F9  = \"$FILE9\"
forrest@1: set F10 = \"$FILE10\"
forrest@1: set GRIDQ   = \"$GRID\"
forrest@1: set BGCQ    = \"$BGC\"
forrest@1: set ENERGYQ = \"$ENERGY\"
forrest@1: set DIRSQ   = \"$DIR_S\"
forrest@1: set DIROQ   = \"$DIR_O\"
forrest@1: set DIRCQ   = \"$DIR_SCRIPTS\"
forrest@1: 
forrest@1: set COMPAREQ = \"$COMPARE\"
forrest@1: set MODELN   = \"model1\"
forrest@1: 
forrest@1: 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"
forrest@1:  
forrest@1: echo "`date` - Running diagnostics for model ${MODELQ}"
forrest@1: # CLAMP metric processing of model1
forrest@1: echo "`date` - Running 00.initial.ncl"
forrest@1: ncl $INPUT_TEXT $DIR_SCRIPTS/00.initial.ncl
forrest@1: echo "`date` - Running 01.npp.ncl"
forrest@1: ncl $INPUT_TEXT $DIR_SCRIPTS/01.npp.ncl
forrest@1: echo "`date` - Running 01.lai.ncl"
forrest@1: ncl $INPUT_TEXT $DIR_SCRIPTS/02.lai.ncl
forrest@1: 
forrest@1: if ($FILE3 != "") then
forrest@1: echo "`date` - Running 03.co2.ncl"
forrest@1: ncl $INPUT_TEXT $DIR_SCRIPTS/03.co2.ncl
forrest@1: endif
forrest@1: 
forrest@1: echo "`date` - Running 04.biomass.ncl"
forrest@1: ncl $INPUT_TEXT $DIR_SCRIPTS/04.biomass.ncl
forrest@1: echo "`date` - Running 06.fluxnet.ncl"
forrest@1: ncl $INPUT_TEXT $DIR_SCRIPTS/06.fluxnet.ncl
forrest@1: echo "`date` - Running 07.beta.ncl"
forrest@1: ncl $INPUT_TEXT $DIR_SCRIPTS/07.beta.ncl
forrest@1: echo "`date` - Running 08.turnover.ncl"
forrest@1: ncl $INPUT_TEXT $DIR_SCRIPTS/08.turnover.ncl
forrest@1: 
forrest@1: if ($BGC != "casa") then
forrest@1: echo "`date` - Running 09.carbon_sink.ncl"
forrest@1: ncl $INPUT_TEXT $DIR_SCRIPTS/09.carbon_sink.ncl
forrest@1: else
forrest@1: echo "`date` - Running 09x.carbon_sink.ncl"
forrest@1: ncl $INPUT_TEXT $DIR_SCRIPTS/09x.carbon_sink.ncl
forrest@1: endif
forrest@1: 
forrest@1: if ($BGC != "casa") then
forrest@1: echo "`date` - Running 10.fire.ncl"
forrest@1: ncl $INPUT_TEXT $DIR_SCRIPTS/10.fire.ncl
forrest@1: endif
forrest@1: 
forrest@1: echo "`date` - Running 11.ameriflux.ncl"
forrest@1: ncl $INPUT_TEXT $DIR_SCRIPTS/11.ameriflux.ncl
forrest@1: echo "`date` - Running 99.final.ncl"
forrest@1: ncl $INPUT_TEXT $DIR_SCRIPTS/99.final.ncl
forrest@1: 
forrest@1: #*******************************************************
forrest@1: # user modification-(3)
forrest@1: 
forrest@1: # model2
forrest@1: set MODEL2 = i01.43cn
forrest@1: set DIR_M  = /lustre/wolf-ddn/scratch/hof/cases1/i01.43cn_out/
forrest@1: set FILE1  = i01.43cn_out_ANN_climo.nc
forrest@1: set FILE2  = i01.43cn_out_MONS_climo.nc
forrest@1: set FILE3  = 
forrest@1: set FILE4  = $FILE1
forrest@1: set FILE5  = $FILE1
forrest@1: set FILE6  = $FILE1
forrest@1: set FILE7  = i01.43cn_Fire_C_2876-2900_monthly.nc
forrest@1: set FILE8  = i01.43cn_ameriflux_2876-2900_monthly.nc
forrest@1: set FILE9  = $FILE1
forrest@1: set FILE10 = $FILE2
forrest@1: set GRID   = 1.9
forrest@1: set BGC    = cn
forrest@1: set ENERGY = old
forrest@1: #*******************************************************
forrest@1: # create model2 directory by copying templates
forrest@1: if ($FILE3 != "") then
forrest@1:    set TEMPLATE1 = template_1-model
forrest@1: else
forrest@1:    set TEMPLATE1 = template_1-model_noCO2
forrest@1: endif
forrest@1: #cp -r $DIR_SCRIPTS/$TEMPLATE1 $MODEL2
forrest@1: rm -rf $MODEL2
forrest@1: tar xvpf $DIR_SCRIPTS/$TEMPLATE1.tar
forrest@1: mv $TEMPLATE1 $MODEL2
forrest@1: 
forrest@1: # add quote, to be usesd in INPUT_TEXT
forrest@1: set MODELQ  = \"$MODEL2\"
forrest@1: set DIRMQ   = \"$DIR_M\"
forrest@1: set F1  = \"$FILE1\"
forrest@1: set F2  = \"$FILE2\"
forrest@1: set F3  = \"$FILE3\"
forrest@1: set F4  = \"$FILE4\"
forrest@1: set F5  = \"$FILE5\"
forrest@1: set F6  = \"$FILE6\"
forrest@1: set F7  = \"$FILE7\"
forrest@1: set F8  = \"$FILE8\"
forrest@1: set F9  = \"$FILE9\"
forrest@1: set F10 = \"$FILE10\"
forrest@1: set GRIDQ   = \"$GRID\"
forrest@1: set BGCQ    = \"$BGC\"
forrest@1: set ENERGYQ = \"$ENERGY\"
forrest@1: 
forrest@1: set COMPAREQ = \"$COMPARE\"
forrest@1: set MODELN   = \"model2\"
forrest@1: 
forrest@1: 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"
forrest@1:  
forrest@1: echo "`date` - Running diagnostics for model ${MODELQ}"
forrest@1: # CLAMP metric processing of model2
forrest@1: echo "`date` - Running 00.inital.ncl"
forrest@1: ncl $INPUT_TEXT $DIR_SCRIPTS/00.initial.ncl
forrest@1: echo "`date` - Running 01.npp.ncl"
forrest@1: ncl $INPUT_TEXT $DIR_SCRIPTS/01.npp.ncl
forrest@1: echo "`date` - Running 02.lai.ncl"
forrest@1: ncl $INPUT_TEXT $DIR_SCRIPTS/02.lai.ncl
forrest@1: 
forrest@1: if ($FILE3 != "") then
forrest@1: echo "`date` - Running 03.co2.ncl"
forrest@1: ncl $INPUT_TEXT $DIR_SCRIPTS/03.co2.ncl
forrest@1: endif
forrest@1: 
forrest@1: echo "`date` - Running 04.biomass.ncl"
forrest@1: ncl $INPUT_TEXT $DIR_SCRIPTS/04.biomass.ncl
forrest@1: echo "`date` - Running 06.fluxnet.ncl"
forrest@1: ncl $INPUT_TEXT $DIR_SCRIPTS/06.fluxnet.ncl
forrest@1: echo "`date` - Running 07.beta.ncl"
forrest@1: ncl $INPUT_TEXT $DIR_SCRIPTS/07.beta.ncl
forrest@1: echo "`date` - Running 08.turnover.ncl"
forrest@1: ncl $INPUT_TEXT $DIR_SCRIPTS/08.turnover.ncl
forrest@1: 
forrest@1: if ($BGC != "casa") then
forrest@1: echo "`date` - Running 09.carbon_sink.ncl"
forrest@1: ncl $INPUT_TEXT $DIR_SCRIPTS/09.carbon_sink.ncl
forrest@1: else
forrest@1: echo "`date` - Running 09x.carbon_sink.ncl"
forrest@1: ncl $INPUT_TEXT $DIR_SCRIPTS/09x.carbon_sink.ncl
forrest@1: endif
forrest@1: 
forrest@1: if ($BGC != "casa") then
forrest@1: echo "`date` - Running 10.fire.ncl"
forrest@1: ncl $INPUT_TEXT $DIR_SCRIPTS/10.fire.ncl
forrest@1: endif
forrest@1: 
forrest@1: echo "`date` - Running 11.ameriflux.ncl"
forrest@1: ncl $INPUT_TEXT $DIR_SCRIPTS/11.ameriflux.ncl
forrest@1: echo "`date` - Running 99.final.ncl"
forrest@1: ncl $INPUT_TEXT $DIR_SCRIPTS/99.final.ncl
forrest@1: 
forrest@1: # create a tar file from the final output
forrest@1: echo "`date` - tarring up final output into all.tar"
forrest@1: tar cf - $MODEL1 $MODEL2 $COMPARE > all.tar
forrest@1: