Index: /issm/trunk-jpl/src/m/classes/model.m
===================================================================
--- /issm/trunk-jpl/src/m/classes/model.m	(revision 21498)
+++ /issm/trunk-jpl/src/m/classes/model.m	(revision 21499)
@@ -308,8 +308,4 @@
 			end
 
-			%lat long
-			if numel(md.mesh.lat) ==md.mesh.numberofvertices,  md.mesh.lat=project2d(md,md.mesh.lat,1); end
-			if numel(md.mesh.long)==md.mesh.numberofvertices, md.mesh.long=project2d(md,md.mesh.long,1); end
-
 			%Initialize with the 2d mesh
 			mesh=mesh2d();
@@ -319,4 +315,7 @@
 			mesh.numberofelements=md.mesh.numberofelements2d;
 			mesh.elements=md.mesh.elements2d;
+			if numel(md.mesh.lat) ==md.mesh.numberofvertices,  mesh.lat=project2d(md,md.mesh.lat,1); end
+			if numel(md.mesh.long)==md.mesh.numberofvertices,	mesh.long=project2d(md,md.mesh.long,1); end
+			mesh.epsg=md.mesh.epsg;
 			if ~isnan(md.mesh.vertexonboundary), mesh.vertexonboundary=project2d(md,md.mesh.vertexonboundary,1); end
 			if ~isnan(md.mesh.elementconnectivity), mesh.elementconnectivity=project2d(md,md.mesh.elementconnectivity,1); end
Index: /issm/trunk-jpl/src/m/classes/model.py
===================================================================
--- /issm/trunk-jpl/src/m/classes/model.py	(revision 21498)
+++ /issm/trunk-jpl/src/m/classes/model.py	(revision 21499)
@@ -757,10 +757,4 @@
 		md.mask.ice_levelset=project2d(md,md.mask.ice_levelset,1)
 
-		#lat long
-                if isinstance(md.mesh.lat,np.ndarray):
-                    if md.mesh.lat.size==md.mesh.numberofvertices:  md.mesh.lat=project2d(md,md.mesh.lat,1) 
-                if isinstance(md.mesh.long,np.ndarray):
-		    if md.mesh.long.size==md.mesh.numberofvertices: md.mesh.long=project2d(md,md.mesh.long,1) 
-
 		#Initialize with the 2d mesh
 		mesh=mesh2d()
@@ -772,4 +766,9 @@
 		if not np.isnan(md.mesh.vertexonboundary).all(): mesh.vertexonboundary=project2d(md,md.mesh.vertexonboundary,1)
 		if not np.isnan(md.mesh.elementconnectivity).all(): mesh.elementconnectivity=project2d(md,md.mesh.elementconnectivity,1)
+		if isinstance(md.mesh.lat,np.ndarray):
+			if md.mesh.lat.size==md.mesh.numberofvertices:  mesh.lat=project2d(md,md.mesh.lat,1) 
+		if isinstance(md.mesh.long,np.ndarray):
+			if md.mesh.long.size==md.mesh.numberofvertices: mesh.long=project2d(md,md.mesh.long,1) 
+		mesh.epsg=md.mesh.epsg
 		md.mesh=mesh