Changeset 6027 for anuga_work/production/geraldton/build_geraldton.py
- Timestamp:
- Nov 28, 2008, 11:48:45 AM (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_work/production/geraldton/build_geraldton.py
r5789 r6027 92 92 G_coast = Geospatial_data(file_name = coast_in_dir_name) 93 93 print'create Geospatial data7 objects from topographies' 94 G_off = Geospatial_data(file_name = offshore_in_dir_name)94 G_off1 = Geospatial_data(file_name = offshore_in_dir_name) 95 95 print'create Geospatial data8 objects from topographies' 96 G_off 1= Geospatial_data(file_name = offshore_in_dir_name1)96 G_off2 = Geospatial_data(file_name = offshore_in_dir_name1) 97 97 print'create Geospatial data9 objects from topographies' 98 G_off 2= Geospatial_data(file_name = offshore_in_dir_name2)98 G_off3 = Geospatial_data(file_name = offshore_in_dir_name2) 99 99 print'create Geospatial data10 objects from topographies' 100 G_off 3= Geospatial_data(file_name = offshore_in_dir_name3)100 G_off4 = Geospatial_data(file_name = offshore_in_dir_name3) 101 101 102 102 #------------------------------------------------------------------------------- … … 105 105 106 106 print 'clipping port data to area of significance' 107 G_off2_clip = G_off2.clip(project.poly_ CBD_1km, verbose=True)107 G_off2_clip = G_off2.clip(project.poly_buffer_20m, verbose=True) 108 108 print 'add all bathymetry files' 109 109 G_off = G_off1 + G_off2_clip + G_off3 + G_off4 … … 113 113 114 114 print'add all geospatial objects' 115 G = G1_clip + G2 + G 3 + G_off_clip115 G = G1_clip + G2 + G_off_clip 116 116 117 117 print'clip combined geospatial object by bounding polygon'
Note: See TracChangeset
for help on using the changeset viewer.