Index: /issm/trunk-jpl/src/m/contrib/morlighem/ronne/interpBedmachineAntarctica.m
===================================================================
--- /issm/trunk-jpl/src/m/contrib/morlighem/ronne/interpBedmachineAntarctica.m	(revision 23874)
+++ /issm/trunk-jpl/src/m/contrib/morlighem/ronne/interpBedmachineAntarctica.m	(revision 23875)
@@ -27,21 +27,12 @@
 	ncdate='2018-09-26';
 	ncdate='2018-11-14';
-	ncdate='2019-01-24';
 	ncdate='2019-01-30';
-	ncdate='2019-02-15';
-	ncdate='2019-02-21'; % no need of firn correction
-	ncdate='2019-02-22'; % no need of firn correction
-	ncdate='2019-02-24'; % no need of firn correction
+	ncdate='2019-03-12';
+	ncdate='2019-04-02';
+	ncdate='2019-04-15';
 	ncdate='2019-04-15'; % no need of firn correction
 end
 
-date1 = sscanf(ncdate,'%d-%d-%d');
-date2 = datetime(date1(1),date1(2),date1(3));
-
-if date2<datetime(2016,10,24),
-	basename = 'AntarcticaMCdataset';
-else
-  basename = 'BedMachineAntarctica';
-end
+basename = 'BedMachineAntarctica';
 
 %read data
Index: /issm/trunk-jpl/src/m/contrib/morlighem/ronne/interpBedmachineGreenland.m
===================================================================
--- /issm/trunk-jpl/src/m/contrib/morlighem/ronne/interpBedmachineGreenland.m	(revision 23874)
+++ /issm/trunk-jpl/src/m/contrib/morlighem/ronne/interpBedmachineGreenland.m	(revision 23875)
@@ -32,5 +32,5 @@
 	ncdate='2017-09-25'; %BedMachine v3
 	ncdate='2018-06-01';
-	ncdate='2018-08-27';
+	ncdate='2019-04-18';
 end
 
Index: /issm/trunk-jpl/src/m/contrib/morlighem/ronne/interpBedmap2.m
===================================================================
--- /issm/trunk-jpl/src/m/contrib/morlighem/ronne/interpBedmap2.m	(revision 23874)
+++ /issm/trunk-jpl/src/m/contrib/morlighem/ronne/interpBedmap2.m	(revision 23875)
@@ -9,5 +9,5 @@
 %      5. rockmask                     is a mask file showing rock outcrops
 %      6. lakemask_vostok              is a mask file showing the extent of the lake cavity of Lake Vostok
-%      7. bed_uncertainty              is the bed uncertainty grid shown in figure 12 of the manuscript
+%      7. grounded_bed_uncertainty     is the bed uncertainty grid shown in figure 12 of the manuscript
 %      8. thickness_uncertainty_5km    is the thickness uncertainty grid shown in figure 11 of the manuscript
 %      9. coverage                     is a binary grid showing the distribution of ice thickness data used in the grid of ice thickness
@@ -19,8 +19,12 @@
 
 nc = '/home/ModelData/Antarctica/BedMap2/bedmap2_bin/Bedmap2.nc';
-nc = '/Users/mmorligh/temp/Bedmap2.nc';
 
-xdata = double(ncread(nc,'x'));
-ydata = double(ncread(nc,'y'));
+if strcmp(string,'thickness_uncertainty_5km')
+	xdata = double(ncread(nc,'x_5km'));
+	ydata = double(ncread(nc,'y_5km'));
+else
+	xdata = double(ncread(nc,'x'));
+	ydata = double(ncread(nc,'y'));
+end
 
 offset=2;
Index: /issm/trunk-jpl/src/m/contrib/morlighem/ronne/interpChuter2015.m
===================================================================
--- /issm/trunk-jpl/src/m/contrib/morlighem/ronne/interpChuter2015.m	(revision 23875)
+++ /issm/trunk-jpl/src/m/contrib/morlighem/ronne/interpChuter2015.m	(revision 23875)
@@ -0,0 +1,13 @@
+function output = interpChuter2015(X,Y),
+
+ncfile='/home/ModelData/Antarctica/ChuterBamberIceShelfH/ChuterBamber_2015_CS2_ice_equivalent_ice_shelf_thickness_Rignot_gl.nc';
+verbose = 0;
+
+if verbose, disp('   -- Chuter2015: loading coordinates'); end
+xdata = double(ncread(ncfile,'x_dimensions'))';
+ydata = double(ncread(ncfile,'y_dimensions'))';
+
+if verbose, disp(['   -- Chuter2015: loading thickenss']); end
+data  = double(ncread(ncfile,'ice_shelf_thickness'))';
+if verbose, disp(['   -- Chuter2015: interpolating ' string]); end
+output = InterpFromGrid(xdata(1,:),ydata(:,1),data,X,Y);
Index: /issm/trunk-jpl/src/m/contrib/morlighem/ronne/interpDTU19MDT.m
===================================================================
--- /issm/trunk-jpl/src/m/contrib/morlighem/ronne/interpDTU19MDT.m	(revision 23875)
+++ /issm/trunk-jpl/src/m/contrib/morlighem/ronne/interpDTU19MDT.m	(revision 23875)
@@ -0,0 +1,33 @@
+function mdt = interpDTU19MDT(X,Y,varargin);
+
+switch oshostname(),
+	case {'ronne'}
+		rootname='/ronne_2/home/ModelData/Global/DTU19MDT/dtu19mdt.mat';
+	otherwise
+		error('machine not supported yet');
+end
+verbose = 1;
+
+if nargin==3,
+	hemisphere = varargin{1};
+else
+	hemisphere = +1;
+end
+
+if hemisphere==+1,
+	if verbose, disp('   -- DTU19MDT: convert to lat/lon using Greenland projection'); end
+	[LAT,  LON  ] = xy2ll(double(X(:)),double(Y(:)),+1,45,70);
+else
+	if verbose, disp('   -- DTU19MDT: convert to lat/lon using Antarctica projection'); end
+	[LAT,  LON  ] = xy2ll(double(X(:)),double(Y(:)),-1,0,71);
+end
+pos=find(LON<0);
+LON(pos) =360+LON(pos);
+LAT=reshape(LAT,size(X));
+LON=reshape(LON,size(X));
+
+if verbose, disp('   -- DTU19MDT: loading DTU19MDT'); end
+A=load(rootname);
+
+if verbose, disp('   -- DTU19MDT: interpolating'); end
+mdt = InterpFromGrid(A.lon_ext,A.lat_ext,A.mdt_ext,LON,LAT);
Index: /issm/trunk-jpl/src/m/contrib/morlighem/ronne/interpMouginotAnt2017.m
===================================================================
--- /issm/trunk-jpl/src/m/contrib/morlighem/ronne/interpMouginotAnt2017.m	(revision 23874)
+++ /issm/trunk-jpl/src/m/contrib/morlighem/ronne/interpMouginotAnt2017.m	(revision 23875)
@@ -1,35 +1,38 @@
-function [vxout vyout]= interpRignot2012(X,Y),
+function [vxout vyout]= interpMouginotAnt2017(X,Y),
 
-filename = '/home/ModelData/Greenland/VelMouginot/RignotGreenland2012Vel.mat';
+%read data
+switch (oshostname()),
+	case {'ronne'}
+		nc = '/home/ModelData/Antarctica/MouginotVel/vel_nsidc.CF16_2.nc';
+	otherwise
+		error('hostname not supported yet');
+end
 
-
-%Figure out what subset of the matrix should be read
-load(filename,'x','y');
-velfile = matfile(filename);
+xdata = double(ncread(nc,'x'));
+ydata = double(ncread(nc,'y'));
 
 offset=2;
 
 xmin=min(X(:)); xmax=max(X(:));
-posx=find(x<=xmax);
-id1x=max(1,find(x>=xmin,1)-offset);
-id2x=min(numel(x),posx(end)+offset);
+posx=find(xdata<=xmax);
+id1x=max(1,find(xdata>=xmin,1)-offset);
+id2x=min(numel(xdata),posx(end)+offset);
 
 ymin=min(Y(:)); ymax=max(Y(:));
-%posy=find(y>=ymin);
-%id1y=max(1,find(y<=ymax,1)-offset);
-%id2y=min(numel(y),posy(end)+offset);
-posy=find(y<=ymax);
-id1y=max(1,find(y>=ymin,1)-offset);
-id2y=min(numel(y),posy(end)+offset);
+posy=find(ydata>=ymin);
+id1y=max(1,find(ydata<=ymax,1)-offset);
+id2y=min(numel(ydata),posy(end)+offset);
 
-vx = velfile.vx(id1y:id2y,id1x:id2x);
-vy = velfile.vy(id1y:id2y,id1x:id2x);
-x = x(id1x:id2x);
-y = y(id1y:id2y);
+disp(['   -- Mouginot 2017: loading velocities']);
+vxdata = double(ncread(nc,'VX',[id1x id1y],[id2x-id1x+1 id2y-id1y+1],[1 1]))';
+vydata = double(ncread(nc,'VY',[id1x id1y],[id2x-id1x+1 id2y-id1y+1],[1 1]))';
+xdata=xdata(id1x:id2x);
+ydata=ydata(id1y:id2y);
 
-%load(filename);
-vxout = InterpFromGrid(x,y,double(vx),X,Y);
-vyout = InterpFromGrid(x,y,double(vy),X,Y);
+disp(['   -- Mouginot 2017: interpolating ']);
+vxout = InterpFromGrid(xdata,ydata,vxdata,double(X),double(Y));
+vyout = InterpFromGrid(xdata,ydata,vydata,double(X),double(Y));
 
+%return vel if only one output is requested
 if nargout==1,
 	vxout = sqrt(vxout.^2+vyout.^2);
Index: /issm/trunk-jpl/src/m/contrib/morlighem/ronne/interpMouginotAnt2019.m
===================================================================
--- /issm/trunk-jpl/src/m/contrib/morlighem/ronne/interpMouginotAnt2019.m	(revision 23875)
+++ /issm/trunk-jpl/src/m/contrib/morlighem/ronne/interpMouginotAnt2019.m	(revision 23875)
@@ -0,0 +1,39 @@
+function [vxout vyout]= interpMouginotAnt2019(X,Y),
+
+%read data
+switch (oshostname()),
+	case {'ronne'}
+		nc = '/home/ModelData/Antarctica/MouginotVel/v_mix.v13Mar2019.nc';
+	otherwise
+		error('hostname not supported yet');
+end
+
+xdata = double(ncread(nc,'x'));
+ydata = double(ncread(nc,'y'));
+
+offset=2;
+
+xmin=min(X(:)); xmax=max(X(:));
+posx=find(xdata<=xmax);
+id1x=max(1,find(xdata>=xmin,1)-offset);
+id2x=min(numel(xdata),posx(end)+offset);
+
+ymin=min(Y(:)); ymax=max(Y(:));
+posy=find(ydata>=ymin);
+id1y=max(1,find(ydata<=ymax,1)-offset);
+id2y=min(numel(ydata),posy(end)+offset);
+
+disp(['   -- Mouginot 2017: loading velocities']);
+vxdata = double(ncread(nc,'VX',[id1x id1y],[id2x-id1x+1 id2y-id1y+1],[1 1]))';
+vydata = double(ncread(nc,'VY',[id1x id1y],[id2x-id1x+1 id2y-id1y+1],[1 1]))';
+xdata=xdata(id1x:id2x);
+ydata=ydata(id1y:id2y);
+
+disp(['   -- Mouginot 2017: interpolating ']);
+vxout = InterpFromGrid(xdata,ydata,vxdata,double(X),double(Y));
+vyout = InterpFromGrid(xdata,ydata,vydata,double(X),double(Y));
+
+%return vel if only one output is requested
+if nargout==1,
+	vxout = sqrt(vxout.^2+vyout.^2);
+end