Changeset 5784 for anuga_work/production/perth/build_perth.py

Timestamp:

Sep 25, 2008, 1:14:25 PM (16 years ago)

Author:

sexton

Message:

updated commments for perth scripts

File:

• anuga_work/production/perth/build_perth.py (modified) (8 diffs)

anuga_work/production/perth/build_perth.py

@@ -1 +1 @@
 """
-Script for building the elevation data to run tsunami inundation scenario
+Script for building the elevation data to run a tsunami inundation scenario
 for Perth, WA, Australia.
 
@@ -7 +7 @@
 The run_perth.py is reliant on the output of this script.
 
-Ole Nielsen and Duncan Gray, GA - 2005 and Jane Sexton, Nick Bartzis, GA - 2006
 """
 
@@ -50 +49 @@
 # Fine pts file to be clipped to area of interest
 #-------------------------------------------------------------------------------
-print"project.poly_all",project.poly_all
-print"project.combined_dir_name",project.combined_dir_name
+print "project.poly_all", project.poly_all
+print "project.combined_dir_name", project.combined_dir_name
 
 # input elevation directory filenames
@@ -85 +84 @@
 convert_dem_from_ascii2netcdf(island_in_dir_name3, basename_out=island_dir_name3, use_cache=True, verbose=True)
 
-#creates pts file for onshore DEM
+# creates pts file for onshore DEM
 print "creates pts file for onshore DEM"
 dem2pts(onshore_dir_name ,use_cache=True,verbose=True)
 
-#creates pts file for island DEM
+# creates pts file for island DEM
 dem2pts(island_dir_name, use_cache=True, verbose=True)
 dem2pts(island_dir_name1, use_cache=True, verbose=True)
@@ -95 +94 @@
 dem2pts(island_dir_name3, use_cache=True, verbose=True)
 
-#-------------------------------------------------------------------------------
-# Combine datasets into project.combined_dir_name
-#-------------------------------------------------------------------------------
-
+# create onshore pts files
 print'create Geospatial data1 objects from topographies'
 G1 = Geospatial_data(file_name = onshore_dir_name + '.pts')
@@ -109 +105 @@
 print'create Geospatial data5 objects from topographies'
 G5 = Geospatial_data(file_name = island_dir_name3 + '.pts')
+
+# create coastal and offshore pts files
 print'create Geospatial data6 objects from topographies'
 G_coast = Geospatial_data(file_name = coast_in_dir_name)
@@ -120 +118 @@
 G_off3 = Geospatial_data(file_name = offshore_in_dir_name3)
 
+#-------------------------------------------------------------------------------
+# Combine, clip and export dataset 
+#-------------------------------------------------------------------------------
 
 print'add all geospatial objects'
@@ -130 +131 @@
 if access(project.topographies_dir,F_OK) == 0:
     mkdir (project.topographies_dir)
-#G_clipped.export_points_file(project.combined_dir_name + '.pts')
-G_clipped.export_points_file(project.combined_dir_name + '.txt') #Use for comparision in ARC
+G_clipped.export_points_file(project.combined_dir_name + '.pts')
+#G_clipped.export_points_file(project.combined_dir_name + '.txt') #Use for comparision in ARC
 
-import sys
-sys.exit()