#! /bin/csh -f

# magic and hidden files: Another of the wonders of GMT

/bin/rm -rf .gmt*   >& /dev/null

# need lots of sig figs to keep surface from moaning

gmtset D_FORMAT = %.16lg
gmtset ELLIPSOID = Sphere

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

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

# /bin/rm -rf $oDir/cmTile/     >& /dev/null
# /bin/rm -rf $oDir/predTile/   >& /dev/null
# /bin/rm -rf $oDir/grd/        >& /dev/null
/bin/rm -rf $oDir/logs/         >& /dev/null
/bin/rm -rf $oDir/sid/          >& /dev/null
/bin/rm -rf $oDir/jpeg/         >& /dev/null
/bin/rm -rf $oDir/ermapper/     >& /dev/null
/bin/rm -rf $oDir/plots/        >& /dev/null
/bin/rm -rf $oDir/pingMap/      >& /dev/null
    
/bin/mkdir -p  $oDir/tmp
/bin/mkdir -p  $oDir/logs
/bin/mkdir -p  $oDir/grd
/bin/mkdir -p  $oDir/predTile
date        >! $oDir/logs/makeAll.log

set img = '../img/topo.img'
set imgIncInSec = 60

set MOA_public =  /Volumes/publicSonarData/seasat2/data/public
set MOA_private = /Volumes/publicSonarData/seasat2/data/private

# hardcode arctic bathymetry sid and bounds

set minArctic = 69
set maxPred = 70
set arcticSid = 00099
set arcticGrd = '../arctic/IBCAO_ver2_23_GEO_netCDF_1min_trimmed.grd'
set predOverlap = "$arcticGrd $minArctic $maxPred $arcticSid"
echo; echo "Hard coding arctic grid SID:         $arcticSid"
echo "Overlap area of predicted and arctic: -R/-180/180/$minArctic/$maxPred";echo

# the next step is only needed if the SRTM is updated 
# or if you want to copy it over again, or something like that
#
#

# set palsar = '/palsar'
# echo "copying SRTM DEM from "$palsar" to ../srtm/dem"
# 
# if (! -d "$palsar") then
#        echo; echo; echo "ERROR: $palsar is not mounted" ; 
#   echo; echo; echo "ABORTING..."; echo; echo; exit 1
#    endif
#  
# rm -rf ../srtm/dem
# cp -r $palsar/DEM/SRTM$incInSec/ ../srtm/dem
# 
# echo "converting SRTM$incInSec DEM to grd"
# makeAllDemGrd |& cat >> $oDir/logs/makeAll.log
# mask_bad_land_grd.com $oDir
# echo "done with SRTM DEM; ok to unmount $palsar"


# #FIXME: for the time being jjMedian is hacked to avoid dealing with making a real SRTM15...
# echo "making a SRTM land directory"
# mkdir -p $oDir/srtm/grd
# 
echo "converting xyzi tracks to xyzi tiles (eg sort into SRTM tiles)"
makeAllCmTiles.sh $MOA_public $MOA_private $predOverlap $oDir |& cat >> $oDir/logs/makeAll.log
#
# echo "make list of all SID, and search for duplicates"
# makeSidLists $MOA_public $MOA_private $oDir   |& cat >> $oDir/logs/makeAll.log
#
# echo "For debugging, sort all pings by depth, and by lat, lon, and sid"
# makeSorted  $oDir/test |& cat >> $oDir/logs/makeAll.log

echo "convert predicted bathymetry from World of Dave to $incInSec second grd"
makeAllPredTiles $img 60  $arcticGrd $minArctic $maxPred $arcticSid $incInSec $oDir |& cat >> $oDir/logs/makeAll.log

echo "median filter pings with predicted depth, and keep track of sid"
makeAllMedian $incInSec $oDir   |& cat >> $oDir/logs/makeAll.log
echo "add SRTM land data to grd"
makeAllMerged $incInSec $oDir   |& cat >> $oDir/logs/makeAll.log

# echo "find locations with at least one ping"
# makePingMap $arcticSid $incInSec $oDir |& cat >> $oDir/logs/makeAll.log

echo "make postscript plots of topography and hit maps"
# Scale for 8.5 x 11 paper (i.e. 0.25 on n90 tiles)
makeAllPlots 0.25 $incInSec $oDir |& cat >> $oDir/logs/makeAll.log
echo "convert postscript in JPEG"
makeAllJpeg $oDir |& cat >> $oDir/logs/makeAll.log

# echo "make ERMapper files"
#FIXME:
#makeAllErmapper $oDir  |& cat >> $oDir/logs/makeAll.log

date    |& cat >> $oDir/logs/makeAll.log