#!/bin/sh

#------------- Set Directories ------------------------

  WEB=""
  if [ -f /usr/bin/firefox ];then WEB="/usr/bin";fi
  if [ -f /usr/local/bin/firefox ];then WEB="/usr/local/bin";fi

  MDL="${HOME}/hysplit"
  MAP="${MDL}/graphics/arlmap"
  MET="."

  if [ ! -f ASCDATA.CFG ]; then
     echo "-90.0  -180.0"     >ASCDATA.CFG
     echo "1.0     1.0"      >>ASCDATA.CFG
     echo "180     360"      >>ASCDATA.CFG
     echo "2"                >>ASCDATA.CFG
     echo "0.2"              >>ASCDATA.CFG
     echo "'$MDL/bdyfiles/'" >>ASCDATA.CFG
  fi
  echo "### $0 ###"

#--------------------------------------------------------------

  syr=11
  smo=03
  sda=14
  shr=00 
     
  olat=37.4206
  olon=141.0329
  lvl1=1.0
  lvl2=100.0
        
  run=72
  ztop=10000.0
  data=wrf11031400.bin

#------ Create CONTROL file -----------------------------

  echo "$syr $smo $sda $shr    " >CONTROL
  echo "2                      ">>CONTROL
  echo "$olat $olon $lvl1      ">>CONTROL
  echo "$olat $olon $lvl2      ">>CONTROL
  echo "$run                   ">>CONTROL
  echo "0                      ">>CONTROL
  echo "$ztop                  ">>CONTROL
  echo "1                      ">>CONTROL
  echo "$MET/                  ">>CONTROL
  echo "$data                  ">>CONTROL
  echo "1                      ">>CONTROL
  echo "CPAR                   ">>CONTROL
  echo "1.0                    ">>CONTROL
  echo "3.0                    ">>CONTROL
  echo "00 00 00 00 00         ">>CONTROL
  echo "1                      ">>CONTROL
  echo "38.0 140.0             ">>CONTROL
  echo "0.05 0.05              ">>CONTROL
  echo "20.0 30.0              ">>CONTROL
  echo "./                     ">>CONTROL
  echo "tcm.bin                ">>CONTROL
  echo "2                      ">>CONTROL
  echo "0 100                  ">>CONTROL
  echo "00 00 00 00 00         ">>CONTROL
  echo "00 00 00 00 00         ">>CONTROL
  echo "00 03 00               ">>CONTROL
  echo "1                      ">>CONTROL
  echo "1.0 1.0 1.0            ">>CONTROL
  echo "0.001 0.0 0.0 0.0 0.0  ">>CONTROL
  echo "0.0 8.0E-05 8.0E-05    ">>CONTROL
  echo "0.0                    ">>CONTROL
  echo "0.0                    ">>CONTROL

#--------- Create Namelist File ---------------------------

  echo "&SETUP "               >SETUP.CFG
  echo "qcycle = 3, "         >>SETUP.CFG
  echo "ichem = 10, "         >>SETUP.CFG
  echo "delt = 5.0, "         >>SETUP.CFG
  echo "khmax = 24, "         >>SETUP.CFG
  echo "numpar = -2500, "     >>SETUP.CFG
  echo "maxpar = 300000, "    >>SETUP.CFG
  echo "/ "                   >>SETUP.CFG

#------- Run HYSPLIT --------------------------------------

rm -f tcm.bin
${MDL}/exec/hycs_std   

echo "Continue with graphics (y/n)?"
read cmd
if [ "$cmd" == "n" ]; then exit; fi

# -------- Apply emission factors and decay ------

${MDL}/exec/tcmsum -itcm.bin -ofdnpp -t031106 -h11025.8 -c1 -scfactors.txt -pC137

# -------- Time Series at JAEA ---------

${MDL}/exec/c2datem  -ifdnpp.bin -oTokai.txt -mJAEA_C137.txt -c1000.0 -z2

echo "'TITLE&','Cs-137 time series at Tokai-mura&'" >LABELS.CFG
echo "'UNITS&',' mBq&'"                            >>LABELS.CFG
${MDL}/exec/timeplot -iTokai.txt -sJAEA_C137.txt +g1 -p2 -y
mv timeplot.html plot_time.html

# -------- Deposition Plot -------------

echo "'TITLE&','Cs-137 deposition &'" >LABELS.CFG
echo "'UNITS&',' kBq&'"              >>LABELS.CFG
${MDL}/exec/concplot -ifdnpp.bin -oplot_dept -j${MAP} +g1 -z95 -h37.0:140.0 -g0:100 \
      -y0.001 -ukBq -b0 -t0 -k2 -c4 -v1000+500+200+100+50+20+10 -r3

# -------- Simple Dose Plot -------------

echo "'TITLE&','Cs-137 dose rate&'" >LABELS.CFG
echo "'UNITS&',' uR&'"             >>LABELS.CFG
echo "'VOLUM&','/h&'"              >>LABELS.CFG
${MDL}/exec/concplot -ifdnpp.bin -oplot_dose -j${MAP} +g1 -z95 -h37.0:140.0 -g0:100 \
    -k2 -x3.340E-05 -y1.08E-06 -c4 -v0.1+0.05+0.02+0.01+0.005+0.002+0.001+0.0005+0.0002+0.0001 -r2

# -------- Display Graphics -------------

if [[ "$OSTYPE" == "darwin"* ]]; then
   open plot_time.html   
   open plot_dept.html 
   open plot_dose.html 
else
   ${WEB}/firefox plot_time.html & 
   ${WEB}/firefox plot_dept.html &  
   ${WEB}/firefox plot_dose.html &  
fi