Changeset 3276 for production/onslow_2006/run_onslow.py

Timestamp:

Jul 4, 2006, 2:51:19 PM (18 years ago)

Author:

sexton

Message:

update run_onslow

File:

• production/onslow_2006/run_onslow.py (modified) (10 diffs)

production/onslow_2006/run_onslow.py

@@ -18 +18 @@
 import os
 import time
+from shutil import copy
+from os import mkdir, access, F_OK
+import sys
 
 # Related major packages
@@ -25 +28 @@
 from pyvolution.combine_pts import combine_rectangular_points_files 
 from pyvolution.pmesh2domain import pmesh_to_domain_instance
-#from geospatial_data import add_points_files
+from geospatial_data import *
+from pyvolution.util import Screen_Catcher
 
 # Application specific imports
 import project                 # Definition of file names and polygons
-from smf import slump_tsunami  # Function for submarine mudslide 
-
-from shutil import copy
-from os import mkdir, access, F_OK
-
-from geospatial_data import *
-import sys
-from pyvolution.util import Screen_Catcher
-
-#-------------------------------------------------------------------------------
-# Preparation of topographic data
-#
-# Convert ASC 2 DEM 2 PTS using source data and store result in source data
-# Do for coarse and fine data
-# Fine pts file to be clipped to area of interest
-#-------------------------------------------------------------------------------
-
-# filenames
-#coarsedemname = project.coarsedemname
-
-onshore_dem_name = project.onshore_dem_name
-islands_dem_name = project.islands_dem_name
-coast_points = project.coast_dem_name
-offshore_points = project.offshore_dem_name
-
-meshname = project.meshname+'.msh'
-
-source_dir = project.boundarydir
+
+#-------------------------------------------------------------------------------
+# Copy scripts to time stamped output directory and capture screen
+# output to file
+#-------------------------------------------------------------------------------
 
 # creates copy of code in output dir if dir doesn't exist
@@ -74 +55 @@
 sys.stderr = Screen_Catcher(screen_error_name)
 
+print 'USER:    ', project.user
+
+#-------------------------------------------------------------------------------
+# Preparation of topographic data
+#
+# Convert ASC 2 DEM 2 PTS using source data and store result in source data
+# Do for coarse and fine data
+# Fine pts file to be clipped to area of interest
+#-------------------------------------------------------------------------------
+
+# filenames
+onshore_dem_name = project.onshore_dem_name
+islands_dem_name = project.islands_dem_name
+coast_points = project.coast_dem_name
+offshore_points = project.offshore_dem_name
+meshname = project.meshname+'.msh'
+source_dir = project.boundarydir
 
 copied_files = False
@@ -80 +78 @@
 files_used = [onshore_dem_name, offshore_points, coast_points,]
 
-if sys.platform != 'win32':    
-    copied_files = True
-    for name in file_list:
-        copy(name, )
-'''    
-#print' most file', project.MOST_dir + project.boundary_basename+'_ha.nc'
-#if access(project.MOST_dir + project.boundary_basename+'_ha.nc',F_OK) == 1 :
-#    print' most file', project.MOST_dir + project.boundary_basename
-
-'''
 # fine data (clipping the points file to smaller area)
 # creates DEM from asc data
@@ -110 +98 @@
 print'create G1'
 G1 = Geospatial_data(file_name = project.offshore_dem_name + '.xya')
-
 print'create G2'
 G2 = Geospatial_data(file_name = project.onshore_dem_name + '.pts')
-
 print'create G3'
 G3 = Geospatial_data(file_name = project.coast_dem_name + '.xya')
-
 print'create G4'
 G4 = Geospatial_data(file_name = project.islands_dem_name + '.pts')
-
 print'add G1+G2+G3+G4'
 G = G1 + G2 + G3 + G4
-
 print'export G'
 G.export_points_file(project.combined_dem_name + '.pts')
-
 
 #-------------------------------------------------------------------------------                                 
@@ -134 +116 @@
 
 from pmesh.mesh_interface import create_mesh_from_regions
-'''
-# original
-interior_res = 5000
-high_res = 1500
-interior_regions = [[project.poly_onslow, high_res],
-                    [project.poly_thevenard, interior_res],
-                    [project.poly_coast, interior_res]]
-'''
+
 #new
-region_res = 25000
-coast_res = 2500
-onslow_res = 500
+region_res = 200000
+coast_res = 25000
+onslow_res = 5000
 interior_regions = [[project.poly_onslow, onslow_res],
                     [project.poly_coast, coast_res],
@@ -181 +156 @@
 domain.set_quantities_to_be_stored(['stage', 'xmomentum', 'ymomentum'])
 
-#-------------------------------------------------------------------------------
-# Set up scenario (tsunami_source is a callable object used with set_quantity)
-#-------------------------------------------------------------------------------
-'''
-tsunami_source = slump_tsunami(length=30000.0,
-                               depth=400.0,
-                               slope=6.0,
-                               thickness=176.0, 
-                               radius=3330,
-                               dphi=0.23,
-                               x0=project.slump_origin[0], 
-                               y0=project.slump_origin[1], 
-                               alpha=0.0, 
-                               domain=domain)
-
-'''
+
 #-------------------------------------------------------------------------------                                 
 # Setup initial conditions
@@ -253 +213 @@
        'inverted_bathymetry': True},
       #evaluate = True,
-      verbose = True, 
-      dependencies = source_dir + project.boundary_basename + '.sww')
+       verbose = True,
+       dependencies = source_dir + project.boundary_basename + '.sww')
 
 
@@ -270 +230 @@
 
 domain.set_boundary( {'top': Bf, 'topleft': Bf,
-                             'topleft1': Bf, 'bottomleft': Bd,
-                             'bottom': Br, 'bottomright': Br, 'topright': Bd} )
+                      'topleft1': Bf, 'bottomleft': Bd,
+                      'bottom': Br, 'bottomright': Br, 'topright': Bd} )
 
 #-------------------------------------------------------------------------------                                 
@@ -279 +239 @@
 t0 = time.time()
 
-for t in domain.evolve(yieldstep = 450, finaltime = 10800): 
-    domain.write_time()
-    domain.write_boundary_statistics(tags = 'top')     
-
-for t in domain.evolve(yieldstep = 1, finaltime = 17760
-                       ,skip_initial_step = True): 
-    domain.write_time()
-    domain.write_boundary_statistics(tags = 'top')     
-
-for t in domain.evolve(yieldstep = 480, finaltime = 36000
+for t in domain.evolve(yieldstep = 240, finaltime = 7200): 
+    domain.write_time()
+    domain.write_boundary_statistics(tags = 'top')     
+
+for t in domain.evolve(yieldstep = 120, finaltime = 12600
+                       ,skip_initial_step = True): 
+    domain.write_time()
+    domain.write_boundary_statistics(tags = 'top')     
+
+for t in domain.evolve(yieldstep = 60, finaltime = 19800
                        ,skip_initial_step = True): 
     domain.write_time()
     domain.write_boundary_statistics(tags = 'top')     
     
-  
+for t in domain.evolve(yieldstep = 120, finaltime = 25200
+                       ,skip_initial_step = True): 
+    domain.write_time()
+    domain.write_boundary_statistics(tags = 'top')     
+
+for t in domain.evolve(yieldstep = 240, finaltime = 36000
+                       ,skip_initial_step = True): 
+    domain.write_time()
+    domain.write_boundary_statistics(tags = 'top')
+    
 print 'That took %.2f seconds' %(time.time()-t0)