Ignore:
Timestamp:
Mar 20, 2009, 4:02:42 PM (15 years ago)
Author:
kristy
Message:

added datasets for bunbury

File:
1 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/production/busselton/build_busselton.py

    r6062 r6560  
    6262offshore_in_dir_name4 = project.offshore_in_dir_name4
    6363offshore_in_dir_name5 = project.offshore_in_dir_name5
     64offshore_in_dir_name6 = project.offshore_in_dir_name6
     65offshore_in_dir_name7 = project.offshore_in_dir_name7
    6466
    6567# output elevation directory filenames
     
    7375offshore_dir_name4 = project.offshore_dir_name4
    7476offshore_dir_name5 = project.offshore_dir_name5
     77offshore_dir_name6 = project.offshore_dir_name6
     78offshore_dir_name7 = project.offshore_dir_name7
     79
    7580# creates DEM from asc data
    7681print "creates DEMs from asc data"
    7782convert_dem_from_ascii2netcdf(onshore_in_dir_name, basename_out=onshore_dir_name, use_cache=True, verbose=True)
    78 convert_dem_from_ascii2netcdf(offshore_in_dir_name5, basename_out=offshore_dir_name5, use_cache=True, verbose=True)
     83convert_dem_from_ascii2netcdf(offshore_in_dir_name7, basename_out=offshore_dir_name7, use_cache=True, verbose=True)
    7984
    8085# creates pts file for onshore DEM
    8186print "creates pts file for onshore DEM"
    8287dem2pts(onshore_dir_name ,use_cache=True,verbose=True)
    83 dem2pts(offshore_dir_name5 ,use_cache=True,verbose=True)
     88dem2pts(offshore_dir_name7 ,use_cache=True,verbose=True)
    8489
    8590# create onshore pts files
     
    103108G_off4 = Geospatial_data(file_name = offshore_in_dir_name4)
    104109print'create Geospatial data9 objects from topographies'
    105 G_off5 = Geospatial_data(file_name = offshore_dir_name5 + '.pts')
     110G_off5 = Geospatial_data(file_name = offshore_in_dir_name5)
     111print'create Geospatial data9 objects from topographies'
     112G_off6 = Geospatial_data(file_name = offshore_in_dir_name6)
     113print'create Geospatial data9 objects from topographies'
     114G_off7 = Geospatial_data(file_name = offshore_dir_name7 + '.pts')
     115
     116
    106117
    107118
     
    111122
    112123print'add all geospatial objects'
    113 G = G1 + G2 + G3 + G_off + G_off1 + G_off2 + G_off3 + G_off4
     124G = G1 + G2 + G3 + G_off + G_off1 + G_off2 + G_off3 + G_off4 
    114125
    115126print'clip combined geospatial object by bounding polygon'
    116127G_clip = G.clip_outside(project.poly_aoi1)
    117 G_all = G_clip + G_off5
     128G_all = G_clip + G_off5 + G_off6 + G_off7
    118129G_clipped = G_all.clip(project.poly_all)
    119130
Note: See TracChangeset for help on using the changeset viewer.