# Figures for FKt230303 MBARI Mapping AUV compilations of surveys at
#    EMARK, Hydra (Grappe Deux), and Puy des Folles Seamount
#
#----------------------------------------------------------------------
# 1) EMARK (inside corner high at the eastern intersection of Kane Fracture Zone and Mid-Atlantic Ridge)
#
# EMARK massif
#      -R -45.145/-45.0079256/23.4971540/23.6164003
# EMARK Falkor Vent Site
#      -R  -45.012836/-44.966428/23.455438/23.500020

# Make datalists

echo ../20230314m1/datalist.mb-1 -1 > datalist_EMARK_massif_MAUV.mb-1
echo ../20230314m2/datalist.mb-1 -1 >> datalist_EMARK_massif_MAUV.mb-1
echo ../20230315m1/datalist.mb-1 -1 >> datalist_EMARK_massif_MAUV.mb-1
echo ../20230316m1/datalist.mb-1 -1 >> datalist_EMARK_massif_MAUV.mb-1
echo ../20230317m1/datalist.mb-1 -1 >> datalist_EMARK_massif_MAUV.mb-1
echo ../20230317m2/datalist.mb-1 -1 >> datalist_EMARK_massif_MAUV.mb-1
echo ../20230318m1/datalist.mb-1 -1 >> datalist_EMARK_massif_MAUV.mb-1
echo ../20230319m1/datalist.mb-1 -1 >> datalist_EMARK_massif_MAUV.mb-1
echo ../20230319m2/datalist.mb-1 -1 >> datalist_EMARK_massif_MAUV.mb-1

# EMARK vents
echo ../20230321m1/datalist.mb-1 -1 > datalist_EMARK_vents_MAUV.mb-1
echo ../20230321m2/datalist.mb-1 -1 >> datalist_EMARK_vents_MAUV.mb-1
echo ../20230323m1/datalist.mb-1 -1 >> datalist_EMARK_vents_MAUV.mb-1

mbdatalist -Z -I datalist_EMARK_massif_MAUV.mb-1
mbdatalist -Z -I datalist_EMARK_vents_MAUV.mb-1

# Make netCDF bathymetry grids, 1-m resolution, geographic coordinates, in square cells for GIS
mbgrid -I datalist_EMARK_massif_MAUVp.mb-1 \
    -R -45.145/-45.0079256/23.4971540/23.6164003 \
  	-A2 -F5 -N -E1/0! -C5 \
  	-O EMARK_massif_MAUV_Topo1m -V
mbgrdviz -I EMARK_massif_MAUV_Topo1m.grd &

mbgrid -I datalist_EMARK_vents_MAUVp.mb-1 \
  	-A2 -F5 -N -E1/0! -C5 \
  	-O EMARK_vents_MAUV_Topo1m -V
mbgrdviz -I EMARK_vents_MAUV_Topo1m.grd &


# Make ASCII raster bathymetry grids
gmt grdconvert EMARK_massif_MAUV_Topo1m.grd -GEMARK_massif_MAUV_Topo1m.asc=ef+n-99999
gmt grdconvert EMARK_vents_MAUV_Topo1m.grd -GEMARK_vents_MAUV_Topo1m.asc=ef+n-99999

# Make geoTiffs, slope-shaded

mbm_grdtiff -I EMARK_massif_MAUV_Topo1m.grd \
  -G5 -D0/1 -A0.5 -Q \
  -O EMARK_massif_MAUV_Topo1m_slope -V
EMARK_massif_MAUV_Topo1m_slope_tiff.cmd  
# full linear color ramp: 4200 to 700 m depth

mbm_grdtiff -I EMARK_vents_MAUV_Topo1m.grd \
  -G5 -D0/1 -A0.5 -Q \
  -O EMARK_vents_MAUV_Topo1m_slope -V
EMARK_vents_MAUV_Topo1m_slope_tiff.cmd  
# full linear color ramp: 4500 to 3000 m depth


#----------------------------------------------------------------------
# 2) Hydra (also known as Grappe Deux; north of Kane Fracture Zone on the Mid-Atlantic Ridge)
#      -R -45.653453/-45.544935/24.929278/24.977264

# Make datalists
echo ../20230325m1/datalist.mb-1 -1 > datalist_Hydra_MAUV.mb-1
echo ../20230326m1/datalist.mb-1 -1 >> datalist_Hydra_MAUV.mb-1
echo ../20230328m1/datalist.mb-1 -1 >> datalist_Hydra_MAUV.mb-1
echo ../20230329m1/datalist.mb-1 -1 >> datalist_Hydra_MAUV.mb-1
mbdatalist -Z -I datalist_Hydra_MAUV.mb-1

# Make netCDF bathymetry grid, 1-m resolution, geographic coordinates, in square cells for GIS
mbgrid -I datalist_Hydra_MAUVp.mb-1 \
  	-A2 -F5 -N -E1/0! -C5 \
  	-O Hydra_MAUV_Topo1m -V
mbgrdviz -I Hydra_MAUV_Topo1m.grd &

# Make ASCII raster bathymetry grid
gmt grdconvert Hydra_MAUV_Topo1m.grd -GHydra_MAUV_Topo1m.asc=ef+n-99999

# Make geoTiff, slope-shaded

mbm_grdtiff -I Hydra_MAUV_Topo1m.grd \
   -G5 -D0/1 -Q -A0.5 \
   -O Hydra_MAUV_Topo1m_slope -V
Hydra_MAUV_Topo1m_slope_tiff.cmd
# full linear color ramp: 4375 to 2625 m depth

#----------------------------------------------------------------------
# 3A) Puy des Folles Seamount (on the Mid-Atlantic Ridge)
#      -R -45.661427/-45.618509/20.461722/20.534378
#
# Make datalist
echo ../20230310m1/datalist.mb-1 -1 > datalist_PuyDesFolles_MAUV.mb-1
echo ../20230310m2/datalist.mb-1 -1 >> datalist_PuyDesFolles_MAUV.mb-1
echo ../20230402m1/datalist.mb-1 -1 >> datalist_PuyDesFolles_MAUV.mb-1
mbdatalist -Z -I datalist_PuyDesFolles_MAUV.mb-1

# Make netFCDF bathymetry grid, 1-m resolution, geographic coordinates, in square cells for GIS
mbgrid -I datalist_PuyDesFolles_MAUVp.mb-1 \
  	-A2 -F5 -N -E1/0! -C5 \
  	-O PuyDesFolles_MAUV_Topo1m -V

# Make ASCII raster bathymetry grid
gmt grdconvert PuyDesFolles_MAUV_Topo1m.grd -GPuyDesFolles_MAUV_Topo1m.asc=ef+n-99999 -V

# Make geoTiffs, slope-shaded, in several color ramps

mbm_grdtiff -I PuyDesFolles_MAUV_Topo1m.grd \
  -G5 -D0/1 -Q \
  -O PuyDesFolles_MAUV_Topo1m_slope -V
PuyDesFolles_MAUV_Topo1m_slope_tiff.cmd
# full linear color ramp: -2375 to -1750 m depth

mbm_grdtiff -I PuyDesFolles_MAUV_Topo1m.grd \
  -G5 -D0/1 -Q -S \
  -O PuyDesFolles_MAUV_Topo1m_slopeHistEq -V
PuyDesFolles_MAUV_Topo1m_slopeHistEq_tiff.cmd

mbm_grdtiff -I PuyDesFolles_MAUV_Topo1m.grd \
  -G5 -D0/1 -Q -Y -Z -2045/-1895 \
  -O PuyDesFolles_MAUV_Topo1m_slope2045-1895 -V
PuyDesFolles_MAUV_Topo1m_slope2045-1895_tiff.cmd

#----------------------------------------------------------------------
# 3B) Close-ups of three chimney areas flown at 25-m altitude on survey 20230402m1, 
#      Puy des Folles Seamount (on the Mid-Atlantic Ridge)


# 3B.1) PdF "S500 NNW chimney field" at 50 cm resolution
#	 -R -45.6506673/-45.6442437/20.5161780/20.5211826

mbgrid -I ../20230402m1/datalistp.mb-1 \
	-R -45.6506673/-45.6442437/20.5161780/20.5211826 \
	-K PuyDesFolles_MAUV_Topo1m.grd \
  	-A2 -F5 -N -E0.5/0! -C5 \
  	-O PuyDesFolles_NNWchimneys_MAUV_Topo50cm -V
mbgrdviz -I PuyDesFolles_NNWchimneys_MAUV_Topo50cm.grd &

gmt grdconvert PuyDesFolles_NNWchimneys_MAUV_Topo50cm.grd -GPuyDesFolles_NNWchimneys_MAUV_Topo50cm.asc=ef+n-99999 -V

mbm_grdtiff -I PuyDesFolles_NNWchimneys_MAUV_Topo50cm.grd \
  -G5 -D0/1 -Q -Y -Z -1980/-1946 \
  -O PuyDesFolles_NNWchimneys_MAUV_Topo50cm_slope1980-1946 -V
PuyDesFolles_NNWchimneys_MAUV_Topo50cm_slope1980-1946_tiff.cmd


# 3B.2) PdF "Birthday Field" NE chimney field
# 	-R -45.6401853/-45.6327580/20.5101646/20.5158347

mbgrid -I ../20230402m1/datalistp.mb-1 \
	-R -45.6401853/-45.6327580/20.5101646/20.5158347 \
	-K PuyDesFolles_MAUV_Topo1m.grd \
  	-A2 -F5 -N -E0.5/0! -C5 \
  	-O PuyDesFolles_Birthday_MAUV_Topo50cm -V
mbgrdviz -I PuyDesFolles_Birthday_MAUV_Topo50cm.grd &

gmt grdconvert PuyDesFolles_Birthday_MAUV_Topo50cm.grd -GPuyDesFolles_Birthday_MAUV_Topo50cm.asc=ef+n-99999 -V

mbm_grdtiff -I PuyDesFolles_Birthday_MAUV_Topo50cm.grd \
  -G5 -D0/1 -Q -Y -Z -1979/-1909 \
  -O PuyDesFolles_Birthday_MAUV_Topo50cm_slope1979-1909 -V
PuyDesFolles_Birthday_MAUV_Topo50cm_slope1979-1909_tiff.cmd

# 3B.3 Southern chimneys (unobserved)
# 	-R -45.6443388/-45.6374397/20.4871093/20.4939988

mbgrid -I ../20230402m1/datalistp.mb-1 \
	-R -45.6443388/-45.6374397/20.4871093/20.4939988 \
	-K PuyDesFolles_MAUV_Topo1m.grd \
  	-A2 -F5 -N -E0.5/0! -C5 \
  	-O PuyDesFolles_Schimneys_MAUV_Topo50cm -V
mbgrdviz -I PuyDesFolles_Schimneys_MAUV_Topo50cm.grd &

gmt grdconvert PuyDesFolles_Schimneys_MAUV_Topo50cm.grd -GPuyDesFolles_Schimneys_MAUV_Topo50cm.asc=ef+n-99999 -V

mbm_grdtiff -I PuyDesFolles_Schimneys_MAUV_Topo50cm.grd \
  -G5 -D0/1 -Q -Y -Z -1956/-1896 \
  -O PuyDesFolles_Schimneys_MAUV_Topo50cm_slope1956-1896 -V
PuyDesFolles_Schimneys_MAUV_Topo50cm_slope1956-1896_tiff.cmd


