Index: /issm/trunk-jpl/src/m/interp/PointValues.m
===================================================================
--- /issm/trunk-jpl/src/m/interp/PointValues.m	(revision 17718)
+++ /issm/trunk-jpl/src/m/interp/PointValues.m	(revision 17718)
@@ -0,0 +1,14 @@
+function [data_interp]=PointValues(md,data,xpoint,ypoint)
+%POINTVALUES - compute the value of a field on a single point
+%
+%   This routine gets the value of a given field of the model on points
+%
+%   Usage:
+%      [z,data]=PointValues(md,data,X,Y,Z)
+
+%Get bed and surface for each 2d point, offset to make sure that it is inside the glacier system
+%offset=10^-3;
+%bed=InterpFromMeshToMesh2d(md.mesh.elements2d,md.mesh.x2d,md.mesh.y2d,project2d(md,md.geometry.bed,1),xpoint,ypoint)+offset;
+%surface=InterpFromMeshToMesh2d(md.mesh.elements2d,md.mesh.x2d,md.mesh.y2d,project2d(md,md.geometry.surface,1),xpoint,ypoint)-offset;
+
+data_interp=InterpFromMeshToMesh2d(md.mesh.elements,md.mesh.x,md.mesh.y,data,xpoint,ypoint);
Index: /issm/trunk-jpl/src/m/plot/plot_manager.m
===================================================================
--- /issm/trunk-jpl/src/m/plot/plot_manager.m	(revision 17717)
+++ /issm/trunk-jpl/src/m/plot/plot_manager.m	(revision 17718)
@@ -112,4 +112,7 @@
 			plot_sarpwr(md,options,subplotwidth,i)
 			return
+		case 'time_dependant' ,
+			plot_vstime(md,options,nlines,ncols,i)
+			return
 		case 'icefront'
 			plot_icefront(md,options,subplotwidth,i,data)
Index: /issm/trunk-jpl/src/m/plot/plot_profile.m
===================================================================
--- /issm/trunk-jpl/src/m/plot/plot_profile.m	(revision 17717)
+++ /issm/trunk-jpl/src/m/plot/plot_profile.m	(revision 17718)
@@ -31,5 +31,5 @@
 	[datai datatype]=processdata(md,data(:,i),options);
 
-	%resolution
+	%resolution[z,data_interp]=ProfileValues(md,datai,xprof,yprof,resolution);
 	if exist(options,'resolution'),
 		resolution=getfieldvalue(options,'resolution');
Index: /issm/trunk-jpl/src/m/plot/plot_vstime.m
===================================================================
--- /issm/trunk-jpl/src/m/plot/plot_vstime.m	(revision 17718)
+++ /issm/trunk-jpl/src/m/plot/plot_vstime.m	(revision 17718)
@@ -0,0 +1,55 @@
+function plot_vstime(md,options,nlines,ncols,i)
+
+%plot mesh
+subplot(nlines,ncols,i); 
+
+%getting the variable structure
+datastruct  = getfieldvalue(options,'Input');
+structnames = strsplit(datastruct,'.');
+
+%getting the position of the to be ploted point
+location=getfieldvalue(options,'position');
+if ~isnumeric(location) | numel(location)~=2,
+	error('location provided not supported (should be [x y])');
+end
+xpoint=location(1);
+ypoint=location(2);
+
+%gathering a few time related parameters
+FirstTime=md.timestepping.start_time;
+LastTime=md.timestepping.final_time;
+Dt=md.timestepping.time_step;
+OutputDt=md.settings.output_frequency*Dt;
+%Constructiion of the result structure(one step above variable in
+%the structure tree)
+solstruct = md;
+for i=2:numel(structnames)-1,
+	solstruct=solstruct.(char(structnames(i)));
+end
+
+timesteps   = numel(solstruct);
+%now build a table with the variable needed in the end plot
+Value(1:timesteps)=NaN*ones(timesteps,1);
+time(1:timesteps)=NaN*ones(timesteps,1);
+
+for t=1:timesteps,
+	
+	data=solstruct(t).(char(structnames(numel(structnames))));
+
+	[TimeData datatype]=processdata(md,data,options);
+	clear data
+	Value(t)=PointValues(md,TimeData,xpoint,ypoint);
+	time(t)=FirstTime+(t-1)*OutputDt+Dt;
+end
+
+plot(time,Value);
+clear Value time
+
+options=addfielddefault(options,'title',md.miscellaneous.name);
+options=addfielddefault(options,'ylabel',char(structnames(numel(structnames))));
+options=addfielddefault(options,'xlabel','Time');
+options=addfielddefault(options,'colorbar',0);
+options=addfielddefault(options,'axis','auto');
+options=addfielddefault(options,'xlim',[FirstTime LastTime]);
+options=addfielddefault(options,'view',2);
+applyoptions(md,[],options);