Changeset 5653

Timestamp:

Aug 13, 2008, 5:35:07 PM (16 years ago)

Author:

sexton

Message:

updates for gong study

Location:

anuga_work/development/gong_2008

Files:

• export_results.py (added)
• get_timeseries.py (added)
• run_gong_slide.py (modified) (8 diffs)

anuga_work/development/gong_2008/run_gong_slide.py

@@ -37 +37 @@
 #-------------------------------------------------------------------------------
 
-# creates copy of code in output dir 
-copy_code_files(project_slide.outputtimedir,__file__,dirname(project_slide.__file__)+sep+ project_slide.__name__+'.py' )
-myid = 0
-numprocs = 1
-start_screen_catcher(project_slide.outputtimedir, myid, numprocs)
-
-print 'USER:    ', project_slide.user
 
 #-------------------------------------------------------------------------------
@@ -50 +43 @@
 # Convert ASC 2 DEM 2 PTS using source data and store result in source data
 #-------------------------------------------------------------------------------
-
+"""
 # filenames
 on_offshore10_dem_name = project_slide.on_offshore10_dem_name
 nsw_dem_name = project_slide.nsw_dem_name
-meshname = project_slide.meshname+'.msh'
+
 
 # creates DEM from asc data
@@ -88 +81 @@
 G.export_points_file(project_slide.combined_dem_name + '.pts')
 #G.export_points_file(project_slide.combined_dem_name + '.xya')
-
+"""
 #----------------------------------------------------------------------------
 # Create the triangular mesh based on overall clipping polygon with a tagged
@@ -96 +89 @@
 
 from anuga.pmesh.mesh_interface import create_mesh_from_regions
-remainder_res = 100000
-local_res = 25000
-gong_res = 500
+meshname = project_slide.meshname+'.msh'
+remainder_res = 100000*100
+local_res = 25000*1000
+gong_res = 500*1000
 interior_regions = [[project_slide.poly_local, local_res],
                     [project_slide.poly_gong, gong_res],
@@ -171 +165 @@
 thickness = t0
 alpha = a0
-gamma0 = 1.8
+gamma0 = 1.85
 gamma = gamma0
 m0 = 1.
@@ -190 +184 @@
 
 # scaling for Double Gaussian function
-scale0 = 100.
+scale0 = 100. # Bridgette's fiddle 
 scale = scale0 + 10.
 
-for i in range(10):
+for i in range(1):
     scale = scale - 10.   
-    mydir = 'testing_' + str(int(scale))
+    mydir = project_slide.outputdir+'testing_' + str(int(scale))
     print 'dir_comment', mydir
     start_screen_catcher(mydir, 0, 1)
+    # creates copy of code in output dir 
+    copy_code_files(mydir,__file__,dirname(project_slide.__file__)+sep+ project_slide.__name__+'.py' )
+    start_screen_catcher(mydir, 0, 1)
+
+    print 'USER:    ', project_slide.user
+
     tsunami_source = slide_tsunami(length=length,
                                    width=width,
@@ -214 +214 @@
                                    kappad=kappad,
                                    scale=scale,
-                                   domain=domain)
+                                   domain=domain,
+                                   verbose=True)
 
     print 'hello', scale, tsunami_source.wavelength, tsunami_source.a3D
-
+    
     #-------------------------------------------------------------------------------                                 
     # Setup initial conditions
@@ -242 +243 @@
 
     print 'finished'
+