Changeset 2732

Timestamp:

Apr 20, 2006, 3:46:19 PM (19 years ago)

Author:

sexton

Message:

updates on the black screen of death

Location:

production/sydney_2006

Files:

• The Black Screen of Death.doc (modified) (previous)
• run_sydney_smf_polyingo.py (modified) (8 diffs)

production/sydney_2006/run_sydney_smf_polyingo.py

@@ -39 +39 @@
 # Fine pts file to be clipped to area of interest
 #-------------------------------------------------------------------------------
-
 # filenames
 coarsedemname = project.coarsedemname
 finedemname = project.finedemname
-meshname = project.meshname4+'.msh'
+
 
 # coarse data
@@ -72 +71 @@
                                  project.combineddemname + '.pts')
 
-
 #-------------------------------------------------------------------------------                                 
 # Create the triangular mesh based on overall clipping polygon with a tagged
@@ -81 +79 @@
 from pmesh.mesh_interface import create_mesh_from_regions
 
-interior_res1 = 5000
-interior_res2 = 315
+meshname = project.meshname4+'.msh'
+
+interior_res1 = 25000
+interior_res2 = 1315
 print 'interior resolution', interior_res1, interior_res2
 
@@ -123 +123 @@
 #interior_regions = []
 #interior_regions = [[project.harbour_polygon_2, interior_res1],
-                    [project.botanybay_polygon_2, interior_res1]]
+#                    [project.botanybay_polygon_2, interior_res1]]
+
+# finer set
+#interior_regions = [[project.newpoly1, interior_res1],
+#                    #[project.parrariver, interior_res1], #remove for second run
+#                    [project.south1, interior_res1],
+#                    [project.finepolymanly, interior_res2],
+#                    [project.finepolyquay, interior_res2]]
 
 #FIXME: Fix caching of this one once the mesh_interface is ready
 from caching import cache
 
+#_ = cache(create_mesh_from_regions,
+#          project.diffpolygonall,
+#          {'boundary_tags': {'bottom': [0], 
+#                             'right1': [1], 'right0': [2],
+#                             'right2': [3], 'top': [4], 'left1': [5],
+#                             'left2': [6], 'left3': [7]},
+#           'maximum_triangle_area': 100000,
+#           'filename': meshname,           
+#           'interior_regions': interior_regions},
+#          verbose = True)
+
 _ = cache(create_mesh_from_regions,
-          project.diffpolygonall,
-          {'boundary_tags': {'bottom': [0], 
-                             'right1': [1], 'right0': [2],
-                             'right2': [3], 'top': [4], 'left1': [5],
-                             'left2': [6], 'left3': [7]},
-           'maximum_triangle_area': 1000000,
+          project.diffpolygonall_test,
+          {'boundary_tags': {'bottom': [0], 'right': [1],
+                             'top': [2], 'left': [3]},
+           'maximum_triangle_area': 100000,
            'filename': meshname,           
            'interior_regions': interior_regions},
@@ -150 +166 @@
 print 'Number of triangles = ', len(domain)
 print 'The extent is ', domain.get_extent()
-
+print domain.statistics()
 
 domain.set_name(project.basename4)
 domain.set_datadir(project.outputdir)
 domain.set_quantities_to_be_stored(['stage', 'xmomentum', 'ymomentum'])
+
+
 
 
@@ -178 +196 @@
 #-------------------------------------------------------------------------------
 
-from math import sin
-def f(x,y):
-    return x*x
-
 G = Geospatial_data(project.test_pts, project.test_elev) #points, attributes
-domain.set_quantity('stage', tsunami_source)
+domain.set_quantity('stage', 0) #tsunami_source)
 domain.set_quantity('friction', 0.0) # supplied by Benfield, initial value 0.03
-domain.set_quantity('elevation', G,
+domain.set_quantity('elevation', alpha = 10.0,#G, alpha = 550.0, 
                     #filename = project.combineddemname + '.pts',
-                    #filename = project.coarsedemname + '.pts',
-                    use_cache = True,
+                    filename = project.coarsedemname + '.pts',
+                    use_cache = False,
                     verbose = True)
 
-
+from pyvolution.least_squares import fit_to_mesh_file, DEFAULT_ALPHA 
+
+#Add elevation data to the mesh
+#fit_to_mesh_file(mesh_file, point_file, mesh_output_file, DEFAULT_ALPHA,
+#                 verbose=True, expand_search=True, precrop=True)
+#pointname = project.coarsedemname + '.pts'
+#mesh_elevname = project.meshelevname
+#cache(fit_to_mesh_file,(meshname,
+#                 pointname,
+#                 mesh_elevname,
+#                 DEFAULT_ALPHA),
+#      {'verbose': True,
+#       'expand_search': True,
+#       'precrop': True}
+#      ,dependencies = [meshname, pointname]
+#      #,evaluate = True      
+#      ,verbose = False
+#      )
 
 #-------------------------------------------------------------------------------                                 
@@ -200 +231 @@
 
 Br = Reflective_boundary(domain)
-Bw = Time_boundary(domain=domain,
-                   f=lambda t: [(6<t<66)*6, 0, 0])
-domain.set_boundary( {'bottom': Br, 'right1': Br, 'right0': Br,
-                      'right2': Br, 'top': Br, 'left1': Br,
-                      'left2': Br, 'left3': Br} )
+
+#domain.set_boundary( {'bottom': Br, 'right1': Br, 'right0': Br,
+#                      'right2': Br, 'top': Br, 'left1': Br,
+#                      'left2': Br, 'left3': Br} )
+domain.set_boundary( {'bottom': Br, 'right': Br, 'left': Br, 'top': Br} )
 #domain.set_boundary( {'bottom': Br, 'right1': Br, 'right0': Br,
 #                      'right2': Br, 'top': Br, 'left1': Br,
@@ -220 +251 @@
     
 # original
-for t in domain.evolve(yieldstep = 60, finaltime = 10000):
+for t in domain.evolve(yieldstep = 1, finaltime = 10):
     domain.write_time()
     domain.write_boundary_statistics(tags = 'bottom')
-    #stagestep = domain.get_quantity('stage')
-    stagestep = domain.quantities['stage']
-    elevstep = domain.quantities['elevation']
-    depthstep = stagestep - elevstep
-    s = '%.2f, %.2f\n' %(t, depthstep.get_integral())
+    stagestep = domain.get_quantity('stage') 
+    s = '%.2f, %.2f\n' %(t, stagestep.get_integral())
     fid.write(s)
 
 print 'That took %.2f seconds' %(time.time()-t0)
 
-import sys; sys.exit()
-
-for t in domain.evolve(yieldstep = 60, finaltime = 1000):
-    domain.write_time()
-    domain.write_boundary_statistics(tags = 'ocean') #quantities = 'stage')
-
-for t in domain.evolve(yieldstep = 10, finaltime = 10000,
-                       skip_initial_step = True):
-    domain.write_time()
-    domain.write_boundary_statistics(tags = 'ocean') #quantities = 'stage')
-    
-for t in domain.evolve(yieldstep = 60, finaltime = 14400,
-                       skip_initial_step = True):
-    domain.write_time()
-    domain.write_boundary_statistics(tags = 'ocean') #quantities = 'stage'
-
-
-
+
+