Changeset 11249

Timestamp:

01/30/12 09:36:04 (13 years ago)

Author:

seroussi

Message:

appropriate name of plot is plot_icefront not plot_pressureload

Location:

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

Files:

• plot_BC.m (modified) (1 diff)
• plot_icefront.m (moved) (moved from issm/trunk-jpl/src/m/model/plot/plot_pressureload.m ) (7 diffs)
• plot_manager.m (modified) (1 diff)
• plotdoc.m (modified) (1 diff)

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

@@ -2 +2 @@
 
 %plot neuman
-plot_pressureload(md,options,width,i,data)
+plot_icefront(md,options,width,i,data)
 
 hold on

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

@@ -1 @@
-function plot_pressureload(md,options,width,i,data);
-%PLOT_PRESSURELOAD - plot segment on neumann BC
+function plot_icefront(md,options,width,i,data);
+%PLOT_ICEFRONT - plot segment on neumann BC
 %
 %   Usage:
-%      plot_pressureload(md,options,width,i);
+%      plot_icefront(md,options,width,i);
 %
 %   See also: PLOTMODEL
@@ -12 +12 @@
 %process mesh and data
 [x y z elements is2d isplanet]=processmesh(md,[],options);
-pressureload=md.diagnostic.icefront;
+icefront=md.diagnostic.icefront;
 
 if (md.mesh.dimension==2),
@@ -22 +22 @@
 
         %highlight elements on neumann
-        pos=find(pressureload(:,end)==1);
-        pos=pressureload(pos,end-1);
+        pos=find(icefront(:,end)==1);
+        pos=icefront(pos,end-1);
         A=elements(pos,1); B=elements(pos,2); C=elements(pos,3); 
         h2=patch( 'Faces', [A B C], 'Vertices', [x y z],'FaceVertexCData', [1 1 1],'FaceColor','blue','EdgeColor','black');
-        pos=find(pressureload(:,end)==0);
-        pos=pressureload(pos,end-1);
+        pos=find(icefront(:,end)==0);
+        pos=icefront(pos,end-1);
         A=elements(pos,1); B=elements(pos,2); C=elements(pos,3); 
         h3=patch( 'Faces', [A B C], 'Vertices', [x y z],'FaceVertexCData', [1 1 1],'FaceColor','green','EdgeColor','black');
 
-        pos=find(pressureload(:,end)==2);
-        pos=pressureload(pos,end-1);
+        pos=find(icefront(:,end)==2);
+        pos=icefront(pos,end-1);
         A=elements(pos,1); B=elements(pos,2); C=elements(pos,3); 
         h3bis=patch( 'Faces', [A B C], 'Vertices', [x y z],'FaceVertexCData', [1 1 1],'FaceColor','cyan','EdgeColor','black');
 
         %display arrows pointing outward
-        xstart=mean(x(pressureload(:,1:end-2)),2);
-        ystart=mean(y(pressureload(:,1:end-2)),2);
-        length=sqrt((x(pressureload(:,1))-x(pressureload(:,2))).^2 + (y(pressureload(:,1))-y(pressureload(:,2))).^2 );
-        normal(:,1)=cos(atan2((x(pressureload(:,1))-x(pressureload(:,2))) , (y(pressureload(:,2))-y(pressureload(:,1)))));
-        normal(:,2)=sin(atan2((x(pressureload(:,1))-x(pressureload(:,2))) , (y(pressureload(:,2))-y(pressureload(:,1)))));
+        xstart=mean(x(icefront(:,1:end-2)),2);
+        ystart=mean(y(icefront(:,1:end-2)),2);
+        length=sqrt((x(icefront(:,1))-x(icefront(:,2))).^2 + (y(icefront(:,1))-y(icefront(:,2))).^2 );
+        normal(:,1)=cos(atan2((x(icefront(:,1))-x(icefront(:,2))) , (y(icefront(:,2))-y(icefront(:,1)))));
+        normal(:,2)=sin(atan2((x(icefront(:,1))-x(icefront(:,2))) , (y(icefront(:,2))-y(icefront(:,1)))));
         xend=xstart+length.*normal(:,1);
         yend=ystart+length.*normal(:,2);
@@ -58 +58 @@
 
         %highlight elements on neumann
-        pos=find(pressureload(:,end)==1);
-        pos=pressureload(pos,end-1);
+        pos=find(icefront(:,end)==1);
+        pos=icefront(pos,end-1);
         A=elements(pos,1); B=elements(pos,2); C=elements(pos,3); D=elements(pos,4); E=elements(pos,5); F=elements(pos,6);
         h2=patch( 'Faces', [A B C],'Vertices', [x y z],'FaceVertexCData', [1 1 1],'FaceColor','blue','EdgeColor','black');
@@ -66 +66 @@
         patch( 'Faces', [B E F C],'Vertices', [x y z],'FaceVertexCData', [1 1 1],'FaceColor','blue','EdgeColor','black');
         patch( 'Faces', [C A D F],'Vertices', [x y z],'FaceVertexCData', [1 1 1],'FaceColor','blue','EdgeColor','black');
-        pos=find(pressureload(:,end)==0);
-        pos=pressureload(pos,end-1);
+        pos=find(icefront(:,end)==0);
+        pos=icefront(pos,end-1);
         A=elements(pos,1); B=elements(pos,2); C=elements(pos,3); D=elements(pos,4); E=elements(pos,5); F=elements(pos,6);
         h3=patch( 'Faces', [A B C],'Vertices', [x y z],'FaceVertexCData', [1 1 1],'FaceColor','green','EdgeColor','black');
@@ -74 +74 @@
         patch( 'Faces', [B E F C],'Vertices', [x y z],'FaceVertexCData', [1 1 1],'FaceColor','green','EdgeColor','black');
         patch( 'Faces', [C A D F],'Vertices', [x y z],'FaceVertexCData', [1 1 1],'FaceColor','green','EdgeColor','black');
-        pos=find(pressureload(:,end)==2);
-        pos=pressureload(pos,end-1);
+        pos=find(icefront(:,end)==2);
+        pos=icefront(pos,end-1);
         A=elements(pos,1); B=elements(pos,2); C=elements(pos,3); D=elements(pos,4); E=elements(pos,5); F=elements(pos,6);
         h3bis=patch( 'Faces', [A B C],'Vertices', [x y z],'FaceVertexCData', [1 1 1],'FaceColor','cyan','EdgeColor','black');
@@ -84 +84 @@
 
         %display arrows pointing outward
-        xstart=mean(x(pressureload(:,1:end-2)),2);
-        ystart=mean(y(pressureload(:,1:end-2)),2);
-        zstart=mean(z(pressureload(:,1:end-2)),2);
-        length=sqrt((x(pressureload(:,1))-x(pressureload(:,2))).^2 + (y(pressureload(:,1))-y(pressureload(:,2))).^2 );
-        normal(:,1)=cos(atan2((x(pressureload(:,1))-x(pressureload(:,2))) , (y(pressureload(:,2))-y(pressureload(:,1)))));
-        normal(:,2)=sin(atan2((x(pressureload(:,1))-x(pressureload(:,2))) , (y(pressureload(:,2))-y(pressureload(:,1)))));
+        xstart=mean(x(icefront(:,1:end-2)),2);
+        ystart=mean(y(icefront(:,1:end-2)),2);
+        zstart=mean(z(icefront(:,1:end-2)),2);
+        length=sqrt((x(icefront(:,1))-x(icefront(:,2))).^2 + (y(icefront(:,1))-y(icefront(:,2))).^2 );
+        normal(:,1)=cos(atan2((x(icefront(:,1))-x(icefront(:,2))) , (y(icefront(:,2))-y(icefront(:,1)))));
+        normal(:,2)=sin(atan2((x(icefront(:,1))-x(icefront(:,2))) , (y(icefront(:,2))-y(icefront(:,1)))));
         xend=xstart+length.*normal(:,1);
         yend=ystart+length.*normal(:,2);

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

@@ -112 +112 @@
                         plot_sarpwr(md,options,subplotwidth,i)
                         return
-                case 'pressureload'
-                        plot_pressureload(md,options,subplotwidth,i,data)
+                case 'icefront'
+                        plot_icefront(md,options,subplotwidth,i,data)
                         return
                 case 'segments'

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

@@ -15 +15 @@
 disp('                  - ''basal_dragx'' or ''basal_dragy'' : plot a component of the basal drag on the bed (in kPa)');
 disp('                  - ''boundaries'': this will draw all the segment boundaries to the model, including rifts.');
+disp('                  - ''icefront'': this will show segments that are used to define the icefront of the model (Neumann boundary conditions).');
 disp('                  - ''BC'': this will draw all the boundary conditions (Dirichlet and Neumann).');
 disp('                  - ''deviatoricstress_tensor'': plot the components of the deviatoric stress tensor (tauxx,tauyy,tauzz,tauxy,tauxz,tauyz) if computed');