Changeset 4311 for anuga_work/production/dampier_2006/run_dampier_simple.py

Timestamp:

Mar 22, 2007, 8:49:47 AM (17 years ago)

Author:

nick

Message:

updates to dampier

File:

• anuga_work/production/dampier_2006/run_dampier_simple.py (modified) (10 diffs)

anuga_work/production/dampier_2006/run_dampier_simple.py

@@ -36 +36 @@
 from anuga.caching import myhash
 # Application specific imports
-#import project                 # Definition of file names and polygons
+import project                 # Definition of file names and polygons
 
 #------------------------------------------------------------------------------
@@ -43 +43 @@
 #------------------------------------------------------------------------------
 
-#start_screen_catcher(project.output_run_time_dir, myid, numprocs)
+start_screen_catcher(project.output_run_time_dir, myid, numprocs)
 print "Processor Name:",get_processor_name()
 
@@ -52 +52 @@
 #boundaries_dir_name = project.boundaries_dir_name
 
-#tide = project.tide
+tide = project.tide
 
 # creates copy of code in output dir
-'''
 if myid == 0:
     copy_code_files(project.output_run_time_dir,__file__, 
                  dirname(project.__file__)+sep+ project.__name__+'.py' )
 barrier()
-'''
-#print 'USER: ', project.user
-#print 'min triangles', project.trigs_min,
+
+print 'USER: ', project.user
+print 'min triangles', project.trigs_min,
 print 'Note: This is generally about 20% less than the final amount'
 
@@ -117 +116 @@
                              verbose=True)
 '''                            
-all = [[469000,7760000],[470000,7758000],[468000,7758000]]
-#all = [[470000,7760000],[469000,7758000],[468000,7760000]]
+#all = [[469000,7760000],[470000,7758000],[468000,7758000]]
+all = [[470000,7760000],[469000,7758000],[468000,7760000]]
 create_mesh_from_regions(all,
-#                             boundary_tags={'ocean': [0,1,2]},
-                             boundary_tags={'ocean': [ 2],'side': [0, 1]},
-#                             boundary_tags={'ocean': [ 2],'ref_side': [0, 1]},
-                             maximum_triangle_area=100000,
-                             filename=meshes_dir_name,
-                             use_cache=True,
-                             verbose=True)
-import sys
-sys.exit()
+                             boundary_tags={'ocean': [0,1,2]},
+#                             boundary_tags={'ocean': [ 2],'side': [0, 1]},
+                             maximum_triangle_area=50000,
+                             filename=meshes_dir_name,
+                             use_cache=False,
+                             verbose=True)
+
 #-------------------------------------------------------------------------
 # Setup computational domain
@@ -134 +131 @@
 print 'Setup computational domain'
 domain = Domain( meshes_dir_name, verbose=True)
-
 '''
 from caching import cache
@@ -199 +195 @@
 print 'Start Set quantity'
 
-domain.set_quantity('elevation', 
+#domain.set_quantity('elevation', 
 #                    filename = project.combined_dir_name + '.pts',
 # MUST USE TXT FILES FOR CACHING TO WORK!
-                    filename = project.combined_smallest_dir_name + '.txt',
-                    use_cache = True,
-                    verbose = True,
-                    alpha = 0.1)
+#                    filename = project.combined_smallest_dir_name + '.txt',
+#                    use_cache = True,
+#                    verbose = True,
+#                    alpha = 0.1)
+domain.set_quantity('elevation', -42.3)
+
 print 'Finished Set quantity'
 
@@ -252 +250 @@
 Bd = Dirichlet_boundary([tide,0,0])
 
-Bo = Dirichlet_boundary([tide+4.0,0,0])
-
 print'set_boundary'
 ##domain.set_boundary({'back': Br,
@@ -259 +255 @@
 ##                     'ocean': Bf})
 domain.set_boundary({#'back': Br,
-                     'side': Br,
-                     'ref_side': Br,
-                    'ocean': Bd
+                     'side': Bd,
+                    'ocean': Bf
                     })
 print'finish set boundary'
@@ -271 +266 @@
 t0 = time.time()
 
-for t in domain.evolve(yieldstep = 10, finaltime = 3400):
+for t in domain.evolve(yieldstep = 60, finaltime = 3400):
     domain.write_time()
     domain.write_boundary_statistics(tags = 'ocean')
-
-    if allclose(t, 120):
-        domain.set_boundary({'side': Br, 'ocean': Bo})
-
-    if allclose(t, 1020):
-        domain.set_boundary({'side': Br, 'ocean': Bd})
 
 #for t in domain.evolve(yieldstep = 120, finaltime = 9000):
@@ -305 +294 @@
 
 
-
-
-
-
-
-
-
-
-