Changeset 3145

Timestamp:

Jun 13, 2006, 2:38:59 PM (19 years ago)

Author:

duncan

Message:

minor streamlining

Location:

development/dam_2006

Files:

• create_mesh.py (modified) (2 diffs)
• run_dam.py (modified) (4 diffs)

development/dam_2006/create_mesh.py

@@ -49 +49 @@
                             'wall',
                             'wall',
-                            'wave',
+                            'wall', #back of the dam
                             'wall',
                             'wall',
@@ -60 +60 @@
     dam.setTag("dam")
     #m.generateMesh('pzq28.0za0.01a')
-    m.generateMesh('pzq28.0za0.0001a')
+    #m.generate_mesh('pzq28.0za0.0001)
+    m.generate_mesh(maximum_triangle_area=0.0001)
 
     m.export_mesh_file(project.mesh_filename)

development/dam_2006/run_dam.py

@@ -31 +31 @@
 import create_mesh
 
+
+#-------------------------------------------------------------------------------                                 
+# Setup archiving of simulations
+#-------------------------------------------------------------------------------   
+
 copy (project.codedirname, project.outputtimedir + 'project.py')
 copy (project.codedir + 'run_dam.py', project.outputtimedir + 'run_dam.py')
@@ -40 +45 @@
 screen_error_name = project.outputtimedir + "screen_error.txt"
 
-#used to catch screen output to file
-#sys.stdout = Screen_Catcher(screen_output_name)
-#sys.stderr = Screen_Catcher(screen_error_name)
-
-'''
-copied_files = False
-
-# files to be used
-files_used = [onshore_dem_name, offshore_points,]
-
-if sys.platform != 'win32':    
-    copied_files = True
-    for name in file_list:
-        copy(name, )
-''' 
 #-------------------------------------------------------------------------------                                 
 # Create the triangular mesh
 #-------------------------------------------------------------------------------
-#from pmesh.pmesh import mesh
 
 create_mesh.generate() # this creates the mesh
@@ -96 +85 @@
 Bd = Dirichlet_boundary([tide,0,0])
 
-
-# 7 min square wave starting at 0 min, 6m high
-Bw = Time_boundary(domain = domain,
-                   f=lambda t: [(0<t<480)*1, 0, 0])
-
 #domain.set_boundary( {'wall': Br, 'wave': Bw} )
 domain.set_boundary( {'wall': Br, 'wave': Br} )
@@ -110 +94 @@
 t0 = time.time()
 
-for t in domain.evolve(yieldstep = 0.02, finaltime = 10):
+for t in domain.evolve(yieldstep = 0.1, finaltime = 5):
     domain.write_time()