Changeset 3669 for anuga_work/production/hobart_2006/run_hobart.py

Timestamp:

Sep 26, 2006, 11:30:09 AM (17 years ago)

Author:

sexton

Message:

reflecting PATH change in project scripts

File:

• anuga_work/production/hobart_2006/run_hobart.py (modified) (5 diffs)

anuga_work/production/hobart_2006/run_hobart.py

@@ -69 +69 @@
 
 # creates DEM from asc data - 12.5m
-convert_dem_from_ascii2netcdf(onshore_dem_name, use_cache=True, verbose=True)
+#convert_dem_from_ascii2netcdf(onshore_dem_name, use_cache=True, verbose=True)
 
 #creates pts file for onshore DEM - 12.5
-dem2pts(onshore_dem_name,
-        easting_min=project.eastingmin,
-        easting_max=project.eastingmax,
-        northing_min=project.northingmin,
-        northing_max= project.northingmax,
-        use_cache=True, verbose=True)
+#dem2pts(onshore_dem_name,
+#        easting_min=project.eastingmin,
+#        easting_max=project.eastingmax,
+#        northing_min=project.northingmin,
+#        northing_max= project.northingmax,
+#        use_cache=True, verbose=True)
 
 # create DEM from asc data - 25m data
@@ -83 +83 @@
 
 #creates pts file for onshore DEM - 25
-dem2pts(onshore_dem_name_25, use_cache=True, verbose=True)
-
-combine_rectangular_points_files(project.onshore_dem_name + '.pts',
-                                 project.onshore_dem_name_25 + '.pts',
-                                 project.all_onshore_dem_name + '.pts')
+dem2pts(onshore_dem_name_25,
+        easting_min=project.eastingmin25,
+        easting_max=project.eastingmax25,
+        northing_min=project.northingmin25,
+        northing_max= project.northingmax25,
+        use_cache=True, verbose=True)
+
+#combine_rectangular_points_files(project.onshore_dem_name + '.pts',
+#                                 project.onshore_dem_name_25 + '.pts',
+#                                 project.all_onshore_dem_name + '.pts')
 
 print 'adding data sets'
@@ -110 +115 @@
     Geospatial_data(file_name = project.offshore_dem_name_aho15 + '.xya')+\
     Geospatial_data(file_name = project.offshore_dem_name_aho16 + '.xya')+\
-    Geospatial_data(file_name = project.all_onshore_dem_name + '.pts')
+    Geospatial_data(file_name = project.onshore_dem_name_25 + '.pts')
 G.export_points_file(project.combined_dem_name + '.pts')
-
+#Geospatial_data(file_name = project.all_onshore_dem_name + '.pts')
 #----------------------------------------------------------------------------
 # Create the triangular mesh based on overall clipping polygon with a tagged
@@ -134 +139 @@
 _ = cache(create_mesh_from_regions,
           project.polyAll,
-          {'boundary_tags': {'bottom': [0], 'right': [1],
-                             'top': [2], 'left': [3]},
-           'maximum_triangle_area': 100000,
+          {'boundary_tags': {'bottom': [0], 'bright': [1],
+                             'topr': [2], 'top': [3], 'left': [4]},
+           'maximum_triangle_area': 250000,
            'filename': meshname},           
            #'interior_regions': interior_regions},
@@ -230 +235 @@
 
 # for testing
-domain.set_boundary( {'top': Bd, 'left': Bd,
-                      'bottom': Bd, 'right': Bw} )
+domain.set_boundary( {'topr': Bd, 'left': Bd, 'top': Bd,
+                      'bottom': Bd, 'bright': Bw} )
 
 #-------------------------------------------------------------------------------