Merge pull request #182 from sbrus89/CONUS_mesh_fix

Minor clean up of NAEC60to30cr8 mesh This PR cleans up some resolution artifacts around the Bering Sea region.
E3SM-Project · Mar 8, 2019 · fe1019e · fe1019e
2 parents 7cc3fc0 + e884f4d
commit fe1019e
Show file tree

Hide file tree

Showing 2 changed files with 7 additions and 8 deletions.
diff --git a/testing_and_setup/compass/ocean/global_ocean/NAEC60to30cr8/init/define_base_mesh.py b/testing_and_setup/compass/ocean/global_ocean/NAEC60to30cr8/init/define_base_mesh.py
@@ -16,6 +16,7 @@ def cellWidthVsLatLon():
 
     #params["plot_box"] = np.array([-80,-30,30,90]) 
     params["plot_box"] = ct.Entire_Globe
+    params["plot_option"] = False
 
     print "***Gulf Coast***"
     params["region_box"] = ct.US_Gulf_Coast
@@ -49,8 +50,8 @@ def cellWidthVsLatLon():
 
     print "***Aleutian Islands (West)***"
     params["region_box"] = ct.Aleutian_Islands_W
-    params["lon_wrap"]  = True
     params["n_longest"] = 100
+    params["trans_start"] = 200.0*km
     cell_width,lon,lat = ct.coastal_refined_mesh(params,cell_width,lon,lat)
 
     print "****Aleutian Islands (East)***"
@@ -87,7 +88,6 @@ def cellWidthVsLatLon():
     params["restrict_box"] = ct.Empty
     params["trans_width"] = 600.0*km
     params["trans_start"] = 400.0*km
-    params["lon_wrap"] = False
     cell_width,lon,lat = ct.coastal_refined_mesh(params,cell_width,lon,lat)
 
     print "***Greenland***"

diff --git a/testing_and_setup/compass/ocean/jigsaw_to_MPAS/coastal_tools.py b/testing_and_setup/compass/ocean/jigsaw_to_MPAS/coastal_tools.py
@@ -85,13 +85,13 @@
 Bering_Sea_E = {"include": [np.array([-180.0, -168.5, 56.00, 64.0])],
                 "exclude": []}
 
-Bering_Sea_W = {"include": [np.array([161.83, 180.0, 51.13, 66.63])],
+Bering_Sea_W = {"include": [np.array([161.90, 180.0, 57.85, 66.63])],
                 "exclude": []}
 
 Aleutian_Islands_E = {"include": [np.array([-180.0, -151.79, 49.54, 58.83])],
                       "exclude": [np.array([-173.16, -164.37, 55.81, 57.55])]}
 
-Aleutian_Islands_W = {"include": [np.array([171.8, 180.0, 49.54, 58.83])],
+Aleutian_Islands_W = {"include": [np.array([164.9, 180.0, 50.77, 56.31])],
                       "exclude": [np.array([178.5, 179.5, 51.25, 51.75])]}
 
 
@@ -161,10 +161,10 @@
                                              [-63.47, 53.76]])]}
 Bering_Sea_restrict = {"include": [],
                        "exclude": [np.array([[143.46, 51.79],
-                                             [156.55, 51.97],
-                                             [165.97, 62.99],
+                                             [155.65, 50.13],
+                                             [166.23, 63.78],
                                              [148.63, 62.30]]),
-                                   np.array([[161.14, 68.41],
+                                   np.array([[154.36, 68.22],
                                              [-173.80, 65.94],
                                              [-161.81, 72.02],
                                              [163.64, 73.70]])]}
@@ -195,7 +195,6 @@
     "mesh_type": 'EC',
     "dx_max_global": 30.0 * km,
     "dx_min_global": 10.0 * km,
-    "lon_wrap": False,
 
     # Coastal mesh parameters
     "dx_min_coastal": 10.0 * km,