Contents

function [xGL,yGL,GLgeo]=PlotGroundingLines(CtrlVar,MUA,GF,GLgeo,xGL,yGL,varargin)

Plots grounding lines

[xGL,yGL,GLgeo]=PlotGroundingLines(CtrlVar,MUA,GF,GLgeo,xGL,yGL,varargin)

The only required inputs are:

MUA
GF

Other fields can be left empty. However, if the grounding line needs to be plotted repeatedly for same MUA and GF, entering GLgeo, xGL and yGL, obtained as outputs from a previouis call, will speed things up.

varargin is passed over the the plot function and can be any input accepted by the matlap plot function.

Examples:

Plot grounding lines with minimum of required input:

load('MUA-PIG-TWG-Example.mat','MUA','BCs','GF','CtrlVar')
Tarea=TriAreaFE(MUA.coordinates,MUA.connectivity); Tlength=sqrt(2*Tarea) ;
figure ;
[xGL,yGL,GLgeo]=PlotGroundingLines(CtrlVar,MUA,GF) ;

Plot grounding lines in red over the computational mesh in black:

load('MUA-PIG-TWG-Example.mat','MUA','BCs','GF','CtrlVar')
figure
CtrlVar.WhenPlottingMesh_PlotMeshBoundaryCoordinatesToo=0;
CtrlVar.PlotIndividualGLs=1 ;
PlotMuaMesh(CtrlVar,MUA,[],'k') ;
hold on
PlotGroundingLines(CtrlVar,MUA,GF,[],[],[],'color','r','LineWidth',2);

Plot grounding lines twice using outputs from first call in the second call:

  load('MUA-PIG-TWG-Example.mat','MUA','BCs','GF','CtrlVar')
  GLgeo=[] ; xGL=[], yGL=[];
  figure
  tic; [xGL,yGL,GLgeo]=PlotGroundingLines(CtrlVar,MUA,GF,GLgeo,xGL,yGL); toc
  tic; [xGL,yGL,GLgeo]=PlotGroundingLines(CtrlVar,MUA,GF,GLgeo,xGL,yGL); toc

To plot grounding lines for higher-order elements with sub-element resolution set CtrlVar.GLsubdivide=1. Doing so causes higher-order elements to be split into sub-elements allowing for more accurate representation of the grounding line within an element.

If xGL and yGL are empty on input, they are calculated from MUA and GLgeo.

If GLgeo, xGL and yGL are all empty on input, then GLgeo is first calculated from from GF and MUA, and then xGL and yGL are calculated from GLgeo.

If GLgeo, xGL and yGL are all empty, then MUA and GF can not be empty.

If xGL and yGL are not empty, then all other fields can be left empty.

Once xGL and yGL have been calculated, one can also simply plot the grounding lines as:

plot(xGL,yGL)

To plot individual grounding lines in different colours set

CtrlVar.PlotIndividualGLs=1

To just calculate the grounding line, but not plot it set

CtrlVar.PlotGLs=0;
[xGL,yGL,GLgeo]=PlotGroundingLines(CtrlVar,MUA,GF,GLgeo,xGL,yGL);

See also: PlotMuaBoundary, EleBasedGL

narginchk(3,inf)

if isempty(CtrlVar)
    CtrlVar.XYscale=1;
    CtrlVar.PlotIndividualGLs=0;
    CtrlVar.PlotGLs=1;
end

if ~isfield(CtrlVar,'PlotXYscale') ; CtrlVar.PlotXYscale=1 ; end
if ~isfield(CtrlVar,'PlotGLs') ; CtrlVar.PlotGLs=1 ; end
if ~isfield(CtrlVar,'PlotIndividualGLs') ; CtrlVar.PlotIndividualGLs=0 ; end


if nargin<4 || isempty(GLgeo)

    GLgeo=GLgeometry(MUA.connectivity,MUA.coordinates,GF,CtrlVar);

end

if nargin<6 || ( isempty(xGL) || isempty(yGL))
    xa=GLgeo(:,3) ;  xb=GLgeo(:,4) ; ya=GLgeo(:,5) ;  yb=GLgeo(:,6) ;
    [xGL,yGL]=LineUpEdges2([],xa,xb,ya,yb);
end



if CtrlVar.PlotGLs

    if ~CtrlVar.PlotIndividualGLs

        plot(xGL/CtrlVar.PlotXYscale,yGL/CtrlVar.PlotXYscale,varargin{:}) ;
        ax=gca; ax.DataAspectRatio=[1 1 1];

    else

        i=0;
        I=find(isnan(xGL)) ;
        I=[1;I(:)];
        col=['b','r','c','g','k','m'];
        %col=['b','r'];
        for ii=1:numel(I)-1
            i=i+1;

            plot(xGL(I(ii):I(ii+1))/CtrlVar.PlotXYscale,yGL(I(ii):I(ii+1))/CtrlVar.PlotXYscale,'color',col(i),varargin{:}) ;
            axis equal ; hold on ;
            if i==numel(col) ; i=0 ; end
        end
        ax=gca; ax.DataAspectRatio=[1 1 1];
    end
end

end

Published with MATLAB® R2017a