Changeset 4856
- Timestamp:
- Nov 27, 2007, 4:18:39 PM (17 years ago)
- Files:
-
- 1 added
- 33 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_core/documentation/user_manual/demos/cairns/GetTimeseries.py
r4087 r4856 21 21 22 22 # nominate directory location of sww file with associated attribute 23 production_dirs = {'slide': 'Slide'} 23 production_dirs = {#'slide': 'Slide', 24 'fixed_wave': 'Fixed Wave'} 24 25 25 26 # Generate figures … … 36 37 reportname = '', 37 38 plot_quantity = ['stage', 'speed'], 39 generate_fig = False, 38 40 surface = False, 39 41 time_min = None, -
anuga_core/documentation/user_manual/demos/cairns/project.py
r4854 r4856 7 7 import sys 8 8 from time import localtime, strftime, gmtime 9 from anuga.utilities.polygon import read_polygon, plot_polygons _points, \9 from anuga.utilities.polygon import read_polygon, plot_polygons, \ 10 10 polygon_area, is_inside_polygon 11 11 -
anuga_core/documentation/user_manual/demos/cairns/runcairns.py
r4853 r4856 39 39 # Define scenario as either slide or fixed_wave. 40 40 #------------------------------------------------------------------------------ 41 scenario = 'slide'42 #scenario = 'fixed_wave'41 #scenario = 'slide' 42 scenario = 'fixed_wave' 43 43 44 44 if os.access(scenario, os.F_OK) == 0: -
anuga_core/source/anuga/abstract_2d_finite_volumes/test_neighbour_mesh.py
r4851 r4856 973 973 974 974 975 from anuga.utilities.polygon import plot_polygons _points975 from anuga.utilities.polygon import plot_polygons 976 976 977 977 # First do the continuous version of mesh -
anuga_core/source/anuga/fit_interpolate/interpolate.py
r4853 r4856 558 558 559 559 if verbose is True: 560 from anuga.utilities.polygon import plot_polygons _points560 from anuga.utilities.polygon import plot_polygons 561 561 #out_interp_pts = take(interpolation_points,[indices]) 562 562 title = 'Interpolation points fall outside specified mesh' 563 plot_polygons _points([mesh_boundary_polygon,interpolation_points,out_interp_pts],563 plot_polygons([mesh_boundary_polygon,interpolation_points,out_interp_pts], 564 564 ['line','point','outside'],figname='points_boundary_out',label=title,verbose=verbose) 565 565 … … 574 574 # Plot boundary and interpolation points 575 575 if verbose is True: 576 from anuga.utilities.polygon import plot_polygons _points576 from anuga.utilities.polygon import plot_polygons 577 577 title = 'Interpolation function: Polygon and interpolation points' 578 plot_polygons _points([mesh_boundary_polygon,interpolation_points],578 plot_polygons([mesh_boundary_polygon,interpolation_points], 579 579 ['line','point'],figname='points_boundary',label=title,verbose=verbose) 580 580 -
anuga_core/source/anuga/utilities/polygon.py
r4852 r4856 341 341 return abs(poly_area/2) 342 342 343 def plot_polygons _points(polygons_points, style=None,343 def plot_polygons(polygons_points, style=None, 344 344 figname=None, label=None, verbose=False): 345 345 -
anuga_work/production/broome_2006/get_timeseries.py
r4581 r4856 18 18 19 19 # bulli event 20 production_dirs = {'2007070 2_050939_run_final_4.9_exmouth_nbartzis': ''}20 production_dirs = {'20070703_061553_run_final_0_dampier_nbartzis': ''} 21 21 22 # 2007070 2_050939_run_final_4.9_exmouth_nbartzisHAT exmouth23 # 20070 615_062730_run_final_0_exmouth_nbartzis# MSL exmouth24 # 20070 618_064014_run_final_4.9_onslow_nbartzis# HAT onslow25 # 20070 622_004008_run_final_0_onslow_nbartzis# MSL onslow26 # 20070 615_063137_run_final_4.9_dampier_nbartzis# HAT dampier27 # 20070 621_002922_run_final_0_dampier_nbartzis# MSL dampier22 # 20070703_053716_run_final_4.9_exmouth_nbartzis HAT exmouth 23 # 20070703_053608_run_final_0_exmouth_nbartzis # MSL exmouth 24 # 20070703_053822_run_final_4.9_onslow_nbartzis # HAT onslow 25 # 20070703_061501_run_final_0_onslow_nbartzis # MSL onslow 26 # 20070703_062312_run_final_4.9_dampier_nbartzis # HAT dampier 27 # 20070703_061553_run_final_0_dampier_nbartzis # MSL dampier 28 28 29 29 is_parallel = True -
anuga_work/production/broome_2006/project.py
r4853 r4856 7 7 import sys 8 8 from time import localtime, strftime, gmtime 9 from anuga.utilities.polygon import read_polygon, plot_polygons _points, is_inside_polygon, number_mesh_triangles9 from anuga.utilities.polygon import read_polygon, plot_polygons, is_inside_polygon, number_mesh_triangles 10 10 #from anuga.coordinate_transforms.redfearn import degminsec2decimal_degrees, convert_points_from_latlon_to_utm 11 11 from anuga.utilities.system_tools import get_user_name … … 126 126 # Domain definitions 127 127 ############################### 128 from anuga.utilities.polygon import read_polygon, plot_polygons _points, polygon_area, is_inside_polygon128 from anuga.utilities.polygon import read_polygon, plot_polygons, polygon_area, is_inside_polygon 129 129 130 130 poly_all = read_polygon(polygons_dir+'extent_small.csv') -
anuga_work/production/broome_2006/project_urs.py
r4853 r4856 6 6 from os import sep, environ, getenv, getcwd,umask 7 7 from os.path import expanduser, basename 8 from anuga.utilities.polygon import read_polygon, plot_polygons _points, polygon_area, is_inside_polygon, number_mesh_triangles8 from anuga.utilities.polygon import read_polygon, plot_polygons, polygon_area, is_inside_polygon, number_mesh_triangles 9 9 import sys 10 10 from anuga.coordinate_transforms.redfearn import degminsec2decimal_degrees … … 152 152 # Domain definitions 153 153 ############################### 154 from anuga.utilities.polygon import read_polygon, plot_polygons _points, polygon_area, is_inside_polygon154 from anuga.utilities.polygon import read_polygon, plot_polygons, polygon_area, is_inside_polygon 155 155 156 156 poly_all = read_polygon(polygons_dir+'extent_small.csv') … … 160 160 #p2 = read_polygon(polygons_dir+'broome_north_coast_inside_extent.csv') 161 161 #label=None 162 #plot_polygons _points([poly_all,point,p2],['line','point','outside'],figname='test',label=label)163 #plot_polygons _points([poly_all,point],figname='test')162 #plot_polygons([poly_all,point,p2],['line','point','outside'],figname='test',label=label) 163 #plot_polygons([poly_all,point],figname='test') 164 164 #poly_all = read_polygon(polygons_dir+'extent.csv') 165 165 res_poly_all = 150000*res_factor -
anuga_work/production/busselton_2006/project.py
r4853 r4856 7 7 import sys 8 8 from time import localtime, strftime, gmtime 9 from anuga.utilities.polygon import read_polygon, plot_polygons _points, is_inside_polygon, number_mesh_triangles9 from anuga.utilities.polygon import read_polygon, plot_polygons, is_inside_polygon, number_mesh_triangles 10 10 #from anuga.coordinate_transforms.redfearn import degminsec2decimal_degrees, convert_points_from_latlon_to_utm 11 11 from anuga.utilities.system_tools import get_user_name … … 129 129 res_busselton2 = 500 130 130 131 #plot_polygons _points([polyAll,poly_broome1,poly_broome2,poly_broome3],figname='boundingpoly2',verbose=False)131 #plot_polygons([polyAll,poly_broome1,poly_broome2,poly_broome3],figname='boundingpoly2',verbose=False) 132 132 133 133 interior_regions = [[poly_busselton1,res_busselton1],[poly_busselton2,res_busselton2]] -
anuga_work/production/dampier_2006/project.py
r4853 r4856 6 6 import sys 7 7 from time import localtime, strftime, gmtime 8 from anuga.utilities.polygon import read_polygon, plot_polygons _points, polygon_area, is_inside_polygon, number_mesh_triangles8 from anuga.utilities.polygon import read_polygon, plot_polygons, polygon_area, is_inside_polygon, number_mesh_triangles 9 9 from anuga.coordinate_transforms.redfearn import degminsec2decimal_degrees, convert_from_latlon_to_utm 10 10 from anuga.utilities.system_tools import get_user_name … … 237 237 238 238 239 #plot_polygons _points([poly_dampier,poly_karratha,poly_karratha_town,poly_delambre,239 #plot_polygons([poly_dampier,poly_karratha,poly_karratha_town,poly_delambre, 240 240 # poly_coast,poly_NWislands,poly_island0,poly_island1,poly_island2, 241 241 # poly_island3,poly_island4,poly_island5,poly_island6, -
anuga_work/production/dampier_2006/project_urs.py
r4853 r4856 4 4 from os import sep, environ, getenv, getcwd,umask 5 5 from os.path import expanduser, basename 6 from anuga.utilities.polygon import read_polygon, plot_polygons _points, polygon_area, is_inside_polygon, number_mesh_triangles6 from anuga.utilities.polygon import read_polygon, plot_polygons, polygon_area, is_inside_polygon, number_mesh_triangles 7 7 import sys 8 8 from anuga.coordinate_transforms.redfearn import degminsec2decimal_degrees … … 243 243 244 244 245 #plot_polygons _points([poly_dampier,poly_karratha,poly_karratha_town,poly_delambre,245 #plot_polygons([poly_dampier,poly_karratha,poly_karratha_town,poly_delambre, 246 246 # poly_coast,poly_NWislands,poly_island0,poly_island1,poly_island2, 247 247 # poly_island3,poly_island4,poly_island5,poly_island6, -
anuga_work/production/dampier_2006/report/anuga_setup.tex
r4479 r4856 38 38 39 39 40 \begin{figure}[h bt]40 \begin{figure}[h] 41 41 42 \centerline{ \includegraphics[scale=0.5]{../report_figures/dampier_model.jpg}} 43 42 %\centerline{ \includegraphics[scale=0.5]{../report_figures/dampier_model.jpg}} 43 \centerline{ \includegraphics[scale=0.3]{../report_figures/model.jpg}} 44 44 45 \caption{Study area for Dampier scenario highlighting four regions of increased refinement. 45 Region 1: Surrounds onshore gas infrastructure at Dampier with a cell area of 500 m$^2$ (lateral accuracy 30 m) - not shown here.46 Region 1: Surrounds onshore gas infrastructure at Dampier with a cell area of 500 m$^2$ (lateral accuracy 30 m). 46 47 Region 2: Surrounds the gas pipelines with a cell area of 2000 m$^2$ (lateral accuracy 100 m). 47 48 Region 3: Surrounds the coastal regions with a cell area of 10000 m$^2$ (lateral accuracy m). … … 70 71 results are therefore likely to be over estimates. 71 72 72 The following table summarises the modelling parameters; 73 74 75 \begin{table} 76 \begin{center} 77 \caption{Parameters used in ANUGA for the Dampier scenario.} 78 \begin{tabular}{|l|l|l|}\hline 79 Mesh & & \\ \hline 80 & resolution in Region 1 & 500 m$^2$ \\ \hline 81 & resolution in Region 2 & 2000 m$^2$ \\ \hline 82 & resolution in Region 3 & 2000 m$^2$ \\ \hline 83 & remaining resolution & 100 000 m$^2$ \\ \hline 84 Model parameters & & \\ \hline 85 & friction & 0 \\ \hline 86 & minimum stored height & 0.1 m \\ \hline 87 \end{tabular} 88 \end{center} 89 \end{table} 73 %The following table summarises the modelling parameters; 74 % 75 %\begin{table} 76 %\begin{center} 77 %\caption{Parameters used in ANUGA for the Dampier scenario.} 78 %\begin{tabular}{|l|l|l|}\hline 79 %Mesh & & \\ \hline 80 %& resolution in Region 1 & 500 m$^2$ \\ \hline 81 %& resolution in Region 2 & 2000 m$^2$ \\ \hline 82 %& resolution in Region 3 & 2000 m$^2$ \\ \hline 83 %& remaining resolution & 100 000 m$^2$ \\ \hline 84 %Model parameters & & \\ \hline 85 %& friction & 0 \\ \hline 86 %& minimum stored height & 0.1 m \\ \hline 87 %\end{tabular} 88 %\end{center} 89 %\end{table} -
anuga_work/production/dampier_2006/report/mw9_map.tex
r4479 r4856 1 1 \begin{sidewaysfigure} 2 2 %\begin{figure} 3 \centerline{ \includegraphics[scale=0.4]{../report_figures/max_inundation_high_tide_v 2.jpg}}3 \centerline{ \includegraphics[scale=0.4]{../report_figures/max_inundation_high_tide_v3.jpg}} 4 4 \caption{Maximum inundation map for the HAT scenario for Dampier region. Data: WA DLI, DPI and AHO.} 5 5 \label{fig:HAT_max_inundation} -
anuga_work/production/exmouth_2006/export_results.py
r4479 r4856 5 5 from os import sep 6 6 7 time_dir = '' # HAT 1 in Dampier 10000 yr 8 #time_dir = '' # MSL 1 in Dampier 10000 yr 9 #time_dir = '' # HAT 1 in Broome 10000 yr 10 #time_dir = '' # MSL 1 in Broome 10000 yr 11 #time_dir = '' # HAT 1 in Pt Hedland 10000 yr 12 #time_dir = '' # MSL 1 in Pt Hedland 10000 yr 13 #time_dir = '' # HAT 1 in Onslow 10000 yr 14 #time_dir = '' # MSL 1 in Onslow 10000 yr 15 #time_dir = '' # HAT 1 in Exmouth 10000 yr 16 #time_dir = '' # MSL 1 in Exmouth 10000 yr 7 #time_dir = '20070618_035101_run_final_1.4_dampier_nbartzis' # HAT 1 in Dampier 10000 yr 8 #time_dir = '20070618_035114_run_final_0_dampier_nbartzis' # MSL 1 in Dampier 10000 yr 9 #time_dir = '20070619_042043_run_final_1.4_onslow_nbartzis' # HAT 1 in Onslow 10000 yr 10 #time_dir = '20070619_042140_run_final_0_onslow_nbartzis' # MSL 1 in Onslow 10000 yr 11 #time_dir = '20070619_042542_run_final_1.4_exmouth_nbartzis' # HAT 1 in Exmouth 10000 yr 12 time_dir = '20070619_042417_run_final_0_exmouth_nbartzis' # MSL 1 in Exmouth 10000 yr 17 13 cellsize = 25 18 14 timestep = None … … 37 33 38 34 if which_var == 2: # Depth 39 outname = name + '_depth '35 outname = name + '_depth_north' 40 36 quantityname = 'stage-elevation' 41 37 42 38 if which_var == 3: # Speed 43 outname = name + '_speed '39 outname = name + '_speed_north' 44 40 quantityname = '(xmomentum**2 + ymomentum**2)**0.5/(stage-elevation+1.e-6/(stage-elevation))' #Speed 45 41 -
anuga_work/production/exmouth_2006/project.py
r4853 r4856 5 5 from os import sep, environ, getenv, getcwd,umask 6 6 from os.path import expanduser, basename 7 from anuga.utilities.polygon import read_polygon, plot_polygons _points, polygon_area, is_inside_polygon, number_mesh_triangles7 from anuga.utilities.polygon import read_polygon, plot_polygons, polygon_area, is_inside_polygon, number_mesh_triangles 8 8 import sys 9 9 from anuga.coordinate_transforms.redfearn import degminsec2decimal_degrees … … 154 154 # Domain definitions 155 155 ############################### 156 from anuga.utilities.polygon import read_polygon, plot_polygons _points, polygon_area, is_inside_polygon156 from anuga.utilities.polygon import read_polygon, plot_polygons, polygon_area, is_inside_polygon 157 157 158 158 poly_all = read_polygon(polygons_dir+'extent_more_points.csv') … … 222 222 ################################################################### 223 223 224 # exporting asc grid 224 # exporting asc grid - Exmouth town area 225 225 eastingmin = 200663 226 226 eastingmax = 207500 227 227 northingmax = 7574000 228 228 northingmin = 7564500 229 230 # exporting asc grid - carparks etc on NW peninsula 231 e_min_area = 192700 232 e_max_area = 203755 233 n_max_area = 7587330 234 n_min_area = 7579460 -
anuga_work/production/galle_2007/project.py
r4853 r4856 8 8 import sys 9 9 from time import localtime, strftime, gmtime 10 from anuga.utilities.polygon import read_polygon, plot_polygons _points, is_inside_polygon, \10 from anuga.utilities.polygon import read_polygon, plot_polygons, is_inside_polygon, \ 11 11 number_mesh_triangles, polygon_area 12 12 from anuga.utilities.system_tools import get_user_name … … 133 133 [poly_canal, canal_res]] 134 134 135 #plot_polygons _points([bounding_polygon,poly_intermediate,\135 #plot_polygons([bounding_polygon,poly_intermediate,\ 136 136 # poly_shallow,poly_coast,poly_canal],\ 137 137 # figname='boundingpoly',verbose=False) -
anuga_work/production/gold_coast_2007/project.py
r4853 r4856 7 7 import sys 8 8 from time import localtime, strftime, gmtime 9 from anuga.utilities.polygon import read_polygon, plot_polygons _points, polygon_area, is_inside_polygon9 from anuga.utilities.polygon import read_polygon, plot_polygons, polygon_area, is_inside_polygon 10 10 11 11 if sys.platform == 'win32': -
anuga_work/production/hobart_2006/project.py
r4853 r4856 7 7 import sys 8 8 from time import localtime, strftime, gmtime 9 from anuga.utilities.polygon import read_polygon, plot_polygons _points, polygon_area9 from anuga.utilities.polygon import read_polygon, plot_polygons, polygon_area 10 10 from anuga.coordinate_transforms.redfearn import degminsec2decimal_degrees, convert_from_latlon_to_utm 11 11 … … 149 149 # Mark run - morning Fri 29 Sep; surrounds -100m and 20mish elevation 150 150 polyAll = read_polygon(polygondir+'new_extent_2.csv') 151 #plot_polygons _points([polyAll, polyAll2],figname='boundingpoly',verbose=False)151 #plot_polygons([polyAll, polyAll2],figname='boundingpoly',verbose=False) 152 152 #print 'Area of bounding polygon', polygon_area(polyAll) 153 153 polyAll_refine = read_polygon(polygondir+'extent_refine.csv') -
anuga_work/production/newcastle_2006/project_slide.py
r4853 r4856 7 7 import sys 8 8 from time import localtime, strftime, gmtime 9 from anuga.utilities.polygon import read_polygon, plot_polygons _points, polygon_area, is_inside_polygon9 from anuga.utilities.polygon import read_polygon, plot_polygons, polygon_area, is_inside_polygon 10 10 11 11 if sys.platform == 'win32': … … 98 98 poly_entrance_c0 = read_polygon(polygondir+'entrance.csv') 99 99 100 #plot_polygons _points([polyAll,poly_local,poly_newcastle],figname='fig',verbose=True)100 #plot_polygons([polyAll,poly_local,poly_newcastle],figname='fig',verbose=True) 101 101 ################################################################### 102 102 # Clipping regions for export to asc and regions for clipping data -
anuga_work/production/onslow_2006/export_results.py
r4581 r4856 1 import project , os1 import project_grad, os 2 2 import sys 3 3 … … 11 11 #time_dir = '' # HAT 1 in Pt Hedland 10000 yr 12 12 #time_dir = '' # MSL 1 in Pt Hedland 10000 yr 13 #time_dir = '20070 530_082920_run_final_1.5_onslow_nbartzis' # HAT 1 in Onslow 10000 yr13 #time_dir = '20070608_060316_run_final_1.5_onslow_nbartzis' # HAT 1 in Onslow 10000 yr 14 14 #time_dir = '20070531_002753_run_final_0.0_onslow_nbartzis' # MSL 1 in Onslow 10000 yr 15 time_dir = '20070608_062811_run_final_1.5_exmouth_nbartzis' # HAT 1 in Exmouth 10000 yr15 #time_dir = '20070608_062811_run_final_1.5_exmouth_nbartzis' # HAT 1 in Exmouth 10000 yr 16 16 #time_dir = '20070605_234905_run_final_0.0_exmouth_nbartzis' # MSL 1 in Exmouth 10000 yr 17 17 #time_dir = '20060704_063005' # 2006 simulation HAT Mw9 18 #time_dir = '20070619_042140_run_final_1.5_exmouth_nbartzis' # exmouth HAT longer time 19 #time_dir = '20070904_235118_run_final_1.5_jsexton_exmouth_original' # revised model original data 20 #time_dir = '20070905_052311_run_final_1.5_jsexton_exmouth_revised' # revised model original + survey data 21 time_dir = '20070906_061418_run_final_1.5_jsexton_exmouth_original' # refined again 22 #time_dir = '20070907_055936_run_final_1.5_jsexton_exmouth_revised' # refined again 18 23 cellsize = 25 19 24 #cellsize = 150 20 25 timestep = None 21 directory = project.outputdir 22 #name = directory + time_dir + sep + project.scenario_name 23 name = directory+time_dir+sep+project.scenario_name 26 directory = project_grad.output_dir 27 #name = directory + time_dir + sep + project_grad.scenario_name 28 name = directory+time_dir+sep+project_grad.scenario_name 29 #name = directory+time_dir+sep+'onslow_time_35760_0' 24 30 #name = directory+time_dir+sep+'source' 25 31 26 32 from anuga.shallow_water.data_manager import convert_dem_from_ascii2netcdf, dem2pts 27 33 28 elevgrid = directory+time_dir+sep+'onslow_elevation_town'34 #elevgrid = directory+time_dir+sep+'onslow_elevation_town' 29 35 30 convert_dem_from_ascii2netcdf(elevgrid, use_cache=True, verbose=True)36 #convert_dem_from_ascii2netcdf(elevgrid, use_cache=True, verbose=True) 31 37 32 dem2pts(elevgrid,run_up_elev=True,run_up_loc=True)38 #dem2pts(elevgrid,run_up_elev=True,run_up_loc=True) 33 39 34 """35 40 is_parallel = False 36 41 #is_parallel = True … … 40 45 #var = [2,3,4] # depth and speed 41 46 #var = [2] # depth 42 var = [ 0]47 var = [4] 43 48 44 49 for which_var in var: … … 48 53 49 54 if which_var == 1: # Absolute Momentum 50 outname = name + '_momentum _i1'55 outname = name + '_momentum' 51 56 quantityname = '(xmomentum**2 + ymomentum**2)**0.5' 52 57 … … 73 78 timestep = timestep, 74 79 cellsize = cellsize, 75 easting_min = project .e_min_area,76 easting_max = project .e_max_area,77 northing_min = project .n_min_area,78 northing_max = project .n_max_area,80 easting_min = project_grad.e_min_area, 81 easting_max = project_grad.e_max_area, 82 northing_min = project_grad.n_min_area, 83 northing_max = project_grad.n_max_area, 79 84 reduction = max, 80 85 verbose = True, … … 86 91 timestep = timestep, 87 92 cellsize = cellsize, 88 easting_min = project .e_min_area,89 easting_max = project .e_max_area,90 northing_min = project .n_min_area,91 northing_max = project .n_max_area,93 easting_min = project_grad.e_min_area, 94 easting_max = project_grad.e_max_area, 95 northing_min = project_grad.n_min_area, 96 northing_max = project_grad.n_max_area, 92 97 reduction = max, 93 98 verbose = True, 94 99 format = 'asc') 95 100 96 """ -
anuga_work/production/onslow_2006/project.py
r4581 r4856 142 142 n_max_area = 7620000 143 143 144 # old onslow for Gordon Monday 27 July 2007 145 e_min_area = 290828 146 e_max_area = 301021 147 n_min_area = 7596160 148 n_max_area = 7603100 149 144 150 # for test grid for ITEM ppt 145 151 #e_min_area = 207600 -
anuga_work/production/onslow_2006/project_grad.py
r4853 r4856 6 6 from os import sep, environ, getenv, getcwd, umask 7 7 from os.path import expanduser, basename, join 8 from anuga.utilities.polygon import read_polygon, plot_polygons _points, polygon_area, is_inside_polygon, number_mesh_triangles8 from anuga.utilities.polygon import read_polygon, plot_polygons, polygon_area, is_inside_polygon, number_mesh_triangles 9 9 import sys 10 10 from anuga.coordinate_transforms.redfearn import degminsec2decimal_degrees … … 41 41 #setup='final' 42 42 setup = 'trial' 43 which_data = 'original'44 #which_data = 'revised'43 #which_data = 'original' 44 which_data = 'revised' 45 45 boundary_event = 'onslow' 46 46 … … 99 99 100 100 # final topo name 101 combined_name ='onslow_combined_elevation_grad '101 combined_name ='onslow_combined_elevation_grad_%s' %which_data 102 102 combined_small_name = 'onslow_combined_elevation_grad' 103 103 … … 143 143 n_min_area = 7600000 144 144 n_max_area = 7610000 145 # for old onslow town 146 ##e_min_area = 300000 # Eastings min 147 ##e_max_area = 310000 148 ##n_min_area = 7600000 149 ##n_max_area = 7610000 145 150 146 151 #FIXME (Ole): It is bad to read in polygons in project.py. If a file does not exist … … 155 160 'inner_region': 5000, 156 161 'region_50m': 50000, 157 'onslow_v2': 500} 162 'onslow_v2': 500, 163 'beach_left': 200, 164 'beach_right': 200} 158 165 159 166 # Domain -
anuga_work/production/onslow_2006/project_urs.py
r4853 r4856 6 6 import sys 7 7 from time import localtime, strftime, gmtime 8 from anuga.utilities.polygon import read_polygon, plot_polygons _points, polygon_area, is_inside_polygon, number_mesh_triangles8 from anuga.utilities.polygon import read_polygon, plot_polygons, polygon_area, is_inside_polygon, number_mesh_triangles 9 9 from anuga.coordinate_transforms.redfearn import degminsec2decimal_degrees, convert_from_latlon_to_utm 10 10 from anuga.utilities.system_tools import get_user_name, get_host_name … … 39 39 finaltime=25000 40 40 setup='final' 41 source='dampier' 41 #source='exmouth' 42 source='onslow' 43 #source='dampier' 42 44 43 45 if setup =='trial': -
anuga_work/production/perth_2006/project.py
r4853 r4856 7 7 import sys 8 8 from time import localtime, strftime, gmtime 9 from anuga.utilities.polygon import read_polygon, plot_polygons _points, is_inside_polygon, number_mesh_triangles9 from anuga.utilities.polygon import read_polygon, plot_polygons, is_inside_polygon, number_mesh_triangles 10 10 #from anuga.coordinate_transforms.redfearn import degminsec2decimal_degrees, convert_points_from_latlon_to_utm 11 11 from anuga.utilities.system_tools import get_user_name … … 118 118 # Domain definitions 119 119 ############################### 120 from anuga.utilities.polygon import read_polygon, plot_polygons _points, polygon_area, is_inside_polygon120 from anuga.utilities.polygon import read_polygon, plot_polygons, polygon_area, is_inside_polygon 121 121 122 122 poly_all = read_polygon(polygons_dir+'bounding_poly.csv') -
anuga_work/production/pt_hedland_2006/project_urs.py
r4853 r4856 5 5 from os import sep, environ, getenv, getcwd,umask 6 6 from os.path import expanduser, basename 7 from anuga.utilities.polygon import read_polygon, plot_polygons _points, polygon_area, is_inside_polygon, number_mesh_triangles7 from anuga.utilities.polygon import read_polygon, plot_polygons, polygon_area, is_inside_polygon, number_mesh_triangles 8 8 import sys 9 9 from anuga.coordinate_transforms.redfearn import degminsec2decimal_degrees -
anuga_work/production/shark_bay_2007/export_results.py
r4624 r4856 6 6 7 7 #time_dir = '20070711_034317_run_final_0.0_shark_bay_onielsen_10000_tide0.0_coastpolys' # July 2006 event; small domain 8 time_dir = '20070710_023824_run_final_-0.85_shark_bay_onielsen_july2006_tide-0.85_coastpolys' # July 2006; larger domain 9 cellsize = 20 8 #time_dir = '20070710_023824_run_final_-0.85_shark_bay_onielsen_july2006_tide-0.85_coastpolys' # July 2006; larger domain 9 time_dir = '20070717_021240_run_final_0.79_shark_bay_onielsen_experimental_A0.5_T1200' 10 cellsize = 10 10 11 timestep = None 11 12 directory = project.output_dir … … 17 18 print 'output dir:', name 18 19 19 var = [0,2,3] # stage, depth and speed20 #var = [4]20 #var = [0,2,3] # stage, depth and speed 21 var = [4] 21 22 22 23 for which_var in var: -
anuga_work/production/shark_bay_2007/get_timeseries.py
r4624 r4856 13 13 #'20070710_001018_run_final_0.85_shark_bay_onielsen_10000_tide0.85_coastpolys': 'HAT10000'} 14 14 #'20070710_001110_run_final_0.85_shark_bay_onielsen_july2006_tide0.85_coastpolys': 'HAT06', 15 # '20070710_023824_run_final_-0.85_shark_bay_onielsen_july2006_tide-0.85_coastpolys': 'LAT06',16 '20070712_001432_run_final_0.0_shark_bay_onielsen_experimental_A0.5_T100': 'T100'}15 # 20070710_023824_run_final_-0.85_shark_bay_onielsen_july2006_tide-0.85_coastpolys': 'LAT06', 16 #'20070712_001432_run_final_0.0_shark_bay_onielsen_experimental_A0.5_T100': 'T100'} 17 17 #'20070712_002232_run_final_0.0_shark_bay_onielsen_experimental_A0.5_T200': 'T200'} 18 18 #'20070712_002254_run_final_0.0_shark_bay_onielsen_experimental_A0.5_T400': 'T400'} … … 34 34 #'20070712_003626_run_final_0.0_shark_bay_onielsen_experimental_A0.5_T12': 'T12'} 35 35 #'20070712_003638_run_final_0.0_shark_bay_onielsen_experimental_A0.5_T6': 'T6'} 36 #'20070717_021240_run_final_0.79_shark_bay_onielsen_experimental_A0.5_T1200': 'T1200new'} 37 #'20070717_021545_run_final_0.79_shark_bay_onielsen_experimental_A0.5_T1300': 'T1300new'} 38 #'20070717_021601_run_final_0.79_shark_bay_onielsen_experimental_A0.5_T1400': 'T1400new'} 39 #'20070717_021615_run_final_0.79_shark_bay_onielsen_experimental_A0.5_T1500': 'T1500new'} 40 #'20070717_021646_run_final_0.79_shark_bay_onielsen_experimental_A0.5_T1100': 'T1100new'} 41 #'20070717_021656_run_final_0.79_shark_bay_onielsen_experimental_A0.5_T1000': 'T1000new'} 42 #'20070717_021707_run_final_0.79_shark_bay_onielsen_experimental_A0.5_T900': 'T900new'} 43 #'20070717_021717_run_final_0.79_shark_bay_onielsen_experimental_A0.5_T800': 'T800new'} 44 #'20070717_021750_run_final_0.79_shark_bay_onielsen_experimental_A0.5_T700': 'T700new'} 45 #'20070717_021801_run_final_0.79_shark_bay_onielsen_experimental_A0.5_T600': 'T600new'} 46 #'20070717_021814_run_final_0.79_shark_bay_onielsen_experimental_A0.5_T300': 'T300new'} 47 #'20070717_021831_run_final_0.79_shark_bay_onielsen_experimental_A0.5_T3000': 'T3000new'} 48 #'20070716_025846_run_final_0.79_shark_bay_onielsen_july2006_X1_fixedbathy': '2006new'} 49 #'20070713_033049_run_final_0.0_shark_bay_onielsen_july2006_X10_fixedbathy': '2006new10'} 50 #'20070716_011226_run_final_0.0_shark_bay_onielsen_july2006_X50_fixedbathy': '2006new50'} 51 '20070716_024711_run_final_0.79_shark_bay_onielsen_10000_X1_fixedbathy': '10000new'} 36 52 #'': ''} 37 53 … … 49 65 report = False, 50 66 plot_quantity = ['stage', 'speed'], 51 generate_fig = True,67 generate_fig = False, 52 68 time_min = None, 53 69 time_max = None, -
anuga_work/production/shark_bay_2007/project.py
r4853 r4856 7 7 from os import sep, environ, getenv, getcwd, umask 8 8 from os.path import expanduser, basename, join 9 from anuga.utilities.polygon import read_polygon, plot_polygons _points, polygon_area, is_inside_polygon, number_mesh_triangles9 from anuga.utilities.polygon import read_polygon, plot_polygons, polygon_area, is_inside_polygon, number_mesh_triangles 10 10 import sys 11 11 from anuga.coordinate_transforms.redfearn import degminsec2decimal_degrees … … 145 145 146 146 #gauges 147 gauge_name = 'gaugesv 2.csv'147 gauge_name = 'gaugesv3.csv' 148 148 #gauge_name_test = 'gauge_checking_test.csv' 149 149 gauges_dir = anuga_dir+sep+'gauges'+sep -
anuga_work/production/sydney_2006/project_slide.py
r4853 r4856 7 7 import sys 8 8 from time import localtime, strftime, gmtime 9 from anuga.utilities.polygon import read_polygon, plot_polygons _points, polygon_area, is_inside_polygon9 from anuga.utilities.polygon import read_polygon, plot_polygons, polygon_area, is_inside_polygon 10 10 11 11 if sys.platform == 'win32': -
anuga_work/production/wollongong_2006/project_kembla.py
r4853 r4856 7 7 import sys 8 8 from time import localtime, strftime, gmtime 9 from anuga.utilities.polygon import read_polygon, plot_polygons _points, polygon_area, is_inside_polygon9 from anuga.utilities.polygon import read_polygon, plot_polygons, polygon_area, is_inside_polygon 10 10 11 11 #Making assumptions about the location of scenario data -
anuga_work/production/wollongong_2006/project_slide.py
r4853 r4856 7 7 import sys 8 8 from time import localtime, strftime, gmtime 9 from anuga.utilities.polygon import read_polygon, plot_polygons _points, polygon_area, is_inside_polygon9 from anuga.utilities.polygon import read_polygon, plot_polygons, polygon_area, is_inside_polygon 10 10 11 11 if sys.platform == 'win32': -
anuga_work/production/wollongong_2006/run_kembla.py
r4696 r4856 41 41 meshname = project_kembla.meshname+'.msh' 42 42 43 # creates DEM from asc data44 convert_dem_from_ascii2netcdf(dem_name, use_cache=True, verbose=True)45 46 # creates pts file for onshore DEM47 dem2pts(dem_name, use_cache=True, verbose=True)48 49 print 'create onshore'50 G = Geospatial_data(file_name = project_kembla.dem_name + '.pts')51 52 print 'export points'53 G.export_points_file(project_kembla.combined_dem_name + '.pts')43 ### creates DEM from asc data 44 ##convert_dem_from_ascii2netcdf(dem_name, use_cache=True, verbose=True) 45 ## 46 ###creates pts file for onshore DEM 47 ##dem2pts(dem_name, use_cache=True, verbose=True) 48 ## 49 ##print 'create onshore' 50 ##G = Geospatial_data(file_name = project_kembla.dem_name + '.pts') 51 ## 52 ##print 'export points' 53 ##G.export_points_file(project_kembla.combined_dem_name + '.pts') 54 54 #G.export_points_file(project_kembla.combined_dem_name + '.xya') 55 55 … … 114 114 f=lambda t: [(60<t<660)*6, 0, 0]) 115 115 116 domain.set_boundary( {'top': Bd, 116 domain.set_boundary( {'top': Bd,'right': Bw, 'bottom': Bd, 'left': Bd} ) 117 117 118 118 … … 125 125 from anuga.abstract_2d_finite_volumes.quantity import Quantity 126 126 127 for t in domain.evolve(yieldstep = 30, finaltime = 5000):127 for t in domain.evolve(yieldstep = 10, finaltime = 1000): 128 128 domain.write_time() 129 domain.write_boundary_statistics(tags = ' ')129 domain.write_boundary_statistics(tags = 'right') 130 130 stagestep = domain.get_quantity('stage') 131 131
Note: See TracChangeset
for help on using the changeset viewer.