Changeset 4072 for anuga_work/production/sydney_2006

Timestamp:

Dec 11, 2006, 6:45:38 PM (18 years ago)

Author:

sexton

Message:

slide modelling updates (based on inundation extents found)

Location:

anuga_work/production/sydney_2006

Files:

• export_results.py (modified) (2 diffs)
• project_slide.py (modified) (2 diffs)
• run_sydney_slide.py (modified) (1 diff)

anuga_work/production/sydney_2006/export_results.py

@@ -1 @@
-import project, os
+import project_slide, os
+import sys
 
-from anuga.pyvolution.data_manager import sww2dem
-from anuga.pyvolution.ermapper_grids import read_ermapper_grid
+from anuga.shallow_water.data_manager import sww2dem
+from os import sep
 
-name = project.outputname
+time_dir = '20061206_061014' 
+directory = project_slide.outputdir
+name = directory + time_dir + sep + project_slide.basename
 
-#print 'Which variable do you want to export?'
-#which_var = int(raw_input('Stage = 0, Absolute Momentum = 1, Depth = 2, Speed = 3  '  ))
+is_parallel = False
+if is_parallel == True: nodes = 4
+print 'output dir:', name
+
 which_var = 2
-
 if which_var == 0:  # Stage
     outname = name + '_stage'
@@ -15 +19 @@
 
 if which_var == 1:  # Absolute Momentum
-    outname = name + '_momentum'
-    quantityname = '(xmomentum**2 + ymomentum**2)**0.5'  #Absolute momentum
+    outname = name + '_momentum_i1'
+    quantityname = '(xmomentum**2 + ymomentum**2)**0.5'  
 
 if which_var == 2:  # Depth
     outname = name + '_depth'
-    quantityname = 'stage-elevation'  #Depth
+    quantityname = 'stage-elevation'  
 
 if which_var == 3:  # Speed
-    outname = name + '_speed'
+    outname = name + '_speed_i0'
     #quantityname = '((xmomentum/(stage-elevation))**2 + (ymomentum/(stage-elevation))**2)**0.5'  #Speed
-    #quantityname = 'xmomentum/(stage-elevation)'  #Speed
     quantityname = '(xmomentum**2 + ymomentum**2)**0.5/(stage-elevation+1.e-30)'  #Speed
 
-    
-sww2dem(name, basename_out = outname,
-            quantity = quantityname,
-            cellsize = 25,       # Trevor would like this at 25
-            # define region for viz purposes
-            easting_min = project.eminviz,
-            easting_max = project.emaxviz,
-            northing_min = project.nminviz,
-            northing_max = project.nmaxviz,        
-            reduction = max, #this is because we want max quantityname
-            verbose = True,
-            format = 'asc')
+if which_var == 4:  # Elevation
+    outname = name + '_elevation'
+    quantityname = 'elevation'  #Elevation
 
-#Check
-
-#data = read_ermapper_grid(name)
-#print 'Values from %s are in [%f, %f]' %(name, min(data.flat), max(data.flat))
+if is_parallel == True:
+    for i in range(0,nodes):
+        namei = name + '_P%d_%d' %(i,nodes)
+        outnamei = outname + '_P%d_%d' %(i,nodes)
+        print 'start sww2dem for sww file %d' %(i)
+        sww2dem(namei, basename_out = outnamei,
+                    quantity = quantityname,
+                    timestep = 1,
+                    cellsize = 100,      
+                    #easting_min = project.e_min_area,
+                    #easting_max = project.e_max_area,
+                    #northing_min = project.n_min_area,
+                    #northing_max = project.n_max_area,        
+                    reduction = max, 
+                    verbose = True,
+                    format = 'asc')
+else:
+    print 'start sww2dem'
+    sww2dem(name, basename_out = outname,
+                quantity = quantityname,
+                cellsize = 100,      
+                easting_min = project_slide.eastingmin,
+                easting_max = project_slide.eastingmax,
+                northing_min = project_slide.northingmin,
+                northing_max = project_slide.northingmax,        
+                reduction = max, 
+                verbose = True,
+                format = 'asc')

anuga_work/production/sydney_2006/project_slide.py

@@ -93 +93 @@
 
 # sydney digitized polygons
-poly_coast = read_polygon(polygondir+'coast.csv')
+poly_coast = read_polygon(polygondir+'coast2.csv')
 poly_syd1 = read_polygon(polygondir+'local.csv')
 
@@ -129 +129 @@
 
 # exporting asc grid
-#eastingmin = 
-#eastingmax = 
-#northingmin = 
-#northingmax = 
+eastingmin = 337500
+eastingmax = 344060
+northingmin = 6248890
+northingmax = 6274550
 
 ###################################################################

anuga_work/production/sydney_2006/run_sydney_slide.py

@@ -113 +113 @@
 
 from anuga.pmesh.mesh_interface import create_mesh_from_regions
-remainder_res = 150000.
+remainder_res = 250000.
 local_res = 50000.
 coast_res = 500.