#!/bin/sh

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

  MDL="${HOME}/hysplit4"
  MAP="${MDL}/graphics/arlmap"
  MET="${HOME}/TutorialX/files"

  OUT="."
  cd $OUT

  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=15
  shr=12 
     
  olat=37.4206
  olon=141.0329
  lvl1=1.0
  lvl2=100.0
        
  run=24
  ztop=10000.0
  data=wrf11031400.bin

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

  echo "$MET/             " >input.txt
  echo "$data             ">>input.txt
  echo "01                ">>input.txt
  echo "TPP1 01 1000.0    ">>input.txt
  echo "$olat $olon       ">>input.txt
  echo "1                 ">>input.txt
  echo "rain.txt          ">>input.txt
  echo "1                 ">>input.txt
  echo "9999              ">>input.txt
  ${MDL}/exec/xtrct_stn -i  <input.txt

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

  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 "8.33333                ">>CONTROL
  echo "12.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 "$OUT/                  ">>CONTROL
  echo "fdnpp.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 12 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 "&SETUP "                        >SETUP.CFG
  echo "delt = 5.0, "                  >>SETUP.CFG
  echo "initd = 0, "                   >>SETUP.CFG
  echo "khmax = 24, "                  >>SETUP.CFG
  echo "numpar = -2500, "              >>SETUP.CFG
  echo "maxpar = 300000, "             >>SETUP.CFG
  echo "krand = 1, "                   >>SETUP.CFG
  echo "vscales = 5.0, "               >>SETUP.CFG
  echo "vscaleu = 200.0, "             >>SETUP.CFG
  echo "kmix0 = 150, "                 >>SETUP.CFG
  echo "rhb = 100, "                   >>SETUP.CFG
  echo "rht = 80, "                    >>SETUP.CFG
  echo "/ "                            >>SETUP.CFG

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

  cp SETUP.CFG   SETUP.ver0
  cp CONTROL   CONTROL.ver0
  echo "0.0 3.5E+05 5.0E-05    ">>CONTROL.ver0
  echo "0.0                    ">>CONTROL.ver0
  echo "0.0                    ">>CONTROL.ver0
  ${MDL}/exec/hycs_std ver0
  mv fdnpp.bin fdnpp0.bin
  ${MDL}/exec/concplot -ifdnpp0.bin -oplot_dept0 -j${MAP} -z95 -h37.0:140.0 -g0:500 \
                       -y1.0E+10 -b0 -t0 -k2 -c4 -v100+50+20+10+5+2+1 

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

  cp SETUP.CFG   SETUP.ver1
  cp CONTROL   CONTROL.ver1
  echo "0.0 4.0E+04 5.0E-05    ">>CONTROL.ver1
  echo "0.0                    ">>CONTROL.ver1
  echo "0.0                    ">>CONTROL.ver1
  ${MDL}/exec/hycs_std ver1
  mv fdnpp.bin fdnpp1.bin
  ${MDL}/exec/concplot -ifdnpp1.bin -oplot_dept1 -j${MAP} -z95 -h37.0:140.0 -g0:500 \
                       -y1.0E+10 -b0 -t0 -k2 -c4 -v100+50+20+10+5+2+1 

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

  cp SETUP.CFG   SETUP.ver2
  cp CONTROL   CONTROL.ver2
  echo "0.0 4.0E+04 5.0E-05    ">>CONTROL.ver2
  echo "0.0                    ">>CONTROL.ver2
  echo "0.0                    ">>CONTROL.ver2
  ${MDL}/exec/hycs_std ver2
  mv fdnpp.bin fdnpp2.bin
  ${MDL}/exec/concplot -ifdnpp2.bin -oplot_dept2 -j${MAP} -z95 -h37.0:140.0 -g0:500 \
                       -y1.0E+10 -b0 -t0 -k2 -c4 -v100+50+20+10+5+2+1 

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

  cp SETUP.CFG   SETUP.ver3
  cp CONTROL   CONTROL.ver3
  echo "0.0 1.0E+06 1.0E-06    ">>CONTROL.ver3
  echo "0.0                    ">>CONTROL.ver3
  echo "0.0                    ">>CONTROL.ver3
  ${MDL}/exec/hycs_std ver3
  mv fdnpp.bin fdnpp3.bin
  ${MDL}/exec/concplot -ifdnpp3.bin -oplot_dept3 -j${MAP} -z95 -h37.0:140.0 -g0:500 \
                       -y1.0E+10 -b0 -t0 -k2 -c4 -v100+50+20+10+5+2+1 

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

  cp SETUP.CFG   SETUP.ver4
  cp CONTROL   CONTROL.ver4
  echo "0.0 8.0E-05 8.0E-05    ">>CONTROL.ver4
  echo "0.0                    ">>CONTROL.ver4
  echo "0.0                    ">>CONTROL.ver4
  ${MDL}/exec/hycs_std ver4
  mv fdnpp.bin fdnpp4.bin
  ${MDL}/exec/concplot -ifdnpp4.bin -oplot_dept4 -j${MAP} -z95 -h37.0:140.0 -g0:500 \
                       -y1.0E+10 -b0 -t0 -k2 -c4 -v100+50+20+10+5+2+1