Changeset 17720
- Timestamp:
- 04/14/14 14:10:49 (11 years ago)
- Location:
- issm/trunk-jpl/src/m/classes
- Files:
-
- 32 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/m/classes/SMB.m
r17079 r17720 9 9 end 10 10 methods 11 function createxml(obj,fid) % {{{ 12 fprintf(fid, '\n\n'); 13 fprintf(fid, '%s\n', '<!-- surfaceforcings(SMB) -->'); 14 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n','<parameter key ="mass_balance" type="',class(obj.mass_balance),'" default="',obj.mass_balance,'">',' <section name="surfaceforcings(SMB)" />',' <help> surface mass balance [m/yr ice eq] </help>','</parameter>'); 15 16 end % }}} 11 17 function obj = SMB(varargin) % {{{ 12 18 switch nargin -
issm/trunk-jpl/src/m/classes/autodiff.m
r16764 r17720 20 20 %}}} 21 21 methods 22 function createxml(obj,fid) % {{{ 23 fprintf(fid, '<!-- autodiff -->\n'); 24 25 % automatic differentiation parameters 26 fprintf(fid,'%s\n%s\n%s\n','<frame key="1" label="automatic differentiation parameters">','<section name="autodiff" />'); 27 fprintf(fid,'%s%s%s\n%s\n%s%s%s\n%s\n%s\n',' <parameter key ="isautodiff" type="',class(obj.isautodiff),'" optional="false">',' <section name="autodiff" />',' <option value="',convert2str(obj.isautodiff),'" type="string" default="true"></option>',' <help> indicates if the automatic differentiation is activated </help>',' </parameter>'); 28 29 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="dependents" type="',class(obj.dependents),'" default="',convert2str(obj.dependents),'">',' <section name="autodiff" />',' <help> list of dependent variables </help>',' </parameter>'); 30 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="independents" type="',class(obj.independents),'" default="',convert2str(obj.independents),'">',' <section name="autodiff" />',' <help> list of independent variables </help>',' </parameter>'); 31 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="driver" type="',class(obj.driver),'" default="',convert2str(obj.driver),'">',' <section name="autodiff" />',' <help> ADOLC driver (''fos_forward'' or ''fov_forward'') </help>',' </parameter>'); 32 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="obufsize" type="',class(obj.obufsize),'" default="',convert2str(obj.obufsize),'">',' <section name="autodiff" />',' <help> Number of operations per buffer (==OBUFSIZE in usrparms.h) </help>',' </parameter>'); 33 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="lbufsize" type="',class(obj.lbufsize),'" default="',convert2str(obj.lbufsize),'">',' <section name="autodiff" />',' <help> Number of locations per buffer (==LBUFSIZE in usrparms.h) </help>',' </parameter>'); 34 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="cbufsize" type="',class(obj.cbufsize),'" default="',convert2str(obj.cbufsize),'">',' <section name="autodiff" />',' <help> Number of values per buffer (==CBUFSIZE in usrparms.h) </help>',' </parameter>'); 35 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="tbufsize" type="',class(obj.tbufsize),'" default="',convert2str(obj.tbufsize),'">',' <section name="autodiff" />',' <help> Number of taylors per buffer (&lt;=TBUFSIZE in usrparms.h) </help>',' </parameter>'); 36 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="gcTriggerRatio" type="',class(obj.gcTriggerRatio),'" default="',convert2str(obj.gcTriggerRatio),'">',' <section name="autodiff" />',' <help> free location block sorting/consolidation triggered if the ratio between allocated and used locations exceeds gcTriggerRatio </help>',' </parameter>'); 37 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="gcTriggerRatio" type="',class(obj.gcTriggerRatio),'" default="',convert2str(obj.gcTriggerRatio),'">',' <section name="autodiff" />',' <help> free location block sorting/consolidation triggered if the allocated locations exceed gcTriggerMaxSize </help>',' </parameter>'); 38 39 fprintf(fid,'%s\n%s\n','</frame>'); 40 41 end % }}} 22 42 function obj = autodiff(varargin) % {{{ 23 43 switch nargin -
issm/trunk-jpl/src/m/classes/balancethickness.m
r16764 r17720 11 11 end 12 12 methods 13 function createxml(obj,fid) % {{{ 14 fprintf(fid, '\n\n'); 15 fprintf(fid, '%s\n', '<!-- balancethickness -->'); 16 17 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n', '<parameter key ="spcthickness" type="', class(obj.spcthickness),'" default="', convert2str(obj.spcthickness),'">', ' <section name="balancethickness" />',' <help> thickness constraints (NaN means no constraint) [m] </help>','</parameter>'); 18 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n', '<parameter key ="thickening_rate" type="', class(obj.thickening_rate),'" default="', convert2str(obj.thickening_rate),'">', ' <section name="balancethickness" />',' <help> ice thickening rate used in the mass conservation (dh/dt) [m/yr] </help>','</parameter>'); 19 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n', '<parameter key ="stabilization" type="', class(obj.stabilization),'" default="', convert2str(obj.stabilization),'">', ' <section name="balancethickness" />',' <help> 0: None, 1: SU, 2: SSAs artificial diffusivity, 3:DG </help>','</parameter>'); 20 21 end % }}} 13 22 function obj = balancethickness(varargin) % {{{ 14 23 switch nargin -
issm/trunk-jpl/src/m/classes/basalforcings.m
r17071 r17720 11 11 end 12 12 methods 13 function createxml(obj,fid) % {{{ 14 fprintf(fid, '\n\n'); 15 fprintf(fid, '%s\n', '<!-- basalforcings -->'); 16 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n', '<parameter key ="melting_rate" type="', class(obj.melting_rate),'" default="', num2str(obj.melting_rate),'">', ' <section name="basalforcings" />',' <help> basal melting rate (positive if melting) [m/yr] </help>','</parameter>'); 17 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n', '<parameter key ="melting_rate_correction" type="', class(obj.melting_rate_correction),'" default="', num2str(obj.melting_rate_correction),'">', ' <section name="basalforcings" />',' <help> additional melting applied to compensate for dh/dt [m/yr] </help>','</parameter>'); 18 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n', '<parameter key ="geothermalflux" type="', class(obj.geothermalflux),'" default="', num2str(obj.geothermalflux),'">', ' <section name="basalforcings" />',' <help> geothermal heat flux [W/m^2] </help>','</parameter>'); 19 20 end % }}} 13 21 function obj = basalforcings(varargin) % {{{ 14 22 switch nargin -
issm/trunk-jpl/src/m/classes/clusters/generic.m
r17052 r17720 23 23 end 24 24 methods 25 function createxml(obj,fid) % {{{ 26 fprintf(fid, '\n\n'); 27 fprintf(fid, '%s\n', '<!-- generic -->'); 28 29 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n', '<parameter key ="name" type="',class(obj.name),'" default="',convert2str(obj.name),'">', ' <section name="cluster" />',' <help> N/A </help>','</parameter>'); 30 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n', '<parameter key ="login" type="',class(obj.login),'" default="',convert2str(obj.login),'">', ' <section name="cluster" />',' <help> N/A </help>','</parameter>'); 31 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n', '<parameter key ="np" type="',class(obj.np),'" default="',convert2str(obj.np),'">', ' <section name="cluster" />',' <help> N/A </help>','</parameter>'); 32 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n', '<parameter key ="port" type="',class(obj.port),'" default="',convert2str(obj.port),'">', ' <section name="cluster" />',' <help> N/A </help>','</parameter>'); 33 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n', '<parameter key ="codepath" type="',class(obj.codepath),'" default="',convert2str(obj.codepath),'">', ' <section name="cluster" />',' <help> N/A </help>','</parameter>'); 34 35 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n', '<parameter key ="executionpath" type="',class(obj.executionpath),'" default="',convert2str(obj.executionpath),'">', ' <section name="cluster" />',' <help> N/A </help>','</parameter>'); 36 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n', '<parameter key ="etcpath" type="',class(obj.etcpath),'" default="',convert2str(obj.etcpath),'">', ' <section name="cluster" />',' <help> N/A </help>','</parameter>'); 37 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n', '<parameter key ="valgrind" type="',class(obj.valgrind),'" default="',convert2str(obj.valgrind),'">', ' <section name="cluster" />',' <help> N/A </help>','</parameter>'); 38 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n', '<parameter key ="valgrindlib" type="',class(obj.valgrindlib),'" default="',convert2str(obj.valgrindlib),'">', ' <section name="cluster" />',' <help> N/A </help>','</parameter>'); 39 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n', '<parameter key ="valgrindsup" type="',class(obj.valgrindsup),'" default="',convert2str(obj.valgrindsup),'">', ' <section name="cluster" />',' <help> N/A </help>','</parameter>'); 40 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n', '<parameter key ="verbose" type="',class(obj.verbose),'" default="',convert2str(obj.verbose),'">', ' <section name="cluster" />',' <help> N/A </help>','</parameter>'); 41 42 end % }}} 25 43 function cluster=generic(varargin) % {{{ 26 44 -
issm/trunk-jpl/src/m/classes/constants.m
r16764 r17720 11 11 end 12 12 methods 13 function createxml(obj,fid) % {{{ 14 fprintf(fid, '\n\n'); 15 fprintf(fid, '%s\n', '<!-- constants -->'); 16 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n','<parameter key ="g" type="',class(obj.g),'" default="',num2str(obj.g),'">',' <section name="constants" />',' <help> gravitational acceleration [m/s^2] </help>','</parameter>'); 17 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n','<parameter key ="yts" type="',class(obj.yts),'" default="',num2str(obj.yts),'">',' <section name="constants" />',' <help> number of seconds in a year [s/yr] </help>','</parameter>'); 18 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n','<parameter key ="referencetemperature" type="',class(obj.referencetemperature),'" default="',num2str(obj.referencetemperature),'">',' <section name="constants" />',' <help> reference temperature used in the enthalpy model [K] </help>','</parameter>'); 19 end % }}} 13 20 function obj = constants(varargin) % {{{ 14 21 switch nargin -
issm/trunk-jpl/src/m/classes/damage.m
r17674 r17720 30 30 end 31 31 methods 32 function createxml(obj,fid) % {{{ 33 fprintf(fid, '\n\n'); 34 fprintf(fid, '%s\n', '<!-- damage -->'); 35 fprintf(fid, '%s\n', '<!-- Note: this class depends on different input of law -->'); 36 37 %fprintf(fid,'%s%s%s%s\n%s\n%s\n%s\n', '<parameter key ="law" type="logical"', '" default="', num2str(obj.law),'">', ' <section name="damage" />',' <help> damage law (string) from {"undamaged","pralong"} </help>','</parameter>'); 38 % drop-down 39 fprintf(fid,'%s%s%s%s%s\n\t%s\n','<parameter key ="law" type="','alternative','" optional="','false','">','<section name="damage" />'); 40 41 % law = 'undamage' 42 fprintf(fid,'\t%s%s%s%s%s\n\t\t%s\n','<option value="undamage" type="','string','" default="','true','">','<help> law = undamage </help>'); 43 % footer for option 44 fprintf(fid,'\t%s\n%s\n','</option>'); 45 46 % law = 'pralong' 47 fprintf(fid,'\t%s%s%s%s%s\n\t\t%s\n','<option value="pralong" type="','string','" default="','false','">','<help> law = pralong </help>'); 48 49 fprintf(fid,'\t\t%s%s%s%s%s\n\t\t\t%s\n\t\t%s\n', '<parameter key ="stress_threshold" type="',class(obj.stress_threshold),'" default="',num2str(obj.stress_threshold),'">','<help> damage stress threshold [Pa] </help>','</parameter>'); 50 fprintf(fid,'\t\t%s%s%s%s%s\n\t\t\t%s\n\t\t%s\n', '<parameter key ="c1" type="', class(obj.c1),'" default="', num2str(obj.c1),'">', '<help> damage parameter 1 </help>','</parameter>'); 51 fprintf(fid,'\t\t%s%s%s%s%s\n\t\t\t%s\n\t\t%s\n', '<parameter key ="c2" type="', class(obj.c2),'" default="', num2str(obj.c2),'">','<help> damage parameter 2 </help>','</parameter>'); 52 fprintf(fid,'\t\t%s%s%s%s%s\n\t\t\t%s\n\t\t%s\n', '<parameter key ="c3" type="', class(obj.c3),'" default="', num2str(obj.c3),'">','<help> damage parameter 3 [W/m^2] </help>','</parameter>'); 53 fprintf(fid,'\t\t%s%s%s%s%s\n\t\t\t%s\n\t\t%s\n', '<parameter key ="c4" type="', class(obj.c4),'" default="', num2str(obj.c4),'">','<help> damage parameter 4 </help>','</parameter>'); 54 fprintf(fid,'\t\t%s%s%s%s%s\n\t\t\t%s\n\t\t%s\n', '<parameter key ="healing" type="', class(obj.healing),'" default="', num2str(obj.healing),'">','<help> damage healing parameter 1 </help>','</parameter>'); 55 fprintf(fid,'\t\t%s%s%s%s%s\n\t\t\t%s\n\t\t%s\n', '<parameter key ="equiv_stress" type="', class(obj.equiv_stress),'" default="',convert2str(obj.equiv_stress),'">','<help> 0: von Mises </help>','</parameter>'); 56 fprintf(fid,'\t\t%s%s%s%s%s\n\t\t\t%s\n\t\t%s\n', '<parameter key ="requested_outputs" type="', class(obj.requested_outputs),'" default="',convert2str(obj.requested_outputs),'">','<help> additional outputs requested </help>','</parameter>'); 57 58 59 % footer for option 60 fprintf(fid,'\t%s\n%s\n','</option>'); 61 62 63 % footer for drop-down 64 fprintf(fid,'\t%s\n%s\n%s','<help> damage law (string) from {"undamaged","pralong"} </help>','</parameter>'); 65 66 67 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n', '<parameter key ="D" type="', class(obj.D),'" default="', num2str(obj.D),'">', ' <section name="damage" />',' <help> damage tensor (scalar) </help>','</parameter>'); 68 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n', '<parameter key ="law" type="', class(obj.law),'" default="', num2str(obj.law),'">', ' <section name="damage" />',' <help> damage law (string) from {"undamaged","pralong"} </help>','</parameter>'); 69 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n', '<parameter key ="spcdamage" type="', class(obj.spcdamage),'" default="', num2str(obj.spcdamage),'">', ' <section name="damage" />',' <help> damage constraints (NaN means no constraint) </help>','</parameter>'); 70 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n', '<parameter key ="max_damage" type="', class(obj.max_damage),'" default="', num2str(obj.max_damage),'">', ' <section name="damage" />',' <help> maximum possible damage (0&lt;=max_damage&lt;1) </help>','</parameter>'); 71 72 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n', '<parameter key ="stabilization" type="', class(obj.stabilization),'" default="', num2str(obj.stabilization),'">', ' <section name="damage" />',' <help> 0: no, 1: artificial_diffusivity, 2: SUPG </help>','</parameter>'); 73 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n', '<parameter key ="maxiter" type="', class(obj.maxiter),'" default="', num2str(obj.maxiter),'">', ' <section name="damage" />',' <help> maximum number of non linear iterations </help>','</parameter>'); 74 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n', '<parameter key ="penalty_lock" type="', class(obj.penalty_lock),'" default="', num2str(obj.penalty_lock),'">', ' <section name="damage" />',' <help> stabilize unstable damage constraints that keep zigzagging after n iteration (default is 0, no stabilization) </help>','</parameter>'); 75 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n', '<parameter key ="penalty_threshold" type="', class(obj.penalty_threshold),'" default="', num2str(obj.penalty_threshold),'">', ' <section name="damage" />',' <help> threshold to declare convergence of damage evolution solution (default is 0) </help>','</parameter>'); 76 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n', '<parameter key ="penalty_factor" type="', class(obj.penalty_factor),'" default="', num2str(obj.penalty_factor),'">', ' <section name="damage" />',' <help> scaling exponent (default is 3) </help>','</parameter>'); 77 78 end % }}} 79 32 80 function obj = damage(varargin) % {{{ 33 81 switch nargin -
issm/trunk-jpl/src/m/classes/debug.m
r15131 r17720 11 11 end 12 12 methods 13 function createxml(obj,fid) % {{{ 14 fprintf(fid, '\n\n'); 15 fprintf(fid, '%s\n', '<!-- Debug -->'); 16 fprintf(fid,'%s%s%s\n%s\n%s%s%s\n%s\n%s\n', '<parameter key ="valgrind" type="', class(obj.valgrind),'" optional="false">', ' <section name="debug" />',' <option value="',convert2str(obj.valgrind),'" type="string" default="true"></option>',' <help> use Valgrind to debug (0 or 1) </help>','</parameter>'); 17 fprintf(fid,'%s%s%s\n%s\n%s%s%s\n%s\n%s\n', '<parameter key ="gprof" type="', class(obj.gprof),'" optional="false">', ' <section name="debug" />',' <option value="',convert2str(obj.gprof),'" type="string" default="true"></option>',' <help> use gnu-profiler to find out where the time is spent </help>','</parameter>'); 18 fprintf(fid,'%s%s%s\n%s\n%s%s%s\n%s\n', '<parameter key ="profiling" type="', class(obj.profiling),'" optional="false">', ' <section name="debug" />',' <option value="',convert2str(obj.profiling),'" type="string" default="true"></option>',' <help> enables profiling (memory, flops, time) </help>','</parameter>'); 19 end % }}} 13 20 function obj = debug(varargin) % {{{ 14 21 switch nargin -
issm/trunk-jpl/src/m/classes/flaim.m
r16764 r17720 21 21 end 22 22 methods 23 function createxml(obj,fid) % {{{ 24 fprintf(fid, '<!-- flaim -->\n'); 25 26 % Input 27 fprintf(fid,'%s\n%s\n%s\n','<frame key="1" label="Input:">','<section name="flaim" />'); 28 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="targets" type="',class(obj.targets),'" default="',convert2str(obj.targets),'">',' <section name="flaim" />',' <help> name of kml output targets file </help>',' </parameter>'); 29 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="tracks" type="',class(obj.tracks),'" default="',convert2str(obj.tracks),'">',' <section name="flaim" />',' <help> name of kml input tracks file </help>',' </parameter>'); 30 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="flightreqs" type="',class(obj.flightreqs),'" default="',convert2str(obj.flightreqs),'">',' <section name="flaim" />',' <help> structure of kml flight requirements (not used yet) </help>',' </parameter>'); 31 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="criterion" type="',class(obj.criterion),'" default="',convert2str(obj.criterion),'">',' <section name="flaim" />',' <help> element or nodal criterion for flight path evaluation (metric) </help>',' </parameter>'); 32 fprintf(fid,'%s\n%s\n','</frame>'); 33 34 % Arguments 35 fprintf(fid,'%s\n%s\n%s\n','<frame key="2" label="Arguments:">','<section name="flaim" />'); 36 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="gridsatequator" type="',class(obj.gridsatequator),'" default="',convert2str(obj.gridsatequator),'">',' <section name="flaim" />',' <help> number of grids at equator (determines resolution) </help>',' </parameter>'); 37 fprintf(fid,'%s%s%s\n%s\n%s%s%s\n%s\n%s\n',' <parameter key ="usevalueordering" type="',class(obj.usevalueordering),'" optional="false">',' <section name="flaim" />',' <option value="',convert2str(obj.usevalueordering),'" type="string" default="true"></option>',' <help> flag to consider target values for flight path evaluation </help>',' </parameter>'); 38 fprintf(fid,'%s%s%s\n%s\n%s%s%s\n%s\n%s\n',' <parameter key ="split_antimeridian" type="',class(obj.split_antimeridian),'" optional="false">',' <section name="flaim" />',' <option value="',convert2str(obj.split_antimeridian),'" type="string" default="true"></option>',' <help> flag to split polygons on the antimeridian </help>',' </parameter>'); 39 fprintf(fid,'%s\n%s\n','</frame>'); 40 41 % Optimization 42 fprintf(fid,'%s\n%s\n%s\n','<frame key="3" label="Optimization:">','<section name="flaim" />'); 43 fprintf(fid,'%s%s%s\n%s\n%s%s%s\n%s\n%s\n',' <parameter key ="path_optimize" type="',class(obj.path_optimize),'" optional="false">',' <section name="flaim" />',' <option value="',convert2str(obj.path_optimize),'" type="string" default="true"></option>',' <help> optimize? (default false) </help>',' </parameter>'); 44 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="opt_ndir" type="',class(obj.opt_ndir),'" default="',convert2str(obj.opt_ndir),'">',' <section name="flaim" />',' <help> number of directions to test when moving a point. If this value = 1, a random direction is tested. A value > 1 results in directions equally spaced from [0, 2*PI] being tested. For example, 4 would result in directions [0, PI/2, PI, 3PI/2] </help>',' </parameter>'); 45 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="opt_dist" type="',class(obj.opt_dist),'" default="',convert2str(obj.opt_dist),'">',' <section name="flaim" />',' <help> specifies the distance in km (default 25) to move a randomly selected path point on each iteration </help>',' </parameter>'); 46 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="opt_niter" type="',class(obj.opt_niter),'" default="',convert2str(obj.opt_niter),'">',' <section name="flaim" />',' <help> number of iterations (default 30,000) to run for flightplan optimization, i.e. the number of times to randomly select a point and move it. </help>',' </parameter>'); 47 fprintf(fid,'%s\n%s\n','</frame>'); 48 49 % Output 50 fprintf(fid,'%s\n%s\n%s\n','<frame key="4" label="Output:">','<section name="flaim" />'); 51 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="solution" type="',class(obj.solution),'" default="',convert2str(obj.solution),'">',' <section name="flaim" />',' <help> name of kml solution file </help>',' </parameter>'); 52 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="quality" type="',class(obj.quality),'" default="',convert2str(obj.quality),'">',' <section name="flaim" />',' <help> quality of kml solution </help>',' </parameter>'); 53 fprintf(fid,'%s\n%s\n','</frame>'); 54 55 end % }}} 23 56 function obj = flaim(varargin) % {{{ 24 57 switch nargin -
issm/trunk-jpl/src/m/classes/flowequation.m
r17686 r17720 59 59 end 60 60 methods 61 function createxml(obj,fid) % {{{ 62 fprintf(fid, '\n\n'); 63 fprintf(fid, '%s\n', '<!-- flowequation -->'); 64 fprintf(fid,'%s\n%s\n%s\n','<frame key="1" label="Flow equation parameters">','<section name="flowequation" />'); 65 66 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n', '<parameter key ="isSIA" type="', class(obj.isSIA),'" default="', convert2str(obj.isSIA),'">', ' <section name="flowequation" />',' <help> is the Shallow Ice Approximation (SIA) used ? </help>','</parameter>'); 67 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n', '<parameter key ="isSSA" type="', class(obj.isSSA),'" default="', convert2str(obj.isSSA),'">', ' <section name="flowequation" />',' <help> is the Shelfy-Stream Approximation (SSA) used ? </help>','</parameter>'); 68 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n', '<parameter key ="isL1L2" type="', class(obj.isL1L2),'" default="', convert2str(obj.isL1L2),'">', ' <section name="flowequation" />',' <help> is the L1L2 approximation used ? </help>','</parameter>'); 69 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n', '<parameter key ="isHO" type="', class(obj.isHO),'" default="', convert2str(obj.isHO),'">', ' <section name="flowequation" />',' <help> is the Higher-Order (HO) approximation used ? </help>','</parameter>'); 70 71 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n', '<parameter key ="isFS" type="', class(obj.isFS),'" default="', convert2str(obj.isFS),'">', ' <section name="flowequation" />',' <help> are the Full-FS (FS) equations used ? </help>','</parameter>'); 72 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n', '<parameter key ="fe_SSA" type="', class(obj.fe_SSA),'" default="', convert2str(obj.fe_SSA),'">', ' <section name="flowequation" />',' <help> Finite Element for SSA "P1", "P1bubble" "P1bubblecondensed" "P2" </help>','</parameter>'); 73 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n', '<parameter key ="fe_HO" type="', class(obj.fe_HO),'" default="', convert2str(obj.fe_HO),'">', ' <section name="flowequation" />',' <help> Finite Element for HO "P1" "P1bubble" "P1bubblecondensed" "P1xP2" "P2xP1" "P2" </help>','</parameter>'); 74 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n', '<parameter key ="fe_FS" type="', class(obj.fe_FS),'" default="', convert2str(obj.fe_FS),'">', ' <section name="flowequation" />',' <help> Finite Element for FS "P1P1" (debugging only) "P1P1GLS" "MINIcondensed" "MINI" "TaylorHood" </help>','</parameter>'); 75 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n', '<parameter key ="vertex_equation" type="', class(obj.vertex_equation),'" default="', convert2str(obj.vertex_equation),'">', ' <section name="flowequation" />',' <help> flow equation for each vertex </help>','</parameter>'); 76 77 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n', '<parameter key ="element_equation" type="', class(obj.element_equation),'" default="', convert2str(obj.element_equation),'">', ' <section name="flowequation" />',' <help> flow equation for each element </help>','</parameter>'); 78 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n', '<parameter key ="borderSSA" type="', class(obj.borderSSA),'" default="', convert2str(obj.borderSSA),'">', ' <section name="flowequation" />',' <help> vertices on SSAs border (for tiling) </help>','</parameter>'); 79 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n', '<parameter key ="borderHO" type="', class(obj.borderHO),'" default="', convert2str(obj.borderHO),'">', ' <section name="flowequation" />',' <help> vertices on HOs border (for tiling) </help>','</parameter>'); 80 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n', '<parameter key ="borderFS" type="', class(obj.borderFS),'" default="', convert2str(obj.borderFS),'">', ' <section name="flowequation" />',' <help> vertices on FS border (for tiling) </help>','</parameter>'); 81 82 fprintf(fid,'%s\n%s\n','</frame>'); 83 end % }}} 61 84 function obj = flowequation(varargin) % {{{ 62 85 switch nargin -
issm/trunk-jpl/src/m/classes/friction.m
r17051 r17720 11 11 end 12 12 methods 13 function createxml(obj,fid) % {{{ 14 fprintf(fid, '\n\n'); 15 fprintf(fid, '%s\n', '<!-- Friction: Sigma= drag^2 * Neff ^r * u ^s, with Neff=rho_ice*g*thickness+rho_water*g*bed, r=q/p and s=1/p -->'); 16 fprintf(fid,'%s\n%s\n%s\n','<frame key="1" label="Friction: Sigma= drag^2 * Neff ^r * u ^s, with Neff=rho_ice*g*thickness+rho_water*g*bed, r=q/p and s=1/p">','<section name="friction" />'); 17 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n', '<parameter key ="coefficient" type="', class(obj.coefficient),'" default="', convert2str(obj.coefficient),'">', ' <section name="friction" />',' <help> friction coefficient [SI] </help>','</parameter>'); 18 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n', '<parameter key ="p" type="', class(obj.p),'" default="', convert2str(obj.p),'">', ' <section name="friction" />',' <help> p exponent </help>','</parameter>'); 19 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n', '<parameter key ="q" type="', class(obj.q),'" default="', convert2str(obj.q),'">', ' <section name="friction" />',' <help> q exponent </help>','</parameter>'); 20 fprintf(fid,'%s\n%s\n','</frame>'); 21 end % }}} 13 22 function obj = friction(varargin) % {{{ 14 23 switch nargin -
issm/trunk-jpl/src/m/classes/geometry.m
r17590 r17720 33 33 end 34 34 methods 35 function createxml(obj,fid) % {{{ 36 fprintf(fid, '\n\n'); 37 fprintf(fid, '%s\n', '<!-- geometry -->'); 38 fprintf(fid,'%s\n%s\n%s\n','<frame key="1" label="Geometry parameters">','<section name="geometry" />'); 39 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n','<parameter key ="surface" type="','path','" optional="','false','">',' <section name="geometry" />',' <help> surface elevation [m] </help>','</parameter>'); 40 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n','<parameter key ="thickness" type="','path','" optional="','false','">',' <section name="geometry" />',' <help> ice thickness [m] </help>','</parameter>'); 41 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n','<parameter key ="bed" type="','path','" optional="','false','">',' <section name="geometry" />',' <help> bed elevation [m] </help>','</parameter>'); 42 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n','<parameter key ="bathymetry" type="',class(obj.bathymetry),'" default="',convert2str(obj.bathymetry),'">',' <section name="geometry" />',' <help> bathymetry elevation [m] </help>','</parameter>'); 43 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n','<parameter key ="hydrostatic_ratio" type="',class(obj.hydrostatic_ratio),'" default="',convert2str(obj.hydrostatic_ratio),'">',' <section name="geometry" />',' <help> coefficient for ice shelves; thickness correction: hydrostatic_ratio H_obs+ (1-hydrostatic_ratio) H_hydro </help>','</parameter>'); 44 fprintf(fid,'%s\n%s\n','</frame>'); 45 end % }}} 35 46 function obj = geometry(varargin) % {{{ 36 47 switch nargin -
issm/trunk-jpl/src/m/classes/gia.m
r16764 r17720 11 11 end 12 12 methods 13 function createxml(obj,fid) % {{{ 14 fprintf(fid, '<!-- gia -->\n'); 15 16 % gia solution parameters 17 fprintf(fid,'%s\n%s\n%s\n','<frame key="1" label="gia parameters">','<section name="gia" />'); 18 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="mantle_viscosity" type="',class(obj.mantle_viscosity),'" default="',convert2str(obj.mantle_viscosity),'">',' <section name="gia" />',' <help> mantle viscosity[Pa s] </help>',' </parameter>'); 19 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="lithosphere_thickness" type="',class(obj.lithosphere_thickness),'" default="',convert2str(obj.lithosphere_thickness),'">',' <section name="gia" />',' <help> lithosphere thickness (km) </help>',' </parameter>'); 20 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="cross_section_shape" type="',class(obj.cross_section_shape),'" default="',convert2str(obj.cross_section_shape),'">',' <section name="gia" />',' <help> 1: square-edged (default). 2: elliptical. See iedge in GiaDeflectionCore </help>',' </parameter>'); 21 fprintf(fid,'%s\n%s\n','</frame>'); 22 23 end % }}} 13 24 function obj = gia(varargin) % {{{ 14 25 switch nargin -
issm/trunk-jpl/src/m/classes/groundingline.m
r17590 r17720 10 10 end 11 11 methods 12 function createxml(obj,fid) % {{{ 13 fprintf(fid, '\n\n'); 14 fprintf(fid, '%s\n', '<!-- groundingline -->'); 15 16 % Convergence criteria 17 fprintf(fid,'%s\n%s\n%s\n','<frame key="1" label="Grounding line migration parameters">','<section name="groundingline" />'); 18 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n', '<parameter key ="migration" type="', class(obj.migration),'" default="', convert2str(obj.migration),'">', ' <section name="groundingline" />',' <help> type of grounding line migration: "SoftMigration","AgressiveMigration" or "None" </help>','</parameter>'); 19 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n', '<parameter key ="melting_rate" type="', class(obj.melting_rate),'" default="', convert2str(obj.melting_rate),'">', ' <section name="groundingline" />',' <help> melting rate applied when previously grounded parts start floating [m/yr] </help>','</parameter>'); 20 fprintf(fid,'%s\n%s\n','</frame>'); 21 end % }}} 12 22 function obj = groundingline(varargin) % {{{ 13 23 switch nargin -
issm/trunk-jpl/src/m/classes/hydrologyshreve.m
r16764 r17720 15 15 end 16 16 methods 17 function createxml(obj,fid) % {{{ 18 fprintf(fid, '\n\n'); 19 fprintf(fid, '%s\n', '<!-- Hydrology -->'); 20 21 % Convergence criteria 22 fprintf(fid,'%s\n%s\n%s\n','<frame key="1" label="Hydrologyshreve solution parameters">','<section name="hydrologyshreve" />'); 23 24 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n', '<parameter key ="spcwatercolumn" type="', class(obj.spcwatercolumn),'" default="', convert2str(obj.spcwatercolumn),'">', ' <section name="hydrologyshreve" />',' <help> water thickness constraints (NaN means no constraint) [m] </help>','</parameter>'); 25 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n', '<parameter key ="n" type="', class(obj.n),'" default="', convert2str(obj.n),'">', ' <section name="hydrologyshreve" />',' <help> Manning roughness coefficient </help>','</parameter>'); 26 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n', '<parameter key ="CR" type="', class(obj.CR),'" default="', convert2str(obj.CR),'">', ' <section name="hydrologyshreve" />',' <help> tortuosity parameter </help>','</parameter>'); 27 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n', '<parameter key ="p" type="', class(obj.p),'" default="', convert2str(obj.p),'">', ' <section name="hydrologyshreve" />',' <help> dimensionless exponent in Manning velocity formula </help>','</parameter>'); 28 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n', '<parameter key ="q" type="', class(obj.q),'" default="', convert2str(obj.q),'">', ' <section name="hydrologyshreve" />',' <help> dimensionless exponent in Manning velocity formula </help>','</parameter>'); 29 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n', '<parameter key ="kn" type="', class(obj.kn),'" default="', convert2str(obj.kn),'">', ' <section name="hydrologyshreve" />',' <help> parameter in effective pressure formula </help>','</parameter>'); 30 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n', '<parameter key ="stabilization" type="', class(obj.stabilization),'" default="', convert2str(obj.stabilization),'">', ' <section name="hydrologyshreve" />',' <help> artificial diffusivity (default is 1). can be more than 1 to increase diffusivity. </help>','</parameter>'); 31 fprintf(fid,'%s\n%s\n','</frame>'); 32 end % }}} 17 33 function obj = hydrologyshreve(varargin) % {{{ 18 34 switch nargin -
issm/trunk-jpl/src/m/classes/initialization.m
r17686 r17720 19 19 end 20 20 methods 21 function createxml(obj,fid) % {{{ 22 fprintf(fid, '\n\n'); 23 fprintf(fid, '%s\n', '<!-- initialization -->'); 24 fprintf(fid,'%s\n%s\n%s\n','<frame key="1" label="Initial field values">','<section name="initialization" />'); 25 26 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n','<parameter key ="vx" type="',class(obj.vx),'" default="',obj.vx,'">',' <section name="initialization" />',' <help> x component of velocity [m/yr] </help>','</parameter>'); 27 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n','<parameter key ="vy" type="',class(obj.vy),'" default="',obj.vy,'">',' <section name="initialization" />',' <help> y component of velocity [m/yr] </help>','</parameter>'); 28 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n','<parameter key ="vz" type="',class(obj.vz),'" default="',obj.vz,'">',' <section name="initialization" />',' <help> z component of velocity [m/yr] </help>','</parameter>'); 29 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n','<parameter key ="vel" type="',class(obj.vel),'" default="',obj.vel,'">',' <section name="initialization" />',' <help> velocity norm [m/yr] </help>','</parameter>'); 30 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n','<parameter key ="pressure" type="',class(obj.pressure),'" default="',obj.pressure,'">',' <section name="initialization" />',' <help> pressure field [Pa] </help>','</parameter>'); 31 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n','<parameter key ="temperature" type="',class(obj.temperature),'" default="',obj.temperature,'">',' <section name="initialization" />',' <help> fraction of water in the ice </help>','</parameter>'); 32 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n','<parameter key ="waterfraction" type="',class(obj.waterfraction),'" default="',obj.waterfraction,'">',' <section name="initialization" />',' <help> ice thickness [m] </help>','</parameter>'); 33 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n','<parameter key ="sediment_head" type="',class(obj.sediment_head),'" default="',obj.sediment_head,'">',' <section name="initialization" />',' <help> sediment water head of subglacial system [m] </help>','</parameter>'); 34 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n','<parameter key ="epl_head" type="',class(obj.epl_head),'" default="',obj.epl_head,'">',' <section name="initialization" />',' <help> epl water head of subglacial system [m] </help>','</parameter>'); 35 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n','<parameter key ="watercolumn" type="',class(obj.watercolumn),'" default="',obj.watercolumn,'">',' <section name="initialization" />',' <help> thickness of subglacial water [m] </help>','</parameter>'); 36 fprintf(fid,'%s\n%s\n','</frame>'); 37 end % }}} 21 38 function obj = initialization(varargin) % {{{ 22 39 switch nargin -
issm/trunk-jpl/src/m/classes/inversion.m
r16764 r17720 26 26 end 27 27 methods 28 function createxml(obj,fid) % {{{ 29 fprintf(fid, '<!-- inversion -->\n'); 30 31 % inversion parameters 32 fprintf(fid,'%s\n%s\n%s\n','<frame key="1" label="inversion parameters">','<section name="inversion" />'); 33 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="iscontrol" type="',class(obj.iscontrol),'" default="',convert2str(obj.iscontrol),'">',' <section name="inversion" />',' <help> is inversion activated? </help>',' </parameter>'); 34 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="incomplete_adjoint" type="',class(obj.incomplete_adjoint),'" default="',convert2str(obj.incomplete_adjoint),'">',' <section name="inversion" />',' <help> 1: linear viscosity, 0: non-linear viscosity </help>',' </parameter>'); 35 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="control_parameters" type="',class(obj.control_parameters),'" default="',convert2str(obj.control_parameters),'">',' <section name="inversion" />',' <help> ex: {''FrictionCoefficient''}, or {''MaterialsRheologyBbar''} </help>',' </parameter>'); 36 37 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="nsteps" type="',class(obj.nsteps),'" default="',convert2str(obj.nsteps),'">',' <section name="inversion" />',' <help> number of optimization searches </help>',' </parameter>'); 38 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="cost_functions" type="',class(obj.cost_functions),'" default="',convert2str(obj.cost_functions),'">',' <section name="inversion" />',' <help> indicate the type of response for each optimization step </help>',' </parameter>'); 39 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="cost_functions_coefficients" type="',class(obj.cost_functions_coefficients),'" default="',convert2str(obj.cost_functions_coefficients),'">',' <section name="inversion" />',' <help> cost_functions_coefficients applied to the misfit of each vertex and for each control_parameter </help>',' </parameter>'); 40 41 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="cost_function_threshold" type="',class(obj.cost_function_threshold),'" default="',convert2str(obj.cost_function_threshold),'">',' <section name="inversion" />',' <help> misfit convergence criterion. Default is 1%, NaN if not applied </help>',' </parameter>'); 42 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="maxiter_per_step" type="',class(obj.maxiter_per_step),'" default="',convert2str(obj.maxiter_per_step),'">',' <section name="inversion" />',' <help> maximum iterations during each optimization step </help>',' </parameter>'); 43 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="gradient_scaling" type="',class(obj.gradient_scaling),'" default="',convert2str(obj.gradient_scaling),'">',' <section name="inversion" />',' <help> scaling factor on gradient direction during optimization, for each optimization step </help>',' </parameter>'); 44 45 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="step_threshold" type="',class(obj.step_threshold),'" default="',convert2str(obj.step_threshold),'">',' <section name="inversion" />',' <help> decrease threshold for misfit, default is 30% </help>',' </parameter>'); 46 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="min_parameters" type="',class(obj.min_parameters),'" default="',convert2str(obj.min_parameters),'">',' <section name="inversion" />',' <help> absolute minimum acceptable value of the inversed parameter on each vertex </help>',' </parameter>'); 47 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="max_parameters" type="',class(obj.max_parameters),'" default="',convert2str(obj.max_parameters),'">',' <section name="inversion" />',' <help> absolute maximum acceptable value of the inversed parameter on each vertex </help>',' </parameter>'); 48 49 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="vx_obs" type="',class(obj.vx_obs),'" default="',convert2str(obj.vx_obs),'">',' <section name="inversion" />',' <help> observed velocity x component [m/yr] </help>',' </parameter>'); 50 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="vy_obs" type="',class(obj.vy_obs),'" default="',convert2str(obj.vy_obs),'">',' <section name="inversion" />',' <help> observed velocity y component [m/yr] </help>',' </parameter>'); 51 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="vel_obs" type="',class(obj.vel_obs),'" default="',convert2str(obj.vel_obs),'">',' <section name="inversion" />',' <help> observed velocity magnitude [m/yr] </help>',' </parameter>'); 52 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="thickness_obs" type="',class(obj.thickness_obs),'" default="',convert2str(obj.thickness_obs),'">',' <section name="inversion" />',' <help> observed thickness [m]) </help>',' </parameter>'); 53 54 fprintf(fid,'%s\n%s\n','</frame>'); 55 56 fprintf(fid,'%s\n%s\n%s\n','<frame key="2" label="Available cost functions">','<section name="inversion" />'); 57 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="SurfaceAbsVelMisfit" type="','string','" default="','101','">',' <section name="inversion" />',' <help> </help>',' </parameter>'); 58 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="SurfaceRelVelMisfit" type="','string','" default="','102','">',' <section name="inversion" />',' <help> </help>',' </parameter>'); 59 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="SurfaceLogVelMisfit" type="','string','" default="','103','">',' <section name="inversion" />',' <help> </help>',' </parameter>'); 60 61 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="SurfaceLogVxVyMisfit" type="','string','" default="','104','">',' <section name="inversion" />',' <help> </help>',' </parameter>'); 62 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="SurfaceAverageVelMisfit" type="','string','" default="','105','">',' <section name="inversion" />',' <help> </help>',' </parameter>'); 63 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="ThicknessAbsMisfit" type="','string','" default="','106','">',' <section name="inversion" />',' <help> </help>',' </parameter>'); 64 65 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="DragCoefficientAbsGradient" type="','string','" default="','107','">',' <section name="inversion" />',' <help> </help>',' </parameter>'); 66 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="RheologyBbarAbsGradient" type="','string','" default="','108','">',' <section name="inversion" />',' <help> </help>',' </parameter>'); 67 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="ThicknessAbsGradient" type="','string','" default="','109','">',' <section name="inversion" />',' <help> </help>',' </parameter>'); 68 69 fprintf(fid,'%s\n%s\n','</frame>'); 70 71 end % }}} 28 72 function obj = inversion(varargin) % {{{ 29 73 switch nargin … … 32 76 case 1 33 77 if isa(varargin{1},'taoinversion'), 34 disp('converting taoinversion to inversion');78 disp('converting taoinversion'); 35 79 in=varargin{1}; 36 80 obj=inversion(); 37 81 obj.iscontrol = in.iscontrol; 38 obj.tao = 1; 82 obj.tao = 1;rsion to inve 39 83 obj.incomplete_adjoint = in.incomplete_adjoint; 40 84 obj.control_parameters = in.control_parameters; -
issm/trunk-jpl/src/m/classes/mask.m
r17260 r17720 22 22 end 23 23 methods 24 function createxml(obj,fid) % {{{ 25 fprintf(fid, '\n\n'); 26 fprintf(fid, '%s\n', '<!-- mask -->'); 27 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n','<parameter key ="groundedice_levelset" type="',class(obj.groundedice_levelset),'" default="',obj.groundedice_levelset,'">',' <section name="mask" />',' <help> is ice grounded ? grounded ice if > 0, grounding line position if = 0, floating ice if &lt; 0 </help>','</parameter>'); 28 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n','<parameter key ="ice_levelset" type="',class(obj.ice_levelset),'" default="',obj.ice_levelset,'">',' <section name="mask" />',' <help> presence of ice if > 0, icefront position if = 0, no ice if &lt; 0 </help>','</parameter>'); 29 end % }}} 24 30 function obj = mask(varargin) % {{{ 25 31 switch nargin -
issm/trunk-jpl/src/m/classes/masstransport.m
r17555 r17720 43 43 end 44 44 methods 45 function createxml(obj,fid) % {{{ 46 fprintf(fid, '<!-- masstransport -->\n'); 47 48 % Masstransport solution parameters 49 fprintf(fid,'%s\n%s\n%s\n','<frame key="1" label="Masstransport solution parameters">','<section name="masstransport" />'); 50 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="spcthickness" type="',class(obj.spcthickness),'" default="',convert2str(obj.spcthickness),'">',' <section name="masstransport" />',' <help> thickness constraints (NaN means no constraint) [m] </help>',' </parameter>'); 51 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="calvingrate" type="',class(obj.calvingrate),'" default="',convert2str(obj.calvingrate),'">',' <section name="masstransport" />',' <help> calving rate at given location [m/a] </help>',' </parameter>'); 52 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="isfreesurface" type="',class(obj.isfreesurface),'" default="',convert2str(obj.isfreesurface),'">',' <section name="masstransport" />',' <help> do we use free surfaces (FS only) are mass conservation </help>',' </parameter>'); 53 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="min_thickness" type="',class(obj.min_thickness),'" default="',convert2str(obj.min_thickness),'">',' <section name="masstransport" />',' <help> minimum ice thickness allowed [m] </help>',' </parameter>'); 54 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="hydrostatic_adjustment" type="',class(obj.hydrostatic_adjustment),'" default="',convert2str(obj.hydrostatic_adjustment),'">',' <section name="masstransport" />',' <help> adjustment of ice shelves surface and bed elevations: ''Incremental'' or ''Absolute'' </help>',' </parameter>'); 55 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="stabilization" type="',class(obj.stabilization),'" default="',convert2str(obj.stabilization),'">',' <section name="masstransport" />',' <help> 0: no, 1: artificial_diffusivity, 2: streamline upwinding, 3: discontinuous Galerkin </help>',' </parameter>'); 56 fprintf(fid,'%s\n%s\n','</frame>'); 57 58 %Penalty options 59 fprintf(fid,'%s\n%s\n%s\n','<frame key="2" label="Penalty options">','<section name="masstransport" />'); 60 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="vertex_pairing" type="',class(obj.vertex_pairing),'" default="',convert2str(obj.vertex_pairing),'">',' <section name="masstransport" />',' <help> offset used by penalties: penalty = Kmax*10^offset </help>',' </parameter>'); 61 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="penalty_factor" type="',class(obj.penalty_factor),'" default="',convert2str(obj.penalty_factor),'">',' <section name="masstransport" />',' <help> pairs of vertices that are penalized </help>',' </parameter>'); 62 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="requested_outputs" type="',class(obj.requested_outputs),'" default="',convert2str(obj.requested_outputs),'">',' <section name="masstransport" />',' <help> additional outputs requested </help>',' </parameter>'); 63 fprintf(fid,'%s\n%s\n','</frame>'); 64 65 end % }}} 45 66 function obj = masstransport(varargin) % {{{ 46 67 switch nargin -
issm/trunk-jpl/src/m/classes/matice.m
r17456 r17720 30 30 end 31 31 methods 32 function createxml(obj,fid) % {{{ 33 fprintf(fid, '\n\n'); 34 fprintf(fid, '<!-- materials -->\n'); 35 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n','<parameter key ="rho_ice" type="',class(obj.rho_ice),'" default="',convert2str(obj.rho_ice),'">',' <section name="materials" />',' <help> ice density [kg/m^3] </help>','</parameter>'); 36 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n','<parameter key ="rho_water" type="',class(obj.rho_water),'" default="',convert2str(obj.rho_water),'">',' <section name="materials" />',' <help> ocean water density [kg/m^3] </help>','</parameter>'); 37 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n','<parameter key ="rho_freshwater" type="',class(obj.rho_freshwater),'" default="',convert2str(obj.rho_freshwater),'">',' <section name="materials" />',' <help> fresh water density [kg/m^3] </help>','</parameter>'); 38 39 40 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n','<parameter key ="mu_water" type="',class(obj.mu_water),'" default="',convert2str(obj.mu_water),'">',' <section name="materials" />',' <help> water viscosity [N s/m^2] </help>','</parameter>'); 41 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n','<parameter key ="heatcapacity" type="',class(obj.heatcapacity),'" default="',convert2str(obj.heatcapacity),'">',' <section name="materials" />',' <help> heat capacity [J/kg/K] </help>','</parameter>'); 42 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n','<parameter key ="latentheat" type="',class(obj.latentheat),'" default="',convert2str(obj.latentheat),'">',' <section name="materials" />',' <help> latent heat of fusion [J/kg] </help>','</parameter>'); 43 44 45 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n','<parameter key ="thermalconductivity" type="',class(obj.thermalconductivity),'" default="',convert2str(obj.thermalconductivity),'">',' <section name="materials" />',' <help> ice thermal conductivity [W/m/K] </help>','</parameter>'); 46 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n','<parameter key ="temperateiceconductivity" type="',class(obj.temperateiceconductivity),'" default="',convert2str(obj.temperateiceconductivity),'">',' <section name="materials" />',' <help> temperate ice thermal conductivity [W/m/K] </help>','</parameter>'); 47 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n','<parameter key ="meltingpoint" type="',class(obj.meltingpoint),'" default="',convert2str(obj.meltingpoint),'">',' <section name="materials" />',' <help> melting point of ice at 1atm in K </help>','</parameter>'); 48 49 50 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n','<parameter key ="beta" type="',class(obj.beta),'" default="',convert2str(obj.beta),'">',' <section name="materials" />',' <help> rate of change of melting point with pressure [K/Pa] </help>','</parameter>'); 51 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n','<parameter key ="mixed_layer_capacity" type="',class(obj.mixed_layer_capacity),'" default="',convert2str(obj.mixed_layer_capacity),'">',' <section name="materials" />',' <help> mixed layer capacity [W/kg/K] </help>','</parameter>'); 52 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n','<parameter key ="thermal_exchange_velocity" type="',class(obj.thermal_exchange_velocity),'" default="',convert2str(obj.thermal_exchange_velocity),'">',' <section name="materials" />',' <help> thermal exchange velocity [m/s] </help>','</parameter>'); 53 54 55 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n','<parameter key ="rheology_B" type="',class(obj.rheology_B),'" default="',convert2str(obj.rheology_B),'">',' <section name="materials" />',' <help> flow law parameter [Pa/s^(1/n)] </help>','</parameter>'); 56 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n','<parameter key ="rheology_n" type="',class(obj.rheology_n),'" default="',convert2str(obj.rheology_n),'">',' <section name="materials" />',' <help> Glens flow law exponent </help>','</parameter>'); 57 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n','<parameter key ="rheology_law" type="',class(obj.rheology_law),'" default="',convert2str(obj.rheology_law),'">',' <section name="materials" />',' <help> law for the temperature dependance of the rheology: "None", "Paterson", "Arrhenius" or "LliboutryDuval" </help>','</parameter>'); 58 59 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n','<parameter key ="lithosphere_shear_modulus" type="',class(obj.lithosphere_shear_modulus),'" default="',convert2str(obj.lithosphere_shear_modulus),'">',' <section name="materials" />',' <help> Lithosphere shear modulus [Pa] </help>','</parameter>'); 60 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n','<parameter key ="lithosphere_density" type="',class(obj.lithosphere_density),'" default="',convert2str(obj.lithosphere_density),'">',' <section name="materials" />',' <help> Lithosphere density [g/cm^-3] </help>','</parameter>'); 61 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n','<parameter key ="mantle_shear_modulus" type="',class(obj.mantle_shear_modulus),'" default="',convert2str(obj.mantle_shear_modulus),'">',' <section name="materials" />',' <help> Mantle shear modulus [Pa] </help>','</parameter>'); 62 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n','<parameter key ="mantle_density" type="',class(obj.mantle_density),'" default="',convert2str(obj.mantle_density),'">',' <section name="materials" />',' <help> Mantle density [g/cm^-3] </help>','</parameter>'); 63 64 65 end % }}} 32 66 function obj = matice(varargin) % {{{ 33 67 switch nargin -
issm/trunk-jpl/src/m/classes/mesh2d.m
r17703 r17720 104 104 fielddisplay(obj,'hemisphere','Indicate hemisphere ''n'' or ''s'' '); 105 105 end % }}} 106 function createxml(obj,fid) % {{{ 107 fprintf(fid, '<!-- mesh -->\n'); 108 109 % drop-down 110 fprintf(fid,'%s%s%s%s%s\n%s\n','<parameter key ="Dimension" type="','alternative','" optional="','false','">',' <section name="mesh" />'); 111 112 % dimension=0 113 fprintf(fid,'%s%s%s%s%s\n%s\n',' <option value="0" type="','string','" default="','true','">',' <help> dimension is 0 </help>'); 114 % footer for dim=0 115 fprintf(fid,'%s\n%s\n',' </option>'); 116 117 % dimension=3 118 fprintf(fid,'%s%s%s\n%s\n',' <option value="3" type="','string','">',' <help> dimension is 3 </help>'); 119 120 fprintf(fid,' <!-- Elements and vertices of the original 2d mesh -->\n'); 121 %fprintf(fid,'%s\n%s\n%s\n',' <frame key="dim=3" label="Elements and vertices of the original 2d mesh">',' <section name="mesh" />'); 122 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n',' <parameter key ="x2d" type="',class(obj.x2d),'" default="',convert2str(obj.x2d),'">',' <help> vertices x coordinate [m] </help>',' </parameter>'); 123 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n',' <parameter key ="y2d" type="',class(obj.y2d),'" default="',convert2str(obj.y2d),'">',' <help> vertices y coordinate [m] </help>',' </parameter>'); 124 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n',' <parameter key ="elements2d" type="',class(obj.elements2d),'" default="',convert2str(obj.elements2d),'">',' <help> vertex indices of the mesh elements </help>',' </parameter>'); 125 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n',' <parameter key ="numberofvertices2d" type="',class(obj.numberofvertices2d),'" default="',convert2str(obj.numberofvertices2d),'">',' <help> number of vertices </help>',' </parameter>'); 126 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n',' <parameter key ="numberofelements2d" type="',class(obj.numberofelements2d),'" default="',convert2str(obj.numberofelements2d),'">',' <help> number of elements </help>',' </parameter>'); 127 %fprintf(fid,'%s\n%s\n',' </frame>'); 128 % footer for dim=3 129 fprintf(fid,'%s\n%s\n',' </option>'); 130 131 % footer for drop-down 132 fprintf(fid,'%s\n%s\n%s',' <help> select a mesh dimension 0(default) or 3 </help>','</parameter>'); 133 134 135 % elements and vertices 136 fprintf(fid,'%s\n%s\n%s\n','<frame key="1" label="Elements and vertices">','<section name="mesh" />'); 137 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="numberofelements" type="',class(obj.numberofelements),'" default="',convert2str(obj.numberofelements),'">',' <section name="mesh" />',' <help> number of elements </help>',' </parameter>'); 138 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="numberofvertices" type="',class(obj.numberofvertices),'" default="',convert2str(obj.numberofvertices),'">',' <section name="mesh" />',' <help> number of verticies </help>',' </parameter>'); 139 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="elements" type="',class(obj.elements),'" default="',convert2str(obj.elements),'">',' <section name="mesh" />',' <help> vertex indices of the mesh elements </help>',' </parameter>'); 140 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="x" type="',class(obj.x),'" default="',convert2str(obj.x),'">',' <section name="mesh" />',' <help> vertices x coordinate [m] </help>',' </parameter>'); 141 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="y" type="',class(obj.y),'" default="',convert2str(obj.y),'">',' <section name="mesh" />',' <help> vertices y coordinate [m] </help>',' </parameter>'); 142 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="z" type="',class(obj.z),'" default="',convert2str(obj.z),'">',' <section name="mesh" />',' <help> vertices z coordinate [m] </help>',' </parameter>'); 143 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="edges" type="',class(obj.edges),'" default="',convert2str(obj.edges),'">',' <section name="mesh" />',' <help> edges of the 2d mesh (vertex1 vertex2 element1 element2) </help>',' </parameter>'); 144 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="numberofedges" type="',class(obj.numberofedges),'" default="',convert2str(obj.numberofedges),'">',' <section name="mesh" />',' <help> number of edges of the 2d mesh </help>',' </parameter>'); 145 fprintf(fid,'%s\n%s\n','</frame>'); 146 147 % properties 148 fprintf(fid,'%s\n%s\n%s\n','<frame key="2" label="Properties">','<section name="mesh" />'); 149 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="selected dimension" type="',class(obj.dimension),'" default="',convert2str(obj.dimension),'">',' <section name="mesh" />',' <help> mesh dimension </help>',' </parameter>'); 150 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="numberoflayers" type="',class(obj.numberoflayers),'" default="',convert2str(obj.numberoflayers),'">',' <section name="mesh" />',' <help> number of extrusion layers </help>',' </parameter>'); 151 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="vertexonbed" type="',class(obj.vertexonbed),'" default="',convert2str(obj.vertexonbed),'">',' <section name="mesh" />',' <help> lower vertices flags list </help>',' </parameter>'); 152 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="elementonbed" type="',class(obj.elementonbed),'" default="',convert2str(obj.elementonbed),'">',' <section name="mesh" />',' <help> lower elements flags list </help>',' </parameter>'); 153 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="vertexonsurface" type="',class(obj.vertexonsurface),'" default="',convert2str(obj.vertexonsurface),'">',' <section name="mesh" />',' <help> upper vertices flags list </help>',' </parameter>'); 154 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="elementonsurface" type="',class(obj.elementonsurface),'" default="',convert2str(obj.elementonsurface),'">',' <section name="mesh" />',' <help> upper element flag list </help>',' </parameter>'); 155 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="uppervertex" type="',class(obj.uppervertex),'" default="',convert2str(obj.uppervertex),'">',' <section name="mesh" />',' <help> upper vertex list (NaN for vertex on the upper surface) </help>',' </parameter>'); 156 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="upperelements" type="',class(obj.upperelements),'" default="',convert2str(obj.upperelements),'">',' <section name="mesh" />',' <help> upper element list (NaN for element on the upper layer) </help>',' </parameter>'); 157 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="lowervertex" type="',class(obj.lowervertex),'" default="',convert2str(obj.lowervertex),'">',' <section name="mesh" />',' <help> lower vertex list (NaN for vertex on the lower surface) </help>',' </parameter>'); 158 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="lowerelements" type="',class(obj.lowerelements),'" default="',convert2str(obj.lowerelements),'">',' <section name="mesh" />',' <help> lower element list (NaN for element on the lower layer </help>',' </parameter>'); 159 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="vertexonboundary" type="',class(obj.vertexonboundary),'" default="',convert2str(obj.vertexonboundary),'">',' <section name="mesh" />',' <help> vertices on the boundary of the domain flag list </help>',' </parameter>'); 160 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="segments" type="',class(obj.segments),'" default="',convert2str(obj.segments),'">',' <section name="mesh" />',' <help> edges on domain boundary (vertex1 vertex2 element) </help>',' </parameter>'); 161 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="segmentmarkers" type="',class(obj.segmentmarkers),'" default="',convert2str(obj.segmentmarkers),'">',' <section name="mesh" />',' <help> number associated to each segment </help>',' </parameter>'); 162 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="vertexconnectivity" type="',class(obj.vertexconnectivity),'" default="',convert2str(obj.vertexconnectivity),'">',' <section name="mesh" />',' <help> list of vertices connected to vertex_i </help>',' </parameter>'); 163 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="elementconnectivity" type="',class(obj.elementconnectivity),'" default="',convert2str(obj.elementconnectivity),'">',' <section name="mesh" />',' <help> list of vertices connected to element_i </help>',' </parameter>'); 164 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="average_vertex_connectivity" type="',class(obj.average_vertex_connectivity),'" default="',convert2str(obj.average_vertex_connectivity),'">',' <section name="mesh" />',' <help> average number of vertices connected to one vertex </help>',' </parameter>'); 165 fprintf(fid,'%s\n%s\n','</frame>'); 166 167 %extracted model 168 fprintf(fid,'%s\n%s\n%s\n','<frame key="3" label="Extracted Model">','<section name="mesh" />'); 169 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="extractedvertices" type="',class(obj.extractedvertices),'" default="',convert2str(obj.extractedvertices),'">',' <section name="mesh" />',' <help> vertices extracted from the model </help>',' </parameter>'); 170 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="extractedelements" type="',class(obj.extractedelements),'" default="',convert2str(obj.extractedelements),'">',' <section name="mesh" />',' <help> elements extracted from the model </help>',' </parameter>'); 171 fprintf(fid,'%s\n%s\n','</frame>'); 172 173 %projection 174 fprintf(fid,'%s\n%s\n%s\n','<frame key="4" label="Projection">','<section name="mesh" />'); 175 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="lat" type="',class(obj.lat),'" default="',convert2str(obj.lat),'">',' <section name="mesh" />',' <help> vertices latitude [degrees] </help>',' </parameter>'); 176 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="long" type="',class(obj.long),'" default="',convert2str(obj.long),'">',' <section name="mesh" />',' <help> verticies longitude [degrees] </help>',' </parameter>'); 177 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n',' <parameter key ="hemisphere" type="',class(obj.hemisphere),'" default="',convert2str(obj.hemisphere),'">',' <section name="mesh" />',' <help> Indicate hemisphere ''n'' or ''s'' </help>',' </parameter>'); 178 fprintf(fid,'%s\n%s\n','</frame>'); 179 180 end % }}} 106 181 function marshall(obj,md,fid) % {{{ 107 182 WriteData(fid,'enum',DomainTypeEnum(),'data',StringToEnum(['Domain' domaintype(obj)]),'format','Integer'); -
issm/trunk-jpl/src/m/classes/miscellaneous.m
r16764 r17720 11 11 end 12 12 methods 13 function createxml(obj,fid) % {{{ 14 fprintf(fid, '<!-- miscellaneous -->\n'); 15 16 % miscellaneous solution parameters 17 fprintf(fid,'%s\n%s\n%s\n','<frame key="1" label="miscellaneous parameters">','<section name="miscellaneous" />'); 18 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="notes" type="',class(obj.notes),'" default="',convert2str(obj.notes),'">',' <section name="miscellaneous" />',' <help> notes in a cell of strings </help>',' </parameter>'); 19 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="name" type="',class(obj.name),'" default="',convert2str(obj.name),'">',' <section name="miscellaneous" />',' <help> model name </help>',' </parameter>'); 20 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="dummy" type="',class(obj.dummy),'" default="',convert2str(obj.dummy),'">',' <section name="miscellaneous" />',' <help> empty field to store some data </help>',' </parameter>'); 21 fprintf(fid,'%s\n%s\n','</frame>'); 22 23 end % }}} 13 24 function obj = miscellaneous(varargin) % {{{ 14 25 switch nargin -
issm/trunk-jpl/src/m/classes/outputdefinition.m
r16764 r17720 9 9 end 10 10 methods 11 function createxml(obj,fid) % {{{ 12 fprintf(fid, '<!-- outputdefinition -->\n'); 13 14 % outputdefinition solution parameters 15 fprintf(fid,'%s\n%s\n%s\n','<frame key="1" label="outputdefinition">','<section name="outputdefinition" />'); 16 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="definitions" type="',class(obj.definitions),'" default="',convert2str(obj.definitions),'">',' <section name="outputdefinition" />',' <help> list of potential outputs that can be requested, but which need additional data to be defined </help>',' </parameter>'); 17 fprintf(fid,'%s\n%s\n','</frame>'); 18 19 end % }}} 11 20 function obj = outputdefinition(varargin) % {{{ 12 21 switch nargin -
issm/trunk-jpl/src/m/classes/qmu.m
r16304 r17720 24 24 end 25 25 methods 26 function createxml(obj,fid) % {{{ 27 fprintf(fid, '<!-- qmu -->\n'); 28 29 % qmu solution parameters 30 fprintf(fid,'%s\n%s\n%s\n','<frame key="1" label="qmu parameters">','<section name="qmu" />'); 31 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="isdakota" type="',class(obj.isdakota),'" default="',convert2str(obj.isdakota),'">',' <section name="qmu" />',' <help> is qmu analysis activated? </help>',' </parameter>'); 32 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="variables" type="',class(obj.variables),'" default="',convert2str(obj.variables),'">',' <section name="qmu" />',' <help> (arrays of each variable class) </help>',' </parameter>'); 33 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="responses" type="',class(obj.responses),'" default="',convert2str(obj.responses),'">',' <section name="qmu" />',' <help> (arrays of each response class) </help>',' </parameter>'); 34 35 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="numberofresponses" type="',class(obj.numberofresponses),'" default="',convert2str(obj.numberofresponses),'">',' <section name="qmu" />',' <help> number of responses </help>',' </parameter>'); 36 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="params" type="',class(obj.params),'" default="',convert2str(obj.params),'">',' <section name="qmu" />',' <help> (array of method-independent parameters) </help>',' </parameter>'); 37 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="results" type="',class(obj.results),'" default="',convert2str(obj.results),'">',' <section name="qmu" />',' <help> (information from dakota files) </help>',' </parameter>'); 38 39 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="partition" type="',class(obj.partition),'" default="',convert2str(obj.partition),'">',' <section name="qmu" />',' <help> user provided mesh partitioning, defaults to metis if not specified </help>',' </parameter>'); 40 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="numberofpartitions" type="',class(obj.numberofpartitions),'" default="',convert2str(obj.numberofpartitions),'">',' <section name="qmu" />',' <help> number of partitions for semi-discrete qmu </help>',' </parameter>'); 41 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="variabledescriptors" type="',class(obj.variabledescriptors),'" default="',convert2str(obj.variabledescriptors),'">',' <section name="qmu" />',' <help> </help>',' </parameter>'); 42 43 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="responsedescriptors" type="',class(obj.responsedescriptors),'" default="',convert2str(obj.responsedescriptors),'">',' <section name="qmu" />',' <help> use an enthalpy formulation to include temperate ice (default is 0) </help>',' </parameter>'); 44 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="method" type="',class(obj.method),'" default="',convert2str(obj.method),'">',' <section name="qmu" />',' <help> array of dakota_method class </help>',' </parameter>'); 45 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="mass_flux_profile_directory" type="',class(obj.mass_flux_profile_directory),'" default="',convert2str(obj.mass_flux_profile_directory),'">',' <section name="qmu" />',' <help> directory for mass flux profiles </help>',' </parameter>'); 46 47 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="mass_flux_profiles" type="',class(obj.mass_flux_profiles),'" default="',convert2str(obj.mass_flux_profiles),'">',' <section name="qmu" />',' <help> </help>',' </parameter>'); 48 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="mass_flux_segments" type="',class(obj.mass_flux_segments),'" default="',convert2str(obj.mass_flux_segments),'">',' <section name="qmu" />',' <help> </help>',' </parameter>'); 49 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="adjacency" type="',class(obj.adjacency),'" default="',convert2str(obj.adjacency),'">',' <section name="qmu" />',' <help> additional outputs requested </help>',' </parameter>'); 50 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="vertex_weight" type="',class(obj.vertex_weight),'" default="',convert2str(obj.vertex_weight),'">',' <section name="qmu" />',' <help> weight applied to each mesh vertex </help>',' </parameter>'); 51 52 fprintf(fid,'%s\n%s\n','</frame>'); 53 54 end % }}} 26 55 function obj = qmu(varargin) % {{{ 27 56 switch nargin -
issm/trunk-jpl/src/m/classes/rifts.m
r17686 r17720 10 10 end 11 11 methods 12 function createxml(obj,fid) % {{{ 13 fprintf(fid, '\n\n'); 14 fprintf(fid, '%s\n', '<!-- rifts -->'); 15 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n', '<parameter key ="riftstruct" type="', class(obj.riftstruct),'" default="', obj.riftstruct,'">', ' <section name="rifts" />',' <help> structure containing all rift information (vertices coordinates, segments, type of melange, ...) </help>','</parameter>'); 16 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n', '<parameter key ="riftproperties" type="', class(obj.riftproperties),'" default="', obj.riftproperties,'">', ' <section name="rifts" />',' <help> N/A </help>','</parameter>'); 17 18 end % }}} 12 19 function obj = rifts(varargin) % {{{ 13 20 switch nargin -
issm/trunk-jpl/src/m/classes/settings.m
r17100 r17720 19 19 end 20 20 methods 21 function createxml(obj,fid) % {{{ 22 fprintf(fid, '\n\n'); 23 fprintf(fid, '%s\n', '<!-- settings -->'); 24 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n', '<parameter key ="results_on_nodes" type="', class(obj.results_on_nodes),'" default="', convert2str(obj.results_on_nodes),'">', ' <section name="settings" />',' <help> results are output for all the nodes of each element </help>','</parameter>'); 25 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n', '<parameter key ="io_gather" type="', class(obj.io_gather),'" default="', convert2str(obj.io_gather),'">', ' <section name="settings" />',' <help> I/O gathering strategy for result outputs (default 1) </help>','</parameter>'); 26 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n', '<parameter key ="lowmem" type="', class(obj.lowmem),'" default="', convert2str(obj.lowmem),'">', ' <section name="settings" />',' <help> is the memory limited ? (0 or 1) </help>','</parameter>'); 27 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n', '<parameter key ="output_frequency" type="', class(obj.output_frequency),'" default="', convert2str(obj.output_frequency),'">', ' <section name="settings" />',' <help> frequency at which results are saved in all solutions with multiple time_steps </help>','</parameter>'); 28 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n', '<parameter key ="waitonlock" type="', class(obj.waitonlock),'" default="', convert2str(obj.waitonlock),'">', ' <section name="settings" />',' <help> maximum number of minutes to wait for batch results (NaN to deactivate) </help>','</parameter>'); 29 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n', '<parameter key ="upload_server" type="', class(obj.upload_server),'" default="', convert2str(obj.upload_server),'">', ' <section name="settings" />',' <help> server hostname where model should be uploaded </help>','</parameter>'); 30 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n', '<parameter key ="upload_path" type="', class(obj.upload_path),'" default="', convert2str(obj.upload_path),'">', ' <section name="settings" />',' <help> path on server where model should be uploaded </help>','</parameter>'); 31 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n', '<parameter key ="upload_login" type="', class(obj.upload_login),'" default="', convert2str(obj.upload_login),'">', ' <section name="settings" />',' <help> server login </help>','</parameter>'); 32 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n', '<parameter key ="upload_port" type="', class(obj.upload_port),'" default="', convert2str(obj.upload_port),'">', ' <section name="settings" />',' <help> port login (default is 0) </help>','</parameter>'); 33 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n', '<parameter key ="upload_filename" type="', class(obj.upload_filename),'" default="', convert2str(obj.upload_filename),'">', ' <section name="settings" />',' <help> unique id generated when uploading the file to server </help>','</parameter>'); 34 35 end % }}} 21 36 function obj = settings(varargin) % {{{ 22 37 switch nargin -
issm/trunk-jpl/src/m/classes/steadystate.m
r16764 r17720 11 11 end 12 12 methods 13 function createxml(obj,fid) % {{{ 14 fprintf(fid, '<!-- steadystate -->\n'); 15 16 % steadystate solution parameters 17 fprintf(fid,'%s\n%s\n%s\n','<frame key="1" label="steadystate solution parameters">','<section name="steadystate" />'); 18 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="reltol" type="',class(obj.reltol),'" default="',convert2str(obj.reltol),'">',' <section name="steadystate" />',' <help> relative tolerance criterion [K] </help>',' </parameter>'); 19 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="maxiter" type="',class(obj.maxiter),'" default="',convert2str(obj.maxiter),'">',' <section name="steadystate" />',' <help> maximum number of iterations </help>',' </parameter>'); 20 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="requested_outputs" type="',class(obj.requested_outputs),'" default="',convert2str(obj.requested_outputs),'">',' <section name="steadystate" />',' <help> additional requested outputs </help>',' </parameter>'); 21 fprintf(fid,'%s\n%s\n','</frame>'); 22 23 end % }}} 13 24 function obj = steadystate(varargin) % {{{ 14 25 switch nargin -
issm/trunk-jpl/src/m/classes/stressbalance.m
r17686 r17720 26 26 end 27 27 methods 28 function createxml(obj,fid) % {{{ 29 fprintf(fid, '\n\n'); 30 fprintf(fid, '%s\n', '<!-- StressBalance solution parameters: -->'); 31 32 % Convergence criteria 33 fprintf(fid,'%s\n%s\n%s\n','<frame key="1" label="Convergence criteria">','<section name="stressbalance" />'); 34 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n', '<parameter key ="restol" type="', class(obj.restol),'" default="', convert2str(obj.restol),'">', ' <section name="stressbalance" />',' <help> mechanical equilibrium residual convergence criterion </help>','</parameter>'); 35 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n', '<parameter key ="reltol" type="', class(obj.reltol),'" default="', convert2str(obj.reltol),'">', ' <section name="stressbalance" />',' <help> velocity relative convergence criterion, NaN: not applied </help>','</parameter>'); 36 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n', '<parameter key ="abstol" type="', class(obj.abstol),'" default="', convert2str(obj.abstol),'">', ' <section name="stressbalance" />',' <help> velocity absolute convergence criterion, NaN: not applied </help>','</parameter>'); 37 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n', '<parameter key ="isnewton" type="', class(obj.isnewton),'" default="', convert2str(obj.isnewton),'">', ' <section name="stressbalance" />',' <help> 0: Picards fixed point, 1: Newtons method, 2: hybrid </help>','</parameter>'); 38 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n', '<parameter key ="maxiter" type="', class(obj.maxiter),'" default="', convert2str(obj.maxiter),'">', ' <section name="stressbalance" />',' <help> maximum number of nonlinear iterations </help>','</parameter>'); 39 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n', '<parameter key ="viscosity_overshoot" type="', class(obj.viscosity_overshoot),'" default="', convert2str(obj.viscosity_overshoot),'">', ' <section name="stressbalance" />',' <help> over-shooting constant new=new+C*(new-old) </help>','</parameter>'); 40 fprintf(fid,'%s\n%s\n','</frame>'); 41 42 % boundary conditions 43 fprintf(fid,'%s\n%s\n%s\n','<frame key="2" label="Boundary conditions">','<section name="stressbalance" />'); 44 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n', '<parameter key ="spcvx" type="', class(obj.spcvx),'" default="', convert2str(obj.spcvx),'">', ' <section name="stressbalance" />',' <help> x-axis velocity constraint (NaN means no constraint) [m/yr] </help>','</parameter>'); 45 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n', '<parameter key ="spcvy" type="', class(obj.spcvy),'" default="', convert2str(obj.spcvy),'">', ' <section name="stressbalance" />',' <help> y-axis velocity constraint (NaN means no constraint) [m/yr] </help>','</parameter>'); 46 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n', '<parameter key ="spcvz" type="', class(obj.spcvz),'" default="', convert2str(obj.spcvz),'">', ' <section name="stressbalance" />',' <help> z-axis velocity constraint (NaN means no constraint) [m/yr] </help>','</parameter>'); 47 fprintf(fid,'%s\n%s\n','</frame>'); 48 49 % Rift options 50 fprintf(fid,'%s\n%s\n%s\n','<frame key="3" label="Rift options">','<section name="stressbalance" />'); 51 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n', '<parameter key ="rift_penalty_threshold" type="', class(obj.rift_penalty_threshold),'" default="', convert2str(obj.rift_penalty_threshold),'">', ' <section name="stressbalance" />',' <help> threshold for instability of mechanical constraints </help>','</parameter>'); 52 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n', '<parameter key ="rift_penalty_lock" type="', class(obj.rift_penalty_lock),'" default="', convert2str(obj.rift_penalty_lock),'">', ' <section name="stressbalance" />',' <help> number of iterations before rift penalties are locked </help>','</parameter>'); 53 fprintf(fid,'%s\n%s\n','</frame>'); 54 55 %others 56 fprintf(fid,'%s\n%s\n%s\n','<frame key="4" label="Other">','<section name="stressbalance" />'); 57 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n', '<parameter key ="shelf_dampening" type="', class(obj.shelf_dampening),'" default="', convert2str(obj.shelf_dampening),'">', ' <section name="stressbalance" />',' <help> use dampening for floating ice ? Only for FS model </help>','</parameter>'); 58 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n', '<parameter key ="FSreconditioning" type="', class(obj.FSreconditioning),'" default="', convert2str(obj.FSreconditioning),'">', ' <section name="stressbalance" />',' <help> multiplier for incompressibility equation. Only for FS model </help>','</parameter>'); 59 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n', '<parameter key ="referential" type="', class(obj.referential),'" default="', convert2str(obj.referential),'">', ' <section name="stressbalance" />',' <help> local referential </help>','</parameter>'); 60 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n', '<parameter key ="loadingforce" type="', class(obj.loadingforce),'" default="', convert2str(obj.loadingforce),'">', ' <section name="stressbalance" />',' <help> loading force applied on each point [N/m^3] </help>','</parameter>'); 61 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n', '<parameter key ="requested_outputs" type="', class(obj.requested_outputs),'" default="', convert2str(obj.requested_outputs),'">', ' <section name="stressbalance" />',' <help> additional outputs requested </help>','</parameter>'); 62 63 64 fprintf(fid,'%s\n%s\n','</frame>'); 65 end % }}} 28 66 function obj = stressbalance(varargin) % {{{ 29 67 switch nargin -
issm/trunk-jpl/src/m/classes/thermal.m
r17686 r17720 17 17 end 18 18 methods 19 function createxml(obj,fid) % {{{ 20 fprintf(fid, '<!-- thermal -->\n'); 21 22 % thermal solution parameters 23 fprintf(fid,'%s\n%s\n%s\n','<frame key="1" label="Thermal solution parameters">','<section name="thermal" />'); 24 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="spctemperature" type="',class(obj.spctemperature),'" default="',convert2str(obj.spctemperature),'">',' <section name="thermal" />',' <help> temperature constraints (NaN means no constraint) [K] </help>',' </parameter>'); 25 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="penalty_threshold" type="',class(obj.penalty_threshold),'" default="',convert2str(obj.penalty_threshold),'">',' <section name="thermal" />',' <help> 0: no, 1: artificial_diffusivity, 2: SUPG </help>',' </parameter>'); 26 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="stabilization" type="',class(obj.stabilization),'" default="',convert2str(obj.stabilization),'">',' <section name="thermal" />',' <help> maximum number of non linear iterations </help>',' </parameter>'); 27 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="maxiter" type="',class(obj.maxiter),'" default="',convert2str(obj.maxiter),'">',' <section name="thermal" />',' <help> stabilize unstable thermal constraints that keep zigzagging after n iteration (default is 0, no stabilization) </help>',' </parameter>'); 28 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="penalty_lock" type="',class(obj.penalty_lock),'" default="',convert2str(obj.penalty_lock),'">',' <section name="thermal" />',' <help> threshold to declare convergence of thermal solution (default is 0) </help>',' </parameter>'); 29 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="penalty_factor" type="',class(obj.penalty_factor),'" default="',convert2str(obj.penalty_factor),'">',' <section name="thermal" />',' <help> scaling exponent (default is 3) </help>',' </parameter>'); 30 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="isenthalpy" type="',class(obj.isenthalpy),'" default="',convert2str(obj.isenthalpy),'">',' <section name="thermal" />',' <help> use an enthalpy formulation to include temperate ice (default is 0) </help>',' </parameter>'); 31 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="isdynamicbasalspc" type="',class(obj.isdynamicbasalspc),'" default="',convert2str(obj.isdynamicbasalspc),'">',' <section name="thermal" />',' <help> enable dynamic setting of basal forcing. required for enthalpy formulation (default is 0) </help>',' </parameter>'); 32 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="requested_outputs" type="',class(obj.requested_outputs),'" default="',convert2str(obj.requested_outputs),'">',' <section name="thermal" />',' <help> additional outputs requested </help>',' </parameter>'); 33 fprintf(fid,'%s\n%s\n','</frame>'); 34 35 end % }}} 19 36 function obj = thermal(varargin) % {{{ 20 37 switch nargin -
issm/trunk-jpl/src/m/classes/timestepping.m
r16764 r17720 13 13 end 14 14 methods 15 function createxml(obj,fid) % {{{ 16 fprintf(fid, '\n\n'); 17 fprintf(fid, '%s\n', '<!-- Timestepping -->'); 18 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n', '<parameter key ="start_time" type="', class(obj.start_time),'" default="', convert2str(obj.start_time),'">', ' <section name="timestepping" />',' <help> simulation starting time [yr] </help>','</parameter>'); 19 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n', '<parameter key ="final_time" type="', class(obj.final_time),'" default="', convert2str(obj.final_time),'">', ' <section name="timestepping" />',' <help> final time to stop the simulation [yr] </help>','</parameter>'); 20 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n', '<parameter key ="time_step" type="', class(obj.time_step),'" default="', convert2str(obj.time_step),'">', ' <section name="timestepping" />',' <help> length of time steps [yr] </help>','</parameter>'); 21 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n', '<parameter key ="time_adapt" type="', class(obj.time_adapt),'" default="', convert2str(obj.time_adapt),'">', ' <section name="timestepping" />',' <help> use cfl condition to define time step ? (0 or 1) </help>','</parameter>'); 22 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n', '<parameter key ="cfl_coefficient" type="', class(obj.cfl_coefficient),'" default="', convert2str(obj.cfl_coefficient),'">', ' <section name="timestepping" />',' <help> coefficient applied to cfl condition </help>','</parameter>'); 23 24 end % }}} 25 15 26 function obj = timestepping(varargin) % {{{ 16 27 switch nargin -
issm/trunk-jpl/src/m/classes/transient.m
r17644 r17720 16 16 end 17 17 methods 18 function createxml(obj,fid) % {{{ 19 fprintf(fid, '<!-- transient -->\n'); 20 21 % transient solution parameters 22 fprintf(fid,'%s\n%s\n%s\n','<frame key="1" label="transient solution parameters">','<section name="transient" />'); 23 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="ismasstransport" type="',class(obj.ismasstransport),'" default="',convert2str(obj.ismasstransport),'">',' <section name="transient" />',' <help> indicates whether a masstransport solution is used in the transient </help>',' </parameter>'); 24 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="isstressbalance" type="',class(obj.isstressbalance),'" default="',convert2str(obj.isstressbalance),'">',' <section name="transient" />',' <help> indicates whether a stressbalance solution is used in the transient </help>',' </parameter>'); 25 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="isthermal" type="',class(obj.isthermal),'" default="',convert2str(obj.isthermal),'">',' <section name="transient" />',' <help> indicates whether a thermal solution is used in the transient </help>',' </parameter>'); 26 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="isgroundingline" type="',class(obj.isgroundingline),'" default="',convert2str(obj.isgroundingline),'">',' <section name="transient" />',' <help> indicates whether a groundingline migration is used in the transient </help>',' </parameter>'); 27 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="isgia" type="',class(obj.isgia),'" default="',convert2str(obj.isgia),'">',' <section name="transient" />',' <help> indicates whether a postglacial rebound model is used in the transient </help>',' </parameter>'); 28 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="isdamage" type="',class(obj.isdamage),'" default="',convert2str(obj.isdamage),'">',' <section name="transient" />',' <help> indicates whether damage evolution is used in the transient </help>',' </parameter>'); 29 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="islevelset" type="',class(obj.islevelset),'" default="',convert2str(obj.islevelset),'">',' <section name="transient" />',' <help> LEVEL SET DESCRIPTION... </help>',' </parameter>'); 30 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="requested_outputs" type="',class(obj.requested_outputs),'" default="',convert2str(obj.requested_outputs),'">',' <section name="transient" />',' <help> list of additional outputs requested </help>',' </parameter>'); 31 fprintf(fid,'%s\n%s\n','</frame>'); 32 33 end % }}} 18 34 function obj = transient(varargin) % {{{ 19 35 switch nargin -
issm/trunk-jpl/src/m/classes/verbose.m
r15427 r17720 38 38 %}}} 39 39 methods 40 function createxml(obj,fid) % {{{ 41 fprintf(fid, '\n\n'); 42 fprintf(fid, '%s\n', '<!-- verbose -->'); 43 44 fprintf(fid,'%s%s%s\n%s\n%s%s%s\n%s\n%s\n', '<parameter key ="mprocessor" type="', class(obj.mprocessor),'" optional="false">', ' <section name="verbose" />',' <option value="',convert2str(obj.mprocessor),'" type="string" default="true"></option>',' <help> </help>','</parameter>'); 45 fprintf(fid,'%s%s%s\n%s\n%s%s%s\n%s\n%s\n', '<parameter key ="module" type="', class(obj.module),'" optional="false">', ' <section name="verbose" />',' <option value="',convert2str(obj.module),'" type="string" default="true"></option>',' <help> </help>','</parameter>'); 46 fprintf(fid,'%s%s%s\n%s\n%s%s%s\n%s\n%s\n', '<parameter key ="solution" type="', class(obj.solution),'" optional="false">', ' <section name="verbose" />',' <option value="',convert2str(obj.solution),'" type="string" default="true"></option>',' <help> </help>','</parameter>'); 47 fprintf(fid,'%s%s%s\n%s\n%s%s%s\n%s\n%s\n', '<parameter key ="solver" type="', class(obj.solver),'" optional="false">', ' <section name="verbose" />',' <option value="',convert2str(obj.solver),'" type="string" default="true"></option>',' <help> </help>','</parameter>'); 48 fprintf(fid,'%s%s%s\n%s\n%s%s%s\n%s\n%s\n', '<parameter key ="convergence" type="', class(obj.convergence),'" optional="false">', ' <section name="verbose" />',' <option value="',convert2str(obj.convergence),'" type="string" default="true"></option>',' <help> </help>','</parameter>'); 49 fprintf(fid,'%s%s%s\n%s\n%s%s%s\n%s\n%s\n', '<parameter key ="control" type="', class(obj.control),'" optional="false">', ' <section name="verbose" />',' <option value="',convert2str(obj.control),'" type="string" default="true"></option>',' <help> </help>','</parameter>'); 50 fprintf(fid,'%s%s%s\n%s\n%s%s%s\n%s\n%s\n', '<parameter key ="qmu" type="', class(obj.qmu),'" optional="false">', ' <section name="verbose" />',' <option value="',convert2str(obj.qmu),'" type="string" default="true"></option>',' <help> </help>','</parameter>'); 51 fprintf(fid,'%s%s%s\n%s\n%s%s%s\n%s\n', '<parameter key ="autodiff" type="', class(obj.autodiff),'" optional="false">', ' <section name="verbose" />',' <option value="',convert2str(obj.autodiff),'" type="string" default="true"></option>',' <help> </help>','</parameter>'); 52 53 54 end % }}} 40 55 function verbose=verbose(varargin) % {{{ 41 56
Note:
See TracChangeset
for help on using the changeset viewer.