Changeset 6359 for anuga_work/production/australia_ph2/ceduna/build_elevation.py

Timestamp:

Feb 18, 2009, 9:14:00 AM (14 years ago)

Author:

myall

Message:

phase 2 scripts for ceduna, eucla and esperance still not working...

File:

• anuga_work/production/australia_ph2/ceduna/build_elevation.py (modified) (1 diff)

anuga_work/production/australia_ph2/ceduna/build_elevation.py

@@ -50 +50 @@
 # Create Geospatial data from ASCII files
 geospatial_data = {}
-for filename in 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)
+
+# Create Geospatial data from TXT files
+for filename in project.point_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',
+    geospatial_data[filename] = Geospatial_data(file_name=absolute_filename,
                                                 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)