Changeset 5811


Ignore:
Timestamp:
Oct 2, 2008, 4:07:08 PM (15 years ago)
Author:
rwilson
Message:

Doesn't work, yet.

File:
1 edited

Legend:

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

    r5810 r5811  
    9292#creates pts file for onshore DEM
    9393print "creates pts file for onshore DEM"
    94 dem2pts(offshore_dir_name ,use_cache=True,verbose=True)
     94dem2pts(onshore_dir_name ,use_cache=True,verbose=True)
    9595dem2pts(offshore_dir_name1 ,use_cache=True,verbose=True)
    9696dem2pts(offshore_dir_name2 ,use_cache=True,verbose=True)
     
    120120G_off = Geospatial_data(file_name = offshore_in_dir_name)
    121121print'create Geospatial data8 objects from topographies'
    122 G_off1 = Geospatial_data(file_name = offshore_in_dir_name1)
     122G_off1 = Geospatial_data(file_name = offshore_in_dir_name1 + '.pts')
    123123print'create Geospatial data9 objects from topographies'
    124 G_off2 = Geospatial_data(file_name = offshore_in_dir_name2)
     124G_off2 = Geospatial_data(file_name = offshore_in_dir_name2 + '.pts')
    125125print'create Geospatial data10 objects from topographies'
    126 G_off3 = Geospatial_data(file_name = offshore_in_dir_name3)
     126G_off3 = Geospatial_data(file_name = offshore_in_dir_name3 + '.pts')
    127127print'create Geospatial data10 objects from topographies'
    128 G_off3 = Geospatial_data(file_name = offshore_in_dir_name4)
     128G_off4 = Geospatial_data(file_name = offshore_in_dir_name4 + '.pts')
    129129print'create Geospatial data10 objects from topographies'
    130 G_off3 = Geospatial_data(file_name = offshore_in_dir_name5)
     130G_off5 = Geospatial_data(file_name = offshore_in_dir_name5 + '.pts')
    131131print'create Geospatial data10 objects from topographies'
    132 G_off3 = Geospatial_data(file_name = offshore_in_dir_name6)
     132G_off6 = Geospatial_data(file_name = offshore_in_dir_name6 + '.pts')
    133133print'add all geospatial objects'
    134 G = G1 + G_coast + G_off + G_off1 + G_off2 + G_off3
     134G = G1 + G_coast + G_off + G_off1 + G_off2 + G_off3  + G_off4 + G_off5 + G_off6
    135135print'clip combined geospatial object by bounding polygon'
    136136G_clipped = G.clip(project.poly_all)
Note: See TracChangeset for help on using the changeset viewer.