#! /bin/csh -f
#
#   script to mask some bad data from the SRTM grids.
#
if ($#argv != 2) then
	echo "usage: $0 incInSec oDir"
	echo "  example: $0 6 .."
	exit
endif

set incInSec = $1; shift
set oDir = $1; shift


cd $oDir/srtm/grd
ln -s ~/Desktop/SRTM$incInSec/com/path* .
#
grdmask path1.xy -Gpath1_mask.grd -N1/1/NaN -F `grdinfo -I- e140n40.grd` `grdinfo -I e140n40.grd`
mv e140n40.grd e140n40.grd.orig
grdmath path1_mask.grd e140n40.grd.orig MUL = e140n40.grd
#
grdmask path2.xy -Gpath2_mask.grd -N1/1/NaN -F `grdinfo -I- e140n90.grd` `grdinfo -I e140n90.grd`
mv e140n90.grd e140n90.grd.orig
grdmath path2_mask.grd e140n90.grd.orig MUL = e140n90.grd
#
grdmask path3.xy -Gpath3_mask.grd -N1/1/NaN -F `grdinfo -I- e100n40.grd` `grdinfo -I e100n40.grd`
mv e100n40.grd e100n40.grd.orig
grdmath path3_mask.grd e100n40.grd.orig MUL = e100n40.grd
#
rm path*.grd
cd ~/Desktop/SRTM$incInSec/com