Index: /issm/trunk-jpl/examples/IceBridge/Greenland.par
===================================================================
--- /issm/trunk-jpl/examples/IceBridge/Greenland.par	(revision 20745)
+++ /issm/trunk-jpl/examples/IceBridge/Greenland.par	(revision 20746)
@@ -26,4 +26,27 @@
 md.geometry.thickness(pos0)=1;
 md.geometry.surface=md.geometry.thickness+md.geometry.base;
+
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+%Reading IceBridge data for Jakobshavn
+disp('      reading IceBridge Jakobshavn bedrock');
+fid  = fopen('../Data/Jakobshavn_2008_2011_Composite_XYZGrid.txt');
+titles = fgets(fid); 
+data = fscanf(fid,'%g,%g,%g,%g,%g',[5 266400])';
+fclose(fid);
+
+[xi,yi]= ll2xy(md.mesh.lat,md.mesh.long,+1,45,70);
+bed  = flipud(reshape(data(:,5),[360 740])); bed(find(bed==-9999))=NaN;
+bedy = flipud(reshape(data(:,1),[360 740]));
+bedx = flipud(reshape(data(:,2),[360 740]));
+
+%Insert Icebridge bed and recalculate thickness
+bed_jks=InterpFromGridToMesh(bedx(1,:)',bedy(:,1),bed,xi,yi,NaN);
+in=ContourToMesh(md.mesh.elements,md.mesh.x,md.mesh.y,...
+	   'Jak_grounded.exp','node',1);
+bed_jks(~in)=NaN;
+pos=find(~isnan(bed_jks));
+md.geometry.base(pos)=bed_jks(pos);
+md.geometry.thickness=md.geometry.surface-md.geometry.base;
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 
 disp('   Interpolating velocities ');
@@ -58,6 +81,6 @@
 disp('   Set other boundary conditions');
 md.mask.ice_levelset(md.mesh.vertexonboundary==1)=0;
+md.basalforcings.groundedice_melting_rate = zeros(md.mesh.numberofvertices,1);
 md.basalforcings.floatingice_melting_rate = zeros(md.mesh.numberofvertices,1);
-md.basalforcings.groundedice_melting_rate = zeros(md.mesh.numberofvertices,1);
 md.thermal.spctemperature     = [md.initialization.temperature;1]; %impose observed temperature on surface
 md.masstransport.spcthickness    = NaN*ones(md.mesh.numberofvertices,1);
@@ -70,8 +93,7 @@
 
 disp('   Single point constraints');
-%Initialize single point constraint arrarys
+%Initialize single point constraint arrays
 md.stressbalance.referential=NaN*ones(md.mesh.numberofvertices,6);
 md.stressbalance.spcvx = NaN*ones(md.mesh.numberofvertices,1);
 md.stressbalance.spcvy = NaN*ones(md.mesh.numberofvertices,1);
 md.stressbalance.spcvz = NaN*ones(md.mesh.numberofvertices,1);
-
Index: /issm/trunk-jpl/examples/IceBridge/Greenland_solution.par
===================================================================
--- /issm/trunk-jpl/examples/IceBridge/Greenland_solution.par	(revision 20745)
+++ /issm/trunk-jpl/examples/IceBridge/Greenland_solution.par	(revision 20746)
@@ -1,5 +1,5 @@
-%Name and hemisphere
+%Name and Coordinate system
 md.miscellaneous.name='SeaRISEgreenland';
-
+md.mesh.epsg=3413;
 
 disp('   Loading SeaRISE data from NetCDF');
@@ -31,5 +31,6 @@
 disp('      reading IceBridge Jakobshavn bedrock');
 fid  = fopen('../Data/Jakobshavn_2008_2011_Composite_XYZGrid.txt');
-titles = fgets(fid); data = fscanf(fid,'%g,%g,%g,%g,%g',[5 266400])';
+titles = fgets(fid); 
+data = fscanf(fid,'%g,%g,%g,%g,%g',[5 266400])';
 fclose(fid);
 
@@ -44,5 +45,5 @@
 surf_jks=InterpFromGridToMesh(bedx(1,:)',bedy(:,1),surf,xi,yi,NaN);
 in=ContourToMesh(md.mesh.elements,md.mesh.x,md.mesh.y,...
-	'Jak_grounded.exp','node',1);
+	   'Jak_grounded.exp','node',1);
 bed_jks(~in)=NaN;
 surf_jks(~in)=NaN;
@@ -66,6 +67,6 @@
 
 disp('   Interpolating surface mass balance');
-md.surfaceforcings.mass_balance=InterpFromGridToMesh(x1,y1,smb,md.mesh.x,md.mesh.y,0);
-md.surfaceforcings.mass_balance=md.surfaceforcings.mass_balance*md.materials.rho_water/md.materials.rho_ice;
+md.smb.mass_balance=InterpFromGridToMesh(x1,y1,smb,md.mesh.x,md.mesh.y,0);
+md.smb.mass_balance=md.smb.mass_balance*md.materials.rho_water/md.materials.rho_ice;
 
 disp('   Construct basal friction parameters');
@@ -84,5 +85,6 @@
 disp('   Set other boundary conditions');
 md.mask.ice_levelset(md.mesh.vertexonboundary==1)=0;
-md.basalforcings.melting_rate = zeros(md.mesh.numberofvertices,1);
+md.basalforcings.groundedice_melting_rate = zeros(md.mesh.numberofvertices,1);
+md.basalforcings.floatingice_melting_rate = zeros(md.mesh.numberofvertices,1);
 md.thermal.spctemperature     = [md.initialization.temperature;1]; %impose observed temperature on surface
 md.masstransport.spcthickness    = NaN*ones(md.mesh.numberofvertices,1);
@@ -94,14 +96,8 @@
 md.initialization.pressure=md.materials.rho_ice*md.constants.g*md.geometry.thickness;
 
-disp('   Single point constraint for continental model');
-%Set at least one vertex to velocity 0 so as to not get a singular problem (point on the wet peninsula)
+disp('   Single point constraints');
+%Initialize single point constraint arrays
 md.stressbalance.referential=NaN*ones(md.mesh.numberofvertices,6);
 md.stressbalance.spcvx = NaN*ones(md.mesh.numberofvertices,1);
 md.stressbalance.spcvy = NaN*ones(md.mesh.numberofvertices,1);
 md.stressbalance.spcvz = NaN*ones(md.mesh.numberofvertices,1);
-location = 1.0e+06 *[.32011 -2.2039];
-[dist pos]=min(sqrt((md.mesh.x - location(1)).^2 + (md.mesh.y - location(2)).^2));
-md.stressbalance.spcvx(pos) = 0;
-md.stressbalance.spcvy(pos) = 0;
-md.stressbalance.spcvz(pos) = 0;
-