#! /bin/csh -f

# ====================== COMBINE_v5.3.X Run Script =================== 
# Usage: run.combine.uncoupled.csh >&! combine.log &                                
#
# To report problems or request help with this script/program:     
#             http://www.epa.gov/cmaq    (EPA CMAQ Website)
#             http://www.cmascenter.org  (CMAS Website)
# ===================================================================  

# ==================================================================
#> Runtime Environment Options
# ==================================================================

#> Choose compiler and set up CMAQ environment with correct
#> libraries using config.cmaq. Options: intel | gcc | pgi
 setenv compiler intel

 #> Source the config.cmaq file to set the build environment
 cd ../../..
 source ./config_cmaq.csh
       
#> Set General Parameters for Configuring the Simulation
 set VRSN1     = cb6r3_ae7_aq             #> Code Version
 set VRSN      = v532              #> Code Version
 set PROC      = mpi               #> serial or mpi
 set MECH      = cb6r3_ae7_aq      #> Mechanism ID
 set APPL      = cn36              #> Application Name (e.g. Gridname)
                                                      
#> Define RUNID as any combination of parameters above or others. By default,
#> this information will be collected into this one string, $RUNID, for easy
#> referencing in output binaries and log files as well as in other scripts.
 setenv RUNID  ${VRSN1}_${compilerString}_${APPL}

#> Set the build directory if this was not set above 
#> (this is where the CMAQ executable is located by default).
 if ( ! $?BINDIR ) then
  setenv BINDIR $CMAQ_HOME/POST/combine/scripts/BLD_combine_${VRSN}_${compilerString}
 endif

#> Set the name of the executable.
 setenv EXEC combine_${VRSN}.exe

#> Set location of CMAQ repo.  This will be used to point to the correct species definition files.
 setenv REPO_HOME  ${CMAQ_REPO}

#> Set working, input and output directories
 setenv METDIR     /data4/sunzhixu/outdata/mcip        #> Met Output Directory
 setenv CCTMOUTDIR /data4/sunzhixu/outdata/cctm/OB.nodust
#setenv CCTMOUTDIR /data3/lihl/output_CCTM_cb6r3_ae7_aq_intel_cn04_2018.07.SH.O3     
 #> CCTM Output Directory
 setenv POSTDIR    ${CMAQ_DATA}/POST             #> Location where combine file will be written

  if ( ! -e $POSTDIR ) then
	  mkdir $POSTDIR
  endif



# =====================================================================
#> COMBINE Configuration Options
# =====================================================================

#> Set Start and End Days for looping

 set START_DATE = "2020-08-01"     #> beginning date (July 1, 2016)
 set END_DATE   = "2020-08-30"     #> ending date    (July 14, 2016)
 
#> Set location of species definition files for concentration and deposition species.
 #setenv SPEC_CONC $REPO_HOME/POST/combine/scripts/spec_def_files/SpecDef_${MECH}.txt
 #setenv SPEC_DEP  $REPO_HOME/POST/combine/scripts/spec_def_files/SpecDef_Dep_${MECH}.txt
 #setenv SPEC_CONC /data2/sunzhixu/CMAQv5.3.2/POST/combine/scripts/spec_def_files/O3.NO2.NO.txt
 setenv SPEC_CONC /data2/sunzhixu/CMAQv5.3.2/POST/combine/scripts/spec_def_files/SpecDef_cb6r3_ae7_aq.txt
#> Use GENSPEC switch to generate a new specdef file (does not generate output file).
 setenv GENSPEC N


# =====================================================================
#> Begin Loop Through Simulation Days to Create ACONC File
# =====================================================================

#> Set the species definition file for concentration species.
 setenv SPECIES_DEF $SPEC_CONC
 
#> Loop through all days between START_DAY and END_DAY
 set TODAYG = ${START_DATE}
 set TODAYJ = `date -ud "${START_DATE}" +%Y%j` #> Convert YYYY-MM-DD to YYYYJJJ
 set STOP_DAY = `date -ud "${END_DATE}" +%Y%j` #> Convert YYYY-MM-DD to YYYYJJJ

 while ($TODAYJ <= $STOP_DAY )  #>Compare dates in terms of YYYYJJJ

 
  #> Retrieve Calendar day Information
   set YYYY = `date -ud "${TODAYG}" +%Y`
   set YY = `date -ud "${TODAYG}" +%y`
   set MM = `date -ud "${TODAYG}" +%m`
   set DD = `date -ud "${TODAYG}" +%d`
  #> for files that are indexed with Julian day:
   #  set YYYYJJJ = `date -ud "${TODAYG}" +%Y%j` 

  #> Define name of combine output file to save hourly average concentration.
  #> A new file will be created for each month/year.
   setenv OUTFILE /data4/sunzhixu/outdata/cctm/combine/COMBINE_ACONC_${RUNID}_$YYYY.36.08.OBnodust.test.nc

  #> Define name of input files needed for combine program.
  #> File [1]: CMAQ conc/aconc file
  #> File [2]: MCIP METCRO3D file
  #> File [3]: CMAQ APMDIAG file
  #> File [4]: MCIP METCRO2D file
  # setenv INFILE1 $CCTMOUTDIR/CCTM_PA_1_cb6r3_ae7_aq_intel_cn04_${YYYY}_$YYYY$MM$DD.nc
 # setenv INFILE1 $CCTMOUTDIR/CCTM_IRR_1_cb6r3_ae7_aq_intel_cn04_${YYYY}_$YYYY$MM$DD.nc
  setenv INFILE1 $CCTMOUTDIR/CCTM_ACONC_cb6r3_ae7_aq_intel_cn36_${YYYY}_$YYYY$MM$DD.nc  
  #setenv INFILE1 $CCTMOUTDIR/CCTM_ACONC_${RUNID}_$YYYY$MM$DD.nc
  setenv INFILE2 $METDIR/METCRO3D_cn36_$TODAYJ
  #setenv INFILE2 $METDIR/METCRO3D_$YY$MM$DD.nc
  setenv INFILE3 $CCTMOUTDIR/CCTM_APMDIAG_cb6r3_ae7_aq_intel_cn36_${YYYY}_$YYYY$MM$DD.nc
  #setenv INFILE3 $CCTMOUTDIR/CCTM_APMDIAG_${RUNID}_$YYYY$MM$DD.nc
  setenv INFILE4 $METDIR/METCRO2D_cn36_$TODAYJ
  #setenv INFILE4 $METDIR/METCRO2D_$YY$MM$DD.nc

  #> Executable call:
   ${BINDIR}/${EXEC}

  #> Increment both Gregorian and Julian Days
   set TODAYG = `date -ud "${TODAYG}+1days" +%Y-%m-%d` #> Add a day for tomorrow
   set TODAYJ = `date -ud "${TODAYG}" +%Y%j` #> Convert YYYY-MM-DD to YYYYJJJ

 end #Loop to the next Simulation Day


# =====================================================================
#> Begin Loop Through Simulation Days to Create DEP File
# =====================================================================

#> Set the species definition file for concentration species.
 #setenv SPECIES_DEF $SPEC_DEP
 
#> Loop through all days between START_DAY and END_DAY
 #set TODAYG = ${START_DATE}
 #set TODAYJ = `date -ud "${START_DATE}" +%Y%j` #> Convert YYYY-MM-DD to YYYYJJJ
 #set STOP_DAY = `date -ud "${END_DATE}" +%Y%j` #> Convert YYYY-MM-DD to YYYYJJJ

# while ($TODAYJ <= $STOP_DAY )  #>Compare dates in terms of YYYYJJJ
 
  #> Retrieve Calendar day Information
  #set YYYY = `date -ud "${TODAYG}" +%Y`
  #set YY = `date -ud "${TODAYG}" +%y`
  #set MM = `date -ud "${TODAYG}" +%m`
  #set DD = `date -ud "${TODAYG}" +%d`
  #> for files that are indexed with Julian day:
   
#  set YYYYJJJ = `date -ud "${TODAYG}" +%Y%j` 

  #> Define name of combine output file to save hourly total deposition.
  #> A new file will be created for each month/year.
  #setenv OUTFILE ${POSTDIR}/COMBINE_DEP_${RUNID}_$YYYY$MM

  #> Define name of input files needed for combine program.
  #> File [1]: CMAQ DRYDEP file
  #> File [2]: CMAQ WETDEP file
  #> File [3]: MCIP METCRO2D
  #> File [4]: {empty}
   #setenv INFILE1 $CCTMOUTDIR/CCTM_DRYDEP_${RUNID}_$YYYY$MM$DD.nc
   #setenv INFILE2 $CCTMOUTDIR/CCTM_WETDEP1_${RUNID}_$YYYY$MM$DD.nc
   #setenv INFILE3 $METDIR/METCRO2D_$YY$MM$DD.nc
  #setenv INFILE4

  #> Executable call:
   #${BINDIR}/${EXEC}

  #> Increment both Gregorian and Julian Days
   #set TODAYG = `date -ud "${TODAYG}+1days" +%Y-%m-%d` #> Add a day for tomorrow
   #set TODAYJ = `date -ud "${TODAYG}" +%Y%j` #> Convert YYYY-MM-DD to YYYYJJJ

# end #Loop to the next Simulation Day

 
 exit()