Changeset 4058


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

updates for all slide scenarios

Location:
anuga_work/production
Files:
6 edited

Legend:

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

    r4046 r4058  
    2727# onshore data provided by NSW LPI, Krishna merged with existing
    2828# PMD data to make 100m DEM
    29 on_offshore_name = ''
     29nsw100 = 'dem_utm'
    3030
    31 # AHO
    32 offshore_name1 = 'XY100011610'
     31# survey data
     32offshore_name2 = 'surveyAreaB'
     33offshore_name3 = 'surveyAreaC'
     34
     35# AHO data
     36offshore_name4 = '1000003611export'
     37offshore_name5 = '1000003613export'
     38offshore_name6 = '1000003614export'
     39offshore_name7 = '1000003627export'
     40offshore_name8 = '1000003628export'
     41offshore_name9 = 'AHDexport'
    3342
    3443#swollen/ all data output
     
    5463# Necessary if using point datasets, rather than grid
    5564on_offshore_dem_name = datadir + on_offshore_name
    56 offshore_dem_name1 = datadir + offshore_name1
     65offshore_dem_name2 = datadir + offshore_name2
     66offshore_dem_name3 = datadir + offshore_name3
     67offshore_dem_name4 = datadir + offshore_name4
     68offshore_dem_name5 = datadir + offshore_name5
     69offshore_dem_name6 = datadir + offshore_name6
     70offshore_dem_name7 = datadir + offshore_name7
     71offshore_dem_name8 = datadir + offshore_name8
     72offshore_dem_name9 = datadir + offshore_name9
    5773combined_dem_name   = datadir + 'newcastle_combined_elevation'
    5874
     
    8298# Clipping regions for export to asc and regions for clipping data
    8399###################################################################
     100
     101# clipping for nsw grid for conversion
     102eastingmin_nsw = 546000
     103eastingmax_nsw = 350390
     104northingmin_nsw = 6253970
     105northingmax_nsw = 6424600
    84106
    85107# exporting asc grid
  • anuga_work/production/newcastle_2006/run_newcastle_slide.py

    r4036 r4058  
    5252
    5353# filenames
    54 on_offshore_dem_name = project.on_offshore_dem_name
    55 meshname = project.meshname+'.msh'
     54nsw_dem_name = project_slide.nsw_dem_name
     55meshname = project_slide.meshname+'.msh'
    5656
    5757# creates DEM from asc data
    58 convert_dem_from_ascii2netcdf(on_offshore_dem_name, use_cache=True, verbose=True)
     58convert_dem_from_ascii2netcdf(nsw_dem_name, use_cache=True, verbose=True)
    5959
    6060#creates pts file for onshore DEM
    61 dem2pts(on_offshore_dem_name, use_cache=True, verbose=True)
     61dem2pts(nsw_dem_name,
     62        easting_min=project_slide.eastingmin_nsw,
     63        easting_max=project_slide.eastingmax_nsw,
     64        northing_min=project_slide.northingmin_nsw,
     65        northing_max= project_slide.northingmax_nsw,
     66        use_cache=True, verbose=True)
    6267
    6368print 'create offshore'
    64 G1 = Geospatial_data(file_name = project.offshore_dem_name1 + '.xya')+\
    65      Geospatial_data(file_name = project.offshore_dem_name2 + '.xya')+\
    66      Geospatial_data(file_name = project.offshore_dem_name3 + '.xya')+\
    67      Geospatial_data(file_name = project.offshore_dem_name4 + '.xya')+\
    68      Geospatial_data(file_name = project.offshore_dem_name5 + '.xya')+\
    69      Geospatial_data(file_name = project.offshore_dem_name6 + '.xya')+\
    70      Geospatial_data(file_name = project.offshore_dem_name7 + '.xya')+\
    71      Geospatial_data(file_name = project.offshore_dem_name8 + '.xya')+\
    72      Geospatial_data(file_name = project.offshore_dem_name9 + '.xya')+\
    73      Geospatial_data(file_name = project.offshore_dem_name10 + '.xya')+\
    74      Geospatial_data(file_name = project.offshore_dem_name11 + '.xya')+\
    75      Geospatial_data(file_name = project.offshore_dem_name12 + '.xya')+\
    76      Geospatial_data(file_name = project.offshore_dem_name13 + '.xya')+\
    77      Geospatial_data(file_name = project.offshore_dem_name14 + '.xya')+\
    78      Geospatial_data(file_name = project.offshore_dem_name15 + '.xya')+\
    79      Geospatial_data(file_name = project.offshore_dem_name16 + '.xya')+\
    80      Geospatial_data(file_name = project.offshore_dem_name17 + '.xya')+\
    81      Geospatial_data(file_name = project.offshore_dem_name18 + '.xya')+\
    82      Geospatial_data(file_name = project.offshore_dem_name19 + '.xya')+\
    83      Geospatial_data(file_name = project.offshore_dem_name20 + '.xya')+\
    84      Geospatial_data(file_name = project.offshore_dem_name21 + '.xya')+\
    85      Geospatial_data(file_name = project.offshore_dem_name22 + '.xya')+\
    86      Geospatial_data(file_name = project.offshore_interp_dem_name + '.pts')
     69G11 = Geospatial_data(file_name = project_slide.offshore_dem_name2 + '.xya')+\
     70      Geospatial_data(file_name = project_slide.offshore_dem_name3 + '.xya')
     71G12 = Geospatial_data(file_name = project_slide.offshore_dem_name4 + '.xya')+\
     72      Geospatial_data(file_name = project_slide.offshore_dem_name5 + '.xya')+\
     73      Geospatial_data(file_name = project_slide.offshore_dem_name6 + '.xya')+\
     74      Geospatial_data(file_name = project_slide.offshore_dem_name7 + '.xya')+\
     75      Geospatial_data(file_name = project_slide.offshore_dem_name8 + '.xya')+\
     76      Geospatial_data(file_name = project_slide.offshore_dem_name9 + '.xya')
    8777print 'create onshore'
    88 G2 = Geospatial_data(file_name = project.on_offshore_dem_name + '.pts')
     78G4 = Geospatial_data(file_name = project_slide.nsw_dem_name + '.pts')
    8979print 'add'
    90 G = G1 + G2 + G3
     80G = G11.clip(Geospatial_data(project_slide.poly_surveyclip)) +\
     81    G12.clip(Geospatial_data(project_slide.polyAll))
     82    (G4.clip(Geospatial_data(project_slide.polyAll))).clip_outside(Geospatial_data(project_slide.poly_surveyclip))
    9183print 'export points'
    92 G.export_points_file(project.combined_dem_name + '.pts')
    93 G.export_points_file(project.combined_dem_name + '.xya')
     84G.export_points_file(project_slide.combined_dem_name + '.pts')
     85#G.export_points_file(project_slide.combined_dem_name + '.xya')
     86
    9487
    9588#----------------------------------------------------------------------------
  • 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')
  • anuga_work/production/wollongong_2006/project_slide.py

    r4046 r4058  
    2727# onshore data provided by NSW LPI, Krishna merged with existing
    2828# PMD data to make 100m DEM
    29 on_offshore_name = ''
     29on_offshore10_name = 'wollongong_10'
     30nsw100 = 'dem_utm'
    3031
    31 # AHO
    32 offshore_name1 = 'XY100011610'
     32# survey data
     33offshore_name1 = 'surveyAreaA'
     34
     35# AHO data
     36offshore_name4 = '1000003611export'
     37offshore_name5 = '1000003613export'
     38offshore_name6 = '1000003614export'
     39offshore_name7 = '1000003627export'
     40offshore_name8 = '1000003628export'
     41offshore_name9 = 'AHDexport'
    3342
    3443#swollen/ all data output
     
    5564on_offshore_dem_name = datadir + on_offshore_name
    5665offshore_dem_name1 = datadir + offshore_name1
     66offshore_dem_name4 = datadir + offshore_name4
     67offshore_dem_name5 = datadir + offshore_name5
     68offshore_dem_name6 = datadir + offshore_name6
     69offshore_dem_name7 = datadir + offshore_name7
     70offshore_dem_name8 = datadir + offshore_name8
     71offshore_dem_name9 = datadir + offshore_name9
    5772combined_dem_name   = datadir + 'gong_combined_elevation'
    5873
     
    6277
    6378# bounding polygon for study area
    64 polyAll = read_polygon(polygondir+'extentA.csv')
     79polyAll = read_polygon(polygondir+'extentA_up.csv')
    6580
    6681print 'Area of bounding polygon', polygon_area(polyAll)/1000000.0
     
    7085###############################
    7186
    72 # broome digitized polygons
    73 poly_gong1 = read_polygon(polygondir+'.csv')
    74 poly_gong2 = read_polygon(polygondir+'.csv')
    75 poly_gong3 = read_polygon(polygondir+'.csv')
     87# areaA digitized polygons
     88poly_local = read_polygon(polygondir+'local.csv')
     89poly_gong = read_polygon(polygondir+'gong.csv')
    7690
    77 print 'Area of local polygon', polygon_area(poly_gong1)/1000000.0
    78 print 'Area of close polygon', polygon_area(poly_gong2)/1000000.0
    79 print 'Area of coastal polygon', polygon_area(poly_gong3)/1000000.0
     91print 'Area of local polygon', polygon_area(poly_local)/1000000.0
    8092
    8193###################################################################
     
    8395###################################################################
    8496
     97poly_10mclip = read_polygon(polygondir+'10mclip.csv')
     98poly_surveyclip = read_polygon(polygondir+'surveyclip.csv')
     99
     100# clipping for nsw grid for conversion
     101eastingmin_nsw = 272900
     102eastingmax_nsw = 463550
     103northingmin_nsw = 6112950
     104northingmax_nsw = 6269750
     105
    85106# exporting asc grid
    86 eastingmin =
    87 eastingmax =
    88 northingmin =
    89 northingmax =
     107#eastingmin =
     108#eastingmax =
     109#northingmin =
     110#northingmax =
    90111
    91112###################################################################
     
    103124bulli_density = 1.46
    104125bulli_slope = 4.0
    105 bulli_alpha = 126.0 + 90.0
     126bulli_alpha = 126.0 - 90.0
    106127
    107128shovel_depth = 968.0
     
    111132shovel_density = 1.49
    112133shove_slope = 4.0
    113 shovel_alpha = 118.0 + 90.0
     134shovel_alpha = 118.0 - 90.0
    114135
    115136yacaaba_depth = 1119.0
     
    119140yacaaba_density = 1.48
    120141yacaaba_slope = 3.0
    121 yacaaba_alpha = 133.0 + 90.0
     142yacaaba_alpha = 133.0 - 90.0
    122143
  • anuga_work/production/wollongong_2006/run_gong_slide.py

    r4036 r4058  
    5252
    5353# filenames
    54 on_offshore_dem_name = project.on_offshore_dem_name
    55 meshname = project.meshname+'.msh'
     54on_offshore10_dem_name = project_slide.on_offshore10_dem_name
     55nsw_dem_name = project_slide.nsw_dem_name
     56meshname = project_slide.meshname+'.msh'
    5657
    5758# creates DEM from asc data
    58 convert_dem_from_ascii2netcdf(on_offshore_dem_name, use_cache=True, verbose=True)
     59convert_dem_from_ascii2netcdf(on_offshore10_dem_name, use_cache=True, verbose=True)
     60convert_dem_from_ascii2netcdf(nsw_dem_name, use_cache=True, verbose=True)
    5961
    6062#creates pts file for onshore DEM
    61 dem2pts(on_offshore_dem_name, use_cache=True, verbose=True)
     63dem2pts(on_offshore10_dem_name, use_cache=True, verbose=True)
     64dem2pts(nsw_dem_name,
     65        easting_min=project_slide.eastingmin_nsw,
     66        easting_max=project_slide.eastingmax_nsw,
     67        northing_min=project_slide.northingmin_nsw,
     68        northing_max= project_slide.northingmax_nsw,
     69        use_cache=True, verbose=True)
    6270
    6371print 'create offshore'
    64 G1 = Geospatial_data(file_name = project.offshore_dem_name1 + '.xya')+\
    65      Geospatial_data(file_name = project.offshore_dem_name2 + '.xya')+\
    66      Geospatial_data(file_name = project.offshore_dem_name3 + '.xya')+\
    67      Geospatial_data(file_name = project.offshore_dem_name4 + '.xya')+\
    68      Geospatial_data(file_name = project.offshore_dem_name5 + '.xya')+\
    69      Geospatial_data(file_name = project.offshore_dem_name6 + '.xya')+\
    70      Geospatial_data(file_name = project.offshore_dem_name7 + '.xya')+\
    71      Geospatial_data(file_name = project.offshore_dem_name8 + '.xya')+\
    72      Geospatial_data(file_name = project.offshore_dem_name9 + '.xya')+\
    73      Geospatial_data(file_name = project.offshore_dem_name10 + '.xya')+\
    74      Geospatial_data(file_name = project.offshore_dem_name11 + '.xya')+\
    75      Geospatial_data(file_name = project.offshore_dem_name12 + '.xya')+\
    76      Geospatial_data(file_name = project.offshore_dem_name13 + '.xya')+\
    77      Geospatial_data(file_name = project.offshore_dem_name14 + '.xya')+\
    78      Geospatial_data(file_name = project.offshore_dem_name15 + '.xya')+\
    79      Geospatial_data(file_name = project.offshore_dem_name16 + '.xya')+\
    80      Geospatial_data(file_name = project.offshore_dem_name17 + '.xya')+\
    81      Geospatial_data(file_name = project.offshore_dem_name18 + '.xya')+\
    82      Geospatial_data(file_name = project.offshore_dem_name19 + '.xya')+\
    83      Geospatial_data(file_name = project.offshore_dem_name20 + '.xya')+\
    84      Geospatial_data(file_name = project.offshore_dem_name21 + '.xya')+\
    85      Geospatial_data(file_name = project.offshore_dem_name22 + '.xya')+\
    86      Geospatial_data(file_name = project.offshore_interp_dem_name + '.pts')
     72G11 = Geospatial_data(file_name = project_slide.offshore_dem_name1 + '.xya')
     73G12 = Geospatial_data(file_name = project_slide.offshore_dem_name4 + '.xya')+\
     74      Geospatial_data(file_name = project_slide.offshore_dem_name5 + '.xya')+\
     75      Geospatial_data(file_name = project_slide.offshore_dem_name6 + '.xya')+\
     76      Geospatial_data(file_name = project_slide.offshore_dem_name7 + '.xya')+\
     77      Geospatial_data(file_name = project_slide.offshore_dem_name8 + '.xya')+\
     78      Geospatial_data(file_name = project_slide.offshore_dem_name9 + '.xya')
    8779print 'create onshore'
    88 G2 = Geospatial_data(file_name = project.on_offshore_dem_name + '.pts')
     80G2 = Geospatial_data(file_name = project_slide.on_offshore10_dem_name + '.pts')
     81G4 = Geospatial_data(file_name = project_slide.nsw_dem_name + '.pts')
    8982print 'add'
    90 G = G1 + G2 + G3
     83G = G11.clip(Geospatial(project_slide.poly_surveyclip)) +\
     84    G12.clip(Geospatial_data(project_slide.polyAll)) +\
     85    G2.clip(Geospatial_data(project_slide.poly_10mclip)) +\
     86    (G4.clip(Geospatial_data(project_slide.polyAll))).clip_outside(Geospatial_data(project_slide.poly_surveyclip)).clip_outside(Geospatial_data(project_slide.poly_10mclip))
    9187print 'export points'
    92 G.export_points_file(project.combined_dem_name + '.pts')
    93 G.export_points_file(project.combined_dem_name + '.xya')
     88G.export_points_file(project_slide.combined_dem_name + '.pts')
     89#G.export_points_file(project_slide.combined_dem_name + '.xya')
    9490
    9591#----------------------------------------------------------------------------
Note: See TracChangeset for help on using the changeset viewer.