#!/bin/csh

# --------------------------------
# Create Basemap Topgraphy MMI Map 
# --------------------------------

# Set Parameters
# --------------
  set NRUN    = 11C
  set NAME    = MMI.$NRUN     
  set MMIGRD  = base$NRUN.grd
  set REGION  = -123.1/-120.9/037.3/038.7
  set REGION  = -126.5/-119.0/035.0/042.0
  set COORD0  = "-119.01 36.77 126.0"       # lower-left of fd grid
  set SCALE   = 1.0
  set JP      = M05.00
  set JP      = M10.00
  set JP      = M07.50
  set AXIS    = a1.0f1.0/a1.0f1.0WSen
  set MBAR    = fx1.00/0.350/36.0/200.
  set MBAR    = fx2.00/0.700/36.0/200.
  set MBAR    = fx1.50/0.525/36.0/200.
  set TOPOI   = fdtopo.i.grd
  set SPAR    = "-D2.50/-0.8/2.0/0.150h"
  set SPAR    = "-D5.00/-1.6/4.0/0.300h"
  set SPAR    = "-D3.75/-1.2/3.0/0.225h"
  set PPAR    = "-R -JM -B"
  set CPAR    = "-Dfull -Ia/1/100/100/255       -N2/3 -S200/200/255"
  set CPAR    = "-Dfull -Ia/2/100/100/255       -N2/3 -S200/200/255"
  gmtset        COLOR_MODEL HSV
  gmtset        COLOR_MODEL RGB
  gmtset        LABEL_FONT_SIZE 10p
 #gmtset        LABEL_FONT_SIZE 20p
  gmtset        HEADER_FONT_SIZE 24p
 #gmtset        HEADER_FONT_SIZE 48p
  gmtset        PAPER_MEDIA legal
  gmtset        PAPER_MEDIA letter
  gmtset        PAPER_MEDIA archD

# Get MMI Values
# --------------
  if (0) then
  runmmi
  endif

# Create MMI Grd File
# -------------------
  if (0) then
  a2e3d 631 321 < $NAME > tmp.e3d
  imagexyz 1.0 $SCALE 1 < tmp.e3d > tmp.a
  getllv1 $COORD0 < tmp.a > tmp.llv
 #surface tmp.llv -G$MMIGRD -R$REGION -I1.000m -S0.0 -N250 -T0.35i -T0.35b -C
  blockmean tmp.llv -R$REGION  -I0.50m > tmp.xyz
  xyz2grd tmp.xyz -G$MMIGRD    -I1.00m -R$REGION  -N0.0
  endif

# Get Topography
# --------------
  if (0) then
  grdraster 2 -R$REGION  -Gustopo.grd -I0.500m
  grd2xyz ustopo.grd > ustopo.xyz
  blockmean ustopo.xyz -R$REGION   -I1.000m > fd.xyz
  surface fd.xyz -Gfd.grd -R$REGION -I1.000m -S0.0 -N250 -T0.35i -T0.35b -C
  grdgradient fd.grd -A0 -Gfdtopo.i.grd -Nt -M -V >& /dev/null
  endif

# Create Auxiliary Files
# ----------------------
  if (1) then
 #echo "-119.04297   36.80634"  > fdbox.xy 
  echo "-119.01000   36.77000"  > fdbox.xy 
  echo "-123.29400   41.34854" >> fdbox.xy 
  echo "-126.20372   39.65889" >> fdbox.xy 
  echo "-121.88189   35.08035" >> fdbox.xy 
  echo "-119.01000   36.77000" >> fdbox.xy 
  echo "-123.10000   37.30000"  > delta.xy 
  echo "-120.90000   37.30000" >> delta.xy 
  echo "-120.90000   38.80000" >> delta.xy 
  echo "-123.10000   38.80000" >> delta.xy 
  echo "-123.10000   37.30000" >> delta.xy 
  echo "-122.55      37.75   "  >   epi.xy
  endif

# Plot Map
# --------
  if (1) then
  psbasemap -R$REGION  -J$JP -B$AXIS -X2.0 -Y3.0            -P    -K >! base.ps
  grdimage     $MMIGRD -R -JM -Cmmi10.cpt     -I$TOPOI         -O -K >> base.ps
  pscoast $PPAR $CPAR                                          -O -K >> base.ps
  psxy    jennings.xy $PPAR -M -W04/255/000/000                -O -K >> base.ps
 #psxy    jennings.xy $PPAR -M -W08/255/000/000                -O -K >> base.ps
 #psxy    ca_roads.xy $PPAR -M -W0/200                         -O -K >> base.ps
 #psxy    ca_roads.xy $PPAR -M -W4/200                         -O -K >> base.ps
  psxy    fdbox.xy    $PPAR -M -A -W2/0t20_5:20                -O -K >> base.ps
 #psxy    fdbox.xy    $PPAR -M -A -W4/0t20_5:20                -O -K >> base.ps
 #psxy    fdbox.xy    $PPAR -M -A -W2/0                        -O -K >> base.ps
  psxy    delta.xy    $PPAR -M -A -W06/000/000/255             -O -K >> base.ps
 #psxy    delta.xy    $PPAR -M -A -W10/000/000/255             -O -K >> base.ps
 #psxy    epi.xy      $PPAR -Sa0.2 -G255/0/0 -W4/0             -O -K >> base.ps
  psxy    epi.xy      $PPAR -Sa0.3 -G255/0/0 -W4/0             -O -K >> base.ps
 #psxy    epi.xy      $PPAR -Sa0.4 -G255/0/0 -W6/0             -O -K >> base.ps
 #pstext  c1906.txt   $PPAR -Dj-1.05/0.00                      -O -K >> base.ps
 #pstext  c1906.txt   $PPAR -Dj-2.10/0.00                      -O -K >> base.ps
 #pstext  delta.txt   $PPAR -Dj+0.00/0.05                      -O -K >> base.ps
 #pstext  delta.txt   $PPAR -Dj+0.00/0.10                      -O -K >> base.ps
  pstext  cacity.txt  $PPAR                                    -O -K >> base.ps
 #psscale $SPAR -Cmmi10.cpt -B002:"MMI Shaking Intesity":      -O -K >> base.ps
  psbasemap -R -J$JP -B$AXIS -L$MBAR                           -O    >> base.ps
  endif

# Cleanup
# -------
  \rm -f tmp.e3d tmp.a tmp.llv tmp.xyz tmp.grd 
  \rm -f ustopo.xyz ustopo.grd fd.xyz fd.grd
  \rm -f fdbox.xy delta.xy epi.xy 
 #\rm -f $MMIGRD

