Index: /issm/trunk-jpl/src/m/coordsystems/ll2xy.m
===================================================================
--- /issm/trunk-jpl/src/m/coordsystems/ll2xy.m	(revision 26957)
+++ /issm/trunk-jpl/src/m/coordsystems/ll2xy.m	(revision 26958)
@@ -40,4 +40,18 @@
 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;
Index: /issm/trunk-jpl/src/m/coordsystems/xy2ll.m
===================================================================
--- /issm/trunk-jpl/src/m/coordsystems/xy2ll.m	(revision 26957)
+++ /issm/trunk-jpl/src/m/coordsystems/xy2ll.m	(revision 26958)
@@ -34,4 +34,18 @@
 	help xy2ll
 	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