From b7d11c6215c492b5271c4342cd012eb0ffa6dce6 Mon Sep 17 00:00:00 2001 From: "Phillip J. Wolfram" Date: Mon, 9 Sep 2019 17:29:57 -0600 Subject: [PATCH 01/13] Adds Delaware coarse resolution flooding case --- .../build_mesh/config_base_mesh.xml | 17 ++ .../build_mesh/config_culled_mesh.xml | 29 ++++ .../build_mesh/config_driver_init.xml | 8 + .../build_mesh/define_base_mesh.py | 66 ++++++++ .../sandy/NOAA-COOPS_stations.txt | 28 ++++ .../sandy/USGS_stations.txt | 154 ++++++++++++++++++ .../sandy/config_driver.xml | 8 + .../sandy/config_forward.xml | 100 ++++++++++++ .../sandy/config_init.xml | 79 +++++++++ 9 files changed, 489 insertions(+) create mode 100755 testing_and_setup/compass/ocean/hurricane/USDEQU120at30cr10rr2WD/build_mesh/config_base_mesh.xml create mode 100644 testing_and_setup/compass/ocean/hurricane/USDEQU120at30cr10rr2WD/build_mesh/config_culled_mesh.xml create mode 100644 testing_and_setup/compass/ocean/hurricane/USDEQU120at30cr10rr2WD/build_mesh/config_driver_init.xml create mode 100755 testing_and_setup/compass/ocean/hurricane/USDEQU120at30cr10rr2WD/build_mesh/define_base_mesh.py create mode 100644 testing_and_setup/compass/ocean/hurricane/USDEQU120at30cr10rr2WD/sandy/NOAA-COOPS_stations.txt create mode 100644 testing_and_setup/compass/ocean/hurricane/USDEQU120at30cr10rr2WD/sandy/USGS_stations.txt create mode 100644 testing_and_setup/compass/ocean/hurricane/USDEQU120at30cr10rr2WD/sandy/config_driver.xml create mode 100644 testing_and_setup/compass/ocean/hurricane/USDEQU120at30cr10rr2WD/sandy/config_forward.xml create mode 100644 testing_and_setup/compass/ocean/hurricane/USDEQU120at30cr10rr2WD/sandy/config_init.xml 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..357fc37cd0 --- /dev/null +++ b/testing_and_setup/compass/ocean/hurricane/USDEQU120at30cr10rr2WD/build_mesh/config_base_mesh.xml @@ -0,0 +1,17 @@ + + + + + + + + + + jigsaw_to_MPAS.build_mesh + + + 20.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..1cfc56edd3 --- /dev/null +++ b/testing_and_setup/compass/ocean/hurricane/USDEQU120at30cr10rr2WD/build_mesh/config_culled_mesh.xml @@ -0,0 +1,29 @@ + + + + + + + + + + + + geometric_data + + + + + 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..b4d628b894 --- /dev/null +++ b/testing_and_setup/compass/ocean/hurricane/USDEQU120at30cr10rr2WD/build_mesh/define_base_mesh.py @@ -0,0 +1,66 @@ +#!/usr/bin/env python +''' +name: define_base_mesh +authors: Steven Brus, Phillip J. Wolfram + +This function specifies the resolution for a coastal refined mesh for Delaware Bay. +It contains the following resolution resgions: + 1) a QU 120km global background resolution + 2) 10km refinement region from the coast to past the shelf-break from North Carolina to New Hampshire + 3) 5km refinement region from the coast up to the shelf-break from Virginia to past Long Island, New York + 4) 2km refinement region inside Delaware Bay + +''' +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/NOAA-COOPS_stations.txt b/testing_and_setup/compass/ocean/hurricane/USDEQU120at30cr10rr2WD/sandy/NOAA-COOPS_stations.txt new file mode 100644 index 0000000000..64e9d20c8d --- /dev/null +++ b/testing_and_setup/compass/ocean/hurricane/USDEQU120at30cr10rr2WD/sandy/NOAA-COOPS_stations.txt @@ -0,0 +1,28 @@ +-74.1417 40.6367 8519483 Bergen Point West Reach +-73.7650 40.8103 8516945 Kings Point +-71.9594 41.0483 8510560 Montauk +-74.0142 40.7006 8518750 The Battery +-74.4183 39.3550 8534720 Atlantic City +-74.8733 40.0800 8539094 Burlington, Delaware River +-74.9597 38.9678 8536110 Cape May +-75.3750 39.3050 8537121 Ship John Shoal +-75.0430 40.0119 8538886 Tacony-Palmyra Bridge +-75.5883 39.5817 8551762 Delaware City +-75.1192 38.7828 8557380 Lewes +-75.5733 39.5583 8551910 Reedy Point +-76.4816 38.9833 8575512 Annapolis +-76.5783 39.2667 8574680 Baltimore, Fort McHenry, Patapsco River +-76.0387 38.2204 8571421 Bishops Head +-76.0722 38.5742 8571892 Cambridge +-75.8100 39.5267 8573927 Chesapeake City +-75.0917 38.3283 8570283 Ocean City Inlet +-76.4508 38.3172 8577330 Solomons Island +-76.2450 39.2133 8573364 Tolchester Beach +-76.1133 36.9667 8638863 Chesapeake Bay Bridge Tunnel +-75.9884 37.1652 8632200 Kiptopeke +-76.4646 37.9954 8635750 Lewisetta +-76.3017 36.7783 8639348 Money Point +-76.3300 36.9467 8638610 Sewells Point +-75.6858 37.6078 8631044 Wachapreague +-76.2900 37.6161 8636580 Windmill Point +-76.4788 37.2265 8637689 Yorktown USCG Training Center diff --git a/testing_and_setup/compass/ocean/hurricane/USDEQU120at30cr10rr2WD/sandy/USGS_stations.txt b/testing_and_setup/compass/ocean/hurricane/USDEQU120at30cr10rr2WD/sandy/USGS_stations.txt new file mode 100644 index 0000000000..0aa8239ac8 --- /dev/null +++ b/testing_and_setup/compass/ocean/hurricane/USDEQU120at30cr10rr2WD/sandy/USGS_stations.txt @@ -0,0 +1,154 @@ +-73.65944 40.9991 SSS-CT-FFD-001WL Fairfield County +-73.41568 41.09979 SSS-CT-FFD-003WL Fairfield County +-73.36481 41.14762 SSS-CT-FFD-005WL Fairfield County +-73.36998 41.1231 SSS-CT-FFD-006WL Fairfield County +-73.08843 41.25325 SSS-CT-FFD-009WL Fairfield County +-73.10901 41.16316 SSS-CT-FFD-010WL Fairfield County +-73.10901 41.16316 SSS-CT-FFD-012WL Fairfield County +-72.529361 41.269194 SSS-CT-MSX-018WL Middlesex County +-72.352222 41.281111 SSS-CT-MSX-019WL Middlesex County +-72.352222 41.281111 SSS-CT-MSX-020WL Middlesex County +-72.90479 41.27221 SSS-CT-NHV-013WL New Haven County +-72.874936 41.367328 SSS-CT-NHV-014WL New Haven County +-72.663611 41.271778 SSS-CT-NHV-015WL New Haven County +-72.663611 41.271778 SSS-CT-NHV-017WL New Haven County +-72.820642 41.260402 SSS-CT-NHV-018WL New Haven County +-72.90479 41.27221 SSS-CT-NHV-019WL New Haven County +-73.049475 41.211286 SSS-CT-NHV-020WL New Haven County +-71.9846 41.32525 SSS-CT-NLD-015WL New London County +-71.9846 41.32525 SSS-CT-NLD-016WL New London County +-72.1954 41.32147 SSS-CT-NLD-018WL New London County +-72.277583 41.284278 SSS-CT-NLD-019WL New London County +-72.3458 41.3125 SSS-CT-NLD-022WL New London County +-72.277583 41.284278 SSS-CT-NLD-023WL New London County +-72.060917 41.31675 SSS-CT-NLD-025WL New London County +-72.0355 41.334972 SSS-CT-NLD-026WL New London County +-71.967917 41.380806 SSS-CT-NLD-027WL New London County +-71.967667 41.346667 SSS-CT-NLD-029WL New London County +-71.909472 41.344333 SSS-CT-NLD-030WL New London County +-75.397806 39.0585278 SSS-DE-KEN-002WL Kent County +-75.474972 39.325833 SSS-DE-KEN-051WL Kent County +-75.446 39.160528 SSS-DE-KEN-053WL Kent County +-75.609472 39.3092778 SSS-DE-NEW-001WL New Castle County +-75.163639 38.790222 SSS-DE-SUS-008WL Sussex County +-75.084194 38.694528 SSS-DE-SUS-010WL Sussex County +-75.062528 38.513667 SSS-DE-SUS-014WL Sussex County +-75.058139 38.454889 SSS-DE-SUS-015WL Sussex County +-75.161833 38.702694 SSS-DE-SUS-030WL Sussex County +-75.099906 38.625334 SSS-DE-SUS-032WL Sussex County +-75.211972 38.591639 SSS-DE-SUS-033WL Sussex County +-75.197 38.76871 SSS-DE-SUS-057WL Sussex County +-70.333889 43.544722 SSS-ME-CUM-002WL Cumberland County +-70.356111 43.446667 SSS-ME-YOR-001WL York County +-70.381944 43.462222 SSS-ME-YOR-002WL York County +-76.48505 38.976833 SSS-MD-ANN-001WL Anne Arundel County +-76.4763 38.968611 SSS-MD-ANN-003WL Anne Arundel County +-76.60608 39.283694 SSS-MD-BAL-001WL Baltimore County +-76.6058 39.286361 SSS-MD-BAL-003WL Baltimore County +-70.872545 42.816118 424858070522116 Essex County +-70.62355 41.65638 SSS-MA-BAR-023WL Barnstable County +-70.65186 41.60618 SSS-MA-BAR-024WL Barnstable County +-70.54864 41.55283 SSS-MA-BAR-025WL Barnstable County +-70.39557 41.62131 SSS-MA-BAR-026WL Barnstable County +-70.08976 41.66752 SSS-MA-BAR-027WL Barnstable County +-70.27664 41.63949 SSS-MA-BAR-028WL Barnstable County +-69.98068 41.80028 SSS-MA-BAR-031WL Barnstable County +-71.07168 41.51626 SSS-MA-BRI-016WL Bristol County +-70.84321 41.59597 SSS-MA-BRI-017WL Bristol County +-70.918315 42.420961 SSS-MA-ESS-037WL Essex County +-70.88681 42.51951 SSS-MA-ESS-038WL Essex County +-70.615021 42.658992 SSS-MA-ESS-039WL Essex County +-70.827774 42.683689 SSS-MA-ESS-040WL Essex County +-70.820202 42.816934 SSS-MA-ESS-041WL Essex County +-70.7893 42.23905 SSS-MA-NOR-036WL Norfolk County +-70.81283 41.65608 SSS-MA-PLY-018WL Plymouth County +-70.76535 41.71257 SSS-MA-PLY-019WL Plymouth County +-70.62779 41.74369 SSS-MA-PLY-021WL Plymouth County +-70.55507 41.92901 SSS-MA-PLY-032WL Plymouth County +-70.68291 41.98024 SSS-MA-PLY-033WL Plymouth County +-70.64636 42.08254 SSS-MA-PLY-034WL Plymouth County +-70.72573 42.20139 SSS-MA-PLY-035WL Plymouth County +-70.817 42.8982 SSS-NH-ROC-004WL Rockingham County +-70.8233 42.9218 SSS-NH-ROC-005WL Rockingham County +-74.46277778 39.55305556 SSS-NJ-ATL-005WL Atlantic County +-74.6275 39.288333 SSS-NJ-CPM-010WL Cape May County +-74.865556 38.936389 SSS-NJ-CPM-035WL Cape May County +-74.865556 38.936389 SSS-NJ-CPM-035WV Cape May County +-75.04055556 39.395278 SSS-NJ-CUM-020WL Cumberland County +-75.23666667 39.429167 SSS-NJ-CUM-025WL Cumberland County +-74.02317 40.74394 SSS-NJ-HUD-001WL Hudson County +-74.06606 40.79982 SSS-NJ-HUD-002WL Hudson County +-74.24687 40.45911 SSS-NJ-MID-001WL Middlesex County +-74.166309 40.926465 SSS-NJ-PAS-001WL Passaic County +-74.20512 40.64775 SSS-NJ-UNI-001WL Union County +-74.27177 40.59952 SSS-NJ-UNI-002WL Union County +-73.2630556 40.643333 403836073154775 Suffolk County +-73.926474 40.800589 404810735538063 New York County +-73.732617 40.94881 405658073433147 Westchester County +-74.011612 40.579998 SSS-NY-KIN-001WL Kings County +-73.988317 40.704581 SSS-NY-KIN-002WL Kings County +-73.989839 40.676878 SSS-NY-KIN-003WL Kings County +-73.530571 40.877911 SSS-NY-NAS-001WL Nassau County +-73.640682 40.582752 SSS-NY-NAS-004WL Nassau County +-73.458503 40.652383 SSS-NY-NAS-005WL Nassau County +-73.5636111 40.8875 SSS-NY-NAS-006WL Nassau County +-73.4633333 40.857222 SSS-NY-NAS-007WL Nassau County +-73.7101944 40.866222 SSS-NY-NAS-008WL Nassau County +-73.926331 40.877565 SSS-NY-NEW-001WL New York County +-73.858278 40.762294 SSS-NY-QUE-001WL Queens County +-73.836384 40.64533 SSS-NY-QUE-002WL Queens County +-73.828785 40.796509 SSS-NY-QUE-004WL Queens County +-73.822651 40.606152 SSS-NY-QUE-005WL Queens County +-74.059845 40.593884 SSS-NY-RIC-001WL Richmond County +-74.230339 40.501881 SSS-NY-RIC-003WL Richmond County +-74.12768 40.54345 SSS-NY-RIC-004WL Richmond County +-72.558281 41.012585 SSS-NY-SUF-001WL Suffolk County +-72.863197 40.964375 SSS-NY-SUF-002WL Suffolk County +-73.072273 40.946167 SSS-NY-SUF-003WL Suffolk County +-72.750252 40.787121 SSS-NY-SUF-004WL Suffolk County +-72.637741 40.916077 SSS-NY-SUF-005WL Suffolk County +-72.502853 40.848869 SSS-NY-SUF-006WL Suffolk County +-72.502998 40.893311 SSS-NY-SUF-008WL Suffolk County +-72.290298 41.001972 SSS-NY-SUF-009WL Suffolk County +-73.353035 40.900482 SSS-NY-SUF-011WL Suffolk County +-72.470741 40.990695 SSS-NY-SUF-014WL Suffolk County +-72.361443 41.10104 SSS-NY-SUF-015WL Suffolk County +-73.202158 40.634734 SSS-NY-SUF-018WL Suffolk County +-73.264861 40.659318 SSS-NY-SUF-019WL Suffolk County +-73.01338 40.749179 SSS-NY-SUF-021WL Suffolk County +-73.279903 40.685227 SSS-NY-SUF-022WL Suffolk County +-72.189097 40.944283 SSS-NY-SUF-023WL Suffolk County +-71.934378 41.073188 SSS-NY-SUF-024WL Suffolk County +-72.855502 40.746867 SSS-NY-SUF-026WL Suffolk County +-73.1503912 40.747598 SSS-NY-SUF-027WL Suffolk County +-73.719828 40.942755 SSS-NY-WES-001WL Westchester County +-73.78172 40.890385 SSS-NY-WES-003WL Westchester County +-73.73236 40.94759 SSS-NY-WES-100WL Westchester County +-75.37609 39.8328 SSS-PA-DEL-003WL Delaware County +-75.30225 39.85918 SSS-PA-DEL-005WL Delaware County +-75.22887 39.86531 SSS-PA-DEL-006WL Delaware County +-75.00025 40.03425 SSS-PA-PHI-013WL Philadelphia County +-75.20825 39.9338889 SSS-PA-PHI-014WL Philadelphia County +-75.164222 39.88758333 SSS-PA-PHI-016WL Philadelphia County +-71.449851 41.449111 412656071265946 Washington County +-71.2859 41.72614 SSS-RI-BRI-013WL Bristol County +-71.39166 41.68673 SSS-RI-KEN-011WL Kent County +-71.24 41.61965 SSS-RI-NEW-014WL Newport County +-71.19241 41.46496 SSS-RI-NEW-015WL Newport County +-71.85914 41.31029 SSS-RI-WAS-001WL Washington County +-71.60535 41.36532 SSS-RI-WAS-003WL Washington County +-71.76663 41.33482 SSS-RI-WAS-005WL Washington County +-71.64473 41.38102 SSS-RI-WAS-007WL Washington County +-71.51472 41.37726 SSS-RI-WAS-008WL Washington County +-71.41637 41.52809 SSS-RI-WAS-012WL Washington County +-75.40668 37.90318 SSS-VA-ACC-001WL Accomack County +-75.58983 37.73291 SSS-VA-ACC-002WL Accomack County +-75.786703 37.721747 SSS-VA-ACC-003WL Accomack County +-76.3172222222 37.0152777778 SSS-VA-HAM-002WL Hampton City +-76.309926 37.492592 SSS-VA-MAT-001WL Mathews County +-75.841713 37.444888 SSS-VA-NOR-001WL Northampton County +-76.0162777778 37.2649166667 SSS-VA-NOR-003WL Northampton County +-75.95017 37.127509 SSS-VA-NOR-004WL Northampton County +-76.08825 36.9068333333 SSS-VA-VAB-001WL Virginia Beach City +-76.319444 37.1106111 SSS-VA-YOR-003WL York County +-73.1575 40.643161 SSS-NY-SUF-017WL Suffolk County 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..765e605959 --- /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.xml b/testing_and_setup/compass/ocean/hurricane/USDEQU120at30cr10rr2WD/sandy/config_forward.xml new file mode 100644 index 0000000000..bad0d957cf --- /dev/null +++ b/testing_and_setup/compass/ocean/hurricane/USDEQU120at30cr10rr2WD/sandy/config_forward.xml @@ -0,0 +1,100 @@ + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + input.nc + + + input.nc + + + 0002_00:00:00 + + + points.nc + input + initial_only + single_file + forward;analysis + + + + +