Changeset 2624 for production/onslow_2006/run_onslow.py

Timestamp:

Mar 29, 2006, 8:42:34 AM (18 years ago)

Author:

nick

Message:

fixed geospatial_data to not use point_dict in export and updated onslow run to use new geospatial_data methods

File:

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

production/onslow_2006/run_onslow.py

@@ -55 +55 @@
 source_dir = project.boundarydir
 
-# creates copy of code in output dir
+# creates copy of code in output dir if dir doesn't exist
 if access(project.outputdir,F_OK) == 0 :
     mkdir (project.outputdir)
 copy (project.codedirname, project.outputdir + project.codename)
 copy (project.codedir + 'run_onslow.py', project.outputdir + 'run_onslow.py')
-
+print'hello'
+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
 
 '''
@@ -79 +82 @@
 
         '''
+'''
 print'create G1'
 G1 = Geospatial_data(file_name = project.offshore_dem_name + '.xya')
@@ -91 +95 @@
 G.new_export_points_file(project.combined_dem_name + '.pts')
 
-'''
-add_points_files(
-                  project.offshore_dem_name + '.xya',
-                  project.onshore_dem_name + '.pts',
-                  project.combined_dem_name + '.pts')
 '''
 #-------------------------------------------------------------------------------                                 
@@ -106 +105 @@
 
 # original
-interior_res = 5000
+interior_res = 50000
 interior_regions = [[project.poly_onslow, interior_res],
                     [project.poly_thevenard, interior_res],
@@ -119 +118 @@
                              'left': [2], 'bottom': [3],
                              'bottomright': [4], 'topright': [5]},
-           'maximum_triangle_area': 100000,
+           'maximum_triangle_area': 1000000,
            'filename': meshname,           
            'interior_regions': interior_regions},
@@ -167 +166 @@
 domain.set_quantity('friction', 0.0) 
 print 'hi and file',project.combined_dem_name + '.pts'
+
 domain.set_quantity('elevation', 
 #                    0.
@@ -176 +176 @@
                     alpha = 0.1
                     )
+
 print 'hi1'
 
@@ -190 +191 @@
 
 cache(ferret2sww,
-      (source_dir + project.boundary_basename,
+#      (source_dir + project.boundary_basename,
+#       source_dir + project.boundary_basename), 
+      (project.MOST_dir + project.boundary_basename,
        source_dir + project.boundary_basename), 
       {'verbose': True,