Changeset 2902 for production/pt_hedland_2006/run_pt_hedland.py

Timestamp:

May 18, 2006, 10:01:58 AM (19 years ago)

Author:

nick

Message:

update onslow and pt hedland

File:

• production/pt_hedland_2006/run_pt_hedland.py (modified) (14 diffs)

production/pt_hedland_2006/run_pt_hedland.py

@@ -16 +16 @@
 
 # Standard modules
-import os
+import os 
+
+from os import sep
+from os.path import dirname, basename
+
 import time
 
@@ -47 +51 @@
 
 # filenames
+print 'run home', project.home
 onshore_dem_name = project.onshore_dem_name
 offshore_points1 = project.offshore_dem_name1
@@ -54 +59 @@
 
 
+
+
+#import sys; sys.exit()
+
 # creates copy of code in output dir if dir doesn't exist
 if access(project.outputtimedir,F_OK) == 0 :
     mkdir (project.outputtimedir)
-copy (project.codedirname, project.outputtimedir + project.codename)
-copy (project.codedir + 'run_pt_hedland.py', project.outputtimedir + 'run_pt_hedland.py')
-print'output dir', project.outputtimedir
+copy (dirname(project.__file__) +sep+ project.__name__+'.py', project.outputtimedir + project.__name__+'.py')
+copy (__file__, project.outputtimedir + basename(__file__))
+
 
 #normal screen output is stored in 
@@ -67 +76 @@
 #used to catch screen output to file
 sys.stdout = Screen_Catcher(screen_output_name)
-#sys.stderr = Screen_Catcher(screen_output_name)
-#sys.stderr = Screen_Catcher(screen_error_name)
-
-'''
-copied_files = False
-
-# files to be used
-files_used = [onshore_dem_name, offshore_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
-
-
+sys.stderr = Screen_Catcher(screen_error_name)
+
+
+'''
 # fine data (clipping the points file to smaller area)
 # creates DEM from asc data 
@@ -114 +108 @@
 G.export_points_file(project.combined_dem_name + '.pts')
 
-
+'''
 #-------------------------------------------------------------------------------                                 
 # Create the triangular mesh based on overall clipping polygon with a tagged
@@ -123 +117 @@
 from pmesh.mesh_interface import create_mesh_from_regions
 
-region_res = 25000
+region_res = 50000
 coast_res = 2500
-pt_hedland_res = 25000
+pt_hedland_res = 5000
 # derive poly_coast from project.coast_name using alpha_shape
-interior_regions = []#[[project.poly_pt_hedland, pt_hedland_res]]#,
-                    #[project.poly_coast, coast_res],
-                    #[project.poly_region, region_res]]
+interior_regions = [[project.poly_pt_hedland, pt_hedland_res],
+                    [project.poly_region, region_res]]
 
 print 'number of interior regions', len(interior_regions)
@@ -136 +129 @@
 _ = cache(create_mesh_from_regions,
           project.polyAll,
-          {'boundary_tags': {'right': [0], 'bottomright': [1],
-                             'bottomleft': [2], 'left': [3], 'top': [4]},
-           'maximum_triangle_area': 250000,
+#          {'boundary_tags': {'right': [0], 'bottomright': [1],
+#                             'bottomleft': [2], 'left': [3], 'top': [4]},
+          {'boundary_tags': {'top': [0], 'left': [1],
+                             'bottomleft': [2], 'bottomright': [3], 'right': [4]},
+           'maximum_triangle_area': 400000,
            'filename': meshname,           
            'interior_regions': interior_regions},
@@ -147 +142 @@
 # Setup computational domain
 #-------------------------------------------------------------------------------                                 
-
+'''
 domain = pmesh_to_domain_instance(meshname, Domain,
                                   use_cache = False,
                                   verbose = True)
+'''
+
+domain = Domain(meshname, use_cache = False, verbose = True)
 
 print 'Number of triangles = ', len(domain)
@@ -171 +169 @@
 
 domain.set_quantity('elevation', 
-#                    filename = project.onshore_dem_name + '.pts',
                     filename = project.combined_dem_name + '.pts',
-#                    filename = project.offshore_dem_name1 + '.xya',
-#                    filename = project.offshore_dem_name2 + '.xya',
                     use_cache = True,
                     verbose = True,
@@ -181 +176 @@
 
 print 'hi1'
-print 'have sent quantities OK - now exiting'
-import sys; sys.exit()
+#print 'have sent quantities OK - now exiting'
+#import sys; sys.exit()
 
 #-------------------------------------------------------------------------------                                 
@@ -188 +183 @@
 #-------------------------------------------------------------------------------
 print 'start ferret2sww'
+# skipped as results in file SU-AU_clipped is correct for all WA
+
 from pyvolution.data_manager import ferret2sww
 
@@ -199 +196 @@
       (source_dir + project.boundary_basename,
        source_dir + project.boundary_basename), 
-#      (project.MOST_dir + project.boundary_basename,
-#       source_dir + project.boundary_basename), 
       {'verbose': True,
-# note didn't work with the below
-#       'minlat': south - 1,
-#       'maxlat': north + 1,
-#       'minlon': west - 1,
-#       'maxlon': east + 1,
        'minlat': south,
        'maxlat': north,
@@ -228 +218 @@
 Bd = Dirichlet_boundary([tide,0,0])
 
-
 # 7 min square wave starting at 1 min, 6m high
 Bw = Time_boundary(domain = domain,
                    f=lambda t: [(60<t<480)*6, 0, 0])
 
-domain.set_boundary( {'right': Br, 'bottomright': Br,
-                             'bottomleft': Br, 'left': Br, 'top': Bf} )
+#domain.set_boundary( {'right': Br, 'bottomright': Br,
+#                             'bottomleft': Br, 'left': Br, 'top': Bf} )
+domain.set_boundary( {'top': Bf, 'left': Br, 'bottomleft': Br, 'bottomright': Br, 'right': Br} )
                              
 #-------------------------------------------------------------------------------                                 
@@ -242 +232 @@
 t0 = time.time()
 
-for t in domain.evolve(yieldstep = 240, finaltime = 7201): 
-    domain.write_time()
-    domain.write_boundary_statistics(tags = 'top')     
-
-for t in domain.evolve(yieldstep = 120, finaltime = 12601
-                       ,skip_initial_step = True): 
-    domain.write_time()
-    domain.write_boundary_statistics(tags = 'top')     
-
-for t in domain.evolve(yieldstep = 60, finaltime = 19801
+for t in domain.evolve(yieldstep = 240, finaltime = 7200): 
+    domain.write_time()
+    domain.write_boundary_statistics(tags = 'top')     
+
+for t in domain.evolve(yieldstep = 240, finaltime = 12600
+                       ,skip_initial_step = True): 
+    domain.write_time()
+    domain.write_boundary_statistics(tags = 'top')     
+
+for t in domain.evolve(yieldstep = 240, finaltime = 19800
                        ,skip_initial_step = True): 
     domain.write_time()
     domain.write_boundary_statistics(tags = 'top')     
     
-for t in domain.evolve(yieldstep = 120, finaltime = 25201
-                       ,skip_initial_step = True): 
-    domain.write_time()
-    domain.write_boundary_statistics(tags = 'top')     
-
-for t in domain.evolve(yieldstep = 240, finaltime = 36001
+for t in domain.evolve(yieldstep = 240, 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()