diff --git a/src/core_ocean/mode_init/mpas_ocn_init_hurricane.F b/src/core_ocean/mode_init/mpas_ocn_init_hurricane.F
index 6b07d815b8..88b6124c1c 100644
--- a/src/core_ocean/mode_init/mpas_ocn_init_hurricane.F
+++ b/src/core_ocean/mode_init/mpas_ocn_init_hurricane.F
@@ -324,16 +324,8 @@ subroutine ocn_init_setup_hurricane(domain, iErr)!{{{
end if
end do
- ! shift coordinates to accomodate wetting and drying
- do k = 1, nVertLevels
- refBottomDepth(k) = refBottomDepth(k) - globalMinBottomDepth
- end do
-
do iCell = 1, nCellsSolve
- bottomDepth(iCell) = bottomDepth(iCell) - globalMinBottomDepth
- ssh(iCell) = ssh(iCell) + globalMinBottomDepth
-
! make sure depth is thick enough via ssh = TOTAL_DEPTH - bottomDepth
ssh(iCell) = - bottomDepth(iCell) + &
max(ssh(iCell) + bottomDepth(iCell), &
diff --git a/testing_and_setup/compass/ocean/hurricane/USDEQU120at30cr10rr2WD/build_mesh/config_base_mesh.xml b/testing_and_setup/compass/ocean/hurricane/USDEQU120at30cr10rr2WD/build_mesh/config_base_mesh.xml
new file mode 100755
index 0000000000..de249353da
--- /dev/null
+++ b/testing_and_setup/compass/ocean/hurricane/USDEQU120at30cr10rr2WD/build_mesh/config_base_mesh.xml
@@ -0,0 +1,19 @@
+
+
+
+
+
+
+
+
+
+
+
+
+ jigsaw_to_MPAS.build_mesh
+
+
+ 10.0
+
+
+
diff --git a/testing_and_setup/compass/ocean/hurricane/USDEQU120at30cr10rr2WD/build_mesh/config_culled_mesh.xml b/testing_and_setup/compass/ocean/hurricane/USDEQU120at30cr10rr2WD/build_mesh/config_culled_mesh.xml
new file mode 100644
index 0000000000..355c59a252
--- /dev/null
+++ b/testing_and_setup/compass/ocean/hurricane/USDEQU120at30cr10rr2WD/build_mesh/config_culled_mesh.xml
@@ -0,0 +1,30 @@
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+ jigsaw_to_MPAS.inject_bathymetry
+ culled_mesh.nc
+
+
+ --ignore_time
+ -l
+ maxEdges=0
+ allOnCells
+ culled_mesh.nc
+ culled_mesh_vtk
+
+
+
diff --git a/testing_and_setup/compass/ocean/hurricane/USDEQU120at30cr10rr2WD/build_mesh/config_driver_init.xml b/testing_and_setup/compass/ocean/hurricane/USDEQU120at30cr10rr2WD/build_mesh/config_driver_init.xml
new file mode 100644
index 0000000000..5b34eec160
--- /dev/null
+++ b/testing_and_setup/compass/ocean/hurricane/USDEQU120at30cr10rr2WD/build_mesh/config_driver_init.xml
@@ -0,0 +1,8 @@
+
+
+
+
+
+
+
+
diff --git a/testing_and_setup/compass/ocean/hurricane/USDEQU120at30cr10rr2WD/build_mesh/define_base_mesh.py b/testing_and_setup/compass/ocean/hurricane/USDEQU120at30cr10rr2WD/build_mesh/define_base_mesh.py
new file mode 100755
index 0000000000..ffd75cf45d
--- /dev/null
+++ b/testing_and_setup/compass/ocean/hurricane/USDEQU120at30cr10rr2WD/build_mesh/define_base_mesh.py
@@ -0,0 +1,54 @@
+#!/usr/bin/env python
+
+import numpy as np
+import jigsaw_to_MPAS.coastal_tools as ct
+
+def cellWidthVsLatLon():
+ km = 1000.0
+
+ params = ct.default_params
+
+ print("****QU 120 background mesh and enhanced Atlantic (30km)****")
+ params["mesh_type"] = "QU"
+ params["dx_max_global"] = 120.0 * km
+ params["region_box"] = ct.Atlantic
+ params["restrict_box"] = ct.Atlantic_restrict
+ params["plot_box"] = ct.Western_Atlantic
+ params["dx_min_coastal"] = 30.0 * km
+ params["trans_width"] = 5000.0 * km
+ params["trans_start"] = 500.0 * km
+
+ cell_width, lon, lat = ct.coastal_refined_mesh(params)
+
+ print("****Northeast refinement (10km)***")
+ params["region_box"] = ct.Delaware_Bay
+ params["plot_box"] = ct.Western_Atlantic
+ params["dx_min_coastal"] = 10.0 * km
+ params["trans_width"] = 600.0 * km
+ params["trans_start"] = 400.0 * km
+
+ cell_width, lon, lat = ct.coastal_refined_mesh(
+ params, cell_width, lon, lat)
+
+ print("****Delaware regional refinement (5km)****")
+ params["region_box"] = ct.Delaware_Region
+ params["plot_box"] = ct.Delaware
+ params["dx_min_coastal"] = 5.0 * km
+ params["trans_width"] = 175.0 * km
+ params["trans_start"] = 75.0 * km
+
+ cell_width, lon, lat = ct.coastal_refined_mesh(
+ params, cell_width, lon, lat)
+
+ print("****Delaware Bay high-resolution (2km)****")
+ params["region_box"] = ct.Delaware_Bay
+ params["plot_box"] = ct.Delaware
+ params["restrict_box"] = ct.Delaware_restrict
+ params["dx_min_coastal"] = 2.0 * km
+ params["trans_width"] = 100.0 * km
+ params["trans_start"] = 17.0 * km
+
+ cell_width, lon, lat = ct.coastal_refined_mesh(
+ params, cell_width, lon, lat)
+
+ return cell_width / 1000, lon, lat
diff --git a/testing_and_setup/compass/ocean/hurricane/USDEQU120at30cr10rr2WD/sandy/config_driver.xml b/testing_and_setup/compass/ocean/hurricane/USDEQU120at30cr10rr2WD/sandy/config_driver.xml
new file mode 100644
index 0000000000..2bbc9940c7
--- /dev/null
+++ b/testing_and_setup/compass/ocean/hurricane/USDEQU120at30cr10rr2WD/sandy/config_driver.xml
@@ -0,0 +1,8 @@
+
+
+
+
+
+
+
+
diff --git a/testing_and_setup/compass/ocean/hurricane/USDEQU120at30cr10rr2WD/sandy/config_forward_RK4.xml b/testing_and_setup/compass/ocean/hurricane/USDEQU120at30cr10rr2WD/sandy/config_forward_RK4.xml
new file mode 100644
index 0000000000..ca1772b126
--- /dev/null
+++ b/testing_and_setup/compass/ocean/hurricane/USDEQU120at30cr10rr2WD/sandy/config_forward_RK4.xml
@@ -0,0 +1,39 @@
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+ 360
+
+
+
+
diff --git a/testing_and_setup/compass/ocean/hurricane/USDEQU120at30cr10rr2WD/sandy/config_init.xml b/testing_and_setup/compass/ocean/hurricane/USDEQU120at30cr10rr2WD/sandy/config_init.xml
new file mode 100644
index 0000000000..906d4d07d1
--- /dev/null
+++ b/testing_and_setup/compass/ocean/hurricane/USDEQU120at30cr10rr2WD/sandy/config_init.xml
@@ -0,0 +1,42 @@
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+ mesh.nc
+ USGS_stations.txt
+ NOAA-COOPS_stations.txt
+
+
+
+
+
diff --git a/testing_and_setup/compass/ocean/hurricane/USDEQU240at60cr20rr4WD/build_mesh/config_base_mesh.xml b/testing_and_setup/compass/ocean/hurricane/USDEQU240at60cr20rr4WD/build_mesh/config_base_mesh.xml
new file mode 100755
index 0000000000..de249353da
--- /dev/null
+++ b/testing_and_setup/compass/ocean/hurricane/USDEQU240at60cr20rr4WD/build_mesh/config_base_mesh.xml
@@ -0,0 +1,19 @@
+
+
+
+
+
+
+
+
+
+
+
+
+ jigsaw_to_MPAS.build_mesh
+
+
+ 10.0
+
+
+
diff --git a/testing_and_setup/compass/ocean/hurricane/USDEQU240at60cr20rr4WD/build_mesh/config_culled_mesh.xml b/testing_and_setup/compass/ocean/hurricane/USDEQU240at60cr20rr4WD/build_mesh/config_culled_mesh.xml
new file mode 100644
index 0000000000..355c59a252
--- /dev/null
+++ b/testing_and_setup/compass/ocean/hurricane/USDEQU240at60cr20rr4WD/build_mesh/config_culled_mesh.xml
@@ -0,0 +1,30 @@
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+ jigsaw_to_MPAS.inject_bathymetry
+ culled_mesh.nc
+
+
+ --ignore_time
+ -l
+ maxEdges=0
+ allOnCells
+ culled_mesh.nc
+ culled_mesh_vtk
+
+
+
diff --git a/testing_and_setup/compass/ocean/hurricane/USDEQU240at60cr20rr4WD/build_mesh/config_driver_init.xml b/testing_and_setup/compass/ocean/hurricane/USDEQU240at60cr20rr4WD/build_mesh/config_driver_init.xml
new file mode 100644
index 0000000000..5b34eec160
--- /dev/null
+++ b/testing_and_setup/compass/ocean/hurricane/USDEQU240at60cr20rr4WD/build_mesh/config_driver_init.xml
@@ -0,0 +1,8 @@
+
+
+
+
+
+
+
+
diff --git a/testing_and_setup/compass/ocean/hurricane/USDEQU240at60cr20rr4WD/build_mesh/define_base_mesh.py b/testing_and_setup/compass/ocean/hurricane/USDEQU240at60cr20rr4WD/build_mesh/define_base_mesh.py
new file mode 100755
index 0000000000..56b5350cc0
--- /dev/null
+++ b/testing_and_setup/compass/ocean/hurricane/USDEQU240at60cr20rr4WD/build_mesh/define_base_mesh.py
@@ -0,0 +1,54 @@
+#!/usr/bin/env python
+
+import numpy as np
+import jigsaw_to_MPAS.coastal_tools as ct
+
+def cellWidthVsLatLon():
+ km = 1000.0
+
+ params = ct.default_params
+
+ print("****QU 240 background mesh and enhanced Atlantic (60km)****")
+ params["mesh_type"] = "QU"
+ params["dx_max_global"] = 240.0 * km
+ params["region_box"] = ct.Atlantic
+ params["restrict_box"] = ct.Atlantic_restrict
+ params["plot_box"] = ct.Western_Atlantic
+ params["dx_min_coastal"] = 60.0 * km
+ params["trans_width"] = 5000.0 * km
+ params["trans_start"] = 500.0 * km
+
+ cell_width, lon, lat = ct.coastal_refined_mesh(params)
+
+ print("****Northeast refinement (20km)***")
+ params["region_box"] = ct.Delaware_Bay
+ params["plot_box"] = ct.Western_Atlantic
+ params["dx_min_coastal"] = 20.0 * km
+ params["trans_width"] = 600.0 * km
+ params["trans_start"] = 400.0 * km
+
+ cell_width, lon, lat = ct.coastal_refined_mesh(
+ params, cell_width, lon, lat)
+
+ print("****Delaware regional refinement (10km)****")
+ params["region_box"] = ct.Delaware_Region
+ params["plot_box"] = ct.Delaware
+ params["dx_min_coastal"] = 10.0 * km
+ params["trans_width"] = 175.0 * km
+ params["trans_start"] = 75.0 * km
+
+ cell_width, lon, lat = ct.coastal_refined_mesh(
+ params, cell_width, lon, lat)
+
+ print("****Delaware Bay high-resolution (4km)****")
+ params["region_box"] = ct.Delaware_Bay
+ params["plot_box"] = ct.Delaware
+ params["restrict_box"] = ct.Delaware_restrict
+ params["dx_min_coastal"] = 4.0 * km
+ params["trans_width"] = 100.0 * km
+ params["trans_start"] = 17.0 * km
+
+ cell_width, lon, lat = ct.coastal_refined_mesh(
+ params, cell_width, lon, lat)
+
+ return cell_width / 1000, lon, lat
diff --git a/testing_and_setup/compass/ocean/hurricane/USDEQU240at60cr20rr4WD/sandy/config_driver.xml b/testing_and_setup/compass/ocean/hurricane/USDEQU240at60cr20rr4WD/sandy/config_driver.xml
new file mode 100644
index 0000000000..2bbc9940c7
--- /dev/null
+++ b/testing_and_setup/compass/ocean/hurricane/USDEQU240at60cr20rr4WD/sandy/config_driver.xml
@@ -0,0 +1,8 @@
+
+
+
+
+
+
+
+
diff --git a/testing_and_setup/compass/ocean/hurricane/USDEQU240at60cr20rr4WD/sandy/config_forward_RK4.xml b/testing_and_setup/compass/ocean/hurricane/USDEQU240at60cr20rr4WD/sandy/config_forward_RK4.xml
new file mode 100644
index 0000000000..1d63cbafa8
--- /dev/null
+++ b/testing_and_setup/compass/ocean/hurricane/USDEQU240at60cr20rr4WD/sandy/config_forward_RK4.xml
@@ -0,0 +1,39 @@
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+ 72
+
+
+
+
diff --git a/testing_and_setup/compass/ocean/hurricane/USDEQU240at60cr20rr4WD/sandy/config_init.xml b/testing_and_setup/compass/ocean/hurricane/USDEQU240at60cr20rr4WD/sandy/config_init.xml
new file mode 100644
index 0000000000..d637fbedae
--- /dev/null
+++ b/testing_and_setup/compass/ocean/hurricane/USDEQU240at60cr20rr4WD/sandy/config_init.xml
@@ -0,0 +1,37 @@
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+ mesh.nc
+ USGS_stations.txt
+ NOAA-COOPS_stations.txt
+
+
+
+
+
diff --git a/testing_and_setup/compass/ocean/hurricane/USDEQU60at15cr5rr1WD/build_mesh/config_base_mesh.xml b/testing_and_setup/compass/ocean/hurricane/USDEQU60at15cr5rr1WD/build_mesh/config_base_mesh.xml
new file mode 100755
index 0000000000..de249353da
--- /dev/null
+++ b/testing_and_setup/compass/ocean/hurricane/USDEQU60at15cr5rr1WD/build_mesh/config_base_mesh.xml
@@ -0,0 +1,19 @@
+
+
+
+
+
+
+
+
+
+
+
+
+ jigsaw_to_MPAS.build_mesh
+
+
+ 10.0
+
+
+
diff --git a/testing_and_setup/compass/ocean/hurricane/USDEQU60at15cr5rr1WD/build_mesh/config_culled_mesh.xml b/testing_and_setup/compass/ocean/hurricane/USDEQU60at15cr5rr1WD/build_mesh/config_culled_mesh.xml
new file mode 100644
index 0000000000..5c62737241
--- /dev/null
+++ b/testing_and_setup/compass/ocean/hurricane/USDEQU60at15cr5rr1WD/build_mesh/config_culled_mesh.xml
@@ -0,0 +1,32 @@
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+ jigsaw_to_MPAS.inject_bathymetry
+ culled_mesh.nc
+
+
+
+ --ignore_time
+ -l
+ maxEdges=0
+ allOnCells
+ culled_mesh.nc
+ culled_mesh_vtk
+
+
+
+
diff --git a/testing_and_setup/compass/ocean/hurricane/USDEQU60at15cr5rr1WD/build_mesh/config_driver_init.xml b/testing_and_setup/compass/ocean/hurricane/USDEQU60at15cr5rr1WD/build_mesh/config_driver_init.xml
new file mode 100644
index 0000000000..5b34eec160
--- /dev/null
+++ b/testing_and_setup/compass/ocean/hurricane/USDEQU60at15cr5rr1WD/build_mesh/config_driver_init.xml
@@ -0,0 +1,8 @@
+
+
+
+
+
+
+
+
diff --git a/testing_and_setup/compass/ocean/hurricane/USDEQU60at15cr5rr1WD/build_mesh/define_base_mesh.py b/testing_and_setup/compass/ocean/hurricane/USDEQU60at15cr5rr1WD/build_mesh/define_base_mesh.py
new file mode 100755
index 0000000000..dcb941aa2d
--- /dev/null
+++ b/testing_and_setup/compass/ocean/hurricane/USDEQU60at15cr5rr1WD/build_mesh/define_base_mesh.py
@@ -0,0 +1,54 @@
+#!/usr/bin/env python
+
+import numpy as np
+import jigsaw_to_MPAS.coastal_tools as ct
+
+def cellWidthVsLatLon():
+ km = 1000.0
+
+ params = ct.default_params
+
+ print("****QU 60 background mesh and enhanced Atlantic (30km)****")
+ params["mesh_type"] = "QU"
+ params["dx_max_global"] = 60.0 * km
+ params["region_box"] = ct.Atlantic
+ params["restrict_box"] = ct.Atlantic_restrict
+ params["plot_box"] = ct.Western_Atlantic
+ params["dx_min_coastal"] = 15.0 * km
+ params["trans_width"] = 5000.0 * km
+ params["trans_start"] = 500.0 * km
+
+ cell_width, lon, lat = ct.coastal_refined_mesh(params)
+
+ print("****Northeast refinement (5km)***")
+ params["region_box"] = ct.Delaware_Bay
+ params["plot_box"] = ct.Western_Atlantic
+ params["dx_min_coastal"] = 5.0 * km
+ params["trans_width"] = 600.0 * km
+ params["trans_start"] = 400.0 * km
+
+ cell_width, lon, lat = ct.coastal_refined_mesh(
+ params, cell_width, lon, lat)
+
+ print("****Delaware regional refinement (2.5km)****")
+ params["region_box"] = ct.Delaware_Region
+ params["plot_box"] = ct.Delaware
+ params["dx_min_coastal"] = 2.5 * km
+ params["trans_width"] = 175.0 * km
+ params["trans_start"] = 75.0 * km
+
+ cell_width, lon, lat = ct.coastal_refined_mesh(
+ params, cell_width, lon, lat)
+
+ print("****Delaware Bay high-resolution (1km)****")
+ params["region_box"] = ct.Delaware_Bay
+ params["plot_box"] = ct.Delaware
+ params["restrict_box"] = ct.Delaware_restrict
+ params["dx_min_coastal"] = 1.0 * km
+ params["trans_width"] = 100.0 * km
+ params["trans_start"] = 17.0 * km
+
+ cell_width, lon, lat = ct.coastal_refined_mesh(
+ params, cell_width, lon, lat)
+
+ return cell_width / 1000, lon, lat
diff --git a/testing_and_setup/compass/ocean/hurricane/USDEQU60at15cr5rr1WD/sandy/config_driver.xml b/testing_and_setup/compass/ocean/hurricane/USDEQU60at15cr5rr1WD/sandy/config_driver.xml
new file mode 100644
index 0000000000..be1976f5c8
--- /dev/null
+++ b/testing_and_setup/compass/ocean/hurricane/USDEQU60at15cr5rr1WD/sandy/config_driver.xml
@@ -0,0 +1,8 @@
+
+
+
+
+
+
+
+
diff --git a/testing_and_setup/compass/ocean/hurricane/USDEQU60at15cr5rr1WD/sandy/config_forward.xml b/testing_and_setup/compass/ocean/hurricane/USDEQU60at15cr5rr1WD/sandy/config_forward.xml
new file mode 100644
index 0000000000..b84b723e6c
--- /dev/null
+++ b/testing_and_setup/compass/ocean/hurricane/USDEQU60at15cr5rr1WD/sandy/config_forward.xml
@@ -0,0 +1,39 @@
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+ 1080
+
+
+
+
diff --git a/testing_and_setup/compass/ocean/hurricane/USDEQU60at15cr5rr1WD/sandy/config_init.xml b/testing_and_setup/compass/ocean/hurricane/USDEQU60at15cr5rr1WD/sandy/config_init.xml
new file mode 100644
index 0000000000..906d4d07d1
--- /dev/null
+++ b/testing_and_setup/compass/ocean/hurricane/USDEQU60at15cr5rr1WD/sandy/config_init.xml
@@ -0,0 +1,42 @@
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+ mesh.nc
+ USGS_stations.txt
+ NOAA-COOPS_stations.txt
+
+
+
+
+
diff --git a/testing_and_setup/compass/ocean/hurricane/forward_variable_drag.xml b/testing_and_setup/compass/ocean/hurricane/forward_variable_drag.xml
new file mode 100644
index 0000000000..a90c6fbcd8
--- /dev/null
+++ b/testing_and_setup/compass/ocean/hurricane/forward_variable_drag.xml
@@ -0,0 +1,19 @@
+
+
+
+
+
+
+
+
+
+
+ input
+ initial_only
+ forcing.nc
+
+
+
+
+
+
diff --git a/testing_and_setup/compass/ocean/hurricane/forward_wetting_drying_template.xml b/testing_and_setup/compass/ocean/hurricane/forward_wetting_drying_template.xml
new file mode 100644
index 0000000000..84e148250f
--- /dev/null
+++ b/testing_and_setup/compass/ocean/hurricane/forward_wetting_drying_template.xml
@@ -0,0 +1,17 @@
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
diff --git a/testing_and_setup/compass/ocean/hurricane/init_variable_drag.xml b/testing_and_setup/compass/ocean/hurricane/init_variable_drag.xml
new file mode 100644
index 0000000000..e62ecad5e4
--- /dev/null
+++ b/testing_and_setup/compass/ocean/hurricane/init_variable_drag.xml
@@ -0,0 +1,16 @@
+
+
+
+
+
+
+ output
+ forcing.nc
+ truncate
+ 0000_00:00:01
+
+
+
+
+
+
diff --git a/testing_and_setup/compass/ocean/hurricane/init_wetting_drying_template.xml b/testing_and_setup/compass/ocean/hurricane/init_wetting_drying_template.xml
new file mode 100644
index 0000000000..a055e0ad3c
--- /dev/null
+++ b/testing_and_setup/compass/ocean/hurricane/init_wetting_drying_template.xml
@@ -0,0 +1,8 @@
+
+
+
+
+
+
+
+
diff --git a/testing_and_setup/compass/ocean/hurricane/maxwaterlevel_template.xml b/testing_and_setup/compass/ocean/hurricane/maxwaterlevel_template.xml
new file mode 100644
index 0000000000..c40cdacf52
--- /dev/null
+++ b/testing_and_setup/compass/ocean/hurricane/maxwaterlevel_template.xml
@@ -0,0 +1,52 @@
+
+
+
+
+
+
+
+
+
+
+
+
+ output
+ maxWaterLevel.nc
+ 00-01-00_00:00:00
+ 01-01-01_00:00:00
+ truncate
+ pnetcdf
+ timeSeriesStatsCustomAMPKG
+ 00-00-00_01:00:00
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+ input;output
+ restarts/mpaso.rst.am.timeSeriesStatsCustom.$Y-$M-$D_$S.nc
+ output_interval
+ 0001-01-01_00:00:00
+ truncate
+ timeSeriesStatsCustomAMPKG
+ initial_only
+ stream:restart:output_interval
+
+
+