#!/bin/sh
#
# This UNIX shell/Linux script uses the command "wget" to download 2012 base case CAMx output:
#  1) CSV files of time-series data at Texas air qaulity monitoring sites
#  2) CSV files of ozone modeling performance statistics by day/region or by episode/site
#  3) PNG image files of surface-layer hourly ozone, NOx, and isopren concentration plots
# 
# Totla download data size: 1 GB per set
#
# Time-series data types:
#  cellvalu - cell value 
#  bilinear - 2x2-cell bilinear interpolated value
#  03x03avg - 3x3-cell average value
#  03x03max - 3x3-cell maximum value
#  03x03min - 3x3-cell minimum value
#  03x03mdn - 3x3-cell median value


# 2012 179B ozone SIP base case release 4:
MOD=camx650_cb6r4h ; SCN=r6f2_r4f.2012_wrf371_p2ma_i2mSNgqsfc0_h

# 2012 HGB/DFW ozone SIP base case release 3:
#MOD=camx650_cb6r4h ; SCN=r6f1_r4f.2012_wrf371_p2ma_i2mSNgqsfc0_g

# 2012 HGB/DFW ozone SIP base case release 2:
#MOD=camx631_cb6r2h ; SCN=r6a_r4a.2012_wrf371_p2ma_i2mSNgqsfc0_f


# loop over five monthly episodes
for EPI in bc12_12may bc12_12jun bc12_12jul bc12_12aug bc12_12sep
do

 case $EPI in
  bc12_12may) DATE=20120501-20120531;;
  bc12_12jun) DATE=20120601-20120630;;
  bc12_12jul) DATE=20120701-20120731;;
  bc12_12aug) DATE=20120801-20120831;;
  bc12_12sep) DATE=20120901-20120930;;
 esac

 # download time-series CSV files
 URL=https://www.tceq.texas.gov/assets/public/implementation/air/am/modeling/tx/mpe/data/$EPI.$SCN/Time_Series/tx_4km/csv
 wget -x -nH --cut-dirs=8 $URL/ts.${DATE}.${MOD}.tx.${EPI}.${SCN}.tx_4km.{ANWC,ARLA,ARPS,ARUT,ATAS,AUDU,BAYP,BMTC,BOER,BP31,BPLG,BPOS,BUHV,BVES,BYTE,BYWC,CALA,CC02,CCHR,CCHS,CCPM,CCTO,CCWT,CECH,CLEB,CLHS,CLTN,CNR2,COGR,COLE,CRBL,CRSA,DALN,DECT,DENN,DHIC,DISH,DNCG,DONA,DRIP,DRPK,ECES,EJPK,EMTL,FAYT,FHRE,FLMS,FORH,FRIC,FWCB,FWMC,GALV,GRAN,GRAP,GRVL,H03H,HALC,HCHV,HCQA,HLAA,HNWA,HOEA,HROC,HSMA,HTCA,HTMS,HUFM,HUTO,HWAA,INEZ,INHB,ISID,ITLY,JEFC,JIHL,KARN,KATP,KAUF,KDTP,KELC,LABD,LABS,LAC1,LAC2,LAC3,LACS,LAGA,LAR1,LAR2,LGTN,LGVW,LKJK,LPSB,LYNF,MACP,MANF,MARY,MDLO,MDTW,MERC,MEYE,MIPK,MRCG,MSTG,NEWF,NLHS,ODEM,OKPK,PAWC,PECV,PGEL,PIPT,PRKP,REDB,RKWL,RSHR,S40S,S42S,S43S,SAWC,SBFP,SEPA,SHLD,SHWH,SKLF,SMSR,SOLC,SOLE,TAFT,TEXC,TOMB,TXCT,TYLA,UHCC,UHJF,UHSL,UHWL,VCTC,VILT,WALV,WAMA,WHOU,WORA,WTFD}.{1,8}h_{03x03avg,03x03max,03x03mdn,03x03min,bilinear,cellvalu}.csv -o wget.ts.log

 # download stats CSV files
 URL=https://www.tceq.texas.gov/assets/public/implementation/air/am/modeling/tx/mpe/data/$EPI.$SCN/Statistics
 wget -x -nH --cut-dirs=8 $URL/stat_o3.$DATE.$MOD.tx.$EPI.$SCN.tx_4km.{arr,bpa,cc,dfw,hgb,netx,san,vic,wac,tx_4km}.{1,8}h_{bilinear,cellvalu}.csv \
                          $URL/site_stat_o3.$DATE.$MOD.tx.$EPI.$SCN.tx_4km.tx_4km.{1,8}h_{bilinear,cellvalu}.csv -o wget.stats.log

 # download surface layer hourly concentration plot image files (not available for release 3 and 4 yet)
 YM=`echo $DATE | cut -c1-6`
 LD=`echo $DATE | cut -c16-`
 for CHEM in o3 nox isop
 do
  URL=https://www.tceq.texas.gov/assets/public/implementation/air/am/modeling/tx/mpe/data/$EPI.$SCN/Tile_Plot/hourly_$CHEM/tx
  wget -x -nH --cut-dirs=8 $URL/gmap_tp_$CHEM.${YM}0{1,2,3,4,5,6,7,8,9}.{0,1}{0,1,2,3,4,5,6,7,8,9}00.$MOD.tx.$EPI.$SCN.tx_4km.tx.1h.png \
                           $URL/gmap_tp_$CHEM.${YM}0{1,2,3,4,5,6,7,8,9}.2{0,1,2,3}00.$MOD.tx.$EPI.$SCN.tx_4km.tx.1h.png \
                           $URL/gmap_tp_$CHEM.${YM}{1,2}{0,1,2,3,4,5,6,7,8,9}.{0,1}{0,1,2,3,4,5,6,7,8,9}00.$MOD.tx.$EPI.$SCN.tx_4km.tx.1h.png \
                           $URL/gmap_tp_$CHEM.${YM}{1,2}{0,1,2,3,4,5,6,7,8,9}.2{0,1,2,3}00.$MOD.tx.$EPI.$SCN.tx_4km.tx.1h.png \
                           $URL/gmap_tp_$CHEM.${YM}{30,$LD}.{0,1}{0,1,2,3,4,5,6,7,8,9}00.$MOD.tx.$EPI.$SCN.tx_4km.tx.1h.png \
                           $URL/gmap_tp_$CHEM.${YM}{30,$LD}.2{0,1,2,3}00.$MOD.tx.$EPI.$SCN.tx_4km.tx.1h.png \
                           -o wget.plot.$CHEM.log
 done

done