Changeset 2396


Ignore:
Timestamp:
10/08/09 09:55:53 (15 years ago)
Author:
Mathieu Morlighem
Message:

use pairoptions in meshyams

Location:
issm/trunk/src/m/classes/public/mesh
Files:
2 deleted
1 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk/src/m/classes/public/mesh/meshyams.m

    r2296 r2396  
    2222
    2323%recover options
    24 yamsoptions=optionlist2cell(varargin{:});
    25 
    26 %add default options
    27 yamsoptions=process_mesh_options(yamsoptions);
     24options=pairoptions(varargin{:});
     25options=deleteduplicates(options,1);
     26if ~exist(options,'domainoutline'),
     27        error('meshyams error message: no ''domainoutline'' was provided');
     28end
     29if ~exist(options,'velocities'),
     30        error('meshyams error message: no ''velocities'' was provided');
     31end
    2832
    2933%recover some fields
    3034disp('MeshYams Options:')
    31 domainoutline=yamsoptions.domainoutline;   disp(sprintf('   %-15s: ''%s''','DomainOutline',domainoutline));
    32 groundeddomain=yamsoptions.groundeddomain; disp(sprintf('   %-15s: ''%s''','GroundedDomain',groundeddomain));
    33 velocities=yamsoptions.velocities;         disp(sprintf('   %-15s: ''%s''','Velocities',velocities));
    34 resolution=yamsoptions.resolution;         disp(sprintf('   %-15s: %f','Resolution',resolution));
    35 gradation=yamsoptions.gradation;
    36 nsteps=yamsoptions.nsteps;                 disp(sprintf('   %-15s: %i','nsteps',nsteps));
    37 epsilon=yamsoptions.epsilon;               disp(sprintf('   %-15s: %f','epsilon',epsilon));
    38 hmin=yamsoptions.hmin;                     disp(sprintf('   %-15s: %f','hmin',hmin));
    39 hmax=yamsoptions.hmax;                     disp(sprintf('   %-15s: %f\n','hmax',hmax));
     35domainoutline=getfieldvalue(options,'domainoutline');
     36disp(sprintf('   %-15s: ''%s''','DomainOutline',domainoutline));
     37groundeddomain=getfieldvalue(options,'groundeddomain','N/A');
     38disp(sprintf('   %-15s: ''%s''','GroundedDomain',groundeddomain));
     39velocities=getfieldvalue(options,'velocities');
     40disp(sprintf('   %-15s: ''%s''','Velocities',velocities));
     41resolution=getfieldvalue(options,'resolution',5000);
     42disp(sprintf('   %-15s: %f','Resolution',resolution));
     43nsteps=getfieldvalue(options,'nsteps',6);
     44disp(sprintf('   %-15s: %i','nsteps',nsteps));
     45gradation=getfieldvalue(options,'gradation',2*ones(nsteps,1));
     46disp(sprintf('   %-15s: %g','gradation',gradation(1)));
     47epsilon=getfieldvalue(options,'epsilon',3);
     48disp(sprintf('   %-15s: %f','epsilon',epsilon));
     49hmin=getfieldvalue(options,'hmin',500);
     50disp(sprintf('   %-15s: %f','hmin',hmin));
     51hmax=getfieldvalue(options,'hmax',150*10^3);
     52disp(sprintf('   %-15s: %f\n','hmax',hmax));
    4053
    4154%mesh with initial resolution
Note: See TracChangeset for help on using the changeset viewer.