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) return cell_width / 1000, lon, lat
def cellWidthVsLatLon(): km = 1000.0 params = ct.default_params params["mesh_type"] = "QU" params["dx_max_global"] = 300.0 * km ct.Delaware_Region["include"][0] = np.array([-77, -69.8, 37.25, 40.25]) params["region_box"] = ct.Delaware_Region params["plot_box"] = ct.Western_Atlantic params["dx_min_coastal"] = 20.0 * km params["trans_width"] = 600.0 * km params["trans_start"] = 100.0 * km cell_width, lon, lat = ct.coastal_refined_mesh(params) ct.Delaware_Bay["include"][0] = np.array( [-75.61903, -74.7, 38.5, 40.312747]) params["region_box"] = ct.Delaware_Bay params["plot_box"] = ct.Delaware params["restrict_box"] = ct.Delaware_restrict params["dx_min_coastal"] = 3.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
def cellWidthVsLatLon(): km = 1000.0 params = ct.default_params params["mesh_type"] = "QU" params["dx_max_global"] = 120.0 * km 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) 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) 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
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
def cellWidthVsLatLon(): km = 1000.0 params = ct.default_params print("****QU120 background mesh and 10km refinement from NC to NH****") params["mesh_type"] = "QU" params["dx_max_global"] = 120.0 * km 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) print("****5km refinement along coast from VA to NY****") 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("****2km refinement inside Delaware Bay****") 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
def cellWidthVsLatLon(): km = 1000.0 params = ct.default_params SFtoLA = { "include": [np.array([-124.0, -117.5, 34.2, 38.0])], # SF to LA "exclude": [np.array([-122.1, -120.8, 37.7, 39.2])] } # SF Bay Delta WestCoast = np.array([-136.0, -102.0, 22.0, 51]) print("****QU120 background mesh and 300m refinement from SF to LA****") params["mesh_type"] = "QU" params["dx_max_global"] = 120.0 * km params["region_box"] = SFtoLA params["plot_box"] = WestCoast params["dx_min_coastal"] = 3.0 * km params["trans_width"] = 100.0 * km params["trans_start"] = 30.0 * km cell_width, lon, lat = ct.coastal_refined_mesh(params) return cell_width / 1000, lon, lat
def cellWidthVsLatLon(): km = 1000.0 params = ct.default_params params["dx_min_coastal"] = 8.0*km params["trans_width"] = 600.0*km params["trans_start"] = 400.0*km params["mesh_type"] = "EC" #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 params["restrict_box"] = ct.Gulf_restrict cell_width, lon, lat = ct.coastal_refined_mesh(params) print("***East Coast***") params["region_box"] = ct.US_East_Coast params["restrict_box"] = ct.East_Coast_restrict cell_width,lon,lat = ct.coastal_refined_mesh(params,cell_width,lon,lat) print("***Caribbean***") params["region_box"] = ct.Caribbean params["restrict_box"] = ct.Caribbean_restrict params["trans_width"] = 400.0*km params["trans_start"] = 300.0*km params["n_longest"] = 50 cell_width,lon,lat = ct.coastal_refined_mesh(params,cell_width,lon,lat) print("***West Coast***") params["region_box"] = ct.US_West_Coast params["restrict_box"] = ct.Empty params["trans_width"] = 600.0*km params["trans_start"] = 400.0*km params["n_longest"] = 10 cell_width,lon,lat = ct.coastal_refined_mesh(params,cell_width,lon,lat) print("***Alaska***") params["region_box"] = ct.Alaska cell_width,lon,lat = ct.coastal_refined_mesh(params,cell_width,lon,lat) print("***Aleutian Islands (West)***") params["region_box"] = ct.Aleutian_Islands_W 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)***") params["region_box"] = ct.Aleutian_Islands_E cell_width,lon,lat = ct.coastal_refined_mesh(params,cell_width,lon,lat) print("***Bering Sea (East)****") params["region_box"] = ct.Bering_Sea_E params["trans_start"] = 400.0*km params["trans_width"] = 600.0*km params["n_longest"] = 10 cell_width,lon,lat = ct.coastal_refined_mesh(params,cell_width,lon,lat) print("***Bering Sea (West)***") params["region_box"] = ct.Bering_Sea_W params["restrict_box"] = ct.Bering_Sea_restrict params["trans_start"] = 450.0*km params["trans_width"] = 600.0*km params["n_longest"] = 10 cell_width,lon,lat = ct.coastal_refined_mesh(params,cell_width,lon,lat) print("***Hawaii***") params["region_box"] = ct.Hawaii params["restrict_box"] = ct.Empty params["trans_start"] = 100*km params["trans_width"] = 200*km cell_width,lon,lat = ct.coastal_refined_mesh(params,cell_width,lon,lat) print("***Newfoundland***") Newfoundland = {"include":[np.array([-65.0,-50.0,44.0,60.0]), np.array([-65.5,-64.5,61.0,62.0])], "exclude":[]} params["region_box"] = Newfoundland params["restrict_box"] = ct.Empty 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("***Greenland***") params["region_box"] = ct.Greenland params["restrict_box"] = ct.Empty params["trans_width"] = 600.0*km params["trans_start"] = 275.0*km cell_width, lon, lat = ct.coastal_refined_mesh(params,cell_width,lon,lat) print("***Labrador Sea***") params["region_box"] = ct.Empty params["restrict_box"] = ct.Empty params["point_list"] = [np.array([-50.0,55.0])] 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("***Central America (West Coast)***") Central_America = {"include":[np.array([[-110.26,20.69], [-87.84, 8.94 ], [-84.55, 12.03], [-104.26,23.11]]), np.array([[-88.02, 10.47], [-81.53, 6.14], [-81.45, 8.07], [-84.80, 11.51]]), np.array([[-81.92, 7.76], [-76.84, 4.51], [-77.41, 8.22], [-79.23, 9.28]])], "exclude":[]} params["region_box"] = Central_America params["restrict_box"] = ct.Empty params["point_list"] = None 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) return cell_width / km, lon, lat
def cellWidthVsLatLon(): km = 1000.0 params = ct.default_params # Use 1.0 degree for fast scoping, and 0.1 for the final creation: params["ddeg"] = 0.2 #params["ddeg"] = 1.0 #params["ddeg"] = 0.5 params["dx_min_coastal"] = 12.0 * km params["trans_start"] = 400.0 * km params["trans_width"] = 600.0 * km params["n_longest"] = 20 params["mesh_type"] = "EC" params["plot_box"] = ct.Entire_Globe params["plot_option"] = False print("***Gulf Coast***") params["region_box"] = ct.US_Gulf_Coast params["restrict_box"] = ct.Gulf_restrict params["trans_start"] = 400.0 * km params["trans_width"] = 600.0 * km cell_width, lon, lat = ct.coastal_refined_mesh(params) print("***East Coast***") params["region_box"] = ct.US_East_Coast params["restrict_box"] = ct.East_Coast_restrict params["trans_start"] = 400.0 * km params["trans_width"] = 600.0 * km cell_width, lon, lat = ct.coastal_refined_mesh(params, cell_width, lon, lat) print("***Gulf Stream extension***") params["restrict_box"] = ct.Empty params["trans_width"] = 600.0 * km params["region_box"] = { "include": [np.array([-78, -70, 40, 42])], "exclude": [] } params["trans_start"] = 700.0 * km cell_width, lon, lat = ct.coastal_refined_mesh(params, cell_width, lon, lat) params["region_box"] = { "include": [np.array([-74, -50, 42, 55])], "exclude": [] } params["trans_start"] = 1100.0 * km cell_width, lon, lat = ct.coastal_refined_mesh(params, cell_width, lon, lat) print("***Gulf Coast spot***") params["restrict_box"] = ct.Empty params["trans_width"] = 100.0 * km params["region_box"] = { "include": [np.array([-98, -94, 28, 30])], "exclude": [] } params["trans_start"] = 400.0 * km cell_width, lon, lat = ct.coastal_refined_mesh(params, cell_width, lon, lat) print("***Greenland***") params["region_box"] = ct.Greenland params["restrict_box"] = ct.Empty params["trans_width"] = 600.0 * km params["trans_start"] = 275.0 * km cell_width, lon, lat = ct.coastal_refined_mesh(params, cell_width, lon, lat) print("***Greenland-Iceland smooth kink***") params["region_box"] = { "include": [np.array([-45, -40, 59, 62])], "exclude": [] } params["trans_start"] = 700.0 * km cell_width, lon, lat = ct.coastal_refined_mesh(params, cell_width, lon, lat) print("***Hudson Bay***") params["region_box"] = { "include": [np.array([-100, -70, 50, 65])], "exclude": [] } params["restrict_box"] = { "include": [np.array([-100, -70, 50, 65])], "exclude": [] } params["trans_start"] = 600.0 * km params["trans_width"] = 100.0 * km cell_width, lon, lat = ct.coastal_refined_mesh(params, cell_width, lon, lat) print("***Northern Canada***") params["region_box"] = { "include": [np.array([-168, -59, 67, 85])], "exclude": [] } params["restrict_box"] = ct.Empty params["trans_start"] = 275.0 * km params["trans_width"] = 600.0 * km params["n_longest"] = 20 cell_width, lon, lat = ct.coastal_refined_mesh(params, cell_width, lon, lat) print("***Alaska***") params["region_box"] = ct.Alaska params["trans_start"] = 400.0 * km params["trans_width"] = 600.0 * km cell_width, lon, lat = ct.coastal_refined_mesh(params, cell_width, lon, lat) print("***Crook at Alaska/Canada***") params["region_box"] = { "include": [np.array([-144, -127, 68.5, 71])], "exclude": [] } params["restrict_box"] = ct.Empty params["trans_start"] = 550.0 * km params["trans_width"] = 600.0 * km params["n_longest"] = 20 cell_width, lon, lat = ct.coastal_refined_mesh(params, cell_width, lon, lat) print("***Caribbean***") params["region_box"] = ct.Caribbean params["restrict_box"] = ct.Caribbean_restrict params["trans_width"] = 400.0 * km params["trans_start"] = 300.0 * km params["n_longest"] = 50 cell_width, lon, lat = ct.coastal_refined_mesh(params, cell_width, lon, lat) print("***West Coast***") params["region_box"] = ct.US_West_Coast params["restrict_box"] = ct.Empty params["trans_width"] = 600.0 * km params["trans_start"] = 400.0 * km params["n_longest"] = 10 cell_width, lon, lat = ct.coastal_refined_mesh(params, cell_width, lon, lat) print("***Aleutian Islands (West)***") params["region_box"] = ct.Aleutian_Islands_W 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)***") params["region_box"] = ct.Aleutian_Islands_E cell_width, lon, lat = ct.coastal_refined_mesh(params, cell_width, lon, lat) print("***Bering Sea (East)****") params["region_box"] = ct.Bering_Sea_E params["trans_start"] = 400.0 * km params["trans_width"] = 600.0 * km params["n_longest"] = 10 cell_width, lon, lat = ct.coastal_refined_mesh(params, cell_width, lon, lat) print("***Bering Sea (West)***") params["region_box"] = ct.Bering_Sea_W params["restrict_box"] = ct.Bering_Sea_restrict params["trans_start"] = 450.0 * km params["trans_width"] = 600.0 * km params["n_longest"] = 10 cell_width, lon, lat = ct.coastal_refined_mesh(params, cell_width, lon, lat) print("***Newfoundland***") Newfoundland = { "include": [ np.array([-65.0, -50.0, 44.0, 60.0]), np.array([-65.5, -64.5, 61.0, 62.0]) ], "exclude": [] } params["region_box"] = Newfoundland params["restrict_box"] = ct.Empty 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("***Labrador Sea***") params["region_box"] = ct.Empty params["restrict_box"] = ct.Empty params["point_list"] = [np.array([-50.0, 55.0])] 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("***Central America (West Coast)***") Central_America = { "include": [ np.array([[-110.26, 20.69], [-87.84, 8.94], [-84.55, 12.03], [-104.26, 23.11]]), np.array([[-88.02, 10.47], [-81.53, 6.14], [-81.45, 8.07], [-84.80, 11.51]]), np.array([[-81.92, 7.76], [-76.84, 4.51], [-77.41, 8.22], [-79.23, 9.28]]) ], "exclude": [] } params["region_box"] = Central_America params["restrict_box"] = ct.Empty params["point_list"] = None 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) # coast of Maine params["dx_min_coastal"] = 6.0 * km params["trans_start"] = 500.0 * km params["trans_width"] = 300.0 * km params["n_longest"] = 20 params["region_box"] = { "include": [np.array([-70, -68, 43, 45])], "exclude": [] } params["restrict_box"] = { "include": [], "exclude": [np.array([-73, -60, 46, 50])] } params["point_list"] = [[-68.8, 43.8]] cell_width, lon, lat = ct.coastal_refined_mesh(params, cell_width, lon, lat) # coast of Maine params["dx_min_coastal"] = 2.0 * km params["trans_start"] = 150.0 * km params["trans_width"] = 300.0 * km params["n_longest"] = 20 params["region_box"] = { "include": [np.array([-70, -68, 43, 45])], "exclude": [] } params["restrict_box"] = { "include": [], "exclude": [np.array([-73, -50, 46, 50])] } params["point_list"] = [[-68.8, 43.8]] cell_width, lon, lat = ct.coastal_refined_mesh(params, cell_width, lon, lat) return cell_width / km, lon, lat