Changeset 3944 for anuga_core/documentation/user_manual/examples/project.py

Timestamp:

Nov 8, 2006, 5:16:44 PM (17 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:

• anuga_core/documentation/user_manual/examples/project.py (modified) (1 diff)

anuga_core/documentation/user_manual/examples/project.py

@@ -45 +45 @@
 # demo poly 
 j0 = [385000, 6280000]
-j1 = [360000, 6272500]
-j2 = [335000, 6272500]
+j1 = [360000, 6273000]
+j2 = [335000, 6273000]
 j3 = [330000, 6265000]
 j31 = [325000, 6260000]
 j4 = [316000, 6260000]
-j5 = [316000, 6247000]
-j6 = [350000, 6247000]
+j5 = [316000, 6246750]
+j6 = [350000, 6246750]
 j7 = [385000, 6238000]
 
 demopoly = [j0, j1, j2, j3, j31, j4, j5, j6, j7]
 
-from utilities.polygon import read_polygon
-polygonptsfile4 = polygondir + 'poly1'
-polygonptsfile0 = polygondir + 'poly2'
-polygonptsfile1 = polygondir + 'poly3'
-polygonptsfile2 = polygondir + 'poly4'
-polygonptsfile3 = polygondir + 'poly5'
-northern_polygon = read_polygon(polygonptsfile0 + '.csv')
-manly_polygon = read_polygon(polygonptsfile1 + '.csv')
-harbour_polygon = read_polygon(polygonptsfile2 + '.csv')
-southern_polygon = read_polygon(polygonptsfile3 + '.csv')
-top_polygon = read_polygon(polygonptsfile4 + '.csv')
+from anuga.utilities.polygon import read_polygon, plot_polygons
+#polygonptsfile4 = polygondir + 'poly1'
+#polygonptsfile0 = polygondir + 'poly2'
+#polygonptsfile1 = polygondir + 'poly3'
+#polygonptsfile2 = polygondir + 'poly4'
+#polygonptsfile3 = polygondir + 'poly5'
+#northern_polygon = read_polygon(polygonptsfile0 + '.csv')
+#manly_polygon = read_polygon(polygonptsfile1 + '.csv')
+#harbour_polygon = read_polygon(polygonptsfile2 + '.csv')
+#southern_polygon = read_polygon(polygonptsfile3 + '.csv')
+#top_polygon = read_polygon(polygonptsfile4 + '.csv')
+coastal_polygon = read_polygon(polygondir+'coastal'+'.csv')
+shallow_polygon = read_polygon(polygondir+'shallow'+'.csv')
+
+#plot_polygons([demopoly,northern_polygon,manly_polygon,harbour_polygon,southern_polygon,top_polygon],'model_setup',verbose=False)
+plot_polygons([demopoly,coastal_polygon,shallow_polygon],'new_model_setup',verbose=False)
 
 slump_origin = [372500.0, 6255000.0] #Absolute UTM