Changeset 4926 for anuga_work/development/cocos/grid_data.py
 Timestamp:
 Jan 10, 2008, 4:48:13 PM (17 years ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

anuga_work/development/cocos/grid_data.py
r4846 r4926 19 19 copy_code_files(output_dir,__file__) 20 20 21 start_screen_catcher(output_dir) 21 start_screen_catcher(output_dir,verbose=True) 22 print'testing' 22 23 23 24 data_dir='/d/ehn/4/ehn/tsunami/data/bathymetry/Cocos/bathymetry/' … … 78 79 # 4.3m2 79 80 80 res = 10 81 res = 100 81 82 #G = Geospatial_data(file_name = data_dir_name + '.txt') 82 83 #G1=G.clip(poly1) … … 84 85 85 86 86 interior_polys = [[poly2,(grid2_res**2/2)*res],[poly3,(grid3_res**2/2)*res],[poly4,(grid4_res**2/2)*res]] 87 #interior_polys = [[poly2,(grid2_res**2/2)*res],[poly3,(grid3_res**2/2)*res],[poly4,(grid4_res**2/2)*res]] 88 interior_polys = [[poly2,(grid2_res**2)*2],[poly3,(grid3_res**2)*res],[poly4,(grid4_res**2)*res]] 87 89 #interior_polys = [[poly4,4.3*res]] 88 90 … … 94 96 # use_cache=True, 95 97 # verbose=True) 98 max_tri = (grid1_res**2/2)*res 99 print 'max_tri',max_tri 96 100 create_mesh_from_regions(poly1, 97 101 #create_mesh_from_regions(poly3, 98 102 boundary_tags={'side': [0,1,2,3]}, 99 maximum_triangle_area=(grid1_res**2/2)*res, 103 maximum_triangle_area=(grid1_res**2)*res, 104 # maximum_triangle_area=(grid1_res**2)*1, 100 105 # maximum_triangle_area=110*res, 101 106 interior_regions=interior_polys, … … 105 110 106 111 107 value, alpha = find_optimal_smoothing_parameter(data_file=data_dir_name_clip + '.txt', 108 alpha_list=[0.0001,0.001,0.01,0.1,1.0], 109 # alpha_list=[0.01], 110 boundary_poly=poly3, 111 mesh_file=mesh_file, 112 plot_name=None, 113 cache=True, 114 verbose=True) 112 #value, alpha = find_optimal_smoothing_parameter(data_file=data_dir_name_clip + '.txt', 113 ## alpha_list=[0.0001,0.001,0.01,0.1,1.0], 114 # alpha_list=[0.1], 115 # boundary_poly=poly3, 116 # mesh_file=mesh_file, 117 # plot_name=None, 118 # split_factor=0.1, 119 # cache=True, 120 # verbose=True) 115 121 116 print 'Setup computational domain', value, alpha122 #print 'Setup computational domain', value, alpha 117 123 118 124 #domain = cache(Domain, (meshes_dir_name), {'use_cache':True, 'verbose':True}, verbose=True) … … 126 132 use_cache = True, 127 133 verbose = True, 128 alpha = alpha) 134 # alpha = alpha) 135 alpha = 0.01) 129 136 print 'Finished Set quantity' 130 137 131 138 domain.set_name(sww_file) 132 139 domain.set_datadir(output_dir) 133 domain.set_default_order( 2) # Apply second order scheme140 domain.set_default_order(1) # Apply second order scheme 134 141 domain.set_minimum_storable_height(0.01) # Don't store anything less than 1cm 135 142 domain.set_store_vertices_uniquely(False) … … 150 157 timestep = None, 151 158 reduction = max, 159 # cellsize = grid4_res*sqrt(1), 152 160 cellsize = grid4_res*sqrt(res), 153 161 NODATA_value = 9999, … … 165 173 timestep = None, 166 174 reduction = max, 175 # cellsize = grid3_res*sqrt(1), 167 176 cellsize = grid3_res*sqrt(res), 168 177 NODATA_value = 9999, … … 180 189 timestep = None, 181 190 reduction = max, 182 cellsize = grid2_res*sqrt(res), 191 cellsize = grid2_res*sqrt(1), 192 # cellsize = grid2_res*sqrt(res), 183 193 NODATA_value = 9999, 184 194 easting_min = 224424.96, … … 195 205 timestep = None, 196 206 reduction = max, 207 # cellsize = grid1_res*sqrt(1), 197 208 cellsize = grid1_res*sqrt(res), 198 209 NODATA_value = 9999,
Note: See TracChangeset
for help on using the changeset viewer.