Ignore:
Timestamp:
Nov 8, 2006, 5:16:44 PM (16 years ago)
Author:
sexton
Message:

(i) update sww2timeseries so can handle gauges which don't fall within sww domain (ii) script for generating report for dampier based on sww files from parallel setup

File:
1 edited

Legend:

Unmodified
Added
Removed
  • anuga_core/documentation/user_manual/examples/project.py

    r3869 r3944  
    4545# demo poly
    4646j0 = [385000, 6280000]
    47 j1 = [360000, 6272500]
    48 j2 = [335000, 6272500]
     47j1 = [360000, 6273000]
     48j2 = [335000, 6273000]
    4949j3 = [330000, 6265000]
    5050j31 = [325000, 6260000]
    5151j4 = [316000, 6260000]
    52 j5 = [316000, 6247000]
    53 j6 = [350000, 6247000]
     52j5 = [316000, 6246750]
     53j6 = [350000, 6246750]
    5454j7 = [385000, 6238000]
    5555
    5656demopoly = [j0, j1, j2, j3, j31, j4, j5, j6, j7]
    5757
    58 from utilities.polygon import read_polygon
    59 polygonptsfile4 = polygondir + 'poly1'
    60 polygonptsfile0 = polygondir + 'poly2'
    61 polygonptsfile1 = polygondir + 'poly3'
    62 polygonptsfile2 = polygondir + 'poly4'
    63 polygonptsfile3 = polygondir + 'poly5'
    64 northern_polygon = read_polygon(polygonptsfile0 + '.csv')
    65 manly_polygon = read_polygon(polygonptsfile1 + '.csv')
    66 harbour_polygon = read_polygon(polygonptsfile2 + '.csv')
    67 southern_polygon = read_polygon(polygonptsfile3 + '.csv')
    68 top_polygon = read_polygon(polygonptsfile4 + '.csv')
     58from anuga.utilities.polygon import read_polygon, plot_polygons
     59#polygonptsfile4 = polygondir + 'poly1'
     60#polygonptsfile0 = polygondir + 'poly2'
     61#polygonptsfile1 = polygondir + 'poly3'
     62#polygonptsfile2 = polygondir + 'poly4'
     63#polygonptsfile3 = polygondir + 'poly5'
     64#northern_polygon = read_polygon(polygonptsfile0 + '.csv')
     65#manly_polygon = read_polygon(polygonptsfile1 + '.csv')
     66#harbour_polygon = read_polygon(polygonptsfile2 + '.csv')
     67#southern_polygon = read_polygon(polygonptsfile3 + '.csv')
     68#top_polygon = read_polygon(polygonptsfile4 + '.csv')
     69coastal_polygon = read_polygon(polygondir+'coastal'+'.csv')
     70shallow_polygon = read_polygon(polygondir+'shallow'+'.csv')
     71
     72#plot_polygons([demopoly,northern_polygon,manly_polygon,harbour_polygon,southern_polygon,top_polygon],'model_setup',verbose=False)
     73plot_polygons([demopoly,coastal_polygon,shallow_polygon],'new_model_setup',verbose=False)
    6974
    7075slump_origin = [372500.0, 6255000.0] #Absolute UTM
Note: See TracChangeset for help on using the changeset viewer.