Changeset 4856


Ignore:
Timestamp:
Nov 27, 2007, 4:18:39 PM (15 years ago)
Author:
sexton
Message:

rename plot_polygons_points to original name of plot_polygons

Files:
1 added
33 edited

Legend:

Unmodified
Added
Removed
  • anuga_core/documentation/user_manual/demos/cairns/GetTimeseries.py

    r4087 r4856  
    2121
    2222# nominate directory location of sww file with associated attribute
    23 production_dirs = {'slide': 'Slide'}
     23production_dirs = {#'slide': 'Slide',
     24                   'fixed_wave': 'Fixed Wave'}
    2425
    2526# Generate figures
     
    3637                                      reportname = '',
    3738                                      plot_quantity = ['stage', 'speed'],
     39                                      generate_fig = False,
    3840                                      surface = False,
    3941                                      time_min = None,
  • anuga_core/documentation/user_manual/demos/cairns/project.py

    r4854 r4856  
    77import sys
    88from time import localtime, strftime, gmtime
    9 from anuga.utilities.polygon import read_polygon, plot_polygons_points, \
     9from anuga.utilities.polygon import read_polygon, plot_polygons, \
    1010                                    polygon_area, is_inside_polygon
    1111
  • anuga_core/documentation/user_manual/demos/cairns/runcairns.py

    r4853 r4856  
    3939# Define scenario as either slide or fixed_wave.
    4040#------------------------------------------------------------------------------
    41 scenario = 'slide'
    42 #scenario = 'fixed_wave'
     41#scenario = 'slide'
     42scenario = 'fixed_wave'
    4343
    4444if os.access(scenario, os.F_OK) == 0:
  • anuga_core/source/anuga/abstract_2d_finite_volumes/test_neighbour_mesh.py

    r4851 r4856  
    973973
    974974
    975         from anuga.utilities.polygon import plot_polygons_points
     975        from anuga.utilities.polygon import plot_polygons
    976976
    977977        # First do the continuous version of mesh
  • anuga_core/source/anuga/fit_interpolate/interpolate.py

    r4853 r4856  
    558558
    559559                if verbose is True:
    560                     from anuga.utilities.polygon import plot_polygons_points
     560                    from anuga.utilities.polygon import plot_polygons
    561561                    #out_interp_pts = take(interpolation_points,[indices])
    562562                    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],
    564564                                         ['line','point','outside'],figname='points_boundary_out',label=title,verbose=verbose)
    565565
     
    574574            # Plot boundary and interpolation points
    575575            if verbose is True:
    576                 from anuga.utilities.polygon import plot_polygons_points
     576                from anuga.utilities.polygon import plot_polygons
    577577                title = 'Interpolation function: Polygon and interpolation points'
    578                 plot_polygons_points([mesh_boundary_polygon,interpolation_points],
     578                plot_polygons([mesh_boundary_polygon,interpolation_points],
    579579                                     ['line','point'],figname='points_boundary',label=title,verbose=verbose)
    580580
  • anuga_core/source/anuga/utilities/polygon.py

    r4852 r4856  
    341341    return abs(poly_area/2)
    342342
    343 def plot_polygons_points(polygons_points, style=None,
     343def plot_polygons(polygons_points, style=None,
    344344                         figname=None, label=None, verbose=False):
    345345   
  • anuga_work/production/broome_2006/get_timeseries.py

    r4581 r4856  
    1818
    1919# bulli event
    20 production_dirs = {'20070702_050939_run_final_4.9_exmouth_nbartzis': ''}
     20production_dirs = {'20070703_061553_run_final_0_dampier_nbartzis': ''}
    2121
    22 # 20070702_050939_run_final_4.9_exmouth_nbartzis HAT exmouth
    23 # 20070615_062730_run_final_0_exmouth_nbartzis # MSL exmouth
    24 # 20070618_064014_run_final_4.9_onslow_nbartzis # HAT onslow
    25 # 20070622_004008_run_final_0_onslow_nbartzis # MSL onslow
    26 # 20070615_063137_run_final_4.9_dampier_nbartzis # HAT dampier
    27 # 20070621_002922_run_final_0_dampier_nbartzis # MSL dampier
     22# 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
    2828
    2929is_parallel = True
  • anuga_work/production/broome_2006/project.py

    r4853 r4856  
    77import sys
    88from time import localtime, strftime, gmtime
    9 from anuga.utilities.polygon import read_polygon, plot_polygons_points, is_inside_polygon, number_mesh_triangles
     9from anuga.utilities.polygon import read_polygon, plot_polygons, is_inside_polygon, number_mesh_triangles
    1010#from anuga.coordinate_transforms.redfearn import degminsec2decimal_degrees, convert_points_from_latlon_to_utm
    1111from anuga.utilities.system_tools import get_user_name
     
    126126# Domain definitions
    127127###############################
    128 from anuga.utilities.polygon import read_polygon, plot_polygons_points, polygon_area, is_inside_polygon
     128from anuga.utilities.polygon import read_polygon, plot_polygons, polygon_area, is_inside_polygon
    129129
    130130poly_all = read_polygon(polygons_dir+'extent_small.csv')
  • anuga_work/production/broome_2006/project_urs.py

    r4853 r4856  
    66from os import sep, environ, getenv, getcwd,umask
    77from os.path import expanduser, basename
    8 from anuga.utilities.polygon import read_polygon, plot_polygons_points, polygon_area, is_inside_polygon, number_mesh_triangles
     8from anuga.utilities.polygon import read_polygon, plot_polygons, polygon_area, is_inside_polygon, number_mesh_triangles
    99import sys
    1010from anuga.coordinate_transforms.redfearn import degminsec2decimal_degrees
     
    152152# Domain definitions
    153153###############################
    154 from anuga.utilities.polygon import read_polygon, plot_polygons_points, polygon_area, is_inside_polygon
     154from anuga.utilities.polygon import read_polygon, plot_polygons, polygon_area, is_inside_polygon
    155155
    156156poly_all = read_polygon(polygons_dir+'extent_small.csv')
     
    160160#p2 = read_polygon(polygons_dir+'broome_north_coast_inside_extent.csv')
    161161#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')
    164164#poly_all = read_polygon(polygons_dir+'extent.csv')
    165165res_poly_all = 150000*res_factor
  • anuga_work/production/busselton_2006/project.py

    r4853 r4856  
    77import sys
    88from time import localtime, strftime, gmtime
    9 from anuga.utilities.polygon import read_polygon, plot_polygons_points, is_inside_polygon, number_mesh_triangles
     9from anuga.utilities.polygon import read_polygon, plot_polygons, is_inside_polygon, number_mesh_triangles
    1010#from anuga.coordinate_transforms.redfearn import degminsec2decimal_degrees, convert_points_from_latlon_to_utm
    1111from anuga.utilities.system_tools import get_user_name
     
    129129res_busselton2 = 500
    130130
    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)
    132132
    133133interior_regions = [[poly_busselton1,res_busselton1],[poly_busselton2,res_busselton2]]
  • anuga_work/production/dampier_2006/project.py

    r4853 r4856  
    66import sys
    77from time import localtime, strftime, gmtime
    8 from anuga.utilities.polygon import read_polygon, plot_polygons_points, polygon_area, is_inside_polygon, number_mesh_triangles
     8from anuga.utilities.polygon import read_polygon, plot_polygons, polygon_area, is_inside_polygon, number_mesh_triangles
    99from anuga.coordinate_transforms.redfearn import degminsec2decimal_degrees, convert_from_latlon_to_utm
    1010from anuga.utilities.system_tools import get_user_name
     
    237237
    238238
    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,
    240240#                poly_coast,poly_NWislands,poly_island0,poly_island1,poly_island2,
    241241#                poly_island3,poly_island4,poly_island5,poly_island6,
  • anuga_work/production/dampier_2006/project_urs.py

    r4853 r4856  
    44from os import sep, environ, getenv, getcwd,umask
    55from os.path import expanduser, basename
    6 from anuga.utilities.polygon import read_polygon, plot_polygons_points, polygon_area, is_inside_polygon, number_mesh_triangles
     6from anuga.utilities.polygon import read_polygon, plot_polygons, polygon_area, is_inside_polygon, number_mesh_triangles
    77import sys
    88from anuga.coordinate_transforms.redfearn import degminsec2decimal_degrees
     
    243243
    244244
    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,
    246246#                poly_coast,poly_NWislands,poly_island0,poly_island1,poly_island2,
    247247#                poly_island3,poly_island4,poly_island5,poly_island6,
  • anuga_work/production/dampier_2006/report/anuga_setup.tex

    r4479 r4856  
    3838
    3939
    40 \begin{figure}[hbt]
     40\begin{figure}[h]
    4141
    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 
    4445  \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.
     46Region 1: Surrounds onshore gas infrastructure at Dampier with a cell area of 500 m$^2$ (lateral accuracy 30 m).
    4647Region 2: Surrounds the gas pipelines with a cell area of 2000 m$^2$ (lateral accuracy 100 m).
    4748Region 3: Surrounds the coastal regions with a cell area of 10000 m$^2$ (lateral accuracy  m).
     
    7071results are therefore likely to be over estimates.
    7172
    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  
    11\begin{sidewaysfigure}
    22%\begin{figure}
    3 \centerline{ \includegraphics[scale=0.4]{../report_figures/max_inundation_high_tide_v2.jpg}}
     3\centerline{ \includegraphics[scale=0.4]{../report_figures/max_inundation_high_tide_v3.jpg}}
    44\caption{Maximum inundation map for the HAT scenario for Dampier region. Data: WA DLI, DPI and AHO.} 
    55\label{fig:HAT_max_inundation}
  • anuga_work/production/exmouth_2006/export_results.py

    r4479 r4856  
    55from os import sep
    66
    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
     12time_dir = '20070619_042417_run_final_0_exmouth_nbartzis' # MSL 1 in Exmouth 10000 yr
    1713cellsize = 25
    1814timestep = None
     
    3733
    3834    if which_var == 2:  # Depth
    39         outname = name + '_depth'
     35        outname = name + '_depth_north'
    4036        quantityname = 'stage-elevation' 
    4137
    4238    if which_var == 3:  # Speed
    43         outname = name + '_speed'
     39        outname = name + '_speed_north'
    4440        quantityname = '(xmomentum**2 + ymomentum**2)**0.5/(stage-elevation+1.e-6/(stage-elevation))'  #Speed
    4541
  • anuga_work/production/exmouth_2006/project.py

    r4853 r4856  
    55from os import sep, environ, getenv, getcwd,umask
    66from os.path import expanduser, basename
    7 from anuga.utilities.polygon import read_polygon, plot_polygons_points, polygon_area, is_inside_polygon, number_mesh_triangles
     7from anuga.utilities.polygon import read_polygon, plot_polygons, polygon_area, is_inside_polygon, number_mesh_triangles
    88import sys
    99from anuga.coordinate_transforms.redfearn import degminsec2decimal_degrees
     
    154154# Domain definitions
    155155###############################
    156 from anuga.utilities.polygon import read_polygon, plot_polygons_points, polygon_area, is_inside_polygon
     156from anuga.utilities.polygon import read_polygon, plot_polygons, polygon_area, is_inside_polygon
    157157
    158158poly_all = read_polygon(polygons_dir+'extent_more_points.csv')
     
    222222###################################################################
    223223
    224 # exporting asc grid
     224# exporting asc grid - Exmouth town area
    225225eastingmin = 200663
    226226eastingmax = 207500
    227227northingmax = 7574000
    228228northingmin = 7564500
     229
     230# exporting asc grid - carparks etc on NW peninsula
     231e_min_area = 192700
     232e_max_area = 203755
     233n_max_area = 7587330
     234n_min_area = 7579460
  • anuga_work/production/galle_2007/project.py

    r4853 r4856  
    88import sys
    99from time import localtime, strftime, gmtime
    10 from anuga.utilities.polygon import read_polygon, plot_polygons_points, is_inside_polygon, \
     10from anuga.utilities.polygon import read_polygon, plot_polygons, is_inside_polygon, \
    1111                                    number_mesh_triangles, polygon_area
    1212from anuga.utilities.system_tools import get_user_name
     
    133133                    [poly_canal, canal_res]]
    134134
    135 #plot_polygons_points([bounding_polygon,poly_intermediate,\
     135#plot_polygons([bounding_polygon,poly_intermediate,\
    136136#               poly_shallow,poly_coast,poly_canal],\
    137137#              figname='boundingpoly',verbose=False)
  • anuga_work/production/gold_coast_2007/project.py

    r4853 r4856  
    77import sys
    88from time import localtime, strftime, gmtime
    9 from anuga.utilities.polygon import read_polygon, plot_polygons_points, polygon_area, is_inside_polygon
     9from anuga.utilities.polygon import read_polygon, plot_polygons, polygon_area, is_inside_polygon
    1010
    1111if sys.platform == 'win32':
  • anuga_work/production/hobart_2006/project.py

    r4853 r4856  
    77import sys
    88from time import localtime, strftime, gmtime
    9 from anuga.utilities.polygon import read_polygon, plot_polygons_points, polygon_area
     9from anuga.utilities.polygon import read_polygon, plot_polygons, polygon_area
    1010from anuga.coordinate_transforms.redfearn import degminsec2decimal_degrees, convert_from_latlon_to_utm
    1111
     
    149149# Mark run - morning Fri 29 Sep; surrounds -100m and 20mish elevation
    150150polyAll = 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)
    152152#print 'Area of bounding polygon', polygon_area(polyAll)
    153153polyAll_refine = read_polygon(polygondir+'extent_refine.csv')
  • anuga_work/production/newcastle_2006/project_slide.py

    r4853 r4856  
    77import sys
    88from time import localtime, strftime, gmtime
    9 from anuga.utilities.polygon import read_polygon, plot_polygons_points, polygon_area, is_inside_polygon
     9from anuga.utilities.polygon import read_polygon, plot_polygons, polygon_area, is_inside_polygon
    1010
    1111if sys.platform == 'win32':
     
    9898poly_entrance_c0 = read_polygon(polygondir+'entrance.csv')
    9999
    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)
    101101###################################################################
    102102# Clipping regions for export to asc and regions for clipping data
  • anuga_work/production/onslow_2006/export_results.py

    r4581 r4856  
    1 import project, os
     1import project_grad, os
    22import sys
    33
     
    1111#time_dir = '' # HAT 1 in Pt Hedland 10000 yr
    1212#time_dir = '' # MSL 1 in Pt Hedland 10000 yr
    13 #time_dir = '20070530_082920_run_final_1.5_onslow_nbartzis' # HAT 1 in Onslow 10000 yr
     13#time_dir = '20070608_060316_run_final_1.5_onslow_nbartzis' # HAT 1 in Onslow 10000 yr
    1414#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 yr
     15#time_dir = '20070608_062811_run_final_1.5_exmouth_nbartzis' # HAT 1 in Exmouth 10000 yr
    1616#time_dir = '20070605_234905_run_final_0.0_exmouth_nbartzis' # MSL 1 in Exmouth 10000 yr
    1717#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
     21time_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
    1823cellsize = 25
    1924#cellsize = 150
    2025timestep = None
    21 directory = project.outputdir
    22 #name = directory + time_dir + sep + project.scenario_name
    23 name = directory+time_dir+sep+project.scenario_name
     26directory = project_grad.output_dir
     27#name = directory + time_dir + sep + project_grad.scenario_name
     28name = directory+time_dir+sep+project_grad.scenario_name
     29#name = directory+time_dir+sep+'onslow_time_35760_0'
    2430#name = directory+time_dir+sep+'source'
    2531
    2632from anuga.shallow_water.data_manager import convert_dem_from_ascii2netcdf, dem2pts
    2733
    28 elevgrid = directory+time_dir+sep+'onslow_elevation_town'
     34#elevgrid = directory+time_dir+sep+'onslow_elevation_town'
    2935
    30 convert_dem_from_ascii2netcdf(elevgrid, use_cache=True, verbose=True)
     36#convert_dem_from_ascii2netcdf(elevgrid, use_cache=True, verbose=True)
    3137
    32 dem2pts(elevgrid,run_up_elev=True,run_up_loc=True)
     38#dem2pts(elevgrid,run_up_elev=True,run_up_loc=True)
    3339
    34 """
    3540is_parallel = False
    3641#is_parallel = True
     
    4045#var = [2,3,4] # depth and speed
    4146#var = [2] # depth
    42 var = [0]
     47var = [4]
    4348
    4449for which_var in var:
     
    4853
    4954    if which_var == 1:  # Absolute Momentum
    50         outname = name + '_momentum_i1'
     55        outname = name + '_momentum'
    5156        quantityname = '(xmomentum**2 + ymomentum**2)**0.5' 
    5257
     
    7378                        timestep = timestep,
    7479                        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,       
    7984                        reduction = max,
    8085                        verbose = True,
     
    8691                    timestep = timestep,
    8792                    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,       
    9297                    reduction = max,
    9398                    verbose = True,
    9499                    format = 'asc')
    95100
    96 """
  • anuga_work/production/onslow_2006/project.py

    r4581 r4856  
    142142n_max_area = 7620000
    143143
     144# old onslow for Gordon Monday 27 July 2007
     145e_min_area = 290828
     146e_max_area = 301021
     147n_min_area = 7596160
     148n_max_area = 7603100
     149
    144150# for test grid for ITEM ppt
    145151#e_min_area = 207600
  • anuga_work/production/onslow_2006/project_grad.py

    r4853 r4856  
    66from os import sep, environ, getenv, getcwd, umask
    77from os.path import expanduser, basename, join
    8 from anuga.utilities.polygon import read_polygon, plot_polygons_points, polygon_area, is_inside_polygon, number_mesh_triangles
     8from anuga.utilities.polygon import read_polygon, plot_polygons, polygon_area, is_inside_polygon, number_mesh_triangles
    99import sys
    1010from anuga.coordinate_transforms.redfearn import degminsec2decimal_degrees
     
    4141#setup='final'
    4242setup = 'trial'
    43 which_data = 'original'
    44 #which_data = 'revised'
     43#which_data = 'original'
     44which_data = 'revised'
    4545boundary_event = 'onslow'
    4646
     
    9999
    100100# final topo name
    101 combined_name ='onslow_combined_elevation_grad'
     101combined_name ='onslow_combined_elevation_grad_%s' %which_data
    102102combined_small_name = 'onslow_combined_elevation_grad'
    103103
     
    143143n_min_area = 7600000
    144144n_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
    145150               
    146151#FIXME (Ole): It is bad to read in polygons in project.py. If a file does not exist
     
    155160    'inner_region': 5000,
    156161    'region_50m': 50000,
    157     'onslow_v2': 500}
     162    'onslow_v2': 500,
     163    'beach_left': 200,
     164    'beach_right': 200}
    158165
    159166# Domain
  • anuga_work/production/onslow_2006/project_urs.py

    r4853 r4856  
    66import sys
    77from time import localtime, strftime, gmtime
    8 from anuga.utilities.polygon import read_polygon, plot_polygons_points, polygon_area, is_inside_polygon, number_mesh_triangles
     8from anuga.utilities.polygon import read_polygon, plot_polygons, polygon_area, is_inside_polygon, number_mesh_triangles
    99from anuga.coordinate_transforms.redfearn import degminsec2decimal_degrees, convert_from_latlon_to_utm
    1010from anuga.utilities.system_tools import get_user_name, get_host_name
     
    3939finaltime=25000
    4040setup='final'
    41 source='dampier'
     41#source='exmouth'
     42source='onslow'
     43#source='dampier'
    4244
    4345if setup =='trial':
  • anuga_work/production/perth_2006/project.py

    r4853 r4856  
    77import sys
    88from time import localtime, strftime, gmtime
    9 from anuga.utilities.polygon import read_polygon, plot_polygons_points, is_inside_polygon, number_mesh_triangles
     9from anuga.utilities.polygon import read_polygon, plot_polygons, is_inside_polygon, number_mesh_triangles
    1010#from anuga.coordinate_transforms.redfearn import degminsec2decimal_degrees, convert_points_from_latlon_to_utm
    1111from anuga.utilities.system_tools import get_user_name
     
    118118# Domain definitions
    119119###############################
    120 from anuga.utilities.polygon import read_polygon, plot_polygons_points, polygon_area, is_inside_polygon
     120from anuga.utilities.polygon import read_polygon, plot_polygons, polygon_area, is_inside_polygon
    121121
    122122poly_all = read_polygon(polygons_dir+'bounding_poly.csv')
  • anuga_work/production/pt_hedland_2006/project_urs.py

    r4853 r4856  
    55from os import sep, environ, getenv, getcwd,umask
    66from os.path import expanduser, basename
    7 from anuga.utilities.polygon import read_polygon, plot_polygons_points, polygon_area, is_inside_polygon, number_mesh_triangles
     7from anuga.utilities.polygon import read_polygon, plot_polygons, polygon_area, is_inside_polygon, number_mesh_triangles
    88import sys
    99from anuga.coordinate_transforms.redfearn import degminsec2decimal_degrees
  • anuga_work/production/shark_bay_2007/export_results.py

    r4624 r4856  
    66
    77#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
     9time_dir = '20070717_021240_run_final_0.79_shark_bay_onielsen_experimental_A0.5_T1200'
     10cellsize = 10
    1011timestep = None
    1112directory = project.output_dir
     
    1718print 'output dir:', name
    1819
    19 var = [0,2,3] # stage, depth and speed
    20 #var = [4]
     20#var = [0,2,3] # stage, depth and speed
     21var = [4]
    2122
    2223for which_var in var:
  • anuga_work/production/shark_bay_2007/get_timeseries.py

    r4624 r4856  
    1313                   #'20070710_001018_run_final_0.85_shark_bay_onielsen_10000_tide0.85_coastpolys': 'HAT10000'}
    1414                   #'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'}
    1717                    #'20070712_002232_run_final_0.0_shark_bay_onielsen_experimental_A0.5_T200': 'T200'}
    1818                    #'20070712_002254_run_final_0.0_shark_bay_onielsen_experimental_A0.5_T400': 'T400'}
     
    3434                    #'20070712_003626_run_final_0.0_shark_bay_onielsen_experimental_A0.5_T12': 'T12'}
    3535                    #'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'}
    3652                    #'': ''}
    3753
     
    4965                             report = False,
    5066                             plot_quantity = ['stage', 'speed'],
    51                              generate_fig = True,
     67                             generate_fig = False,
    5268                             time_min = None,
    5369                             time_max = None,
  • anuga_work/production/shark_bay_2007/project.py

    r4853 r4856  
    77from os import sep, environ, getenv, getcwd, umask
    88from os.path import expanduser, basename, join
    9 from anuga.utilities.polygon import read_polygon, plot_polygons_points, polygon_area, is_inside_polygon, number_mesh_triangles
     9from anuga.utilities.polygon import read_polygon, plot_polygons, polygon_area, is_inside_polygon, number_mesh_triangles
    1010import sys
    1111from anuga.coordinate_transforms.redfearn import degminsec2decimal_degrees
     
    145145
    146146#gauges
    147 gauge_name = 'gaugesv2.csv'
     147gauge_name = 'gaugesv3.csv'
    148148#gauge_name_test = 'gauge_checking_test.csv'
    149149gauges_dir = anuga_dir+sep+'gauges'+sep
  • anuga_work/production/sydney_2006/project_slide.py

    r4853 r4856  
    77import sys
    88from time import localtime, strftime, gmtime
    9 from anuga.utilities.polygon import read_polygon, plot_polygons_points, polygon_area, is_inside_polygon
     9from anuga.utilities.polygon import read_polygon, plot_polygons, polygon_area, is_inside_polygon
    1010
    1111if sys.platform == 'win32':
  • anuga_work/production/wollongong_2006/project_kembla.py

    r4853 r4856  
    77import sys
    88from time import localtime, strftime, gmtime
    9 from anuga.utilities.polygon import read_polygon, plot_polygons_points, polygon_area, is_inside_polygon
     9from anuga.utilities.polygon import read_polygon, plot_polygons, polygon_area, is_inside_polygon
    1010
    1111#Making assumptions about the location of scenario data
  • anuga_work/production/wollongong_2006/project_slide.py

    r4853 r4856  
    77import sys
    88from time import localtime, strftime, gmtime
    9 from anuga.utilities.polygon import read_polygon, plot_polygons_points, polygon_area, is_inside_polygon
     9from anuga.utilities.polygon import read_polygon, plot_polygons, polygon_area, is_inside_polygon
    1010
    1111if sys.platform == 'win32':
  • anuga_work/production/wollongong_2006/run_kembla.py

    r4696 r4856  
    4141meshname = project_kembla.meshname+'.msh'
    4242
    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')
     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')
    5454#G.export_points_file(project_kembla.combined_dem_name + '.xya')
    5555
     
    114114                   f=lambda t: [(60<t<660)*6, 0, 0])
    115115
    116 domain.set_boundary( {'top': Bd,  'right': Bw, 'bottom': Bd, 'left': Bd} )
     116domain.set_boundary( {'top': Bd,'right': Bw, 'bottom': Bd, 'left': Bd} )
    117117
    118118
     
    125125from anuga.abstract_2d_finite_volumes.quantity import Quantity
    126126
    127 for t in domain.evolve(yieldstep = 30, finaltime = 5000):
     127for t in domain.evolve(yieldstep = 10, finaltime = 1000):
    128128    domain.write_time()
    129     domain.write_boundary_statistics(tags = '')
     129    domain.write_boundary_statistics(tags = 'right')
    130130    stagestep = domain.get_quantity('stage')
    131131   
Note: See TracChangeset for help on using the changeset viewer.