Changeset 6895


Ignore:
Timestamp:
Apr 24, 2009, 10:58:04 AM (14 years ago)
Author:
jgriffin
Message:
 
Location:
anuga_work/production/australia_ph2/bamaga
Files:
4 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/production/australia_ph2/bamaga/build_elevation.py

    r6742 r6895  
    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/australia_ph2/bamaga/build_urs_boundary.py

    r6760 r6895  
    139139                ordering_filename=project.urs_order,
    140140                weights=mux_weights,
    141                 central_meridian=project.central_meridian,
     141                #central_meridian=project.central_meridian,
    142142                zone=project.zone,
    143143                verbose=True)
  • anuga_work/production/australia_ph2/bamaga/project.py

    r6772 r6895  
    1717# Note, the user needs to set up the directory system accordingly
    1818state = 'australia_ph2'
    19 scenario_name = 'bamaga'
     19scenario_name = 'gulf'
    2020
    2121#-------------------------------------------------------------------------------
     
    2626# One or all can be changed each time the run_model script is executed
    2727
    28 central_meridian = 143.5 # Central meridian for projection (optional)
    29 zone = None
    30 import sys
    31 if len(sys.argv) > 1:
    32     event_number = int(sys.argv[1])
    33 else:   
    34     event_number = 31975    # the event number or the mux file name
    35 
    36 event_number_list = [31975, 51390, 63735] # To piggy back multiple events
     28#central_meridian = 143.5 # Central meridian for projection (optional)
     29zone = 54
     30##import sys
     31##if len(sys.argv) > 1:
     32##    event_number = int(sys.argv[1])
     33##else:
     34
     35event_number = 64977    # the event number or the mux file name
     36
     37#event_number_list = [31975, 51390, 63735] # To piggy back multiple events
    3738
    3839tide = 0                # difference between MSL and HAT
     
    4041friction=0.01           # manning's friction coefficient
    4142starttime=0             # start time for simulation
    42 finaltime=60000 # final time for simulation
     43finaltime=80000 # final time for simulation
    4344setup = 'final'         # This can be one of three values
    4445                        #    trial - coarsest mesh, fast
     
    4748
    4849# index is only used when wave = Tb
    49 index = 2000            # index from the PTHA - Y2000 0.257m
     50index = 1914            # index from the PTHA - Y2000 0.257m
    5051wave = 'Tb'             # Bf (sts wave) Tb (index wave)
    5152
     
    7879# Used in build_elevation.py
    7980# Format for ascii grids, as produced in ArcGIS + a projection file
    80 ascii_grid_filenames = ['grid250m_pro'] # 250m grid 2005
     81ascii_grid_filenames = ['gulf_bg'] # 250m grid 2005
    8182
    8283# Format for points is x,y,elevation (with header)
  • anuga_work/production/australia_ph2/bamaga/run_model.py

    r6742 r6895  
    9292# Boundary tags refer to project.landward_boundary
    9393# 4 points equals 5 segments start at N
    94 boundary_tags={'back': range(num_ocean_segments+1,
    95                              num_ocean_segments+num_land_points),
    96                'side': [num_ocean_segments,
    97                         num_ocean_segments+num_land_points],
    98                'ocean': range(num_ocean_segments)}
     94##boundary_tags={'back': range(num_ocean_segments+1,
     95##                             num_ocean_segments+num_land_points),
     96##               'side': [num_ocean_segments,
     97##                        num_ocean_segments+num_land_points],
     98##               'ocean': range(num_ocean_segments)}
     99boundary_tags = {'back':[3],
     100                 'side':[2,4],
     101                 'ocean':[0,1,5]}
     102
     103print 'boundary tags', boundary_tags
    99104
    100105# Build mesh and domain
Note: See TracChangeset for help on using the changeset viewer.