Changeset 3661

Timestamp:

Sep 25, 2006, 2:26:13 PM (18 years ago)

Author:

sexton

Message:

clean up

Location:

anuga_work/production/hobart_2006

Files:

• export_results.py (modified) (3 diffs)
• project.py (modified) (3 diffs)
• run_hobart.py (modified) (5 diffs)

anuga_work/production/hobart_2006/export_results.py

@@ -3 +3 @@
 
 from anuga.shallow_water.data_manager import sww2dem
+#from anuga.pyvolution.ermapper_grids import read_ermapper_grid
+from anuga.abstract_2d_finite_volumes.util import Screen_Catcher
 from os import sep
 
-time_dir = '20060921_083846' #MSL
+time_dir = '20060922_024842' #MSL
 directory = project.outputdir
-name = directory + time_dir + sep + 'source'
+name = directory + time_dir +sep + 'source'
 
+print 'output dir:', name
 #print 'Which variable do you want to export?'
 #which_var = int(raw_input('Stage = 0, Absolute Momentum = 1, Depth = 2, Speed = 3  '  ))
-which_var = 4 # for Hamish/Alex to compare
-#which_var = 2 # for Ingo/Neil to make maps
-
+which_var = 4
+#sys.stderr.write(sys.stdout.data)
 if which_var == 0:  # Stage
     outname = name + '_stage'
@@ -33 +35 @@
 
 if which_var == 4:  # Elevation
-    outname = name + '_elevation'
+    outname = name + '_elevation_node3'
     quantityname = 'elevation'  #Elevation
 
+print 'start sww2dem'
+#sys.stderr.write(sys.stdout.data)
 sww2dem(name, basename_out = outname,
             quantity = quantityname,
-            cellsize = 25,       
+            cellsize = 25,      # Trevor would like this at 25
             # define region for viz purposes
             easting_min = project.e_min_area,
@@ -47 +51 @@
             verbose = True,
             format = 'asc')
+
+#sys.stderr.write(sys.stdout.data)
+
+#Check
+
+#data = read_ermapper_grid(name)
+#print 'Values from %s are in [%f, %f]' %(name, min(data.flat), max(data.flat))

anuga_work/production/hobart_2006/project.py

@@ -15 +15 @@
 #onshore_name = 'hob3_topo' # original
 onshore_name = 'hob5_topo' # 12.5m grid and clipped to 100m elevation or 3000m from coast
-onshore_name_25 = 'hob5_topo_25m' # 25m grid and clipped to 100m elevation or 3000m from coast
+#onshore_name_25 = 'hob5_topo_25m' # 25m grid and clipped to 100m elevation or 3000m from coast
+onshore_name_25 = 'hob6_topo_25m' # 25m grid NOT clipped
 offshore_name_tas1 = 'derwent_2m'
 offshore_name_tas2 = 'derwent_5m'
@@ -121 +122 @@
 
 # region to export for inundation map
-e_min_area = 500000
+e_min_area = 490000
 e_max_area = 580000
-n_min_area = 5200000
+n_min_area = 5160000
 n_max_area = 5275000
 
@@ -142 +143 @@
 polyAll, zone = convert_points_from_latlon_to_utm([d0, d1, d2, d3])
 refzone = zone
-NWcorner = polyAll[0]
-SEcorner = polyAll[2]
+polyAll = [[500000, 5200000],[520000, 5200000],[520000, 5400000],[500000,5400000]]
 
 # region to export for Alex to make bathymetry map
-e_min_area = 480000#NWcorner[0]
-e_max_area = 600000#SEcorner[0]
-n_min_area = 5156000#NWcorner[1]
-n_max_area = 5285000#SEcorner[1]
+e_min_area = 500000 #490000
+e_max_area = 560000 #580000
+n_min_area = 5240000#5160000
+n_max_area = 5260000#5270000
 
 #Interior region - Hobart city area + Glenorchy, Kingston

anuga_work/production/hobart_2006/run_hobart.py

@@ -88 +88 @@
                                  project.onshore_dem_name_25 + '.pts',
                                  project.all_onshore_dem_name + '.pts')
-#print 'local offshore data sets'
-#G1 = Geospatial_data(file_name = project.offshore_dem_name_local1 + '.xya')
-#G2 = Geospatial_data(file_name = project.offshore_dem_name_local2 + '.xya')
-#G3 = Geospatial_data(file_name = project.offshore_dem_name_local3 + '.xya')
-#G4 = Geospatial_data(file_name = project.offshore_dem_name_local4 + '.xya')
-#G11 = G1 + G2 + G3 + G4
-#print 'aho offshore data sets'
-#G5 = Geospatial_data(file_name = project.offshore_dem_name_aho1 + '.xya')
-#G6 = Geospatial_data(file_name = project.offshore_dem_name_aho2 + '.xya')
-#G7 = Geospatial_data(file_name = project.offshore_dem_name_aho3 + '.xya')
-#G8 = Geospatial_data(file_name = project.offshore_dem_name_aho4 + '.xya')
-#G9 = Geospatial_data(file_name = project.offshore_dem_name_aho5 + '.xya')
-#G10 = Geospatial_data(file_name = project.offshore_dem_name_aho6 + '.xya')
-#G11 = Geospatial_data(file_name = project.offshore_dem_name_aho7 + '.xya')
-#G12 = Geospatial_data(file_name = project.offshore_dem_name_aho8 + '.xya')
-#G13 = Geospatial_data(file_name = project.offshore_dem_name_aho9 + '.xya')
-#G14 = Geospatial_data(file_name = project.offshore_dem_name_aho10 + '.xya')
-#G15 = Geospatial_data(file_name = project.offshore_dem_name_aho11 + '.xya')
-#G16 = Geospatial_data(file_name = project.offshore_dem_name_aho12 + '.xya')
-#G17 = Geospatial_data(file_name = project.offshore_dem_name_aho13 + '.xya')
-#G18 = Geospatial_data(file_name = project.offshore_dem_name_aho14 + '.xya')
-#G19 = Geospatial_data(file_name = project.offshore_dem_name_aho15 + '.xya')
-#G20 = Geospatial_data(file_name = project.offshore_dem_name_aho16 + '.xya')
-#G12 = G5+G6+G7+G8+G9+G10+G11+G12+G13+G14+G15+G16+G17+G18+G19+G20
-#print 'local onshore dem'
-#G21 = Geospatial_data(file_name = project.onshore_dem_name + '.pts')
-print 'coast'
-#G22 = Geospatial_data(file_name = project.coast_dem_name + '.xya')
-#print 'adding data sets'
-#G = G11 + G12 + G21 + G22
-#G = G1+G2+G3+G4+G5+G6+G7+G8+G9+G10+G11+G12+G13+G14+G15+G16+G17+G18+G19+G20+G21+G22
-#G = G5+G6+G7+G8+G9+G10+G11+G12+G13+G14+G15+G16+G17+G18+G19+G20+G21+G22
-#G = G21+G22+G5+G6+G7+G8+G9+G10+G11+G12+G13+G14+G15+G16+G17+G18+G19+G20
-#G = G21+G22+G5
-#G = Geospatial_data(file_name = project.onshore_dem_name + '.pts') #+ \
-#    Geospatial_data(file_name = project.coast_dem_name + '.xya') #+ \
-#    Geospatial_data(file_name = project.offshore_dem_name_aho1 + '.xya')
+
+print 'adding data sets'
 G = Geospatial_data(file_name = project.offshore_dem_name_local1 + '.xya')+\
     Geospatial_data(file_name = project.offshore_dem_name_local2 + '.xya')+\
@@ -148 +113 @@
 G.export_points_file(project.combined_dem_name + '.pts')
 
-#-------------------------------------------------------------------------------                                 
+#----------------------------------------------------------------------------
 # Create the triangular mesh based on overall clipping polygon with a tagged
 # boundary and interior regions defined in project.py along with
@@ -171 +136 @@
           {'boundary_tags': {'bottom': [0], 'right': [1],
                              'top': [2], 'left': [3]},
-           'maximum_triangle_area': 250000,
-           'filename': meshname,           
-           'interior_regions': interior_regions},
+           'maximum_triangle_area': 100000,
+           'filename': meshname},           
+           #'interior_regions': interior_regions},
           verbose = True, evaluate=True)
 
@@ -211 +176 @@
 
 #-------------------------------------------------------------------------------                                 
-# Setup boundary conditions (all reflective)
+# Setup boundary conditions
 #-------------------------------------------------------------------------------
 print 'start ferret2sww'
@@ -274 +239 @@
 t0 = time.time()
 
-for t in domain.evolve(yieldstep = 240, finaltime = 7200): 
+for t in domain.evolve(yieldstep = 60, finaltime = 120): 
     domain.write_time()
     domain.write_boundary_statistics(tags = 'bottom')     
 
-for t in domain.evolve(yieldstep = 120, finaltime = 12600
-                       ,skip_initial_step = True): 
-    domain.write_time()
-    domain.write_boundary_statistics(tags = 'bottom')     
+#for t in domain.evolve(yieldstep = 120, finaltime = 12600
+#                       ,skip_initial_step = True): 
+#    domain.write_time()
+#    domain.write_boundary_statistics(tags = 'bottom')     
     
 print 'That took %.2f seconds' %(time.time()-t0)