Changeset 4000

Timestamp:

Nov 16, 2006, 10:47:58 AM (17 years ago)

Author:

sexton

Message:

update GIS documentation

Location:

anuga_core

Files:

• documentation/user_manual/How to develop a polygon for scenario models using ArcGIS.pdf (modified) (previous)
• source/anuga/extras/thinbathy.py (modified) (2 diffs)

anuga_core/source/anuga/extras/thinbathy.py

@@ -15 +15 @@
 
     fid = NetCDFFile(input_file, 'r')
+    print fid.variables
 
+    # for first event
     latitudes = fid.variables['LAT'][::lat_thin]
     longitudes = fid.variables['LON'][::lon_thin]
@@ -21 +23 @@
 
     outfile = NetCDFFile(output_file, 'w')
-    outfile.createDimension('LON', fid.dimensions['LON']/lon_thin) # + 1)
-    outfile.createDimension('LAT', fid.dimensions['LAT']/lat_thin) # + 1)
+    outfile.createDimension('LON', fid.dimensions['LON']/lon_thin)
+    outfile.createDimension('LAT', fid.dimensions['LAT']/lat_thin)
+    
+##    # for second event
+##    latitudes = fid.variables['LAT601_690'][::lat_thin]
+##    longitudes = fid.variables['LON645_705'][::lon_thin]
+##    elevations = fid.variables['ELEVATION'][::lat_thin,::lon_thin]
+##
+##    outfile = NetCDFFile(output_file, 'w')
+##    outfile.createDimension('LON', fid.dimensions['LON645_705']/lon_thin + 1)
+##    outfile.createDimension('LAT', fid.dimensions['LAT601_690']/lat_thin + 1)
 
     import Numeric