Changeset 4358


Ignore:
Timestamp:
Apr 4, 2007, 12:44:04 PM (17 years ago)
Author:
nick
Message:

update to broome

Location:
anuga_work/production/broome_2006
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/production/broome_2006/build_broome.py

    r4290 r4358  
    8181
    8282#creates pts file for island DEM
    83 dem2pts(offshore_in_dir_name1, use_cache=True, verbose=True)
    84 dem2pts(offshore_in_dir_name2, use_cache=True, verbose=True)
     83dem2pts(offshore_dir_name1, use_cache=True, verbose=True)
     84dem2pts(offshore_dir_name2, use_cache=True, verbose=True)
    8585
    8686print'create Geospatial data1 objects from topographies'
     
    9393G_off = Geospatial_data(file_name = offshore_in_dir_name + '.txt')
    9494print'create Geospatial data4 objects from offshore1'
    95 G_off1 = Geospatial_data(file_name = offshore_in_dir_name1 + '.pts')
     95G_off1 = Geospatial_data(file_name = offshore_dir_name1 + '.pts')
    9696print'create Geospatial data4 objects from offshore2'
    97 G_off2 = Geospatial_data(file_name = offshore_in_dir_name2 + '.pts')
     97G_off2 = Geospatial_data(file_name = offshore_dir_name2 + '.pts')
    9898
    9999print'add all geospatial objects'
     
    101101
    102102print'clip combined geospatial object by bounding polygon'
    103 G_clipped = G.clip(project.bounding_polygon)
     103G_clipped = G.clip(project.poly_all)
    104104#FIXME: add a clip function to pts
    105105#print'shape of clipped data', G_clipped.get_data_points().shape
    106106
    107107print'split'
    108 G_main, G_small = G_clipped.split(.10)
     108G_small, G_other = G_clipped.split(.10,verbose=True)
    109109
    110110print'export combined DEM file'
  • anuga_work/production/broome_2006/project.py

    r4290 r4358  
    119119
    120120poly_all = read_polygon(polygons_dir+'extent.csv')
    121 res_poly_all = 500000
     121res_poly_all = 150000
    122122
    123123###############################
     
    125125###############################
    126126
    127 poly_1 = read_polygon(polygons_dir+'broome_north_coast.csv')
    128 res_1 = 50000
     127poly_0 = read_polygon(polygons_dir+'neg20_coast_contour_pts.csv')
     128res_0 = 10000
     129poly_1 = read_polygon(polygons_dir+'broome_north_coast_inside_extent.csv')
     130res_1 = 5000
    129131
    130 poly_2 = read_polygon(polygons_dir+'broome_south_coast.csv')
    131 res_2 = 50000
     132poly_2 = read_polygon(polygons_dir+'broome_south_coast_inside_extent.csv')
     133res_2 = 5000
    132134
    133135poly_3 = read_polygon(polygons_dir+'Broome_town_pts.csv')
    134 res_3 = 10000
     136res_3 = 2000
    135137
    136138poly_4 = read_polygon(polygons_dir+'Broome_inner_town_pts.csv')
    137 res_4 = 1000
     139res_4 = 500
    138140#assert zone == refzone
    139141
    140 interior_regions = [[poly_1,res_1],[poly_2,res_2]
     142interior_regions = [[poly_0,res_0],[poly_1,res_1],[poly_2,res_2]
    141143                     ,[poly_3,res_3],[poly_4,res_4]]
    142144
Note: See TracChangeset for help on using the changeset viewer.