Changeset 6424 for anuga_core/source/anuga/coordinate_transforms/redfearn.py

Timestamp:

Feb 26, 2009, 10:43:04 PM (14 years ago)

Author:

ole

Message:

Experimented with scale factors in redfearn

File:

• anuga_core/source/anuga/coordinate_transforms/redfearn.py (modified) (2 diffs)

anuga_core/source/anuga/coordinate_transforms/redfearn.py

@@ -38 +38 @@
 
 def redfearn(lat, lon, false_easting=None, false_northing=None,
-             zone=None, central_meridian=None):
+             zone=None, central_meridian=None, scale_factor=None):
     """Compute UTM projection using Redfearn's formula
 
@@ -63 +63 @@
     a = 6378137.0                       #Semi major axis
     inverse_flattening = 298.257222101  #1/f
-    K0 = 0.9996                         #Central scale factor    
+    if scale_factor is None:
+        K0 = 0.9996                         #Central scale factor    
+    else:
+        K0 = scale_factor
+    #print 'scale', K0
     zone_width = 6                      #Degrees