Changeset 17931
- Timestamp:
- 05/05/14 14:11:31 (11 years ago)
- Location:
- issm/trunk-jpl/src/m/classes
- Files:
-
- 20 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/m/classes/balancethickness.m
r17720 r17931 17 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 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>'); 19 20 % balancethickness drop-down (1,2, or 3) 21 fprintf(fid,'%s\n%s\n%s\n', '<parameter key ="stabilization" type="alternative" optional="false">',' <section name="balancethickness" />',' <help> 0: None, 1: SU, 2: SSAs artificial diffusivity, 3:DG </help>'); 22 fprintf(fid,'%s\n',' <option value="1" type="string" default="true"> </option>'); 23 fprintf(fid,'%s\n',' <option value="2" type="string" default="false"> </option>'); 24 fprintf(fid,'%s\n%s\n',' <option value="3" type="string" default="false"> </option>','</parameter>'); 20 25 21 26 end % }}} -
issm/trunk-jpl/src/m/classes/damage.m
r17906 r17931 72 72 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>'); 73 73 74 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>'); 74 % stabilization (0,1, or 2) drop-down 75 fprintf(fid,'%s\n%s\n%s\n%s\n', '<parameter key ="stabilization" type="alternative" optional="false">',' <section name="damage" />',' <help> 0: no, 1: artificial_diffusivity, 2: SUPG </help>'); 76 fprintf(fid, '%s\n', ' <option value="0" type="string" default="true"></option>'); 77 fprintf(fid, '%s\n', ' <option value="1" type="string" default="false"></option>'); 78 fprintf(fid, '%s\n', ' <option value="2" type="string" default="false"></option>'); 79 fprintf(fid, '%s\n','</parameter>'); 80 75 81 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>'); 76 82 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>'); -
issm/trunk-jpl/src/m/classes/debug.m
r17720 r17931 14 14 fprintf(fid, '\n\n'); 15 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>'); 16 17 %valgrind drop-down (0 or 1) 18 fprintf(fid,'%s\n%s\n%s%s%s\n%s\n%s\n', '<parameter key ="valgrind" type="alternative" optional="false">',' <section name="debug" />',' <help> use Valgrind to debug (0 or 1) </help>'); 19 fprintf(fid,'%s\n',' <option value="0" type="string" default="true"> </option>'); 20 fprintf(fid,'%s\n%s\n',' <option value="1" type="string" default="false"> </option>','</parameter>'); 21 22 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>'); 23 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 24 end % }}} 20 25 function obj = debug(varargin) % {{{ -
issm/trunk-jpl/src/m/classes/flowequation.m
r17720 r17931 70 70 71 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>'); 72 73 % fe_SSA drop-down (P1, P1bubble, P1bubblecondensed, P2) 74 fprintf(fid,'%s\n%s\n%s\n%s\n', '<parameter key ="fe_SSA" type="alternative" optional="false">',' <section name="flowequation" />',' <help> Finite Element for SSA "P1", "P1bubble" "P1bubblecondensed" "P2" </help>'); 75 fprintf(fid,'%s\n',' <option value="P1" type="string" default="true"> </option>'); 76 fprintf(fid,'%s\n',' <option value="P1bubble" type="string" default="false"> </option>'); 77 fprintf(fid,'%s\n',' <option value="P1bubblecondensed" type="string" default="false"> </option>'); 78 fprintf(fid,'%s\n%s\n',' <option value="P2" type="string" default="false"> </option>','</parameter>'); 79 80 %fe_HO drop-down (P1, P1bubble, P1bubblecondensed, P1xP2, P2xP1, P2) 81 fprintf(fid,'%s\n%s\n%s\n%s\n', '<parameter key ="fe_HO" type="alternative" optional="false">',' <section name="flowequation" />',' <help> Finite Element for HO "P1" "P1bubble" "P1bubblecondensed" "P1xP2" "P2xP1" "P2" </help>'); 82 fprintf(fid,'%s\n',' <option value="P1" type="string" default="true"> </option>'); 83 fprintf(fid,'%s\n',' <option value="P1bubble" type="string" default="false"> </option>'); 84 fprintf(fid,'%s\n',' <option value="P1bubblecondensed" type="string" default="false"> </option>'); 85 fprintf(fid,'%s\n',' <option value="P1xP2" type="string" default="false"> </option>'); 86 fprintf(fid,'%s\n',' <option value="P2xP1" type="string" default="false"> </option>'); 87 fprintf(fid,'%s\n%s\n',' <option value="P2" type="string" default="false"> </option>','</parameter>'); 88 89 % fe_FS drop-down (P1P1, P1P1GLS, MINIcondensed, MINI, TaylowHood) 90 fprintf(fid,'%s\n%s\n%s\n%s\n', '<parameter key ="fe_FS" type="alternative" optional="false">',' <section name="flowequation" />',' <help> Finite Element for FS "P1P1" (debugging only) "P1P1GLS" "MINIcondensed" "MINI" "TaylorHood" </help>'); 91 fprintf(fid,'%s\n',' <option value="P1P1" type="string" default="true"> </option>'); 92 fprintf(fid,'%s\n',' <option value="P1P1GLS" type="string" default="false"> </option>'); 93 fprintf(fid,'%s\n',' <option value="MINIcondensed" type="string" default="false"> </option>'); 94 fprintf(fid,'%s\n',' <option value="MINI" type="string" default="false"> </option>'); 95 fprintf(fid,'%s\n%s\n',' <option value="TaylorHood" type="string" default="false"> </option>','</parameter>'); 96 75 97 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 98 -
issm/trunk-jpl/src/m/classes/geometry.m
r17720 r17931 33 33 end 34 34 methods 35 function createxml(obj,fid) % {{{ 35 %function createxml(obj,fid) 36 function createxml(obj, fid)% {{{ 36 37 fprintf(fid, '\n\n'); 37 38 fprintf(fid, '%s\n', '<!-- geometry -->'); … … 40 41 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 42 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%s\n','<parameter key ="bathymetry" type="',class(obj.bathymetry),'" default="',convert2str(obj.bathymetry),'">',' <section name="geometry" />',' <help> bathymetry elevation [m] </help>','</parameter>'); 43 44 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 45 fprintf(fid,'%s\n%s\n','</frame>'); -
issm/trunk-jpl/src/m/classes/gia.m
r17720 r17931 18 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 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>'); 20 %cross_section_shape drop-down (1 o r 2) 21 fprintf(fid,'%s\n%s\n%s\n%s\n',' <parameter key ="cross_section_shape" type="alternative" optional="false">',' <section name="gia" />',' <help> 1: square-edged (default). 2: elliptical. See iedge in GiaDeflectionCore </help>'); 22 fprintf(fid,'%s\n',' <option value="1" type="string" default="true"> </option>'); 23 fprintf(fid,'%s\n%s\n',' <option value="2" type="string" default="false"> </option>','</parameter>'); 24 21 25 fprintf(fid,'%s\n%s\n','</frame>'); 22 26 -
issm/trunk-jpl/src/m/classes/groundingline.m
r17720 r17931 16 16 % Convergence criteria 17 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>'); 18 19 % migration (SoftMigration, AggressiveMigration, or None)drop-down 20 fprintf(fid,'%s\n%s\n%s\n%s\n','<parameter key ="migration" type="alternative" optional="false">',' <section name="groundingline" />',' <help> type of grounding line migration: "SoftMigration","AgressiveMigration" or "None" </help>'); 21 fprintf(fid,'%s\n',' <option value="SoftMigration" type="string" default="true"> </option>'); 22 fprintf(fid,'%s\n',' <option value="AggressiveMigration" type="string" default="false"> </option>'); 23 fprintf(fid, '%s\n%s\n',' <option value="None" type="string" default="false"></option>','</parameter>'); 24 19 25 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 26 fprintf(fid,'%s\n%s\n','</frame>'); -
issm/trunk-jpl/src/m/classes/inversion.m
r17920 r17931 25 25 end 26 26 methods 27 function createxml(obj,fid) % {{{ 28 fprintf(fid, '<!-- inversion -->\n'); 29 30 % inversion parameters 31 fprintf(fid,'%s\n%s\n%s\n','<frame key="1" label="inversion parameters">','<section name="inversion" />'); 32 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>'); 33 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>'); 34 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>'); 35 36 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>'); 37 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>'); 38 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>'); 39 40 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>'); 41 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>'); 42 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>'); 43 44 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>'); 45 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>'); 46 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>'); 47 48 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>'); 49 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>'); 50 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>'); 51 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>'); 52 53 fprintf(fid,'%s\n%s\n','</frame>'); 54 55 fprintf(fid,'%s\n%s\n%s\n','<frame key="2" label="Available cost functions">','<section name="inversion" />'); 56 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>'); 57 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>'); 58 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>'); 59 60 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>'); 61 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>'); 62 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>'); 63 64 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>'); 65 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>'); 66 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>'); 67 68 fprintf(fid,'%s\n%s\n','</frame>'); 69 70 end % }}} 27 function createxml(obj,fid) % {{{ 28 fprintf(fid, '<!-- inversion -->\n'); 29 30 % inversion parameters 31 fprintf(fid,'%s\n%s\n%s\n','<frame key="1" label="inversion parameters">','<section name="inversion" />'); 32 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>'); 33 34 % incompleteadjoing drop-down (0 or 1) 35 fprintf(fid,'%s\n%s\n%s\n%s\n',' <parameter key ="incomplete_adjoint" type="alternative" optional="false">',' <section name="inversion" />',' <help> 1: linear viscosity, 0: non-linear viscosity </help>'); 36 fprintf(fid,'%s\n',' <option value="0" type="string" default="true"> </option>'); 37 fprintf(fid,'%s\n%s\n',' <option value="1" type="string" default="false"> </option>','</parameter>'); 38 39 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>'); 40 41 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>'); 42 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>'); 43 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>'); 44 45 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>'); 46 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>'); 47 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>'); 48 49 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>'); 50 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>'); 51 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>'); 52 53 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>'); 54 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>'); 55 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>'); 56 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>'); 57 58 fprintf(fid,'%s\n%s\n','</frame>'); 59 60 fprintf(fid,'%s\n%s\n%s\n','<frame key="2" label="Available cost functions">','<section name="inversion" />'); 61 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>'); 62 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>'); 63 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>'); 64 65 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>'); 66 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>'); 67 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>'); 68 69 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>'); 70 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>'); 71 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>'); 72 73 fprintf(fid,'%s\n%s\n','</frame>'); 74 75 end % }}} 71 76 function obj = inversion(varargin) % {{{ 72 77 switch nargin -
issm/trunk-jpl/src/m/classes/masstransport.m
r17720 r17931 52 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 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>'); 54 55 % hydrostatic_adjustment drop-down (incremental or absolute ) 56 fprintf(fid,'%s\n%s\n%s\n%s\n',' <parameter key ="hydrostatic_adjustment" type="alternative" optional ="false">',' <section name="masstransport" />',' <help> adjustment of ice shelves surface and bed elevations: ''Incremental'' or ''Absolute'' </help>'); 57 fprintf(fid,'%s\n',' <option value="Incremental" type="string" default="true"> </option>'); 58 fprintf(fid,'%s\n%s\n',' <option value="Absolute" type="string" default="false"> </option>','</parameter>'); 59 60 %stabilization drop-down (0, 1, 2, or 3) 61 fprintf(fid,'%s\n%s\n%s\n%s\n',' <parameter key ="stabilization" type="alternative" optional="false">',' <section name="masstransport" />',' <help> 0: no, 1: artificial_diffusivity, 2: streamline upwinding, 3: discontinuous Galerkin </help>'); 62 fprintf(fid,'%s\n',' <option value="0" type="string" default="true"> </option>'); 63 fprintf(fid,'%s\n',' <option value="1" type="string" default="false"> </option>'); 64 fprintf(fid,'%s\n',' <option value="2" type="string" default="false"> </option>'); 65 fprintf(fid,'%s\n%s\n',' <option value="3" type="string" default="false"> </option>','</parameter>'); 66 56 67 fprintf(fid,'%s\n%s\n','</frame>'); 57 68 -
issm/trunk-jpl/src/m/classes/matice.m
r17720 r17931 55 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 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>'); 57 58 % rheology_law drop-down 59 fprintf(fid,'%s\n%s\n%s\n%s\n','<parameter key ="rheology_law" type="alternative" optional="false">',' <section name="materials" />',' <help> law for the temperature dependance of the rheology: "None", "Paterson", "Arrhenius" or "LliboutryDuval" </help>'); 60 fprintf(fid,'%s\n',' <option value="None" type="string" default="true"> </option>'); 61 fprintf(fid,'%s\n',' <option value="Paterson" type="string" default="false"> </option>'); 62 fprintf(fid,'%s\n',' <option value="Arrhenius" type="string" default="false"> </option>'); 63 fprintf(fid,'%s\n%s\n',' <option value="LliboutryDuval" type="string" default="false"> </option>','</parameter>'); 64 58 65 59 66 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>'); -
issm/trunk-jpl/src/m/classes/mesh2d.m
r17831 r17931 105 105 end % }}} 106 106 function createxml(obj,fid) % {{{ 107 fprintf(fid, '<!-- mesh-->\n');107 fprintf(fid, '<!-- 2D tria Mesh (horizontal) -->\n'); 108 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>'); 109 %elements and vertices 110 fprintf(fid,'%s\n%s\n%s\n','<frame key="1" label="Elements and vertices">','<section name="mesh" />'); 111 fprintf(fid,'%s%s%s%s%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>'); 112 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n',' <parameter key ="numberofvertices" type="',class(obj.numberofvertices),'" default="',convert2str(obj.numberofvertices),'">',' <section name="mesh" />',' <help> number of vertices </help>',' </parameter>'); 113 fprintf(fid,'%s%s%s%s%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 114 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 115 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 116 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 117 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 118 fprintf(fid,'%s\n%s\n','</frame>'); 146 119 147 120 % properties 148 121 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 122 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 123 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>'); … … 175 138 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 139 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>'); 140 % choice (hemisphere) 'n' or 's' 141 fprintf(fid,'%s\n%s\n%s\n',' <parameter key ="hemisphere" type="alternative" optional="false">',' <section name="mesh" />',' <help> Indicate hemisphere ''n'' or ''s'' </help>'); 142 fprintf(fid,'%s\n',' <option value="n" type="string" default="true"> </option>'); 143 fprintf(fid,'%s\n',' <option value="s" type="string" default="false"> </option>'); 144 fprintf(fid,'%s\n',' </parameter>'); 178 145 fprintf(fid,'%s\n%s\n','</frame>'); 179 146 -
issm/trunk-jpl/src/m/classes/mesh2dvertical.m
r17700 r17931 29 29 end 30 30 methods 31 function createxml(obj,fid) % {{{ 32 fprintf(fid, '<!-- 2d Mesh (Vertical) -->\n'); 33 34 %elements and vertices 35 fprintf(fid,'%s\n%s\n%s\n',' <frame key="1" label="Elements and vertices">',' <section name="mesh" />'); 36 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n',' <parameter key ="numberofelements" type="',class(obj.numberofelements),'" default="',convert2str(obj.numberofelements),'">',' <help> number of elements </help>',' </parameter>'); 37 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n',' <parameter key ="numberofvertices" type="',class(obj.numberofvertices),'" default="',convert2str(obj.numberofvertices),'">',' <help> number of vertices </help>',' </parameter>'); 38 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n',' <parameter key ="elements" type="',class(obj.elements),'" default="',convert2str(obj.elements),'">',' <help> vertex indices of the mesh elements </help>',' </parameter>'); 39 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>'); 40 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>'); 41 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>'); 42 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>'); 43 fprintf(fid,'%s\n%s\n',' </frame>'); 44 45 % properties 46 fprintf(fid,'%s\n%s\n%s\n','<frame key="2" label="Properties">','<section name="mesh" />'); 47 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>'); 48 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="vertexonbase" type="',class(obj.vertexonbase),'" default="',convert2str(obj.vertexonbase),'">',' <section name="mesh" />',' <help> vertices on the bed of the domain flag list </help>',' </parameter>'); 49 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>'); 50 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>'); 51 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>'); 52 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>'); 53 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>'); 54 fprintf(fid,'%s\n%s\n','</frame>'); 55 56 %projection 57 fprintf(fid,'%s\n%s\n%s\n','<frame key="3" label="Projection">','<section name="mesh" />'); 58 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>'); 59 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>'); 60 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>'); 61 fprintf(fid,'%s\n%s\n','</frame>'); 62 63 end % }}} 31 64 function obj = mesh2dvertical(varargin) % {{{ 32 65 switch nargin -
issm/trunk-jpl/src/m/classes/mesh3dprisms.m
r17886 r17931 42 42 end 43 43 methods 44 function createxml(obj,fid) % {{{ 45 fprintf(fid, '<!-- 3D prism Mesh -->\n'); 46 47 % Elements and verticies of the original 2d mesh 48 fprintf(fid,'%s\n%s\n%s\n',' <frame key="1" label="Elements and vertices of the orginal 2d mesh">',' <section name="mesh3dprisms" />'); 49 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>'); 50 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>'); 51 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>'); 52 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="x2d" type="',class(obj.x2d),'" default="',convert2str(obj.x2d),'">',' <section name="mesh3dprisms" />',' <help> vertices x coordinate [m] </help>',' </parameter>'); 53 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="y2d" type="',class(obj.y2d),'" default="',convert2str(obj.y2d),'">',' <section name="mesh3dprisms" />',' <help> vertices y coordinate [m] </help>',' </parameter>'); 54 fprintf(fid,'%s\n%s\n',' </frame>'); 55 56 % Elements and vertices of the extruded 3d mesh 57 fprintf(fid,'%s\n%s\n%s\n',' <frame key="2" label="Elements and vertices of the orginal 3d mesh">',' <section name="mesh3dprisms" />'); 58 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n',' <parameter key ="numberofelements" type="',class(obj.numberofelements),'" default="',convert2str(obj.numberofelements),'">',' <help> number of elements </help>',' </parameter>'); 59 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n',' <parameter key ="numberofvertices" type="',class(obj.numberofvertices),'" default="',convert2str(obj.numberofvertices),'">',' <help> number of vertices </help>',' </parameter>'); 60 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n',' <parameter key ="elements" type="',class(obj.elements),'" default="',convert2str(obj.elements),'">',' <help> vertex indices of the mesh elements </help>',' </parameter>'); 61 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="mesh3dprisms" />',' <help> vertices x coordinate [m] </help>',' </parameter>'); 62 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="mesh3dprisms" />',' <help> vertices y coordinate [m] </help>',' </parameter>'); 63 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="z" type="',class(obj.y),'" default="',convert2str(obj.y),'">',' <section name="mesh3dprisms" />',' <help> vertices z coordinate [m] </help>',' </parameter>'); 64 fprintf(fid,'%s\n%s\n',' </frame>'); 65 66 % properties 67 fprintf(fid,'%s\n%s\n%s\n','<frame key="3" label="Properties">','<section name="mesh3dprisms" />'); 68 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="mesh3dprisms" />',' <help> number of extrusion layers </help>',' </parameter>'); 69 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="vertexonbase" type="',class(obj.vertexonbase),'" default="',convert2str(obj.vertexonbase),'">',' <section name="mesh3dprisms" />',' <help> lower vertices flags list </help>',' </parameter>'); 70 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="elementontbase" type="',class(obj.elementontbase),'" default="',convert2str(obj.elementontbase),'">',' <section name="mesh3dprisms" />',' <help> lower elements flags list </help>',' </parameter>'); 71 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="mesh3dprisms" />',' <help> upper vertices flags list </help>',' </parameter>'); 72 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="elementontsurface" type="',class(obj.elementontsurface),'" default="',convert2str(obj.elementontsurface),'">',' <section name="mesh3dprisms" />',' <help> upper elements flags list </help>',' </parameter>'); 73 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="mesh3dprisms" />',' <help> upper vertex list (NaN for vertex on the upper surface) </help>',' </parameter>'); 74 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="mesh3dprisms" />',' <help> upper element list (NaN for element on the upper layer) </help>',' </parameter>'); 75 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="mesh3dprisms" />',' <help> lower vertex list (NaN for vertex on the lower surface) </help>',' </parameter>'); 76 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="mesh3dprisms" />',' <help> element list (NaN for element on the lower layer) </help>',' </parameter>'); 77 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="mesh3dprisms" />',' <help> vertices on the boundary of the domain flag list </help>',' </parameter>'); 78 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="mesh3dprisms" />',' <help> list of vertices connected to vertex_i </help>',' </parameter>'); 79 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="mesh3dprisms" />',' <help> average number of vertices connected to one vertex </help>',' </parameter>'); 80 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="mesh3dprisms" />',' <help> average number of vertices connected to one vertex </help>',' </parameter>'); 81 fprintf(fid,'%s\n%s\n','</frame>'); 82 83 % Extracted model 84 fprintf(fid,'%s\n%s\n%s\n','<frame key="4" label="Extracted Model">','<section name="mesh3dprisms" />'); 85 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="mesh3dprisms" />',' <help> vertices extracted from the model </help>',' </parameter>'); 86 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="mesh3dprisms" />',' <help> elements extracted from the model </help>',' </parameter>'); 87 fprintf(fid,'%s\n%s\n','</frame>'); 88 89 % Projection 90 fprintf(fid,'%s\n%s\n%s\n','<frame key="5" label="Projection">','<section name="mesh3dprisms" />'); 91 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="mesh3dprisms" />',' <help> vertices latitude [degrees] </help>',' </parameter>'); 92 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="mesh3dprisms" />',' <help> verticies longitude [degrees] </help>',' </parameter>'); 93 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="mesh3dprisms" />',' <help> Indicate hemisphere ''n'' or ''s'' </help>',' </parameter>'); 94 fprintf(fid,'%s\n%s\n','</frame>'); 95 96 end % }}}cd 44 97 function self = mesh3dprisms(varargin) % {{{ 45 98 switch nargin -
issm/trunk-jpl/src/m/classes/mesh3dtetras.m
r17700 r17931 40 40 end 41 41 methods 42 function createxml(obj,fid) % {{{ 43 fprintf(fid, '<!-- 3D Tetra Mesh -->\n'); 44 45 % Elements and verticies of the original 2d mesh 46 fprintf(fid,'%s\n%s\n%s\n',' <frame key="1" label="Elements and vertices of the orginal 2d mesh">',' <section name="mesh3dprisms" />'); 47 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>'); 48 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>'); 49 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>'); 50 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="x2d" type="',class(obj.x2d),'" default="',convert2str(obj.x2d),'">',' <section name="mesh3dprisms" />',' <help> vertices x coordinate [m] </help>',' </parameter>'); 51 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="y2d" type="',class(obj.y2d),'" default="',convert2str(obj.y2d),'">',' <section name="mesh3dprisms" />',' <help> vertices y coordinate [m] </help>',' </parameter>'); 52 fprintf(fid,'%s\n%s\n',' </frame>'); 53 54 % Elements and vertices of the extruded 3d mesh 55 fprintf(fid,'%s\n%s\n%s\n',' <frame key="2" label="Elements and vertices of the orginal 3d mesh">',' <section name="mesh3dprisms" />'); 56 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n',' <parameter key ="numberofelements" type="',class(obj.numberofelements),'" default="',convert2str(obj.numberofelements),'">',' <help> number of elements </help>',' </parameter>'); 57 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n',' <parameter key ="numberofvertices" type="',class(obj.numberofvertices),'" default="',convert2str(obj.numberofvertices),'">',' <help> number of vertices </help>',' </parameter>'); 58 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n',' <parameter key ="elements" type="',class(obj.elements),'" default="',convert2str(obj.elements),'">',' <help> vertex indices of the mesh elements </help>',' </parameter>'); 59 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="mesh3dprisms" />',' <help> vertices x coordinate [m] </help>',' </parameter>'); 60 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="mesh3dprisms" />',' <help> vertices y coordinate [m] </help>',' </parameter>'); 61 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="z" type="',class(obj.y),'" default="',convert2str(obj.y),'">',' <section name="mesh3dprisms" />',' <help> vertices z coordinate [m] </help>',' </parameter>'); 62 fprintf(fid,'%s\n%s\n',' </frame>'); 63 64 % properties 65 fprintf(fid,'%s\n%s\n%s\n','<frame key="3" label="Properties">','<section name="mesh3dprisms" />'); 66 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="mesh3dprisms" />',' <help> number of extrusion layers </help>',' </parameter>'); 67 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="vertexonbase" type="',class(obj.vertexonbase),'" default="',convert2str(obj.vertexonbase),'">',' <section name="mesh3dprisms" />',' <help> lower vertices flags list </help>',' </parameter>'); 68 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="elementontbase" type="',class(obj.elementontbase),'" default="',convert2str(obj.elementontbase),'">',' <section name="mesh3dprisms" />',' <help> lower elements flags list </help>',' </parameter>'); 69 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="mesh3dprisms" />',' <help> upper vertices flags list </help>',' </parameter>'); 70 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n',' <parameter key ="elementontsurface" type="',class(obj.elementontsurface),'" default="',convert2str(obj.elementontsurface),'">',' <section name="mesh3dprisms" />',' <help> upper elements flags list </help>',' </parameter>'); 71 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="mesh3dprisms" />',' <help> upper vertex list (NaN for vertex on the upper surface) </help>',' </parameter>'); 72 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="mesh3dprisms" />',' <help> upper element list (NaN for element on the upper layer) </help>',' </parameter>'); 73 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="mesh3dprisms" />',' <help> lower vertex list (NaN for vertex on the lower surface) </help>',' </parameter>'); 74 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="mesh3dprisms" />',' <help> element list (NaN for element on the lower layer) </help>',' </parameter>'); 75 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="mesh3dprisms" />',' <help> vertices on the boundary of the domain flag list </help>',' </parameter>'); 76 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="mesh3dprisms" />',' <help> list of vertices connected to vertex_i </help>',' </parameter>'); 77 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="mesh3dprisms" />',' <help> average number of vertices connected to one vertex </help>',' </parameter>'); 78 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="mesh3dprisms" />',' <help> average number of vertices connected to one vertex </help>',' </parameter>'); 79 fprintf(fid,'%s\n%s\n','</frame>'); 80 81 % Extracted model 82 fprintf(fid,'%s\n%s\n%s\n','<frame key="4" label="Extracted Model">','<section name="mesh3dprisms" />'); 83 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="mesh3dprisms" />',' <help> vertices extracted from the model </help>',' </parameter>'); 84 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="mesh3dprisms" />',' <help> elements extracted from the model </help>',' </parameter>'); 85 fprintf(fid,'%s\n%s\n','</frame>'); 86 87 % Projection 88 fprintf(fid,'%s\n%s\n%s\n','<frame key="5" label="Projection">','<section name="mesh3dprisms" />'); 89 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="mesh3dprisms" />',' <help> vertices latitude [degrees] </help>',' </parameter>'); 90 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="mesh3dprisms" />',' <help> verticies longitude [degrees] </help>',' </parameter>'); 91 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="mesh3dprisms" />',' <help> Indicate hemisphere ''n'' or ''s'' </help>',' </parameter>'); 92 fprintf(fid,'%s\n%s\n','</frame>'); 93 94 end % }}}cd 42 95 function self = mesh3dtetras(varargin) % {{{ 43 96 switch nargin -
issm/trunk-jpl/src/m/classes/model.m
r17811 r17931 1095 1095 md.verbose = verbose(); 1096 1096 md.settings = settings(); 1097 md.toolkits = toolkits();1097 %md.toolkits = toolkits(); 1098 1098 md.cluster = generic(); 1099 1099 md.balancethickness = balancethickness(); -
issm/trunk-jpl/src/m/classes/settings.m
r17720 r17931 24 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 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>'); 26 27 % lowmen drop-down (0 or 1) 28 fprintf(fid,'%s\n%s\n%s\n%s\n', '<parameter key ="lowmem" type="alternative" optional="true">', ' <section name="settings" />',' <help> is the memory limited ? (0 or 1) </help>'); 29 fprintf(fid,'%s\n',' <option value="0" type="string" default="true"> </option>'); 30 fprintf(fid,'%s\n%s\n',' <option value="1" type="string" default="false"> </option>','</parameter>'); 31 27 32 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 33 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>'); -
issm/trunk-jpl/src/m/classes/stressbalance.m
r17720 r17931 35 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 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>'); 37 38 %is newton drop-down (0,1,or 2) 39 fprintf(fid,'%s\n%s\n%s\n%s\n', '<parameter key ="isnewton" type="alternative" optional="false">', ' <section name="stressbalance" />',' <help> 0: Picards fixed point, 1: Newtons method, 2: hybrid </help>'); 40 fprintf(fid,'%s\n',' <option value="0" type="string" default="true"> </option>'); 41 fprintf(fid,'%s\n',' <option value="1" type="string" default="false"> </option>'); 42 fprintf(fid,'%s\n%s\n',' <option value="2" type="string" default="false"> </option>','</parameter>'); 43 38 44 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 45 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>'); -
issm/trunk-jpl/src/m/classes/thermal.m
r17720 r17931 22 22 % thermal solution parameters 23 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>'); 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 26 % penalty_threshold drop-down (0, 1, or 2) 27 fprintf(fid,'%s\n%s\n%s\n%s\n',' <parameter key ="penalty_threshold" type="alternative" optional="false">',' <section name="thermal" />',' <help> 0: no, 1: artificial_diffusivity, 2: SUPG </help>'); 28 fprintf(fid,'%s\n',' <option value="0" type="string" default="true"> </option>'); 29 fprintf(fid,'%s\n',' <option value="1" type="string" default="false"> </option>'); 30 fprintf(fid,'%s\n%s\n',' <option value="2" type="string" default="false"> </option>','</parameter>'); 31 32 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>'); 33 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>'); 34 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>'); 35 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>'); 36 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>'); 37 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>'); 38 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 39 fprintf(fid,'%s\n%s\n','</frame>'); 34 40 -
issm/trunk-jpl/src/m/classes/timestepping.m
r17720 r17931 19 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 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>'); 21 % time_adapt 0 or 1 drop down 22 fprintf(fid,'%s\n%s\n%s\n%s\n','<parameter key ="time_adapt" type="alternative" optional="false">',' <section name="timestepping" />',' <help> use cfl condition to define time step ? (0 or 1) </help>'); 23 fprintf(fid, '%s\n', ' <option value="0" type="string" default="true"></option>'); 24 fprintf(fid, '%s\n%s\n', ' <option value="1" type="string" default="false"></option>','</parameter>'); 22 25 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 26 -
issm/trunk-jpl/src/m/classes/transient.m
r17748 r17931 27 27 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>'); 28 28 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>'); 29 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 ="isdamageevolution" type="',class(obj.isdamageevolution),'" default="',convert2str(obj.isdamageevolution),'">',' <section name="transient" />',' <help> indicates whether damage evolution is used in the transient </help>',' </parameter>'); 30 30 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>'); 31 31 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>');
Note:
See TracChangeset
for help on using the changeset viewer.