Changeset 16287
- Timestamp:
- 10/02/13 20:30:24 (11 years ago)
- Location:
- issm/trunk-jpl/src/m/classes
- Files:
-
- 2 added
- 1 deleted
- 2 copied
- 1 moved
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/m/classes/spheremesh.m
r16222 r16287 1 % PLANETMESH class definition1 %SPHEREMESH class definition 2 2 % 3 3 % Usage: 4 % planetmesh=planetmesh();4 % spheremesh=spheremesh(); 5 5 6 classdef planetmesh6 classdef spheremesh 7 7 properties (SetAccess=public) 8 8 x = NaN; … … 23 23 end 24 24 methods 25 function obj = planetmesh(varargin) % {{{25 function obj = spheremesh(varargin) % {{{ 26 26 switch nargin 27 27 case 0 … … 42 42 function md = checkconsistency(obj,md,solution,analyses) % {{{ 43 43 44 md = checkfield(md,' planetmesh.x','NaN',1,'size',[md.planetmesh.numberofvertices 1]);45 md = checkfield(md,' planetmesh.y','NaN',1,'size',[md.planetmesh.numberofvertices 1]);46 md = checkfield(md,' planetmesh.z','NaN',1,'size',[md.planetmesh.numberofvertices 1]);47 md = checkfield(md,' planetmesh.r','NaN',1,'size',[md.planetmesh.numberofvertices 1]);48 md = checkfield(md,' planetmesh.theta','NaN',1,'size',[md.planetmesh.numberofvertices 1]);49 md = checkfield(md,' planetmesh.phi','NaN',1,'size',[md.planetmesh.numberofvertices 1]);50 md = checkfield(md,' planetmesh.elements','NaN',1,'>',0,'values',1:md.planetmesh.numberofvertices);51 if(md. planetmesh.dimension==2),52 md = checkfield(md,' planetmesh.elements','size',[md.planetmesh.numberofelements 3]);44 md = checkfield(md,'spheremesh.x','NaN',1,'size',[md.spheremesh.numberofvertices 1]); 45 md = checkfield(md,'spheremesh.y','NaN',1,'size',[md.spheremesh.numberofvertices 1]); 46 md = checkfield(md,'spheremesh.z','NaN',1,'size',[md.spheremesh.numberofvertices 1]); 47 md = checkfield(md,'spheremesh.r','NaN',1,'size',[md.spheremesh.numberofvertices 1]); 48 md = checkfield(md,'spheremesh.theta','NaN',1,'size',[md.spheremesh.numberofvertices 1]); 49 md = checkfield(md,'spheremesh.phi','NaN',1,'size',[md.spheremesh.numberofvertices 1]); 50 md = checkfield(md,'spheremesh.elements','NaN',1,'>',0,'values',1:md.spheremesh.numberofvertices); 51 if(md.spheremesh.dimension==2), 52 md = checkfield(md,'spheremesh.elements','size',[md.spheremesh.numberofelements 3]); 53 53 else 54 md = checkfield(md,' planetmesh.elements','size',[md.planetmesh.numberofelements 6]);54 md = checkfield(md,'spheremesh.elements','size',[md.spheremesh.numberofelements 6]); 55 55 end 56 if any(~ismember(1:md. planetmesh.numberofvertices,sort(unique(md.planetmesh.elements(:)))));57 md = checkmessage(md,'orphan nodes have been found. Check the planetmesh outline');56 if any(~ismember(1:md.spheremesh.numberofvertices,sort(unique(md.spheremesh.elements(:))))); 57 md = checkmessage(md,'orphan nodes have been found. Check the spheremesh outline'); 58 58 end 59 md = checkfield(md,' planetmesh.dimension','values',[2 3]);60 md = checkfield(md,' planetmesh.numberoflayers','>=',0);61 md = checkfield(md,' planetmesh.numberofelements','>',0);62 md = checkfield(md,' planetmesh.numberofvertices','>',0);63 if (md. planetmesh.dimension==2),64 md = checkfield(md,' planetmesh.average_vertex_connectivity','>=',9,'message','''planetmesh.average_vertex_connectivity'' should be at least 9 in 2d');59 md = checkfield(md,'spheremesh.dimension','values',[2 3]); 60 md = checkfield(md,'spheremesh.numberoflayers','>=',0); 61 md = checkfield(md,'spheremesh.numberofelements','>',0); 62 md = checkfield(md,'spheremesh.numberofvertices','>',0); 63 if (md.spheremesh.dimension==2), 64 md = checkfield(md,'spheremesh.average_vertex_connectivity','>=',9,'message','''spheremesh.average_vertex_connectivity'' should be at least 9 in 2d'); 65 65 else 66 md = checkfield(md,' planetmesh.average_vertex_connectivity','>=',24,'message','''planetmesh.average_vertex_connectivity'' should be at least 24 in 3d');66 md = checkfield(md,'spheremesh.average_vertex_connectivity','>=',24,'message','''spheremesh.average_vertex_connectivity'' should be at least 24 in 3d'); 67 67 end 68 md = checkfield(md,' planetmesh.elementconnectivity','size',[md.planetmesh.numberofelements 3],'NaN',1);68 md = checkfield(md,'spheremesh.elementconnectivity','size',[md.spheremesh.numberofelements 3],'NaN',1); 69 69 70 70 %Solution specific checks … … 72 72 case MasstransportSolutionEnum(), 73 73 if md.masstransport.stabilization==3, 74 md = checkfield(md,' planetmesh.dimension','values',2,'message','Discontinuous Galerkin only supported for 2d planetmeshes');74 md = checkfield(md,'spheremesh.dimension','values',2,'message','Discontinuous Galerkin only supported for 2d spheremeshes'); 75 75 end 76 76 case TransientSolutionEnum(), 77 77 if md.transient.ismasstransport & md.masstransport.stabilization==3, 78 md = checkfield(md,' planetmesh.dimension','values',2,'message','Discontinuous Galerkin only supported for 2d planetmeshes');78 md = checkfield(md,'spheremesh.dimension','values',2,'message','Discontinuous Galerkin only supported for 2d spheremeshes'); 79 79 end 80 80 case ThermalSolutionEnum(), 81 md = checkfield(md,' planetmesh.dimension','values',3,'message','thermal solution only supported on 3d planetmeshes');81 md = checkfield(md,'spheremesh.dimension','values',3,'message','thermal solution only supported on 3d spheremeshes'); 82 82 end 83 83 end % }}} … … 97 97 98 98 disp(sprintf('\n Properties:')); 99 fielddisplay(obj,'dimension',' planetmesh dimension (2d or 3d)');99 fielddisplay(obj,'dimension','spheremesh dimension (2d or 3d)'); 100 100 fielddisplay(obj,'numberoflayers','number of extrusion layers'); 101 101
Note:
See TracChangeset
for help on using the changeset viewer.