#!/bin/sh

  MDL="${HOME}/hysplit"
  OUT="${MDL}/working"
  MET="${HOME}/Tutorial/captex"
  cd $OUT

  echo "### $0 ###"

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

#  -D[Input metdata directory name with ending /]
#  -F[input metdata file name]
#  -Y[Map center latitude]
#  -X[Map center longitude]
#  -R[Map radius (deg lat)]
#  -V[Variable name to contour (e.g. TEMP)]
#  -L[Level of variable (sfc=1)]
#  -O[Output time offset (hrs)]
#  -T[Output time interval (hrs)]
#  -C[Color (1/3) or B&W (0/2); 0/1=lines 2/3=nolines]
#  -G[Graphics map file (arlmap) or shapefiles.txt]
#  -M[Maximum contour value (Auto=-1.0)]
#  -I[Increment between contours (Auto=-1.0)]
#  -A[Arcview text output]
  
  
  ${MDL}/exec/contour -g${MDL}/graphics/arlmap -d$MET/ -fcaptex2_gblr.bin -y40.0 -x-80.0 -vHGTS -l4 -o48 -r25.0 -c1
  if [[ "$OSTYPE" == "darwin"* ]]; then
     open contour.ps  
  else
     gs contour.ps  
  fi