Ignore:
Timestamp:
Dec 5, 2006, 5:30:41 PM (18 years ago)
Author:
sexton
Message:

updates for all slide scenarios

Location:
anuga_work/production/sydney_2006
Files:
2 edited

Legend:

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

    r4056 r4058  
    3131nsw100 = 'dem_utm'
    3232
    33 # AHO
     33# survey data
    3434offshore_name1 = 'surveyAreaA'
    3535offshore_name2 = 'surveyAreaB'
    3636offshore_name3 = 'surveyAreaC'
     37
     38# AHO data
     39offshore_name4 = '1000003611export'
     40offshore_name5 = '1000003613export'
     41offshore_name6 = '1000003614export'
     42offshore_name7 = '1000003627export'
     43offshore_name8 = '1000003628export'
     44offshore_name9 = 'AHDexport'
    3745
    3846#swollen/ all data output
     
    6371offshore_dem_name2 = datadir + offshore_name2
    6472offshore_dem_name3 = datadir + offshore_name3
     73offshore_dem_name4 = datadir + offshore_name4
     74offshore_dem_name5 = datadir + offshore_name5
     75offshore_dem_name6 = datadir + offshore_name6
     76offshore_dem_name7 = datadir + offshore_name7
     77offshore_dem_name8 = datadir + offshore_name8
     78offshore_dem_name9 = datadir + offshore_name9
    6579combined_dem_name   = datadir + 'sydney_combined_elevation'
    6680
     
    134148bulli_density = 1.46
    135149bulli_slope = 4.0
    136 bulli_alpha = 126.0 + 90.0
     150bulli_alpha = 126.0 - 90.0
    137151
    138152shovel_depth = 968.0
     
    142156shovel_density = 1.49
    143157shove_slope = 4.0
    144 shovel_alpha = 118.0 + 90.0
     158shovel_alpha = 118.0 - 90.0
    145159
    146160yacaaba_depth = 1119.0
     
    150164yacaaba_density = 1.48
    151165yacaaba_slope = 3.0
    152 yacaaba_alpha = 133.0 + 90.0
     166yacaaba_alpha = 133.0 - 90.0
  • anuga_work/production/sydney_2006/run_sydney_slide.py

    r4056 r4058  
    8383
    8484print 'create offshore'
    85 G1 = Geospatial_data(file_name = project_slide.offshore_dem_name1 + '.xya')+\
    86      Geospatial_data(file_name = project_slide.offshore_dem_name2 + '.xya')+\
    87      Geospatial_data(file_name = project_slide.offshore_dem_name3 + '.xya')+\
     85G11 = Geospatial_data(file_name = project_slide.offshore_dem_name1 + '.xya')+\
     86      Geospatial_data(file_name = project_slide.offshore_dem_name2 + '.xya')+\
     87      Geospatial_data(file_name = project_slide.offshore_dem_name3 + '.xya')
     88G12 = Geospatial_data(file_name = project_slide.offshore_dem_name4 + '.xya')+\
     89      Geospatial_data(file_name = project_slide.offshore_dem_name5 + '.xya')+\
     90      Geospatial_data(file_name = project_slide.offshore_dem_name6 + '.xya')+\
     91      Geospatial_data(file_name = project_slide.offshore_dem_name7 + '.xya')+\
     92      Geospatial_data(file_name = project_slide.offshore_dem_name8 + '.xya')+\
     93      Geospatial_data(file_name = project_slide.offshore_dem_name9 + '.xya')
    8894print 'create onshore'
    8995G2 = Geospatial_data(file_name = project_slide.on_offshore25_dem_name + '.pts')
     
    9197G4 = Geospatial_data(file_name = project_slide.nsw_dem_name + '.pts')
    9298print 'add'
    93 G = G1.clip(Geospatial(project_slide.poly_surveyclip)) +\
     99G = G11.clip(Geospatial_data(project_slide.poly_surveyclip)) +\
     100    G12.clip(Geospatial_data(project_slide.polyAll))
    94101    G2.clip(Geospatial_data(project_slide.poly_25mclip)) +\
    95102    G3.clip(Geospatial_data(project_slide.poly_origsyd)) +\
    96     (G4.clip(Geospatial_data(project_slide.polyAll))).clip_outside(Geospatial_data(project_slide.poly_surveyclip))
     103    (G4.clip(Geospatial_data(project_slide.polyAll))).clip_outside(Geospatial_data(project_slide.poly_surveyclip)).clip_outside(Geospatial_data(project_slide.poly_origsyd))
    97104print 'export points'
    98105G.export_points_file(project_slide.combined_dem_name + '.pts')
Note: See TracChangeset for help on using the changeset viewer.