Index: /issm/trunk-jpl/src/m/contrib/chenggong/dataprocessing/analyzeCalvingFront.m
===================================================================
--- /issm/trunk-jpl/src/m/contrib/chenggong/dataprocessing/analyzeCalvingFront.m	(revision 27965)
+++ /issm/trunk-jpl/src/m/contrib/chenggong/dataprocessing/analyzeCalvingFront.m	(revision 27966)
@@ -5,6 +5,4 @@
 %                   [Xmain, positionx, positiony, time], icemask, calving rate,
 %                   melting rate and velocity magnitude, sigma, thickness
-% Author: Cheng Gong
-% Last modified: 2020-09-25
 
 transient = isfield(md.results,'TransientSolution');
Index: /issm/trunk-jpl/src/m/contrib/chenggong/dataprocessing/averageOverTime.m
===================================================================
--- /issm/trunk-jpl/src/m/contrib/chenggong/dataprocessing/averageOverTime.m	(revision 27965)
+++ /issm/trunk-jpl/src/m/contrib/chenggong/dataprocessing/averageOverTime.m	(revision 27966)
@@ -13,6 +13,4 @@
 % if within the range, do linear interpolation.
 %
-%   Author: Cheng Gong
-%   Last modified: 2020-09-09
 if nargin < 4
     endP = startP;
Index: /issm/trunk-jpl/src/m/contrib/chenggong/dataprocessing/fillInNan.m
===================================================================
--- /issm/trunk-jpl/src/m/contrib/chenggong/dataprocessing/fillInNan.m	(revision 27965)
+++ /issm/trunk-jpl/src/m/contrib/chenggong/dataprocessing/fillInNan.m	(revision 27966)
@@ -3,6 +3,4 @@
 %		data should have the same size as md.mesh.x
 %
-% Author: Cheng Gong
-% Last modified: 2021-12-08
 nanvflag = find(isnan(data));
 NNanv = length(nanvflag);
Index: /issm/trunk-jpl/src/m/contrib/chenggong/dataprocessing/interpZeroPos.m
===================================================================
--- /issm/trunk-jpl/src/m/contrib/chenggong/dataprocessing/interpZeroPos.m	(revision 27965)
+++ /issm/trunk-jpl/src/m/contrib/chenggong/dataprocessing/interpZeroPos.m	(revision 27966)
@@ -14,6 +14,4 @@
 %                x0 = -------------
 %                        y1-y2
-% Author: Cheng Gong
-% Last modified: 2020-10-20
 [m, n] = size(X);
 [my, k] = size(Y);
Index: /issm/trunk-jpl/src/m/contrib/chenggong/dataprocessing/npsd.m
===================================================================
--- /issm/trunk-jpl/src/m/contrib/chenggong/dataprocessing/npsd.m	(revision 27965)
+++ /issm/trunk-jpl/src/m/contrib/chenggong/dataprocessing/npsd.m	(revision 27966)
@@ -7,6 +7,4 @@
 %   freq: frequency(non-negative half)
 %
-%   Author: Cheng Gong
-%   Date: 2021-08-17
 
 % length of data
Index: /issm/trunk-jpl/src/m/contrib/chenggong/dataprocessing/projectToFlowlines.m
===================================================================
--- /issm/trunk-jpl/src/m/contrib/chenggong/dataprocessing/projectToFlowlines.m	(revision 27965)
+++ /issm/trunk-jpl/src/m/contrib/chenggong/dataprocessing/projectToFlowlines.m	(revision 27966)
@@ -6,6 +6,4 @@
 %	fy			-	y coordinates of the flowline
 %
-% Author: Cheng Gong
-% Last modified: 2021-01-27
 
 function valueC = projectToFlowlines(md, pValue, fx, fy)
Index: /issm/trunk-jpl/src/m/contrib/chenggong/dataprocessing/psd.m
===================================================================
--- /issm/trunk-jpl/src/m/contrib/chenggong/dataprocessing/psd.m	(revision 27965)
+++ /issm/trunk-jpl/src/m/contrib/chenggong/dataprocessing/psd.m	(revision 27966)
@@ -7,6 +7,4 @@
 %   freq: frequency(non-negative half)
 %
-%   Author: Cheng Gong
-%   Date: 2021-08-17
 
 % length of data
Index: /issm/trunk-jpl/src/m/contrib/chenggong/dataprocessing/wfDistance.m
===================================================================
--- /issm/trunk-jpl/src/m/contrib/chenggong/dataprocessing/wfDistance.m	(revision 27965)
+++ /issm/trunk-jpl/src/m/contrib/chenggong/dataprocessing/wfDistance.m	(revision 27966)
@@ -9,6 +9,4 @@
 %   w2:  Wasserstein-Fourier Distance, or W2 distance of sx and sy
 %
-%   Author: Cheng Gong
-%   Date: 2021-08-18
 N = length(x);
 Sinv=linspace(0,1,N);
Index: /issm/trunk-jpl/src/m/contrib/chenggong/modeldata/interpFromITSLIVE.m
===================================================================
--- /issm/trunk-jpl/src/m/contrib/chenggong/modeldata/interpFromITSLIVE.m	(revision 27965)
+++ /issm/trunk-jpl/src/m/contrib/chenggong/modeldata/interpFromITSLIVE.m	(revision 27966)
@@ -1,3 +1,3 @@
-function [vx_out, vy_out, vx_err, vy_err] = interpFromITSLIVE(X,Y,Tstart,Tend,varargin)
+function [vx_out, vy_out] = interpFromITSLIVE(X,Y,Tstart,Tend,varargin)
 	%interpFromITSLIVE: 
 	%	Interpolate ITS_LIVE velocity data to the given mesh
@@ -28,9 +28,4 @@
 		dataTime(i) = date2decyear(datenum(tempConv{end}, 'yyyy'));
 	end
-	disp(['  Found ', num2str(Ndata), ' records in ', foldername]);
-	disp(['    from ', datestr(decyear2date(min(dataTime)),'yyyy-mm-dd'), ' to ', datestr(decyear2date(max(dataTime)),'yyyy-mm-dd') ]);
-
-	% 0000 is a ref velocity
-
 	% find all the data files with Tstart<=t<=Tend
 	dataInd = (dataTime>=Tstart) & (dataTime<=Tend);
@@ -63,6 +58,4 @@
 	vx_out = zeros(numel(X)+1, numel(timeToload)); 
 	vy_out = zeros(numel(X)+1, numel(timeToload)); 
-	vx_err = zeros(numel(X)+1, numel(timeToload)); 
-	vy_err = zeros(numel(X)+1, numel(timeToload)); 
 	for i = 1:length(dataToLoad)
 
@@ -70,6 +63,4 @@
 		vx = (ncread(filename,'vx',[id1xh id1yh],[id2xh-id1xh+1 id2yh-id1yh+1],[1 1]));
 		vy = (ncread(filename,'vy',[id1xh id1yh],[id2xh-id1xh+1 id2yh-id1yh+1],[1 1]));
-		vxerr = (ncread(filename,'vx_err',[id1xh id1yh],[id2xh-id1xh+1 id2yh-id1yh+1],[1 1]));
-		vyerr = (ncread(filename,'vy_err',[id1xh id1yh],[id2xh-id1xh+1 id2yh-id1yh+1],[1 1]));
 
 		vx(vx<-32760) = nan;
@@ -79,8 +70,4 @@
 		vy_out(1:end-1,i) = InterpFromGrid(xh, yh, double(vy'), X, Y);
 		vy_out(end, i) = timeToload(i);
-		vx_err(1:end-1,i) = InterpFromGrid(xh, yh, double(vxerr'), X, Y);
-		vx_err(end, i) = timeToload(i);
-		vy_err(1:end-1,i) = InterpFromGrid(xh, yh, double(vyerr'), X, Y);
-		vy_err(end, i) = timeToload(i);
 	end
 end
Index: /issm/trunk-jpl/src/m/contrib/chenggong/modeldata/interpISMIP6Temp.m
===================================================================
--- /issm/trunk-jpl/src/m/contrib/chenggong/modeldata/interpISMIP6Temp.m	(revision 27965)
+++ /issm/trunk-jpl/src/m/contrib/chenggong/modeldata/interpISMIP6Temp.m	(revision 27966)
@@ -4,7 +4,4 @@
 %	X and Y are the coordinates of the mesh
 %
-%  Author: Cheng Gong
-%  Last modified: 2021-12-06
-
 filename = '/totten_1/ModelData/Greenland/ISMIP6/GreenlandISMIP6-Morlighem-2020-10-01.nc';
 
Index: /issm/trunk-jpl/src/m/contrib/chenggong/modeldata/interpMonthlyIceMaskGreene.m
===================================================================
--- /issm/trunk-jpl/src/m/contrib/chenggong/modeldata/interpMonthlyIceMaskGreene.m	(revision 27965)
+++ /issm/trunk-jpl/src/m/contrib/chenggong/modeldata/interpMonthlyIceMaskGreene.m	(revision 27966)
@@ -13,6 +13,4 @@
 %   - optional 5th input argument: path to the data file, by default it is the path on totten
 %
-%   Author: Cheng Gong
-%   Last modified: 2023-04-19
 
 % set icemask=-1 for the region with rocks
Index: /issm/trunk-jpl/src/m/contrib/chenggong/visualization/imageNonUni.m
===================================================================
--- /issm/trunk-jpl/src/m/contrib/chenggong/visualization/imageNonUni.m	(revision 27965)
+++ /issm/trunk-jpl/src/m/contrib/chenggong/visualization/imageNonUni.m	(revision 27966)
@@ -2,6 +2,4 @@
 %imageNonUni - draw an image with non uniform grid
 %
-% Author: Cheng Gong
-% Last modified: 2020-08-24
 x = varargin{1};
 y = varargin{2};
Index: /issm/trunk-jpl/src/m/contrib/chenggong/visualization/plotCompareTransientFlowline.m
===================================================================
--- /issm/trunk-jpl/src/m/contrib/chenggong/visualization/plotCompareTransientFlowline.m	(revision 27965)
+++ /issm/trunk-jpl/src/m/contrib/chenggong/visualization/plotCompareTransientFlowline.m	(revision 27966)
@@ -15,6 +15,4 @@
 %   yl:			ylim value
 %
-%   Author: Cheng Gong
-%   Date: 2021-12-06
 
 N = length(velList);