Changeset 6376 for anuga_work/production/busselton/standardised_version/build_elevation.py

Timestamp:

Feb 23, 2009, 2:17:06 PM (14 years ago)

Author:

myall

Message:

Ready for use, fixed known problems

File:

• anuga_work/production/busselton/standardised_version/build_elevation.py (modified) (1 diff)

anuga_work/production/busselton/standardised_version/build_elevation.py

@@ -50 +50 @@
 # Create Geospatial data from ASCII files
 geospatial_data = {}
-for filename in project.ascii_grid_filenames:
-    absolute_filename = join(project.topographies_folder, filename)
-    convert_dem_from_ascii2netcdf(absolute_filename,
-                                  basename_out=absolute_filename,
-                                  use_cache=True,
-                                  verbose=True)
-    dem2pts(absolute_filename, use_cache=True, verbose=True)
+if not project.ascii_grid_filenames == []:
+    for filename in project.ascii_grid_filenames:
+        absolute_filename = join(project.topographies_folder, filename)
+        convert_dem_from_ascii2netcdf(absolute_filename,
+                                      basename_out=absolute_filename,
+                                      use_cache=True,
+                                      verbose=True)
+        dem2pts(absolute_filename, use_cache=True, verbose=True)
 
-    geospatial_data[filename] = Geospatial_data(file_name=absolute_filename+'.pts',
-                                                verbose=True)
+        geospatial_data[filename] = Geospatial_data(file_name=absolute_filename+'.pts',
+                                                    verbose=True)
 
 # Create Geospatial data from TXT files
-for filename in project.point_filenames:
-    absolute_filename = join(project.topographies_folder, filename)
-    geospatial_data[filename] = Geospatial_data(file_name=absolute_filename,
-                                                verbose=True)
-
+if not project.point_filenames == []:
+    for filename in project.point_filenames:
+        absolute_filename = join(project.topographies_folder, filename)
+        geospatial_data[filename] = Geospatial_data(file_name=absolute_filename,
+                                                    verbose=True)
 
 #-------------------------------------------------------------------------------