Changeset 4056


Ignore:
Timestamp:
Dec 5, 2006, 3:54:51 PM (17 years ago)
Author:
sexton
Message:

update sydney script based on survey data

Location:
anuga_work/production/sydney_2006
Files:
1 added
2 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/production/sydney_2006/project_slide.py

    r4046 r4056  
    3232
    3333# AHO
    34 offshore_name1 = 'XY100011610'
     34offshore_name1 = 'surveyAreaA'
     35offshore_name2 = 'surveyAreaB'
     36offshore_name3 = 'surveyAreaC'
    3537
    3638#swollen/ all data output
     
    5961nsw_dem_name = datadir + nsw100
    6062offshore_dem_name1 = datadir + offshore_name1
     63offshore_dem_name2 = datadir + offshore_name2
     64offshore_dem_name3 = datadir + offshore_name3
    6165combined_dem_name   = datadir + 'sydney_combined_elevation'
    6266
     
    6670
    6771# bounding polygon for study area
    68 polyAll = read_polygon(polygondir+'extent2.csv')
     72polyAll = read_polygon(polygondir+'extent3.csv')
    6973
    7074print 'Area of bounding polygon', polygon_area(polyAll)/1000000.0
     
    98102northingmin100 = 6227588
    99103northingmax100 = 6285975
     104
     105poly_surveyclip = read_polygon(polygondir+'surveyclip.csv')
    100106
    101107# clipping for nsw grid for conversion
  • anuga_work/production/sydney_2006/run_sydney_slide.py

    r4046 r4056  
    8282        use_cache=True, verbose=True)
    8383
    84 '''
    8584print 'create offshore'
    8685G1 = Geospatial_data(file_name = project_slide.offshore_dem_name1 + '.xya')+\
    8786     Geospatial_data(file_name = project_slide.offshore_dem_name2 + '.xya')+\
    8887     Geospatial_data(file_name = project_slide.offshore_dem_name3 + '.xya')+\
    89      Geospatial_data(file_name = project_slide.offshore_dem_name4 + '.xya')+\
    90      Geospatial_data(file_name = project_slide.offshore_dem_name5 + '.xya')+\
    91      Geospatial_data(file_name = project_slide.offshore_dem_name6 + '.xya')+\
    92      Geospatial_data(file_name = project_slide.offshore_dem_name7 + '.xya')+\
    93      Geospatial_data(file_name = project_slide.offshore_dem_name8 + '.xya')+\
    94      Geospatial_data(file_name = project_slide.offshore_dem_name9 + '.xya')+\
    95      Geospatial_data(file_name = project_slide.offshore_dem_name10 + '.xya')+\
    96      Geospatial_data(file_name = project_slide.offshore_dem_name11 + '.xya')+\
    97      Geospatial_data(file_name = project_slide.offshore_dem_name12 + '.xya')+\
    98      Geospatial_data(file_name = project_slide.offshore_dem_name13 + '.xya')+\
    99      Geospatial_data(file_name = project_slide.offshore_dem_name14 + '.xya')+\
    100      Geospatial_data(file_name = project_slide.offshore_dem_name15 + '.xya')+\
    101      Geospatial_data(file_name = project_slide.offshore_dem_name16 + '.xya')+\
    102      Geospatial_data(file_name = project_slide.offshore_dem_name17 + '.xya')+\
    103      Geospatial_data(file_name = project_slide.offshore_dem_name18 + '.xya')+\
    104      Geospatial_data(file_name = project_slide.offshore_dem_name19 + '.xya')+\
    105      Geospatial_data(file_name = project_slide.offshore_dem_name20 + '.xya')+\
    106      Geospatial_data(file_name = project_slide.offshore_dem_name21 + '.xya')+\
    107      Geospatial_data(file_name = project_slide.offshore_dem_name22 + '.xya')+\
    108      Geospatial_data(file_name = project_slide.offshore_interp_dem_name + '.pts')
    109 '''
    11088print 'create onshore'
    11189G2 = Geospatial_data(file_name = project_slide.on_offshore25_dem_name + '.pts')
     
    11391G4 = Geospatial_data(file_name = project_slide.nsw_dem_name + '.pts')
    11492print 'add'
    115 #G5
    116 G = G2.clip(Geospatial_data(project_slide.poly_25mclip)) +\
     93G = G1.clip(Geospatial(project_slide.poly_surveyclip)) +\
     94    G2.clip(Geospatial_data(project_slide.poly_25mclip)) +\
    11795    G3.clip(Geospatial_data(project_slide.poly_origsyd)) +\
    118     (G4.clip(Geospatial_data(project_slide.polyAll))).clip_outside(Geospatial_data(project_slide.poly_origsyd))
     96    (G4.clip(Geospatial_data(project_slide.polyAll))).clip_outside(Geospatial_data(project_slide.poly_surveyclip))
    11997print 'export points'
    12098G.export_points_file(project_slide.combined_dem_name + '.pts')
    121 G.export_points_file(project_slide.combined_dem_name + '.xya')
     99#G.export_points_file(project_slide.combined_dem_name + '.xya')
    122100
    123101#----------------------------------------------------------------------------
     
    128106
    129107from anuga.pmesh.mesh_interface import create_mesh_from_regions
    130 remainder_res = 1000000.
    131 local_res = 50000.
    132 coast_res = 10000.
     108remainder_res = 125000.
     109local_res = 30000.
     110coast_res = 500.
    133111interior_regions = [[project_slide.poly_syd1, local_res],
    134112                    [project_slide.poly_coast, coast_res]]
Note: See TracChangeset for help on using the changeset viewer.