#!/bin/sh

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

  MDL="${HOME}/hysplit4"
  MAP="${MDL}/graphics/arlmap"
  MET="${HOME}/Tutorial/japan"

  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=14
  shr=09
     
  olat=37.4206
  olon=141.0329
  lvl1=1.0
  lvl2=100.0
        
  run=192
  ztop=10000.0
  data=gdas11-22.bin

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

  echo "$syr $smo $sda $shr    " >CONTROL0
  echo "2                      ">>CONTROL0
  echo "$olat $olon $lvl1      ">>CONTROL0
  echo "$olat $olon $lvl2      ">>CONTROL0
  echo "$run                   ">>CONTROL0
  echo "0                      ">>CONTROL0
  echo "$ztop                  ">>CONTROL0
  echo "1                      ">>CONTROL0
  echo "$MET/                  ">>CONTROL0
  echo "$data                  ">>CONTROL0
  echo "1                      ">>CONTROL0
  echo "CPAR                   ">>CONTROL0
  echo "0.0                    ">>CONTROL0
  echo "0.0                    ">>CONTROL0
  echo "00 00 00 00 00         ">>CONTROL0
  echo "2                      ">>CONTROL0
  echo "38.0 140.0             ">>CONTROL0
  echo "0.25 0.25              ">>CONTROL0
  echo "20.0 30.0              ">>CONTROL0
  echo "$OUT/                  ">>CONTROL0
  echo "deposit.bin            ">>CONTROL0
  echo "1                      ">>CONTROL0
  echo "0                      ">>CONTROL0
  echo "00 00 00 00 00         ">>CONTROL0
  echo "00 00 00 00 00         ">>CONTROL0
  echo "00 192 00              ">>CONTROL0
  echo "38.0 140.0             ">>CONTROL0
  echo "0.25 0.25              ">>CONTROL0
  echo "20.0 30.0              ">>CONTROL0
  echo "$OUT/                  ">>CONTROL0
  echo "airconc.bin            ">>CONTROL0
  echo "1                      ">>CONTROL0
  echo "100                    ">>CONTROL0
  echo "00 00 00 00 00         ">>CONTROL0
  echo "00 00 00 00 00         ">>CONTROL0
  echo "00 03 00               ">>CONTROL0
  echo "1                      ">>CONTROL0
  echo "1.0 1.0 1.0            ">>CONTROL0
  echo "0.001 0.0 0.0 0.0 0.0  ">>CONTROL0
  echo "0.0 8.0E-05 8.0E-05    ">>CONTROL0
  echo "0.0                    ">>CONTROL0

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

  echo "&SETUP "                        >SETUP.CFG
  echo "delt = 15, "                   >>SETUP.CFG
  echo "initd = 0, "                   >>SETUP.CFG
  echo "khmax = 12, "                  >>SETUP.CFG
  echo "numpar = -2000, "              >>SETUP.CFG
  echo "maxpar = 1000000, "            >>SETUP.CFG
  echo "efile = '../files/EMITIMES',"  >>SETUP.CFG
  echo "/ "                            >>SETUP.CFG

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

  cp CONTROL0    CONTROL
  echo "0.0   ">>CONTROL

  rm -f deposit.bin airconc.bin

  ${MDL}/exec/hycs_std  

  echo "'TITLE&','Without Resuspension&'"  >LABELS.CFG    
  ${MDL}/exec/concplot -ideposit.bin -oemit_surf1.ps -j${MAP} -c4 -y0.001 -z80 -ukBq -g0:1000 -h37.0:144.0 -v100+50+20+10+5+2+1
  ${MDL}/exec/concplot -iairconc.bin -oemit_surf2.ps -j${MAP} -c4 -x1.0   -z80 -uBq -g0:1000 -h37.0:144.0 -v100+10+1+0.1+0.01+0.001

  ${MDL}/exec/c2datem  -iairconc.bin -oTokai.txt -m../files/JAEA_C137_all.txt -c1000.0
  echo "'TITLE&','Cs-137 time series at Tokai-mura&'" >LABELS.CFG
  echo "'UNITS&',' mBq&'"                            >>LABELS.CFG
  ${MDL}/exec/timeplot -iTokai.txt -s../files/JAEA_C137_all.txt -p2 -y
  mv timeplot.ps emit_surf3.ps
 
# ---------------------------------------------------------

  cp CONTROL0        CONTROL
  echo "1.0E-06   ">>CONTROL

  rm -f deposit.bin airconc.bin

  ${MDL}/exec/hycs_std  

  echo "'TITLE&','With Resuspension&'"  >LABELS.CFG    
  ${MDL}/exec/concplot -ideposit.bin -oemit_surf4.ps -j${MAP} -c4 -y0.001 -z80 -ukBq -g0:1000 -h37.0:144.0 -v100+50+20+10+5+2+1
  ${MDL}/exec/concplot -iairconc.bin -oemit_surf5.ps -j${MAP} -c4 -x1.0   -z80 -uBq -g0:1000 -h37.0:144.0 -v100+10+1+0.1+0.01+0.001

  ${MDL}/exec/c2datem  -iairconc.bin -oTokai.txt -m../files/JAEA_C137_all.txt -c1000.0
  echo "'TITLE&','Cs-137 time series at Tokai-mura&'" >LABELS.CFG
  echo "'UNITS&',' mBq&'"                            >>LABELS.CFG
  ${MDL}/exec/timeplot -iTokai.txt -s../files/JAEA_C137_all.txt -p2 -y
  mv timeplot.ps emit_surf6.ps