source: issm/oecreview/Archive/26740-27031/ISSM-26957-26958.diff@ 27032

Index: ../trunk-jpl/src/m/coordsystems/ll2xy.m
===================================================================
--- ../trunk-jpl/src/m/coordsystems/ll2xy.m     (revision 26957)
+++ ../trunk-jpl/src/m/coordsystems/ll2xy.m     (revision 26958)
@@ -39,6 +39,20 @@
        error('bad usage');
 end
 
+%Move to CoordTransform now...
+if exist('CoordTransform_matlab')==3
+       disp('Calling CoordTransform instead, make sure to change your MATLAB script');
+       if sgn==+1
+               assert(delta==45); assert(slat ==70);
+               [x y]=CoordTransform(lat,lon,'EPSG:4326','EPSG:3413');
+               return;
+       else
+               assert(delta==0); assert(slat ==71);
+               [x y]=CoordTransform(lat,lon,'EPSG:4326','EPSG:3031');
+               return;
+       end
+end
+
 % Conversion constant from degrees to radians
 cde  = 57.29577951;
 % Radius of the earth in meters
Index: ../trunk-jpl/src/m/coordsystems/xy2ll.m
===================================================================
--- ../trunk-jpl/src/m/coordsystems/xy2ll.m     (revision 26957)
+++ ../trunk-jpl/src/m/coordsystems/xy2ll.m     (revision 26958)
@@ -35,6 +35,20 @@
        error('bad usage');
 end
 
+%Move to CoordTransform now...
+if exist('CoordTransform_matlab')==3
+       disp('Calling CoordTransform instead, make sure to change your MATLAB script');
+       if sgn==+1
+               assert(delta==45); assert(slat ==70);
+               [lat lon]=CoordTransform(x, y,'EPSG:3413','EPSG:4326');
+               return;
+       else
+               assert(delta==0); assert(slat ==71);
+               [lat lon]=CoordTransform(x, y,'EPSG:3031','EPSG:4326');
+               return;
+       end
+end
+
 if nargout~=3 & nargout~=2,
        help xy2ll
        error('bad usage');