#! /bin/csh -f

if ($#argv != 1) then
	echo "usage: $0 oDir"
	echo "  example: $0 .."
	exit
endif

set oDir = $1; shift


/bin/rm -rf   $oDir/srtm/grd  >& /dev/null
/bin/mkdir -p $oDir/srtm/grd 

# convert USGS dem format to GMT grd format,
# then convert NOAA costal multibeam to ASCII (yes that isn't a grd)
#
# Note we use the GLAS data for antarctica; mo better than GEBCO
# Note the lake files include elevations below lake level

jjConvertDem2Grd $oDir/srtm/dem/E020N40_lake.DEM   020  060 -010  40 $oDir/srtm/grd/e020n40.grd &
jjConvertDem2Grd $oDir/srtm/dem/E020N90_lake.DEM   020  060  040  90 $oDir/srtm/grd/e020n90.grd &
jjConvertDem2Grd $oDir/srtm/dem/E020S10.DEM   020  060 -060 -10 $oDir/srtm/grd/e020s10.grd &
jjConvertDem2Grd $oDir/srtm/dem/E060N40.DEM   060  100 -010  40 $oDir/srtm/grd/e060n40.grd &
jjConvertDem2Grd $oDir/srtm/dem/E060N90.DEM   060  100  040  90 $oDir/srtm/grd/e060n90.grd &
jjConvertDem2Grd $oDir/srtm/dem/E060S10.DEM   060  100 -060 -10 $oDir/srtm/grd/e060s10.grd &
jjConvertDem2Grd $oDir/srtm/dem/E100N40.DEM   100  140 -010  40 $oDir/srtm/grd/e100n40.grd &
jjConvertDem2Grd $oDir/srtm/dem/E100N90_lake.DEM   100  140  040  90 $oDir/srtm/grd/e100n90.grd &
wait
jjConvertDem2Grd $oDir/srtm/dem/E100S10.DEM   100  140 -060 -10 $oDir/srtm/grd/e100s10.grd &
jjConvertDem2Grd $oDir/srtm/dem/E140N40.DEM   140  180 -010  40 $oDir/srtm/grd/e140n40.grd &
jjConvertDem2Grd $oDir/srtm/dem/E140N90.DEM   140  180  040  90 $oDir/srtm/grd/e140n90.grd &
jjConvertDem2Grd $oDir/srtm/dem/E140S10.DEM   140  180 -060 -10 $oDir/srtm/grd/e140s10.grd &
jjConvertDem2Grd $oDir/srtm/dem/W020N40.DEM  -020  020 -010  40 $oDir/srtm/grd/w020n40.grd &
jjConvertDem2Grd $oDir/srtm/dem/W020N90.DEM  -020  020  040  90 $oDir/srtm/grd/w020n90.grd &
jjConvertDem2Grd $oDir/srtm/dem/W020S10.DEM  -020  020 -060 -10 $oDir/srtm/grd/w020s10.grd &
jjConvertDem2Grd $oDir/srtm/dem/W060N40.DEM  -060 -020 -010  40 $oDir/srtm/grd/w060n40.grd &
wait
jjConvertDem2Grd $oDir/srtm/dem/W060N90.DEM  -060 -020  040  90 $oDir/srtm/grd/w060n90.grd &
jjConvertDem2Grd $oDir/srtm/dem/W060S10.DEM  -060 -020 -060 -10 $oDir/srtm/grd/w060s10.grd &
jjConvertDem2Grd $oDir/srtm/dem/W100N40.DEM  -100 -060 -010  40 $oDir/srtm/grd/w100n40.grd &
jjConvertDem2Grd $oDir/srtm/dem/W100N90_lake.DEM  -100 -060  040  90 $oDir/srtm/grd/w100n90.grd &
jjConvertDem2Grd $oDir/srtm/dem/W100S10.DEM  -100 -060 -060 -10 $oDir/srtm/grd/w100s10.grd &
jjConvertDem2Grd $oDir/srtm/dem/W140N40.DEM  -140 -100 -010  40 $oDir/srtm/grd/w140n40.grd &
jjConvertDem2Grd $oDir/srtm/dem/W140N90.DEM  -140 -100  040  90 $oDir/srtm/grd/w140n90.grd &
jjConvertDem2Grd $oDir/srtm/dem/W140S10.DEM  -140 -100 -060 -10 $oDir/srtm/grd/w140s10.grd &
wait
jjConvertDem2Grd $oDir/srtm/dem/W180N40.DEM  -180 -140 -010  40 $oDir/srtm/grd/w180n40.grd &
jjConvertDem2Grd $oDir/srtm/dem/W180N90.DEM  -180 -140  040  90 $oDir/srtm/grd/w180n90.grd &
jjConvertDem2Grd $oDir/srtm/dem/W180S10.DEM  -180 -140 -060 -10 $oDir/srtm/grd/w180s10.grd &

# Note we use the GLAS data for antarctica; mo better than GEBCO

jjConvertDem2Grd $oDir/srtm/dem/W180S60.glas -180 -120 -090 -60 $oDir/srtm/grd/w180s60.grd &
jjConvertDem2Grd $oDir/srtm/dem/W120S60.glas -120 -060 -090 -60 $oDir/srtm/grd/w120s60.grd &
jjConvertDem2Grd $oDir/srtm/dem/W060S60.glas -060  000 -090 -60 $oDir/srtm/grd/w060s60.grd &
jjConvertDem2Grd $oDir/srtm/dem/W000S60.glas  000  060 -090 -60 $oDir/srtm/grd/w000s60.grd &
jjConvertDem2Grd $oDir/srtm/dem/E060S60.glas  060  120 -090 -60 $oDir/srtm/grd/e060s60.grd &
jjConvertDem2Grd $oDir/srtm/dem/E120S60.glas  120  180 -090 -60 $oDir/srtm/grd/e120s60.grd &
wait

date