Changeset 4689

Timestamp:

Aug 28, 2007, 6:00:20 PM (17 years ago)

Author:

ole

Message:

Steep point back log

Location:

anuga_work/production/shark_bay_2007

Files:

• project.py (modified) (3 diffs)
• run_shark_bay.py (modified) (1 diff)
• run_shark_bay_frequency_sweep.py (modified) (2 diffs)

anuga_work/production/shark_bay_2007/project.py

@@ -37 +37 @@
 # For frequency response study
 amplitude = 0.5
-period = 6
+period = 3000
 
 #momentum_scale=50 # Experiment
@@ -46 +46 @@
 friction=0.01
 finaltime=40000
-starttime=6000 # Action starts around 9000
+starttime=8000 # Action starts around 9000
 setup='final'
 #setup='trial'
@@ -68 +68 @@
     print'final'
     res_factor=1.0
-    time_thinning=4
-    yieldstep=60
+    time_thinning=1
+    yieldstep=10
 
 
 
 #dir_comment='_'+setup+'_'+str(tide)+'_'+str(source)+'_'+\
-#             str(user)+'_'+boundary_event+'_A'+str(amplitude)+'_T'+str(period)
+#             str(user)+'_'+boundary_event+'_X'+str(momentum_scale)+'_fixedbathy'
 
 dir_comment='_'+setup+'_'+str(tide)+'_'+str(source)+'_'+\
-             str(user)+'_'+boundary_event+'_X'+str(momentum_scale)+'_fixedbathy'
+             str(user)+'_'+boundary_event
+
+
+if boundary_event == 'experimental':
+    dir_comment += '_A'+str(amplitude)+'_T'+str(period)
+    
 
 # elevation data filenames

anuga_work/production/shark_bay_2007/run_shark_bay.py

@@ -222 +222 @@
     
 
-    for t in domain.evolve(yieldstep = kwargs['yieldstep'], finaltime = 21600
+    for t in domain.evolve(yieldstep = kwargs['yieldstep'], finaltime = 18000
                        ,skip_initial_step = True): 
         domain.write_time()

anuga_work/production/shark_bay_2007/run_shark_bay_frequency_sweep.py

@@ -106 +106 @@
                                  interior_regions=project.interior_regions,
                                  filename=project.mesh_name+'.msh',
+                                 fail_if_polygons_outside=False,
                                  use_cache=False,
                                  verbose=True)
@@ -215 +216 @@
 
     for t in domain.evolve(yieldstep = 5,
-                           finaltime = kwargs['finaltime']): 
+                           finaltime = 10000): 
         domain.write_time()
         domain.write_boundary_statistics(tags = 'ocean')