#!/bin/bash

	# ASSUMES GMT 5
	
	#	SWITCH GMT VERSION MANUALLY
	gmtswitch /Users/esg006/gmt5  

	

	PARFILE=$1
	TRAXIS=$2
    MASKOUT=$3
   
    source $PARFILE
    
    rm age_mask_ll.xy traxis.xy
    
    RLBOUNDS="-R${WLON}/${ELON}/${SLAT}/${NLAT}"
    
    ln -s /radar3/flex/age/age_mask/${TRAXIS} traxis.xy
	
	gmt gmtselect traxis.xy $RLBOUNDS -V > age_mask_ll.xy 
	echo $RLBOUNDS
	
	case "$MASKOUT" in
	
	# assumes trench axis is running WEST to EAST
	N)	echo "Masking out to the NORTH of the trench"
		tail -n1 age_mask_ll.xy | awk '{print $2}' > tmplat
		read TLAT < tmplat
		echo "${ELON} ${TLAT}" >> age_mask_ll.xy 
		echo "${ELON} ${NLAT}" >> age_mask_ll.xy
		echo "${WLON} ${NLAT}" >> age_mask_ll.xy
		head -n1 age_mask_ll.xy | awk '{print $2}' > tmplat
		read TLAT < tmplat
		echo "${WLON} ${TLAT}" >> age_mask_ll.xy
		head -n1 age_mask_ll.xy >> age_mask_ll.xy
		;;

	# assumes trench axis is running NORTH to SOUTH
	E)	echo "Masking out to the EAST of the trench"
		tail -n1 age_mask_ll.xy | awk '{print $1}' > tmplon
		read TLON < tmplon
		echo "${TLON} ${SLAT}" >> age_mask_ll.xy
		echo "${ELON} ${SLAT}" >> age_mask_ll.xy
		echo "${ELON} ${NLAT}" >> age_mask_ll.xy
		head -n1 age_mask_ll.xy | awk '{print $1}' > tmplon
		read TLON < tmplon
		echo "${TLON} ${NLAT}" >> age_mask_ll.xy
		head -n1 age_mask_ll.xy >> age_mask_ll.xy
		;;

	# assumes trench axis is running EAST to WEST
	S)	echo "Masking out to the SOUTH of the trench"
		tail -n1 age_mask_ll.xy | awk -v LON="$WLON" '{print $1,$LON}' >> age_mask_ll.xy
		echo "${WLON} ${SLAT}" >> age_mask_ll.xy
		echo "${ELON} ${SLAT}" >> age_mask_ll.xy
		head -n1 age_mask_ll.xy | awk -v LON="$ELON" '{print $1,$LON}' >> age_mask_ll.xy
		head -n1 age_mask_ll.xy >> age_mask_ll.xy
		;;

	# assumes trench axis is running SOUTH to NORTH
	W)	echo "Masking out to the WEST of the trench"
		tail -n1 age_mask_ll.xy | awk '{print $1}' > tmplon
		read TLON < tmplon
		echo "${TLON} ${NLAT}" >> age_mask_ll.xy
		echo "${WLON} ${NLAT}" >> age_mask_ll.xy
		echo "${WLON} ${SLAT}" >> age_mask_ll.xy
		head -n1 age_mask_ll.xy | awk '{print $1}' > tmplon
		read TLON < tmplon
		echo "${TLON} ${SLAT}" >> age_mask_ll.xy
 		head -n1 age_mask_ll.xy >> age_mask_ll.xy
		;;
	esac