source: issm/oecreview/Archive/18296-19100/ISSM-18823-18824.diff

Last change on this file was 19102, checked in by Mathieu Morlighem, 10 years ago

NEW: added 18296-19100
File size: 7.6 KB

  • ../trunk-jpl/src/m/classes/m1qn3inversion.m

     
    6868                        if ~IssmConfig('_HAVE_M1QN3_'),
    6969                                md = checkmessage(md,['M1QN3 has not been installed, ISSM needs to be reconfigured and recompiled with M1QN3']);
    7070                        end
    71 
    7271                        num_controls=numel(md.inversion.control_parameters);
    7372                        num_costfunc=size(md.inversion.cost_functions,2);
    7473
    … …  
    7675                        md = checkfield(md,'fieldname','inversion.incomplete_adjoint','values',[0 1]);
    7776                        md = checkfield(md,'fieldname','inversion.control_parameters','cell',1,'values',...
    7877                                {'BalancethicknessThickeningRate' 'FrictionCoefficient' 'MaterialsRheologyBbar' 'DamageDbar',...
    79                                 'Vx' 'Vy' 'Thickness' 'BalancethicknessOmega' 'BalancethicknessApparentMassbalance'});
     78                                'Vx' 'Vy' 'Thickness' 'BalancethicknessOmega' 'BalancethicknessApparentMassbalance','MaterialsRheologyB'});
    8079                        md = checkfield(md,'fieldname','inversion.control_scaling_factors','size',[1 num_controls],'>',0,'NaN',1);
    8180                        md = checkfield(md,'fieldname','inversion.maxsteps','numel',1,'>=',0);
    8281                        md = checkfield(md,'fieldname','inversion.maxiter','numel',1,'>=',0);
    8382                        md = checkfield(md,'fieldname','inversion.dxmin','numel',1,'>',0);
    8483                        md = checkfield(md,'fieldname','inversion.gttol','numel',1,'>',0);
    85                         md = checkfield(md,'fieldname','inversion.cost_functions','size',[1 num_costfunc],'values',[101:105 201 501:506 601:604]);
     84                        md = checkfield(md,'fieldname','inversion.cost_functions','size',[1 num_costfunc],'values',[101:105 201 501:507 601:604]);
    8685                        md = checkfield(md,'fieldname','inversion.cost_functions_coefficients','size',[md.mesh.numberofvertices num_costfunc],'>=',0);
    8786                        md = checkfield(md,'fieldname','inversion.min_parameters','size',[md.mesh.numberofvertices num_controls]);
    8887                        md = checkfield(md,'fieldname','inversion.max_parameters','size',[md.mesh.numberofvertices num_controls]);
    … …  
    180179                        pos=find(obj.cost_functions==504); data(pos)=ThicknessAlongGradientEnum();
    181180                        pos=find(obj.cost_functions==505); data(pos)=ThicknessAcrossGradientEnum();
    182181                        pos=find(obj.cost_functions==506); data(pos)=BalancethicknessMisfitEnum();
     182                        pos=find(obj.cost_functions==507); data(pos)=RheologyBAbsGradientEnum();
    183183                        pos=find(obj.cost_functions==601); data(pos)=SurfaceAbsMisfitEnum();
    184184                        WriteData(fid,'data',data,'enum',InversionCostFunctionsEnum(),'format','DoubleMat','mattype',3);
    185185                        WriteData(fid,'data',num_cost_functions,'enum',InversionNumCostFunctionsEnum(),'format','Integer');

  • ../trunk-jpl/src/m/classes/inversionvalidation.m

     
    5858                        md = checkfield(md,'fieldname','inversion.incomplete_adjoint','values',[0 1]);
    5959                        md = checkfield(md,'fieldname','inversion.control_parameters','cell',1,'values',...
    6060                                {'BalancethicknessThickeningRate' 'FrictionCoefficient' 'MaterialsRheologyBbar' 'DamageDbar' 'Vx' 'Vy' 'Thickness',...
    61                                 'BalancethicknessOmega' 'BalancethicknessApparentMassbalance'});
     61                                'BalancethicknessOmega' 'BalancethicknessApparentMassbalance','MaterialsRheologyB'});
    6262                        md = checkfield(md,'fieldname','inversion.control_scaling_factors','size',[1 num_controls],'>',0,'NaN',1);
    63                         md = checkfield(md,'fieldname','inversion.cost_functions','size',[1 num_costfunc],'values',[101:105 201 501:506 601:604]);
     63                        md = checkfield(md,'fieldname','inversion.cost_functions','size',[1 num_costfunc],'values',[101:105 201 501:507 601:604]);
    6464                        md = checkfield(md,'fieldname','inversion.cost_functions_coefficients','size',[md.mesh.numberofvertices num_costfunc],'>=',0);
    6565                        md = checkfield(md,'fieldname','inversion.min_parameters','size',[md.mesh.numberofvertices num_controls]);
    6666                        md = checkfield(md,'fieldname','inversion.max_parameters','size',[md.mesh.numberofvertices num_controls]);
    … …  
    149149                        pos=find(obj.cost_functions==504); data(pos)=ThicknessAlongGradientEnum();
    150150                        pos=find(obj.cost_functions==505); data(pos)=ThicknessAcrossGradientEnum();
    151151                        pos=find(obj.cost_functions==506); data(pos)=BalancethicknessMisfitEnum();
     152                        pos=find(obj.cost_functions==507); data(pos)=RheologyBAbsGradientEnum();
    152153                        pos=find(obj.cost_functions==601); data(pos)=SurfaceAbsMisfitEnum();
    153154                        WriteData(fid,'data',data,'enum',InversionCostFunctionsEnum(),'format','DoubleMat','mattype',3);
    154155                        WriteData(fid,'data',num_cost_functions,'enum',InversionNumCostFunctionsEnum(),'format','Integer');

  • ../trunk-jpl/src/m/classes/inversion.m

     
    133133                        md = checkfield(md,'fieldname','inversion.incomplete_adjoint','values',[0 1]);
    134134                        md = checkfield(md,'fieldname','inversion.control_parameters','cell',1,'values',...
    135135                                {'BalancethicknessThickeningRate' 'FrictionCoefficient' 'MaterialsRheologyBbar' 'DamageDbar' 'Vx' 'Vy' 'Thickness',...
    136                                 'BalancethicknessOmega' 'BalancethicknessApparentMassbalance'});
     136                                'BalancethicknessOmega' 'BalancethicknessApparentMassbalance','MaterialsRheologyB'});
    137137                        md = checkfield(md,'fieldname','inversion.nsteps','numel',1,'>=',0);
    138138                        md = checkfield(md,'fieldname','inversion.maxiter_per_step','size',[md.inversion.nsteps 1],'>=',0);
    139139                        md = checkfield(md,'fieldname','inversion.step_threshold','size',[md.inversion.nsteps 1]);
    140                         md = checkfield(md,'fieldname','inversion.cost_functions','size',[1 num_costfunc],'values',[101:105 201 501:506 601:604]);
     140                        md = checkfield(md,'fieldname','inversion.cost_functions','size',[1 num_costfunc],'values',[101:105 201 501:507 601:604]);
    141141                        md = checkfield(md,'fieldname','inversion.cost_functions_coefficients','size',[md.mesh.numberofvertices num_costfunc],'>=',0);
    142142                        md = checkfield(md,'fieldname','inversion.gradient_scaling','size',[md.inversion.nsteps num_controls]);
    143143                        md = checkfield(md,'fieldname','inversion.min_parameters','size',[md.mesh.numberofvertices num_controls]);
    … …  
    241241                        pos=find(obj.cost_functions==504); data(pos)=ThicknessAlongGradientEnum();
    242242                        pos=find(obj.cost_functions==505); data(pos)=ThicknessAcrossGradientEnum();
    243243                        pos=find(obj.cost_functions==506); data(pos)=BalancethicknessMisfitEnum();
     244                        pos=find(obj.cost_functions==507); data(pos)=RheologyBAbsGradientEnum();
    244245                        pos=find(obj.cost_functions==601); data(pos)=SurfaceAbsMisfitEnum();
    245246                        WriteData(fid,'data',data,'enum',InversionCostFunctionsEnum(),'format','DoubleMat','mattype',3);
    246247                        WriteData(fid,'data',num_cost_functions,'enum',InversionNumCostFunctionsEnum(),'format','Integer');

  • ../trunk-jpl/src/m/classes/clusters/greenplanet.m

     
    88classdef greenplanet
    99    properties (SetAccess=public) 
    1010                 % {{{
    11                  name='greenplanet'
    12                  login='';
    13                  numnodes=20;
    14                  cpuspernode=8;
    15                  port=8000;
    16                  queue='c6145';
    17                  codepath='';
    18                  executionpath='';
    19                  interactive=0;
    20                  time=24*60;
    21                  memory=2;
     11                 name          = 'greenplanet'
     12                 login         = '';
     13                 numnodes      = 20;
     14                 cpuspernode   = 8;
     15                 port          = 8000;
     16                 queue         = 'c6145';
     17                 codepath      = '';
     18                 executionpath = '';
     19                 interactive   = 0;
     20                 time          = 24*60;
     21                 memory        = 2;
    2222         end
    2323         properties (SetAccess=private)
    2424                 np=20*8;
Note: See TracBrowser for help on using the repository browser.