Changeset 12375

Timestamp:

06/06/12 14:18:43 (13 years ago)

Author:

Mathieu Morlighem

Message:

Fixed plot_profile

Location:

issm/trunk-jpl/src/m/model

Files:

• ProfileValues.m (modified) (1 diff)
• plot/plot_profile.m (modified) (4 diffs)

issm/trunk-jpl/src/m/model/ProfileValues.m

@@ -15 +15 @@
 
 %Some useful parameters
-layers=ceil(mean(md.geometry.thickness)/res_v);
-Z=bed:resolution:surface;
+layers=ceil(mean(md.geometry.thickness)/resolution);
+Z=(bed:resolution:surface)';
 X=xprof*ones(size(Z));
 Y=yprof*ones(size(Z));

issm/trunk-jpl/src/m/model/plot/plot_profile.m

@@ -9 +9 @@
 %process model
 [x_m y_m z_m elements_m is2d isplanet]=processmesh(md,[],options);
-if ~is2d, error('only 3d model supported'); end
+if is2d, error('only 3d model supported'); end
 
 %Get number of curves and generate random colors
@@ -16 +16 @@
 color=eval([ colorm '(numcurves);']);
 options=removefield(options,'colormap',0); %back to default colormap
+
+%Get coordinates
+location=getfieldvalue(options,'profile');
+if ~isnumeric(location) | numel(location)~=2,
+        error('location provided not supported (should be [x y])');
+end
+xprof=location(1);
+yprof=location(2);
 
 %Loop over number of curves
@@ -32 +40 @@
 
         %Compute profile value
-        [z,data_interp]=ProfileValues(md,data,xprof,yprof,resolution)
+        [z,data_interp]=ProfileValues(md,data,xprof,yprof,resolution);
 
         %plot profile
         subplot(nlines,ncols,i)
-        A=elements(:,1); B=elements(:,2); C=elements(:,3);  D=elements(:,4); 
-        plot(data_interp,z,'color',color(i,:),'LineWidth',getfieldvalue(options,'linewidth',1));
+        plot(data_interp,z,'color',color(i,:),'LineWidth',getfieldvalue(options,'linewidth',1),'LineStyle','-');
         hold on;
 end
@@ -45 +52 @@
 options=addfielddefault(options,'colorbar',0);
 options=addfielddefault(options,'ylabel','z');
+options=addfielddefault(options,'view',2);
 applyoptions(md,[],options);