Changeset 6427


Ignore:
Timestamp:
Feb 27, 2009, 9:02:48 AM (15 years ago)
Author:
jgriffin
Message:

updated batemans bay scripts

Location:
anuga_work/production/new_south_wales/batemans_bay
Files:
4 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/production/new_south_wales/batemans_bay/build_elevation.py

    r6386 r6427  
    5959        dem2pts(absolute_filename, use_cache=True, verbose=True)
    6060
    61         geospatial_data[filename] = Geospatial_data(file_name=absolute_filename+'.pts',
     61        G_grid = Geospatial_data(file_name=absolute_filename+'.pts',
    6262                                                    verbose=True)
     63        print 'Clip geospatial object'
     64        geospatial_data[filename] = G_grid.clip(project.bounding_polygon)
    6365
    6466# Create Geospatial data from TXT files
     
    6668    for filename in project.point_filenames:
    6769        absolute_filename = join(project.topographies_folder, filename)
    68         geospatial_data[filename] = Geospatial_data(file_name=absolute_filename,
     70        G_points = Geospatial_data(file_name=absolute_filename,
    6971                                                    verbose=True)
     72        print 'Clip geospatial object'
     73        geospatial_data[filename] = G_points.clip(project.bounding_polygon)
    7074
    7175#-------------------------------------------------------------------------------
     
    8286    G += geospatial_data[key]
    8387       
    84 print 'Clip combined geospatial data'
     88#print 'Clip combined geospatial data'
    8589##G_clip = G.clip_outside(project.poly_aoi1)
    8690##G_all = G_clip + geospatial_data[project.point_filenames[4]]
    8791#G_clipped = G_all.clip(project.poly_all)
    88 G_clip = G.clip(project.bounding_polygon)
     92#G_clip = G.clip(project.bounding_polygon)
    8993
    9094
    9195print 'Export combined DEM file'
    92 G_clip.export_points_file(project.combined_elevation + '.pts')
     96G.export_points_file(project.combined_elevation + '.pts')
    9397print 'Do txt version too'
    9498# Use for comparision in ARC
    95 G_clip.export_points_file(project.combined_elevation + '.txt')
     99G.export_points_file(project.combined_elevation + '.txt')
    96100
  • anuga_work/production/new_south_wales/batemans_bay/export_results_max.py

    r6386 r6427  
    2929##time_dirs = [time_dir1, time_dir2]#, time_dir3]
    3030
    31 time_dir1 = '20090130_165532_run_final_0.6_27255_extend_dt_kvanputt'
     31time_dir1 = '20090225_153333_run_final_0_31821_jgriffin'
    3232##time_dir2 = '20081211_162311_run_final_0_27255_alpha0.1_kvanputt'
    3333##time_dir3 = '20081211_162346_run_final_0_68693_alpha0.1_kvanputt'
     
    3939 
    4040
    41 cellsize = 250
     41cellsize = 20
    4242##cellsize = 5
    4343
  • anuga_work/production/new_south_wales/batemans_bay/project.py

    r6386 r6427  
    1616# this section needs to be updated to reflect the modelled community.
    1717# Note, the user needs to set up the directory system accordingly
    18 state = 'western_australia'
    19 scenario_name = 'busselton'
    20 scenario_folder = 'busselton_tsunami_scenario'
     18state = 'new_south_wales'
     19scenario_name = 'batemans_bay'
     20scenario_folder = 'batemans_bay_tsunami_scenario_2009'
    2121
    2222#-------------------------------------------------------------------------------
     
    2727# One or all can be changed each time the run_model script is executed
    2828tide = 0                # difference between MSL and HAT
    29 event_number = 27197    # the event number or the mux file name
     29event_number = 31821    # the event number or the mux file name
    3030alpha = 0.1             # smoothing parameter for mesh
    3131friction=0.01           # manning's friction coefficient
    3232starttime=0             # start time for simulation
    33 finaltime=80000         # final time for simulation
     33finaltime=1000         # final time for simulation
    3434
    3535setup = 'final'         # This can be one of three values
     
    5858# Used in build_elevation.py
    5959# Format for ascii grids, as produced in ArcGIS + a projection file
    60 ascii_grid_filenames = ['busselton_v2',   # Topo
    61                         'grid_250m_clip'] # Busselton Topo
    62 
     60ascii_grid_filenames = ['1a',# Topographic data
     61                        '1b',
     62                        '2b',
     63                        '2a_3',
     64                        '3b',
     65                        '3a',
     66                        '4a_2',
     67                        '4b',]
     68                   
    6369# Format for point is x,y,elevation (with header)
    64 point_filenames = ['Busselton_Contour0.txt',    # Coastline
    65                    'Busselton_BeachSurvey.txt', # Beach survey
    66                    'Busselton_NavyFinal.txt',   # Bathymetry
    67                    'Busselton_Chart.txt',       # Bathymetry Charts
    68                    'Busselton_Digitised.txt',   # Digitised Fairsheet
    69                    'Busselton_250m.txt',        # 250m
    70                    'Bunbury_TIN.txt',           # Bunbury aoi TIN'd in ArcGIS
    71                    'Busselton_TIN.txt',         # Busselton aoi TIN'd in ArcGIS
    72                    'XYAHD_clip.txt']            # To extend boundary
     70point_filenames = ['Batemans_BBHD_MGA_1995.csv',
     71                   'SD100015374_MGA.csv',
     72                   'SD100016502_MGA.csv',
     73                   'SD100016517_MGA.csv',
     74                   'SD100031996_batemans_mga_clip_written.csv',
     75                   'tomaga_offshore_AHD_MGA_1997.csv']
     76         
    7377
    7478### Add csv header list to all files in point_filenames
     
    8589# INTERIOR REGIONS -  for designing the mesh
    8690# Used in run_model.py
    87 # Format for points easting,northing (no header)                   
    88 interior_regions_data = [['busselton_1km.csv', 500],
    89                          ['bunbury_1km.csv', 500],
    90                          ['busselton_2km.csv', 10000],
    91                          ['bunbury_2km.csv', 10000],
    92                          ['island1.csv', 10000],
    93                          ['island2.csv', 10000],
    94                          ['coast_5km_d20m.csv', 40000]]
     91# Format for points easting,northing (no header)
     92
     93interior_regions_data = [['area_of_interest.csv', 500],
     94                        ['area_of_significance.csv', 2500],
     95                        ['shallow_water.csv', 10000]]
     96
    9597
    9698# LAND - used to set the initial stage/water to be offcoast only
    9799# Used in run_model.py.  Format for points easting,northing (no header)
    98 land_initial_conditions_filename = [['initial_condition_extend.csv', 0],
    99                                     ['initial_condition_marina.csv', 0]]
     100land_initial_conditions_filename = [['initial_condition_extend.csv', 0]]
     101                                 
    100102
    101103# GAUGES - for creating timeseries at a specific point
    102104# Used in get_timeseries.py. 
    103105# Format easting,northing,name,elevation (with header)
    104 gauges_filename = 'gauges.csv'
     106##gauges_filename = 'gauges.csv'
    105107
    106108# BUILDINGS EXPOSURE - for identifying inundated houses
    107109# Used in run_building_inundation.py
    108110# Format latitude,longitude etc (geographic)
    109 building_exposure_filename = 'busselton_res_clip.csv' # from NEXIS
     111##building_exposure_filename = 'busselton_res_clip.csv' # from NEXIS
    110112
    111113# BOUNDING POLYGON - used in build_boundary.py and run_model.py respectively
     
    237239# The absolute pathname for the gauges file
    238240# Used for get_timeseries.py
    239 gauges = join(gauges_folder, gauges_filename)       
     241##gauges = join(gauges_folder, gauges_filename)       
    240242
    241243# The absolute pathname for the building file
    242244# Used for run_building_inundation.py
    243 building_exposure = join(gauges_folder, building_exposure_filename)
     245##building_exposure = join(gauges_folder, building_exposure_filename)
    244246
    245247# full path to where MUX files (or meta-files) live
  • anuga_work/production/new_south_wales/batemans_bay/run_model.py

    r6386 r6427  
    100100                                    interior_regions=project.interior_regions,
    101101                                    mesh_filename=project.meshes,
    102                                     use_cache=True,
     102                                    use_cache=False,
    103103                                    verbose=True)
    104104print domain.statistics()
Note: See TracChangeset for help on using the changeset viewer.