Changeset 9821


Ignore:
Timestamp:
09/16/11 12:15:20 (14 years ago)
Author:
seroussi
Message:

almost finished displays

Location:
issm/trunk/src/m/classes
Files:
15 edited

Legend:

Unmodified
Added
Removed

  • issm/trunk/src/m/classes/inversion.m

    r9820 r9821  
    161161                        fielddisplay(obj,'vel_obs','observed velocity magnitude [m/a]');
    162162                        fielddisplay(obj,'thickness_obs','observed thickness [m]');
    163                         disp('Available responses:');
     163                        disp('Available cost functions:');
    164164                        disp('   101: SurfaceAbsVelMisfit');
    165165                        disp('   102: SurfaceRelVelMisfit');

  • issm/trunk/src/m/classes/mesh.m

    r9783 r9821  
    66classdef mesh
    77        properties (SetAccess=public)
    8                 x                 = modelfield('default',NaN,'marshall',true,'format','DoubleMat','mattype',1);
    9                 y                 = modelfield('default',NaN,'marshall',true,'format','DoubleMat','mattype',1);
    10                 z                 = modelfield('default',NaN,'marshall',true,'format','DoubleMat','mattype',1);
    11                 elements          = modelfield('default',NaN,'marshall',true,'format','DoubleMat','mattype',2);
    12                 dimension         = modelfield('default',0,'marshall',true,'format','Integer');
    13                 numberoflayers    = modelfield('default',0,'marshall',true,'format','Integer');
    14                 numberofelements  = modelfield('default',0,'marshall',true,'format','Integer');
    15                 numberofvertices  = modelfield('default',0,'marshall',true,'format','Integer');
    16                 numberofedges     = modelfield('default',0,'marshall',true,'format','Integer');
     8                x                           = modelfield('default',NaN,'marshall',true,'format','DoubleMat','mattype',1);
     9                y                           = modelfield('default',NaN,'marshall',true,'format','DoubleMat','mattype',1);
     10                z                           = modelfield('default',NaN,'marshall',true,'format','DoubleMat','mattype',1);
     11                elements                    = modelfield('default',NaN,'marshall',true,'format','DoubleMat','mattype',2);
     12                dimension                   = modelfield('default',0,'marshall',true,'format','Integer');
     13                numberoflayers              = modelfield('default',0,'marshall',true,'format','Integer');
     14                numberofelements            = modelfield('default',0,'marshall',true,'format','Integer');
     15                numberofvertices            = modelfield('default',0,'marshall',true,'format','Integer');
     16                numberofedges               = modelfield('default',0,'marshall',true,'format','Integer');
    1717               
    18                 lat        = modelfield('default',NaN,'marshall',false);
    19                 long       = modelfield('default',NaN,'marshall',false);
    20                 hemisphere = modelfield('default',NaN,'marshall',false);
    21 
    22                 elementonbed      = modelfield('default',NaN,'marshall',true,'format','BooleanMat','mattype',2);
    23                 elementonsurface  = modelfield('default',NaN,'marshall',true,'format','BooleanMat','mattype',2);
    24                 vertexonbed       = modelfield('default',NaN,'marshall',true,'format','BooleanMat','mattype',1);
    25                 vertexonsurface   = modelfield('default',NaN,'marshall',true,'format','BooleanMat','mattype',1);
    26                 lowerelements     = modelfield('default',NaN,'marshall',true,'format','DoubleMat','mattype',2);
    27                 lowervertex       = modelfield('default',NaN,'marshall',false);
    28                 upperelements     = modelfield('default',NaN,'marshall',true,'format','DoubleMat','mattype',2);
    29                 uppervertex       = modelfield('default',NaN,'marshall',false);
    30                 vertexonboundary  = modelfield('default',NaN,'marshall',false);
    31 
    32                 edges               = modelfield('default',NaN,'marshall',true,'format','DoubleMat','mattype',3);
    33                 segments            = modelfield('default',NaN,'marshall',false);
    34                 segmentmarkers      = modelfield('default',NaN,'marshall',false);
    35                 vertexconnectivity  = modelfield('default',NaN,'marshall',false);
    36                 elementconnectivity = modelfield('default',NaN,'marshall',true,'format','DoubleMat','mattype',3);
     18                lat                         = modelfield('default',NaN,'marshall',false);
     19                long                        = modelfield('default',NaN,'marshall',false);
     20                hemisphere                  = modelfield('default',NaN,'marshall',false);
     21
     22                elementonbed                = modelfield('default',NaN,'marshall',true,'format','BooleanMat','mattype',2);
     23                elementonsurface            = modelfield('default',NaN,'marshall',true,'format','BooleanMat','mattype',2);
     24                vertexonbed                 = modelfield('default',NaN,'marshall',true,'format','BooleanMat','mattype',1);
     25                vertexonsurface             = modelfield('default',NaN,'marshall',true,'format','BooleanMat','mattype',1);
     26                lowerelements               = modelfield('default',NaN,'marshall',true,'format','DoubleMat','mattype',2);
     27                lowervertex                 = modelfield('default',NaN,'marshall',false);
     28                upperelements               = modelfield('default',NaN,'marshall',true,'format','DoubleMat','mattype',2);
     29                uppervertex                 = modelfield('default',NaN,'marshall',false);
     30                vertexonboundary            = modelfield('default',NaN,'marshall',false);
     31
     32                edges                       = modelfield('default',NaN,'marshall',true,'format','DoubleMat','mattype',3);
     33                segments                    = modelfield('default',NaN,'marshall',false);
     34                segmentmarkers              = modelfield('default',NaN,'marshall',false);
     35                vertexconnectivity          = modelfield('default',NaN,'marshall',false);
     36                elementconnectivity         = modelfield('default',NaN,'marshall',true,'format','DoubleMat','mattype',3);
    3737                average_vertex_connectivity = modelfield('default',0,'marshall',true,'format','Integer');
    3838
    39                 x2d                = modelfield('default',NaN,'marshall',false);
    40                 y2d                = modelfield('default',NaN,'marshall',false);
    41                 elements2d         = modelfield('default',NaN,'marshall',true,'format','DoubleMat','mattype',3);
    42                 numberofvertices2d = modelfield('default',0,'marshall',true,'format','Integer');
    43                 numberofelements2d = modelfield('default',0,'marshall',true,'format','Integer');
    44 
    45                 extractedvertices  = modelfield('default',NaN,'marshall',false);
    46                 extractedelements  = modelfield('default',NaN,'marshall',false);
     39                x2d                         = modelfield('default',NaN,'marshall',false);
     40                y2d                         = modelfield('default',NaN,'marshall',false);
     41                elements2d                  = modelfield('default',NaN,'marshall',true,'format','DoubleMat','mattype',3);
     42                numberofvertices2d          = modelfield('default',0,'marshall',true,'format','Integer');
     43                numberofelements2d          = modelfield('default',0,'marshall',true,'format','Integer');
     44
     45                extractedvertices           = modelfield('default',NaN,'marshall',false);
     46                extractedelements           = modelfield('default',NaN,'marshall',false);
    4747        end
    4848        methods
    … …  
    170170                        fielddisplay(obj,'z','vertices z coordinate');
    171171                        fielddisplay(obj,'edges','edges of the 2d mesh (vertex1 vertex2 element1 element2)');
     172                        fielddisplay(obj,'numberofedges','number of edges of the 2d mesh');
    172173
    173174                        disp(sprintf('\n      Properties:'));
    … …  
    178179                        fielddisplay(obj,'vertexonsurface','upper vertices flags list');
    179180                        fielddisplay(obj,'elementonsurface','upper elements flags list');
     181                        fielddisplay(obj,'uppervertex','upper vertex list (NaN for vertex on the upper surface)');
     182                        fielddisplay(obj,'upperelement','upper element list (NaN for element on the upper layer)');
     183                        fielddisplay(obj,'lowervertex','lower vertex list (NaN for vertex on the lower surface)');
     184                        fielddisplay(obj,'lowerelement','lower element list (NaN for element on the lower layer');
     185                        fielddisplay(obj,'vertexonboundary','vertices on the boundary of the domain flag list');
     186                       
     187                        fielddisplay(obj,'segments','edges on domain boundary (vertex1 vertex2 element)');
     188                        fielddisplay(obj,'segmentsmarkers','number associated to each segment');
     189                        fielddisplay(obj,'vertexconnectivity','list of vertices connected to vertex_i');
     190                        fielddisplay(obj,'elementconnectivity','list of vertices connected to element_i');
     191                        fielddisplay(obj,'average_vertex_connectivity','average number of vertices connected to one vertex');
     192
     193                        disp(sprintf('\n      Extracted model:'));
     194                        fielddisplay(obj,'extractedvertices','vertices extracted from the model');
     195                        fielddisplay(obj,'extractedelements','elements extracted from the model');
     196
     197                        disp(sprintf('\n      Projection:'));
     198                        fielddisplay(obj,'lat','vertices latitude');
     199                        fielddisplay(obj,'long','vertices longitude');
     200                        fielddisplay(obj,'hemisphere','Indicate hemisphere ''n'' or ''s'' ');
    180201                end % }}}
    181202        end

  • issm/trunk/src/m/classes/miscellaneous.m

    r9732 r9821  
    4141                        end
    4242                end % }}}
     43                function disp(obj) % {{{
     44                        disp(sprintf('   miscellaneous parameters:'));
     45
     46                        fielddisplay(obj,'notes','notes in a cell of strings');
     47                        fielddisplay(obj,'name','model name');
     48                        fielddisplay(obj,'dummy','empty field to store some data');
     49
     50                end % }}}
    4351        end
    4452end

  • issm/trunk/src/m/classes/private.m

    r9739 r9821  
    3939
    4040                end % }}}
     41                function disp(obj) % {{{
     42                        disp(sprintf('   private parameters: do not change'));
     43
     44                        fielddisplay(obj,'runtimename','name of the run launched');
     45                        fielddisplay(obj,'bamg','structure with mesh properties construced if bamg is used to mesh the domain');
     46                        fielddisplay(obj,'solution','type of solution launched');
     47
     48                end % }}}
    4149        end
    4250end

  • issm/trunk/src/m/classes/prognostic.m

    r9751 r9821  
    5757                function disp(obj) % {{{
    5858                        disp(sprintf('   Prognostic solution parameters:'));
    59 
    60                         disp(sprintf('\n      transient:'));
     59                        fielddisplay(obj,'spcthickness','thickness constraints (NaN means no constraint)');
     60                        fielddisplay(obj,'hydrostatic_adjustment','adjustment of ice shelves surface and bed elevations: ''Incremental'' or ''Absolute'' ');
    6161                        fielddisplay(obj,'stabilization','0->no, 1->artificial_diffusivity, 3->discontinuous Galerkin');
    6262
    63                         disp(sprintf('\n      boundary conditions:'));
    64                         fielddisplay(obj,'spcthickness','thickness constraints (NaN means no constraint)');
     63                        disp(sprintf('\n      %s','Penalty options:'));
     64                        fielddisplay(obj,'penalty_factor','offset used by penalties: penalty = Kmax*10^offset');
     65                        fielddisplay(obj,'vertex_pairing','pairs of vertices that are penalized');
    6566
    6667                end % }}}

  • issm/trunk/src/m/classes/qmu.m

    r9752 r9821  
    66classdef qmu
    77        properties (SetAccess=public)
     8                isdakota                    = modelfield('default',0,'marshall',true,'format','Boolean');
    89                variables                   = modelfield('default',struct(),'marshall',false);
    910                responses                   = modelfield('default',struct(),'marshall',false);
    … …  
    1112                params                      = modelfield('default',struct(),'marshall',false);
    1213                results                     = modelfield('default',struct(),'marshall',false);
    13                 isdakota                    = modelfield('default',0,'marshall',true,'format','Boolean');
    1414                partition                   = modelfield('default',NaN,'marshall',true,'format','DoubleMat','mattype',2);
    1515                numberofpartitions          = modelfield('default',0,'marshall',true,'format','Integer');
    … …  
    8787                end % }}}
    8888                function disp(obj) % {{{
     89                        disp(sprintf('   qmu parameters:'));
    8990
    9091                        fielddisplay(obj,'isdakota','is qmu analysis activated?');
    91 
    9292                        for i=1:numel(obj.variables)
    9393                                disp(sprintf('         variables%s:  (arrays of each variable class)',...
    … …  
    104104                                end
    105105                        end
    106 
    107 
    108106                        for i=1:numel(obj.responses)
    109107                                disp(sprintf('         responses%s:  (arrays of each response class)',...
    … …  
    120118                                end
    121119                        end
    122 
    123 
    124                         disp(sprintf('         qmu_method:  (array of dakota_method class)'));
     120                        fielddisplay(obj,'numberofresponses','number of responses')
    125121                        for i=1:numel(obj.method);
    126122                                if strcmp(class(obj.method(i)),'dakota_method')
    … …  
    129125                                end
    130126                        end
    131 
    132127                        for i=1:numel(obj.params)
    133128                                disp(sprintf('         qmu_params%s:  (array of method-independent parameters)',...
    … …  
    144139                                end
    145140                        end
    146 
    147141                        for i=1:numel(obj.results)
    148142                                disp(sprintf('         results%s:  (information from dakota files)',...
    … …  
    159153                                end
    160154                        end
     155                        fielddisplay(obj,'partition','user provided mesh partitionition, defaults to metis if not specified')
     156                        fielddisplay(obj,'numberofpartitions','number of partitions for semi-descrete qmu')
     157                        fielddisplay(obj,'variabledescriptors','');
     158                        fielddisplay(obj,'responsedescriptors','');
     159                        fielddisplay(obj,'method','array of dakota_method class');
     160                        fielddisplay(obj,'mass_flux_profile_directory','directory for mass flux profiles');
     161                        fielddisplay(obj,'mass_flux_profiles','list of mass_flux profiles');
     162                        fielddisplay(obj,'mass_flux_segments','');
     163                        fielddisplay(obj,'adjacency','');
     164                        fielddisplay(obj,'vertex_weight','weight applied to each mesh vertex');
    161165
    162                         disp(sprintf('         numberofpartitions   : %i (number of partitions for semi-descrete qmu)',obj.numberofpartitions));
    163                         disp(sprintf('         partition    : [%i] (user provided mesh partitionition, defaults to metis if not specified)',length(obj.partition)));
    164166                end % }}}
    165167        end

  • issm/trunk/src/m/classes/radaroverlay.m

    r9621 r9821  
    3535                        end
    3636                end % }}}
     37                function disp(obj) % {{{
     38                        disp(sprintf('   radaroverlay parameters:'));
     39
     40                        fielddisplay(obj,'pwr','radar power image (matrix)');
     41                        fielddisplay(obj,'x','corresponding x coordinates');
     42                        fielddisplay(obj,'y','corresponding y coordinates');
     43
     44                end % }}}
    3745        end
    3846end

  • issm/trunk/src/m/classes/rifts.m

    r9798 r9821  
    5555                        end
    5656                end % }}}
     57                function disp(obj) % {{{
     58                        disp(sprintf('   rifts parameters:'));
     59
     60                        fielddisplay(obj,'numrifts','number of of rifts');
     61                        fielddisplay(obj,'riftstruct','structure containing all rift information (vertices coordinates, segments, type of melange, ...)');
     62                        fielddisplay(obj,'riftproperties','');
     63
     64                end % }}}
    5765        end
    5866end

  • issm/trunk/src/m/classes/settings.m

    r9749 r9821  
    5959
    6060                end % }}}
     61                function disp(obj) % {{{
     62                        disp(sprintf('   general settings parameters:'));
     63
     64                        fielddisplay(obj,'io_gather','I/O gathering strategy for result outputs (default 1)');
     65                        fielddisplay(obj,'lowmem','is the memory limited ? (0 or 1)');
     66                        fielddisplay(obj,'results_on_vertices','provide results on vertices instead of patch (0 or 1)');
     67                        fielddisplay(obj,'output_frequency','frequency at which results are saved in all solutions with multiple time_steps');
     68                        fielddisplay(obj,'waitonlock','maximum number of minutes to wait for batch results, or return 0');
     69
     70                end % }}}
    6171        end
    6272end

  • issm/trunk/src/m/classes/solver.m

    r9739 r9821  
    99         end
    1010         methods
    11                  function disp(obj) % {{{1
    12                          %  display the object
    13                          disp(sprintf('class ''%s'' object ''%s'' = \n',class(o),inputname(1)));
    14                          %  display the options
    15                          for i=1:size(obj.options,1),
    16                                  analysis=obj.options{i,1};
    17                                  ioptions=obj.options{i,2};
    18                                  string=PetscString(ioptions);
    19                                  disp(sprintf('   %s -> ''%s''',EnumToString(analysis),string));
    20                          end
    21                          disp(sprintf('\n'));
    22                  end
    23                  %}}}
    2411                 function obj=addoptions(obj,analysis,solveroptions) % {{{1
    2512                         %first, find out if analysis has already been supplied
    … …  
    119106                         end
    120107                 end % }}}
     108                 function disp(obj) % {{{1
     109                        disp(sprintf('   solver parameters:'));
     110
     111                        for i=1:size(obj.options,1),
     112                                analysis=obj.options{i,1};
     113                                ioptions=obj.options{i,2};
     114                                string=PetscString(ioptions);
     115                                disp(sprintf('   %s -> ''%s''',EnumToString(analysis),string));
     116                        end
     117                 end
     118                 %}}}
    121119         end
    122120 end

  • issm/trunk/src/m/classes/steadystate.m

    r9754 r9821  
    5454                        end
    5555                end % }}}
     56                function disp(obj) % {{{
     57                        disp(sprintf('   steadystate solution parameters:'));
     58
     59                        fielddisplay(obj,'reltol','relative tolerance criterion');
     60                        fielddisplay(obj,'maxiter','maximum number of iterations');
     61                        fielddisplay(obj,'requested_outputs','additional requested outputs');
     62
     63                end % }}}
    5664        end
    5765end

  • issm/trunk/src/m/classes/surfaceforcings.m

    r9739 r9821  
    3838
    3939                end % }}}
     40                function disp(obj) % {{{
     41                        disp(sprintf('   surface forcings parameters:'));
     42
     43                        fielddisplay(obj,'accumulation_rate','surface accumulation rate [m]');
     44                        fielddisplay(obj,'ablation_rate','surface ablation rate [m]');
     45                        fielddisplay(obj,'mass_balance','surface mass balance [m]');
     46
     47                end % }}}
    4048        end
    4149end

  • issm/trunk/src/m/classes/thermal.m

    r9751 r9821  
    5959                        disp(sprintf('   Thermal solution parameters:'));
    6060
    61                         disp(sprintf('\n      parameters:'));
     61                        fielddisplay(obj,'spctemperature','temperature constraints (NaN means no constraint)');
     62                        fielddisplay(obj,'stabilization','0->no, 1->artificial_diffusivity, 2->SUPG');
     63                        fielddisplay(obj,'maxiter','maximum number of non linear iterations');
    6264                        fielddisplay(obj,'penalty_lock','stabilize unstable thermal constraints that keep zigzagging after n iteration (default is 0, no stabilization)');
    6365                        fielddisplay(obj,'penalty_threshold','threshold to declare convergence of thermal solution (default is 0)');
    64 
    65                         disp(sprintf('\n      boundary conditions:'));
    66                         fielddisplay(obj,'spctemperature','constraints flag list (first column) and values (second column)');
    6766
    6867                end % }}}

  • issm/trunk/src/m/classes/timestepping.m

    r9755 r9821  
    5757                        end
    5858                end % }}}
     59                function disp(obj) % {{{
     60                        disp(sprintf('   timestepping parameters:'));
     61
     62                        fielddisplay(obj,'time_step','length of time steps [yrs]');
     63                        fielddisplay(obj,'final_time','final time to stop the simulation [yrs]');
     64                        fielddisplay(obj,'time_adapt','use cfl condition to define time step ? (0 or 1) ');
     65                        fielddisplay(obj,'cfl_coefficient','coefficient applied to cfl condition');
     66
     67                end % }}}
    5968        end
    6069end

  • issm/trunk/src/m/classes/transient.m

    r9752 r9821  
    4444
    4545                end % }}}
    46                 function disp(obj) % {{{
    47                         disp('Transient parameters on a pure solution basis:');
    48                         fielddisplay(obj,'isprognostic','indicates if a prognostic solution is used in the transient');
    49                         fielddisplay(obj,'isthermal','indicates if a thermal solution is used in the transient');
    50                         fielddisplay(obj,'isdiagnostic','indicates if a diagnostic solution is used in the transient');
    51                         fielddisplay(obj,'isgroundingline','indicates if a groundingline migration is used in the transient');
    52 
    53                 end % }}}
    5446                function checkconsistency(obj,md) % {{{
    5547                        if ~ismember(md.transient.isdiagnostic,[0 1]),
    … …  
    6759
    6860                end % }}}
     61                function disp(obj) % {{{
     62                        disp(sprintf('   balance thickness solution parameters:'));
     63
     64                        fielddisplay(obj,'isprognostic','indicates if a prognostic solution is used in the transient');
     65                        fielddisplay(obj,'isthermal','indicates if a thermal solution is used in the transient');
     66                        fielddisplay(obj,'isdiagnostic','indicates if a diagnostic solution is used in the transient');
     67                        fielddisplay(obj,'isgroundingline','indicates if a groundingline migration is used in the transient');
     68
     69                end % }}}
    6970        end
    7071end
Note: See TracChangeset for help on using the changeset viewer.