Changeset 5786 for anuga_work/production/busselton/build_busselton.py

Timestamp:

Sep 25, 2008, 2:24:38 PM (15 years ago)

Author:

kristy

Message:

Updated all scripts to reflect Perth format. Also added 250m scripts

File:

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

anuga_work/production/busselton/build_busselton.py

@@ -1 @@
-"""Script for running tsunami inundation scenario for Dampier, WA, Australia.
+"""
+Script for building the elevation data to run a tsunami inundation scenario
+for busselton, WA, Australia.
 
-Source data such as elevation and boundary data is assumed to be available in
-directories specified by project.py
-The output sww file is stored in project.output_time_dir 
+Input: elevation data from project.py
+Output: pts file stored in project.topographies_dir 
+The run_busselton.py is reliant on the output of this script.
 
-The scenario is defined by a triangular mesh created from project.polygon,
-the elevation data and a simulated submarine landslide.
-
-Ole Nielsen and Duncan Gray, GA - 2005 and Jane Sexton, Nick Bartzis, GA - 2006
 """
 
@@ -27 +25 @@
 from anuga.shallow_water.data_manager import convert_dem_from_ascii2netcdf, dem2pts
 from anuga.geospatial_data.geospatial_data import *
-from anuga.shallow_water.data_manager import start_screen_catcher, copy_code_files,store_parameters
+from anuga.shallow_water.data_manager import start_screen_catcher, copy_code_files
 
 # Application specific imports
@@ -42 +40 @@
 start_screen_catcher(project.output_build_time_dir)
 
-print 'USER: ', project.user
+print 'USER:    ', project.user
 
 #-------------------------------------------------------------------------------
@@ -51 +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
 
-# topography directory filenames
+# input elevation directory filenames
 onshore_in_dir_name = project.onshore_in_dir_name
 coast_in_dir_name = project.coast_in_dir_name
@@ -65 +63 @@
 offshore_in_dir_name5 = project.offshore_in_dir_name5
 
+# output elevation directory filenames
 onshore_dir_name = project.onshore_dir_name
 coast_dir_name = project.coast_dir_name
@@ -74 +73 @@
 offshore_dir_name4 = project.offshore_dir_name4
 offshore_dir_name5 = project.offshore_dir_name5
-
 # creates DEM from asc data
-print "creates DEMs from ascii data"
+print "creates DEMs from asc data"
 convert_dem_from_ascii2netcdf(onshore_in_dir_name, basename_out=onshore_dir_name, 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,
-#        easting_min=project.eastingmin,
-#        easting_max=project.eastingmax,
-#        northing_min=project.northingmin,
-#        northing_max= project.northingmax,
-        use_cache=True, 
-        verbose=True)
+dem2pts(onshore_dir_name ,use_cache=True,verbose=True)
 
-#creates pts file for island DEM
-#dem2pts(island_dir_name, use_cache=True, verbose=True)
+# create onshore pts files
+print'create Geospatial data1 objects from topographies'
+G1 = Geospatial_data(file_name = onshore_dir_name + '.pts')
 
-print'create Geospatial data1 objects from topographies'
-G1 = Geospatial_data(file_name = onshore_dir_name + '.pts',verbose=True)
-G2 = Geospatial_data(file_name = coast_in_dir_name + '.txt',verbose=True)
-G3 = Geospatial_data(file_name = coast_in_dir_name1 + '.txt',verbose=True)
-G_off = Geospatial_data(file_name = offshore_in_dir_name + '.txt',verbose=True)
-G_off1 = Geospatial_data(file_name = offshore_in_dir_name1 + '.txt',verbose=True)
-G_off2 = Geospatial_data(file_name = offshore_in_dir_name2 + '.txt',verbose=True)
-G_off3 = Geospatial_data(file_name = offshore_in_dir_name3 + '.txt',verbose=True)
-G_off4 = Geospatial_data(file_name = offshore_in_dir_name4 + '.txt',verbose=True)
-G_off5 = Geospatial_data(file_name = offshore_in_dir_name5 + '.txt',verbose=True)
+# create coastal and offshore pts files
+print'create Geospatial data2 objects from topographies'
+G3 = Geospatial_data(file_name = coast_in_dir_name)
+print'create Geospatial data3 objects from topographies'
+G4 = Geospatial_data(file_name = coast_in_dir_name1)
+print'create Geospatial data4 objects from topographies'
+G_off = Geospatial_data(file_name = offshore_in_dir_name)
+print'create Geospatial data5 objects from topographies'
+G_off1 = Geospatial_data(file_name = offshore_in_dir_name1)
+print'create Geospatial data6 objects from topographies'
+G_off2 = Geospatial_data(file_name = offshore_in_dir_name2)
+print'create Geospatial data7 objects from topographies'
+G_off3 = Geospatial_data(file_name = offshore_in_dir_name3)
+print'create Geospatial data8 objects from topographies'
+G_off4 = Geospatial_data(file_name = offshore_in_dir_name4)
+print'create Geospatial data9 objects from topographies'
+G_off5 = Geospatial_data(file_name = offshore_in_dir_name5)
+
+
+#-------------------------------------------------------------------------------
+# Combine, clip and export dataset 
+#-------------------------------------------------------------------------------
+
 print'add all geospatial objects'
 G = G1 + G2 + G3 + G_off + G_off1 + G_off2 + G_off3 + G_off4 + G_off5 
@@ -111 +117 @@
 if access(project.topographies_dir,F_OK) == 0:
     mkdir (project.topographies_dir)
-G_clipped.export_points_file(project.combined_dir_name + '.txt')
+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
 
-print'project.combined_dir_name + .txt',project.combined_dir_name + '.txt'
-
-###-------------------------------------------------------------------------
-### Convert URS to SWW file for boundary conditions 
-###-------------------------------------------------------------------------
-##print 'starting to create boundary conditions'
-##from anuga.shallow_water.data_manager import urs2sww, urs_ungridded2sww
-##
-##boundaries_in_dir_name = project.boundaries_in_dir_name
-##print 'boundaries_in_dir_name',project.boundaries_in_dir_name
-##
-##
-###import sys; sys.exit()
-##
-##urs_ungridded2sww(project.boundaries_in_dir_name, project.boundaries_in_dir_name,
-##                  verbose=True, mint=4000, maxt=80000, zscale=1)
-##