Changeset 8715 for issm/trunk/src/c/modules/Ll2xyx/Ll2xyx.cpp

Timestamp:

06/24/11 14:15:04 (14 years ago)

Author:

jschierm

Message:

Extract defaults from Xy2llx and Ll2xyx into separate functions (and fix reformatting).

File:

• issm/trunk/src/c/modules/Ll2xyx/Ll2xyx.cpp (modified) (3 diffs)

issm/trunk/src/c/modules/Ll2xyx/Ll2xyx.cpp

@@ -22 +22 @@
 %                               -1 : south latitude (default is mer=0  lat=71)
 */
-        double  delta,slat;
-        bool    flag=true;
+        double  central_meridian,standard_parallel;
 
-        /*  Get central_meridian and standard_parallel depending on hemisphere  */
-        if (sgn ==  1) {
-                delta= 45;
-                slat = 70;
-                _printf_(flag,"Info: creating coordinates in polar stereographic (Std Latitude: 70N Meridian: 45).\n");
-        }
-        else if (sgn == -1) {
-                delta= 0;
-                slat = 71;
-                _printf_(flag,"Info: creating coordinates in polar stereographic (Std Latitude: 71S Meridian: 0).\n");
-        }
-        else _error_("Sign should be either +1 or -1.\n");
+        Ll2xydef(&central_meridian,&standard_parallel,sgn);
 
-        return(Ll2xyx(x,y, lat,lon,ncoord, sgn,delta,slat));
+        return(Ll2xyx(x,y,lat,lon,ncoord,sgn,central_meridian,standard_parallel));
 }
 
@@ -59 +47 @@
         int     i,iret=0;
         double  delta,slat;
-        double  cde,re,ex2,ex,pi;
+        double  cde,re,ex2,ex;
         double  latitude,longitude;
         double  T,rho,sl,tc,mc;
@@ -105 +93 @@
         return(iret);
 }
+
+void Ll2xydef(double* pdelta, double* pslat, int sgn){
+/*  This is a cpp conversion of the following:
+%LL2XY - converts lat long to polar stereographic
+%
+%   Converts from geodetic latitude and longitude to Polar 
+%   Stereographic (X,Y) coordinates for the polar regions.
+%   Author: Michael P. Schodlok, December 2003 (map2ll)
+%
+%   Usage:
+%      [x,y] = ll2xy(lat,lon,sgn)
+%      [x,y] = ll2xy(lat,lon,sgn,central_meridian,standard_parallel)
+%
+%      - sgn = Sign of latitude +1 : north latitude (default is mer=45 lat=70)
+%                               -1 : south latitude (default is mer=0  lat=71)
+*/
+        bool    flag=true;
+
+        /*  Get central_meridian and standard_parallel depending on hemisphere  */
+        if (sgn ==  1) {
+                *pdelta= 45;
+                *pslat = 70;
+                _printf_(flag,"Info: creating coordinates in polar stereographic (Std Latitude: 70N Meridian: 45).\n");
+        }
+        else if (sgn == -1) {
+                *pdelta= 0;
+                *pslat = 71;
+                _printf_(flag,"Info: creating coordinates in polar stereographic (Std Latitude: 71S Meridian: 0).\n");
+        }
+        else _error_("Sign should be either +1 or -1.\n");
+
+        return;
+}