#!/bin/bash
# GMT script to create a wrapped interferogram of deformation specifically for 
# AVN and JJ
# last modified Mon Apr  9 14:32:38 EDT 2012
# the Hector Mine region
##############################
#  DECLARE VARIABLES
    BASEGRD=/home/anewman/MGM/HectorMine/DEM/HM.grd
    ILLUM=`basename $BASEGRD`.ill
    InSARGRD=$1.grd

    OUTFILE=`basename $1 .hgt`.ps
    CPTGRAY=Gray.cpt
    CPTInSAR=Interfere.cpt 
    RANGE="-R-116.9/-115.42/34.09/35.31"
    PROJ="-JM7.3i"
    BGN=" $RANGE $PROJ -K -B0.1a.2WeSn -P -Y2.5i -X.95i "
    MID=" $RANGE $PROJ -K -O "
    END=" $RANGE $PROJ -O"
    RESAMPLE=0.00055556 # super high res
    RESAMPLE=0.001
##############################

##############################
# CREATE PLOT
gmtset PLOT_DEGREE_FORMAT D  PSIMAGE_FORMAT bin PAPER_MEDIA letter  BASEMAP_TYPE plain
psbasemap  $BGN > $OUTFILE

# make a gray cpt file for plotting background topography
echo "-10000 200 20000 200 " >$CPTGRAY
ZMIN=`grdinfo $InSARGRD | awk '$2~"z_min:"{printf "%i", $3-0.5}'`
ZMAX=`grdinfo $InSARGRD | awk '$2~"z_min:"{printf "%i", $5+0.5}'`
DZ=` echo $ZMIN $ZMAX | awk '{print ($2-$1)/10}'`
# PLOT BASEMAP
    grdgradient $BASEGRD -G$ILLUM -A225 -Ne0.6  # also needed for InSAR PLOT
#    grdimage $BASEGRD $MID  -I$ILLUM -C$CPTGRAY  >> $OUTFILE
 
# DOWN SAMPLE DATA   
    # downsample illumination gradient
      echo resampling $ILLUM to $ILLUM.rs
       grdsample $ILLUM -G$ILLUM.rs -I$RESAMPLE $RANGE
    # downsample interferogram phase data
      echo resampling $InSARGRD to $InSARGRD.rs
      echo grdsample  $InSARGRD -G$InSARGRD.rs -I$RESAMPLE $RANGE
      grdsample  $InSARGRD -G$InSARGRD.rs -I$RESAMPLE $RANGE
    # Plot INSAR DATA  
       makecpt  -Cinsar2 -Z -D -T-3.1416/3.1416/1.5708  > $CPTInSAR 
       grdimage $InSARGRD.rs $MID  -I$ILLUM.rs -Q -E200  -C$CPTInSAR  >> $OUTFILE 
# 
      psscale -D9.3/-1.5/8/.2h -C$CPTInSAR  -B1.57a3.14 -I -O -K >>$OUTFILE
      pstext  $END -N <<....EOF >> $OUTFILE 
         `echo $RANGE | cut -c 3- | awk -F"/" ' {print ($1+$2)/2,$4+0.12}' ` 18 0 0 6 Hector Mine Earthquake: Final Inteferogram
         `echo $RANGE | cut -c 3- | awk -F"/" ' {print ($1+$2)/2,$4+0.06}' ` 12 0 0 6 Synthetic Aperture Radar Interferometry:99/09/15 - 99/10/20
         `echo $RANGE | cut -c 3- | awk -F"/" ' {print ($1+$2)/2,$3-0.20}' ` 14 0 0 6 LOS Phase Change [28.3 mm/fringe]
         `echo $RANGE | cut -c 3- | awk -F"/" ' {print ($1+$2)/2,$3-0.28}' ` 10 0 2 1 `basename $InSARGRD .grd`
#         `echo $RANGE | cut -c 3- | awk -F"/" ' {print ($1+$2)/2,$3-0.09}' ` 12 0 0 6 "28.3mm/fringe"
....EOF
##############################

#   
# WANT TO CHANGE THE SCALE?
  #perl -pi.bak -e 's/0\.24\ 0\.24\ scale/0\.10\ 0\.10\ scale/g' $OUTFILE

# put hidden stamp in file that will denote its source
  echo " %% created by ${USER} using ${HOST}:${PWD}/$0 $* " >>$OUTFILE
  
gs  $OUTFILE

#pstoimg $OUTFILE

# if  ps2pdf $OUTFILE ; then rm $OUTFILE ; else echo PDF creation failed ; fi

#pstoimg `basename $OUTFILE .ps`.pdf