Changeset 5095

Timestamp:

Mar 3, 2008, 9:30:31 AM (17 years ago)

Author:

duncan

Message:

Running different scenarios

Location:

anuga_work/development/near_shore_PMD

Files:

• create_mesh.py (modified) (5 diffs)
• run_beach.py (modified) (8 diffs)

anuga_work/development/near_shore_PMD/create_mesh.py

@@ -8 +8 @@
 
 
-def generate(mesh_filename, maximum_triangle_area=300000):
+def generate(mesh_filename, maximum_triangle_area=300000, thinner=False):
     """
     """
@@ -14 +14 @@
     xleft   = -100  # Beach
     xright  = 10000  # wave generated here
-    ytop    = 15000
-    ybottom = -15000
+    if thinner:
+        ytop    = 1000
+        ybottom = -1000
+    else:
+        ytop    = 2000
+        ybottom = -2000
+        
 
     point_sw = [xleft, ybottom]
@@ -29 +34 @@
     # local refinement
     xleft   = -80
-    xright  = 500
-    ytop    = 1000
-    ybottom = -1000
+    xright  = 2000
+    if thinner:
+        ytop    = 900
+        ybottom = -900
+    else:
+        ytop    = 1000
+        ybottom = -1000
 
     point_sw = [xleft, ybottom]
@@ -43 +52 @@
                         point_sw]
 
-    interior_regions = [[polygon_2, maximum_triangle_area/10.]]
+    interior_regions = [[polygon_2, maximum_triangle_area/2.]]
     
     m = create_mesh_from_regions(bounding_polygon,
@@ -57 +66 @@
 #-------------------------------------------------------------
 if __name__ == "__main__":
-    generate("aa.tsh")
+    generate("aa.tsh", thinner=True)

anuga_work/development/near_shore_PMD/run_beach.py

@@ -23 +23 @@
 from os import path, sep
 from os.path import dirname  #, basename
-
+from math import sin, pi
 
 # Related major packages
@@ -31 +31 @@
                             Transmissive_Momentum_Set_Stage_boundary
 from anuga.fit_interpolate.interpolate import interpolate_sww2csv
-from anuga.abstract_2d_finite_volumes.util import start_screen_catcher, \
-     copy_code_files, file_function
+from anuga.abstract_2d_finite_volumes.util import copy_code_files, \
+     file_function
+
+from anuga.shallow_water.data_manager import start_screen_catcher
+
 from anuga.abstract_2d_finite_volumes.generic_boundary_conditions\
      import File_boundary_time
@@ -59 +62 @@
             #>>> 0.01888*(1000**0.8)
             #4.7424415826900885
+
+            #>>> -0.01888*10000**0.8
+            #-29.922783473665834
     return z 
 
@@ -65 +71 @@
     if is_trial_run is True:
         outputdir_name += '_test'
-        yieldstep = 0.1
-        finaltime = 15.
-        maximum_triangle_area=300000
+        yieldstep = 10
+        finaltime = 1200.
+        maximum_triangle_area=3000
+        thinner=True
     else:
-        yieldstep = 0.02
-        finaltime = 15.1
-        maximum_triangle_area=30000
+        yieldstep = 5.
+        finaltime = 1200
+        maximum_triangle_area=50
+        thinner=True
         
         
     pro_instance = project.Project(['data','near_shore_PMD'],
-                                   outputdir_name=outputdir_name,
-                                   home='.')
+                                   outputdir_name=outputdir_name)
     print "The output dir is", pro_instance.outputdir
     copy_code_files(pro_instance.outputdir,__file__,
                     dirname(project.__file__) \
                     + sep + project.__name__+'.py')
-    copy (pro_instance.codedir + 'run_dam.py',
-          pro_instance.outputdir + 'run_dam.py')
+    copy (pro_instance.codedir + 'run_beach.py',
+          pro_instance.outputdir + 'run_beach.py')
     copy (pro_instance.codedir + 'create_mesh.py',
           pro_instance.outputdir + 'create_mesh.py')
@@ -95 +102 @@
     # creates copy of code in output dir
     if is_trial_run is False:
-        start_screen_catcher(pro_instance.outputdir, rank, pypar.size())
+        start_screen_catcher(pro_instance.outputdir) #, rank, pypar.size())
 
     print 'USER:    ', pro_instance.user
@@ -105 +112 @@
     #gate_position = 12.0
     create_mesh.generate(mesh_filename,
-                         maximum_triangle_area=maximum_triangle_area)
+                         maximum_triangle_area=maximum_triangle_area,
+                         thinner=thinner)
 
     head,tail = path.split(mesh_filename)
@@ -142 +150 @@
     #Bts = Transmissive_Momentum_Set_Stage_boundary(domain, function)
     Br = Reflective_boundary(domain)
-    Bd = Dirichlet_boundary([1000.,0,0]) 
-    domain.set_boundary( {'wall': Br, 'wave': Bd} )
+    Bd = Dirichlet_boundary([10.,0,0])
+    Bwp = Transmissive_Momentum_Set_Stage_boundary(domain,
+                   lambda t: [(1.2*sin(2*t*pi/10)), 0.0, 0.0])
+    Bws = Transmissive_Momentum_Set_Stage_boundary(domain,
+                   lambda t: [1.2*sin(2*t*pi/10)+1.0*sin(2*t*pi/9.5),
+                              0.0, 0.0])
+    domain.set_boundary( {'wall': Br, 'wave': Bws} )
 
     #-------------------------------------------------------------------------
@@ -180 +193 @@
 #-------------------------------------------------------------
 if __name__ == "__main__":
-    main( is_trial_run = True,
-         outputdir_name='Hinwood_low_stage_low_velocity_draft')
+    main( is_trial_run = False,
+         outputdir_name='Primary_wave_thin_50')