Changeset 16325


Ignore:
Timestamp:
10/08/13 08:08:56 (11 years ago)
Author:
Mathieu Morlighem
Message:

NEW: plotting capabilities for mesh2d vertical

Location:
issm/trunk-jpl/src/m
Files:
1 added
3 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk-jpl/src/m/plot/applyoptions.m

    r16051 r16325  
    4949
    5050%view
    51 if md.mesh.dimension==3 & ~exist(options,'layer'),
     51if strcmp(meshtype(md.mesh),'3D') & ~exist(options,'layer'),
    5252        view(getfieldvalue(options,'view',3));
    5353else
     
    6060        eval(['axis ' getfieldvalue(options,'axis')]);
    6161else
    62         if ((md.mesh.dimension==2) | exist(options,'layer')),
     62        if strcmpi(meshtype(md.mesh),'2Dhorizontal') | exist(options,'layer'),
    6363                axis tight equal;
    6464        else
  • issm/trunk-jpl/src/m/plot/processdata.m

    r13730 r16325  
    1515if (iscell(data) | isempty(data) | length(data)==0 | (length(data)==1 & ~isstruct(data) & isnan(data))),
    1616        error('plotmodel error message: data provided is empty');
     17end
     18
     19%Process NaN if any (do not know before mask is applied)
     20if exist(options,'nan')
     21        data(find(isnan(data)))=getfieldvalue(options,'nan',0);
     22end
     23
     24%special case for mesg 2dvertical
     25if strcmp(meshtype(md.mesh),'2Dvertical'),
     26        [data datatype] = processdata(md.mesh,md,data,options);
     27        return;
    1728end
    1829
     
    5162datasize=size(data);
    5263
    53 %Process NaN if any (do not know before mask is applied)
    54 if exist(options,'nan')
    55         data(find(isnan(data)))=getfieldvalue(options,'nan',0);
    56 end
    5764%non patch processing
    5865if datatype~=4 & datatype~=5,
  • issm/trunk-jpl/src/m/plot/processmesh.m

    r13730 r16325  
    1313if md.mesh.numberofvertices==md.mesh.numberofelements
    1414        error(['plot error message: the number of elements is the same as the number of nodes...']);
     15end
     16
     17%special case for mesg 2dvertical
     18if strcmp(meshtype(md.mesh),'2Dvertical'),
     19        [x y z elements is2d isplanet] = processmesh(md.mesh,options);
     20        return;
    1521end
    1622
Note: See TracChangeset for help on using the changeset viewer.