Changeset 18478 for issm/trunk-jpl/src/m/coordsystems/utm2ll.m

Timestamp:

09/02/14 12:05:17 (11 years ago)

Author:

Mathieu Morlighem

Message:

CHG: new version after my suggestion to the author

File:

• issm/trunk-jpl/src/m/coordsystems/utm2ll.m (modified) (3 diffs)

issm/trunk-jpl/src/m/coordsystems/utm2ll.m

@@ -1 +1 @@
 function [lat,lon]=utm2ll(x,y,f,datum,varargin)
 %UTM2LL UTM to Lat/Lon coordinates precise conversion.
-%   [LAT,LON]=UTM2LL(X,Y,ZONE) converts UTM coordinates X,Y (in meters)
+%       [LAT,LON]=UTM2LL(X,Y,ZONE) converts UTM coordinates X,Y (in meters)
 %       defined in the UTM ZONE (integer) to latitude LAT and longitude LON 
 %       (in degrees). Default datum is WGS84.
@@ -30 +30 @@
 %                  Notes Techniques NT/G 76, janvier 1995.
 %
-%   Author: Francois Beauducel, <beauducel@ipgp.fr>
-%   Created: 2001-08-23
-%   Updated: 2014-04-20
+%       Author: Francois Beauducel, <beauducel@ipgp.fr>
+%       Created: 2001-08-23
+%       Updated: 2014-04-20
+
 
 %       Copyright (c) 2001-2014, François Beauducel, covered by BSD License.
@@ -105 +106 @@
 eps = 1e-11;    % minimum residue for latitude computation
 
-K0 = 0.9996;                                                            % UTM scale factor
-X0 = 500000;                                                            % UTM false East (m)
-Y0 = 1e7*(f < 0);                                                       % UTM false North (m)
-P0 = 0;                                                                         % UTM origin latitude (rad)
-L0 = (6*abs(f) - 183)/D0;                                       % UTM origin longitude (rad)
+K0 = 0.9996;                                    % UTM scale factor
+X0 = 500000;                                    % UTM false East (m)
+Y0 = 1e7*(f < 0);                               % UTM false North (m)
+P0 = 0;                                         % UTM origin latitude (rad)
+L0 = (6*abs(f) - 183)/D0;                       % UTM origin longitude (rad)
 E1 = sqrt((A1^2 - (A1*(1 - 1/F1))^2)/A1^2);     % ellpsoid excentricity
 N = K0*A1;