Changeset 6279 for anuga_work/production/busselton/standardised_version/build_busselton.py

Timestamp:

Feb 5, 2009, 2:52:42 PM (16 years ago)

Author:

kristy

Message:

Updated standardised, not tested yet

File:

• anuga_work/production/busselton/standardised_version/build_busselton.py (modified) (4 diffs)

anuga_work/production/busselton/standardised_version/build_busselton.py

@@ -30 +30 @@
 # output to file
 #------------------------------------------------------------------------------
-copy_code_files(project.output_build_time_dir,__file__, 
+copy_code_files(project.output_build,__file__, 
                 os.path.dirname(project.__file__)+os.sep+\
                 project.__name__+'.py' )
-start_screen_catcher(project.output_build_time_dir)
+start_screen_catcher(project.output_build)
 
 
@@ -44 +44 @@
 #------------------------------------------------------------------------------
 
-print 'project.poly_all', project.poly_all
-print 'project.combined_dir_name', project.combined_dir_name
+print 'project.bounding_polygon', project.bounding_polygon
+print 'project.combined_elevation_basename', project.combined_elevation_basename
 
 # Create Geospatial data from ASCII files
 geospatial_data = {}
 for filename in project.ascii_grid_filenames:
-    absolute_filename = project.topographies_dir + filename
+    absolute_filename = project.topographies_folder + os.sep + filename
     convert_dem_from_ascii2netcdf(absolute_filename,
                                   basename_out=absolute_filename,
@@ -62 +62 @@
 # Create Geospatial data from TXT files
 for filename in project.point_filenames:
-    absolute_filename = project.topographies_dir + filename
+    absolute_filename = project.topographies_folder  + os.sep + filename
     geospatial_data[filename] = Geospatial_data(file_name=absolute_filename,
                                                 verbose=True)
@@ -71 +71 @@
 #-------------------------------------------------------------------------------
 
-print 'Add geospatial objects except', project.offshore_name5
-G = None
-for key in geospatial_data:
-    if key != project.offshore_name5:
-        G += geospatial_data[key]
+print 'Add geospatial objects' # except', project.offshore_name5
+G = geospatial_data
+##None
+##for key in geospatial_data:
+##    if key != project.offshore_name5:
+##        G += geospatial_data[key]
         
 print 'Clip combined geospatial data'
-G_clip = G.clip_outside(project.poly_aoi1)
-G_all = G_clip + geospatial_data[project.offshore_name5]   # G_off5
-G_clipped = G_all.clip(project.poly_all)
+##G_clip = G.clip_outside(project.poly_aoi1)
+##G_all = G_clip + geospatial_data[project.offshore_name5]   # G_off5
+#G_clipped = G_all.clip(project.poly_all)
+G_clip = G.clip(project.bounding_polygon)
+
 
 print 'Export combined DEM file'
-G_clipped.export_points_file(project.combined_dir_name + '.pts')
+G_clip.export_points_file(project.combined_dir_name + '.pts')
 print 'Do txt version too'
 # Use for comparision in ARC
-G_clipped.export_points_file(project.combined_dir_name + '.txt')
+G_clip.export_points_file(project.combined_dir_name + '.txt')