Ignore:
Timestamp:
Feb 3, 2009, 11:51:59 AM (10 years ago)
Author:
ole
Message:

Looped elevation data conversions

File:
1 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/production/busselton/standardised_version/project.py

    r6256 r6261  
    66#------------------------------------------------------------------------------
    77
     8import os
    89from os import sep, environ, getenv, getcwd
    910from os.path import expanduser
     
    5253finaltime=80000         # final time for simulation
    5354
    54 setup='final'  # Final can be replaced with trial or basic.
     55setup='trial'  # Final can be replaced with trial or basic.
    5556               # Either will result in a coarser mesh that will allow a
    5657               # faster, but less accurate, simulation.
     
    5859if setup =='trial':
    5960    print'trial'
    60     res_factor=10
    61     time_thinning=48
     61    res_factor=100
     62    time_thinning=96
    6263    yieldstep=240
    6364if setup =='basic':
     
    102103#------------------------------------------------------------------------------
    103104
    104 # FIXME(Ole): Use lists of elevation data files as in Shark Bay study:
    105 ## elevation data filenames
    106 #ascii_grid_filenames = ['10m_dem_without_survey', '50m_dem_without_10m_dem',
    107 #                        'bathysteeppt', 'bathyleft', 'bathyright']
    108 #point_filenames = ['field_survey_north', 'field_survey_south',
    109 #                   'clipped_bathymetry_final', 'coast_points_final']
    110 
    111 
    112105
    113106# elevation data used in build_busselton.py
     
    140133landward = 'landward_bounding_polygon.csv'
    141134
     135
    142136#------------------------------------------------------------------------------
    143137# Output Elevation Data
     
    164158# Location of input and output data
    165159#------------------------------------------------------------------------------
    166 # where the input data sits
    167 onshore_in_dir_name = topographies_in_dir + onshore_name #topo
    168 coast_in_dir_name = topographies_in_dir + coast_name #coastline
    169 coast_in_dir_name1 = topographies_in_dir + coast_name1 #beach survey
    170 offshore_in_dir_name = topographies_in_dir + offshore_name #bathymetry
    171 offshore_in_dir_name1 = topographies_in_dir + offshore_name1 #bathymetry Charts
    172 offshore_in_dir_name2 = topographies_in_dir + offshore_name2 #Digitised Fairsheet
    173 offshore_in_dir_name3 = topographies_in_dir + offshore_name3 #250m
    174 offshore_in_dir_name4 = topographies_in_dir + offshore_name4 #Bunbury DPI
    175 offshore_in_dir_name5 = topographies_in_dir + offshore_name5 #Busselton Topo
    176 
    177 # where the output data sits
    178 onshore_dir_name = topographies_dir + onshore_name
    179 
    180 coast_dir_name = topographies_dir + coast_name
    181 coast_dir_name1 = topographies_dir + coast_name1
    182 
    183 offshore_dir_name = topographies_dir + offshore_name
    184 offshore_dir_name1 = topographies_dir + offshore_name1
    185 offshore_dir_name2 = topographies_dir + offshore_name2
    186 offshore_dir_name3 = topographies_dir + offshore_name3
    187 offshore_dir_name4 = topographies_dir + offshore_name4
    188 offshore_dir_name5 = topographies_dir + offshore_name5
    189 
    190 # where the combined elevation file sits
     160
     161ascii_grid_filenames = [onshore_name,   # Topo
     162                        offshore_name5] # Busselton Topo
     163
     164point_filenames = [coast_name,     # Coastline
     165                   coast_name1,    # Beach survey
     166                   offshore_name,  # Bathymetry
     167                   offshore_name1, # Bathymetry Charts
     168                   offshore_name2, # Digitised Fairsheet
     169                   offshore_name3, # 250m
     170                   offshore_name4] # Bunbury DPI
     171
     172
     173# Where the combined elevation file sits
    191174combined_dir_name = topographies_dir + combined_name
    192175combined_smaller_name_dir = topographies_dir + combined_smaller_name
    193176
    194 # where the mesh sits (this is created during the run_busselton.py)
     177# Where the mesh sits (this is created during the run_busselton.py)
    195178meshes_dir_name = meshes_dir + scenario_name+'.msh'
    196179
    197 # where the boundary ordering files sit (this is used within build_boundary.py)
     180# Where the boundary ordering files sit (this is used within build_boundary.py)
    198181order_filename_dir = boundaries_dir + order_filename
    199182
    200 # where the landward points of boundary extent sit (this is used within run_busselton.py)
     183# Where the landward points of boundary extent sit (this is used within run_busselton.py)
    201184landward_dir = boundaries_dir + landward
    202185
    203 # where the event sts files sits (this is created during the build_boundary.py)
     186# Where the event sts files sits (this is created during the build_boundary.py)
    204187boundaries_dir_event = boundaries_dir + str(event_number) + sep
    205188boundaries_dir_mux = muxhome
    206 
    207 # where the directory of the output filename sits
     189urs_boundary_name = os.path.join(boundaries_dir_event,
     190                                 scenario_name)
     191
     192
     193# Where the directory of the output filename sits
    208194output_build_time_dir = output_dir+build_time+dir_comment+sep   #used for build_busselton.py
    209195output_run_time_dir = output_dir+run_time+dir_comment+sep       #used for run_busselton.py
    210196output_run_time_dir_name = output_run_time_dir + scenario_name  #Used by post processing
    211197
    212 #w here the directory of the gauges sit
     198# Where the directory of the gauges sit
    213199gauges_dir_name = gauges_dir + gauge_name       #used for get_timeseries.py
    214200building_in_dir_name = gauges_dir + building + '.csv'    #used for run_building_inundation.py
     
    218204#------------------------------------------------------------------------------
    219205
    220 #Land, to set the initial stage/water to be offcoast only
     206# Land, to set the initial stage/water to be offcoast only
    221207poly_mainland = read_polygon(polygons_dir+'initial_condition.csv')
    222208
    223 #Land, to set the initial stage/water to be offcoast only
     209# Land, to set the initial stage/water to be offcoast only
    224210poly_marina = read_polygon(polygons_dir+'initial_condition_marina.csv')
    225211
     
    262248
    263249   
    264 trigs_min = number_mesh_triangles(interior_regions, poly_all, res_poly_all)
     250trigs_min = number_mesh_triangles(interior_regions, poly_all,
     251                                  res_poly_all)
    265252print 'min estimated number of triangles', trigs_min
    266253   
     
    270257#------------------------------------------------------------------------------
    271258
    272 # exporting asc grid for Busselton
     259# ASCII export grid for Busselton
    273260xminBusselton = 340000
    274261xmaxBusselton = 352000
     
    276263ymaxBusselton = 6280000
    277264
    278 # exporting asc grid for Bunbury
     265# ASCII export grid for Bunbury
    279266xminBunbury = 369000
    280267xmaxBunbury = 381000
Note: See TracChangeset for help on using the changeset viewer.