Changeset 8689 for issm/trunk/src/mex/Ll2xy/Ll2xy.cpp

Timestamp:

06/22/11 15:51:08 (14 years ago)

Author:

jschierm

Message:

Updated versions of Xy2ll and Ll2xy modules and x-layers.

File:

• issm/trunk/src/mex/Ll2xy/Ll2xy.cpp (modified) (3 diffs)

issm/trunk/src/mex/Ll2xy/Ll2xy.cpp

@@ -18 +18 @@
         /* output: */
         double  *x=NULL,*y=NULL;
+        int     iret=0;
 
         /*Boot module: */
@@ -38 +39 @@
                 if   (nlat == 1) printf("  lat=%g\n",lat[0]);
                 else             printf("  lat=[%d values]\n",nlat);
+//      for (i=0; i<nlat; i++) printf("  lat[%d]=%g\n",i,lat[i]);
         FetchData(&lon,&nlon,LON_IN);
         if (verbose)
                 if   (nlon == 1) printf("  lon=%g\n",lon[0]);
                 else             printf("  lon=[%d values]\n",nlon);
+//      for (i=0; i<nlon; i++) printf("  lon[%d]=%g\n",i,lon[i]);
         FetchData(&sgn,SGN_IN);
         if (verbose) printf("  sgn=%d\n",sgn);
@@ -71 +74 @@
         if (verbose) printf("Calling core:\n");
         if (options->GetOption("central_meridian") && options->GetOption("standard_parallel"))
-                Ll2xyx(x,y,
-                           lat,lon,ncoord,
-                           sgn,cm,sp);
+                iret=Ll2xyx(x,y,
+                                        lat,lon,ncoord,
+                                        sgn,cm,sp);
         else
-                Ll2xyx(x,y,
-                           lat,lon,ncoord,
-                           sgn);
+                iret=Ll2xyx(x,y,
+                                        lat,lon,ncoord,
+                                        sgn);
+        if (verbose) printf("  iret=%d\n",iret);
 
         /*Write data: */