#!/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=67
  smo=01
  sda=08
  shr=00
     
  olat=50.87
  olon=-1.218
  lvl1=24.0
        
  run=48
  ztop=10000.0
  data=gdas11-22.bin

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

  echo "&SETUP "               >SETUP.CFG
  echo "numpar = -1000, "     >>SETUP.CFG
  echo "maxpar = 50000, "     >>SETUP.CFG
  echo "/ "                   >>SETUP.CFG

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

for data in ecmwf ncar2; do

  echo "$syr $smo $sda $shr    " >CONTROL
  echo "1                      ">>CONTROL
  echo "$olat $olon $lvl1      ">>CONTROL
  echo "$run                   ">>CONTROL
  echo "0                      ">>CONTROL
  echo "$ztop                  ">>CONTROL
  echo "1                      ">>CONTROL
  echo "$MET/                  ">>CONTROL
  echo "FMD_${data}.bin        ">>CONTROL
  echo "1                      ">>CONTROL
  echo "FM_7                   ">>CONTROL
  echo "3.333E+4               ">>CONTROL
  echo "48.0                   ">>CONTROL
  echo "00 00 00 00 00         ">>CONTROL
  echo "1                      ">>CONTROL
  echo "0.0  0.0               ">>CONTROL
  echo "0.01 0.01              ">>CONTROL
  echo "3.0  3.0               ">>CONTROL
  echo "$OUT/                  ">>CONTROL
  echo "cdump_${data}          ">>CONTROL
  echo "1                      ">>CONTROL
  echo "100                    ">>CONTROL
  echo "00 00 00 00 00         ">>CONTROL
  echo "00 00 00 00 00         ">>CONTROL
  echo "00 24 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

  ${MDL}/exec/hycs_std  

  echo "'TITLE&','Foot and Mouth using ${data} data &'" >LABELS.CFG
  echo "'UNITS&','TCID&'"
  ${MDL}/exec/concplot -icdump_${data} -oFMD_${data} -j${MAP} -h50.8:-1.2 -g0:50 -x24.0 -c4 -v0.001+0.0001+0.00001+0.000001+0.0000001

done