Index: /issm/trunk-jpl/src/m/coordsystems/ll2xy.m
===================================================================
--- /issm/trunk-jpl/src/m/coordsystems/ll2xy.m	(revision 27037)
+++ /issm/trunk-jpl/src/m/coordsystems/ll2xy.m	(revision 27038)
@@ -44,7 +44,14 @@
 	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;
+		if delta==45 && slat==70      %BedMachine
+			[x y]=CoordTransform(lat,lon,'EPSG:4326','EPSG:3413'); return;
+		elseif delta==0 && slat==75   %IBCAO
+			[x y]=CoordTransform(lat,lon,'EPSG:4326','EPSG:3996'); return;
+		elseif delta==39 && slat==71  %Bamber
+			Bamber_proj = '+proj=stere +lat_0=90 +lat_ts=71 +lon_0=-39 +k=1 +x_0=0 +y_0=0 +ellps=WGS84 +datum=WGS84 +units=m +no_defs';
+			[x y]=CoordTransform(lat,lon,'EPSG:4326',Bamber_proj); return;
+		else
+			error('not supported yet');
+		end
 	else
 		assert(delta==0); assert(slat ==71);
Index: /issm/trunk-jpl/src/m/coordsystems/xy2ll.m
===================================================================
--- /issm/trunk-jpl/src/m/coordsystems/xy2ll.m	(revision 27037)
+++ /issm/trunk-jpl/src/m/coordsystems/xy2ll.m	(revision 27038)
@@ -40,7 +40,14 @@
 	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;
+		if delta==45 && slat==70      %BedMachine
+			[lat lon]=CoordTransform(x, y,'EPSG:3413','EPSG:4326'); return;
+		elseif delta==0 && slat==75   %IBCAO
+			[lat lon]=CoordTransform(x, y,'EPSG:3996','EPSG:4326'); return;
+		elseif delta==39 && slat==71  %Bamber
+			Bamber_proj = '+proj=stere +lat_0=90 +lat_ts=71 +lon_0=-39 +k=1 +x_0=0 +y_0=0 +ellps=WGS84 +datum=WGS84 +units=m +no_defs';
+			[lat lon]=CoordTransform(x, y, Bamber_proj, 'EPSG:4326'); return;
+		else
+			error('not supported yet');
+		end
 	else
 		assert(delta==0); assert(slat ==71);