Changeset 5646


Ignore:
Timestamp:
Aug 13, 2008, 8:22:25 AM (16 years ago)
Author:
kristy
Message:

updated with more data within the 500m2 mesh area.

File:
1 edited

Legend:

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

    r5525 r5646  
    6262offshore_in_dir_name2 = project.offshore_in_dir_name2
    6363offshore_in_dir_name3 = project.offshore_in_dir_name3
     64offshore_in_dir_name4 = project.offshore_in_dir_name4
     65offshore_in_dir_name5 = project.offshore_in_dir_name5
    6466
    6567onshore_dir_name = project.onshore_dir_name
     
    7072offshore_dir_name2 = project.offshore_dir_name2
    7173offshore_dir_name3 = project.offshore_dir_name3
     74offshore_dir_name4 = project.offshore_dir_name4
     75offshore_dir_name5 = project.offshore_dir_name5
    7276
    7377# creates DEM from asc data
     
    96100G_off2 = Geospatial_data(file_name = offshore_in_dir_name2 + '.txt',verbose=True)
    97101G_off3 = Geospatial_data(file_name = offshore_in_dir_name3 + '.txt',verbose=True)
    98                          
     102G_off4 = Geospatial_data(file_name = offshore_in_dir_name4 + '.txt',verbose=True)
     103G_off5 = Geospatial_data(file_name = offshore_in_dir_name5 + '.txt',verbose=True)
    99104print'add all geospatial objects'
    100 G = G1 + G2 + G3 + G_off + G_off1 + G_off2 + G_off3
     105G = G1 + G2 + G3 + G_off + G_off1 + G_off2 + G_off3 + G_off4 + G_off5
    101106
    102107print'clip combined geospatial object by bounding polygon'
     
    110115print'project.combined_dir_name + .txt',project.combined_dir_name + '.txt'
    111116
    112 #-------------------------------------------------------------------------
    113 # Convert URS to SWW file for boundary conditions
    114 #-------------------------------------------------------------------------
    115 print 'starting to create boundary conditions'
    116 from anuga.shallow_water.data_manager import urs2sww, urs_ungridded2sww
    117 
    118 boundaries_in_dir_name = project.boundaries_in_dir_name
    119 print 'boundaries_in_dir_name',project.boundaries_in_dir_name
    120 
    121 
    122 #import sys; sys.exit()
    123 
    124 urs_ungridded2sww(project.boundaries_in_dir_name, project.boundaries_in_dir_name,
    125                   verbose=True, mint=4000, maxt=80000, zscale=1)
    126 
     117###-------------------------------------------------------------------------
     118### Convert URS to SWW file for boundary conditions
     119###-------------------------------------------------------------------------
     120##print 'starting to create boundary conditions'
     121##from anuga.shallow_water.data_manager import urs2sww, urs_ungridded2sww
     122##
     123##boundaries_in_dir_name = project.boundaries_in_dir_name
     124##print 'boundaries_in_dir_name',project.boundaries_in_dir_name
     125##
     126##
     127###import sys; sys.exit()
     128##
     129##urs_ungridded2sww(project.boundaries_in_dir_name, project.boundaries_in_dir_name,
     130##                  verbose=True, mint=4000, maxt=80000, zscale=1)
     131##
Note: See TracChangeset for help on using the changeset viewer.