Changeset 20708


Ignore:
Timestamp:
06/07/16 11:33:03 (9 years ago)
Author:
Mathieu Morlighem
Message:

CHG: revised tutorial

File:
1 edited

Legend:

Unmodified
Added
Removed

  • issm/trunk-jpl/examples/Inversion/runme.m

    r18275 r20708  
    11step=1;
    22if step==1
    3         %Generate observation
     3        %Generate observations
    44        md = model;
    55        md = triangle(md,'DomainOutline.exp',100000);
    … …  
    88        md = setflowequation(md,'SSA','all');
    99        md.cluster = generic('np',2);
    10 
    11         md = solve(md,StressbalanceSolutionEnum);
    12 
    13         plotmodel(md,'data',md.materials.rheology_B,'caxis',[ 1.3 1.9]*10^8,'figure',1,'gridded#all',1);
    14         plotmodel(md,'data',md.results.StressbalanceSolution.Vel,'figure',2,'gridded#all',1); 
    15 
     10        md = solve(md,StressbalanceSolutionEnum());
     11        plotmodel(md,'axis#all','tight','data',md.materials.rheology_B,'caxis',[ 1.3 1.9]*10^8,'title','"True" B',...
     12                'data',md.results.StressbalanceSolution.Vel,'title','"observed velocities"')
    1613        save model1 md
    1714end
    … …  
    2724        md.inversion.vel_obs = md.results.StressbalanceSolution.Vel;
    2825
    29         md = solve(md,StressbalanceSolutionEnum);
    30 
    31         plotmodel(md,'data',md.materials.rheology_B,'caxis',[ 1.3 1.9]*10^8,'figure',1);
    32         plotmodel(md,'data',md.results.StressbalanceSolution.Vel,'figure',2); 
    33 
     26        md = solve(md,StressbalanceSolutionEnum());
     27        plotmodel(md,'axis#all','tight','data',md.materials.rheology_B,'caxis',[ 1.3 1.9]*10^8,'title','B first guess',...
     28                'data',md.results.StressbalanceSolution.Vel,'title','modeled velocities')
    3429        save model2 md
    3530end
    … …  
    4540        md.inversion.cost_functions = 101;
    4641        md.inversion.cost_functions_coefficients = ones(md.mesh.numberofvertices,1);
    47         md.inversion.min_parameters    = paterson(273)*ones(md.mesh.numberofvertices,1);
    48         md.inversion.max_parameters    = paterson(200)*ones(md.mesh.numberofvertices,1);
     42        md.inversion.min_parameters    = cuffey(273)*ones(md.mesh.numberofvertices,1);
     43        md.inversion.max_parameters    = cuffey(200)*ones(md.mesh.numberofvertices,1);
    4944
    5045        %Go solve!
    5146        md.verbose=verbose(0);
    52         md=solve(md,StressbalanceSolutionEnum);
    53 
    54         plotmodel(md,'data',md.results.StressbalanceSolution.MaterialsRheologyBbar,'caxis',[ 1.3 1.9]*10^8,'figure',1);
    55         plotmodel(md,'data',md.results.StressbalanceSolution.Vel,'figure',2); 
     47        md=solve(md,StressbalanceSolutionEnum());
     48        plotmodel(md,'axis#all','tight','data',md.results.StressbalanceSolution.MaterialsRheologyBbar,'caxis',[ 1.3 1.9]*10^8,'title','inferred B',...
     49                'data',md.results.StressbalanceSolution.Vel,'title','modeled velocities')
    5650end
    5751if step==4
    … …  
    6761        md.inversion.cost_functions_coefficients      = ones(md.mesh.numberofvertices,1);
    6862        md.inversion.cost_functions_coefficients(:,2) = 10^-16*ones(md.mesh.numberofvertices,1);
    69         md.inversion.min_parameters    = paterson(273)*ones(md.mesh.numberofvertices,1);
    70         md.inversion.max_parameters    = paterson(200)*ones(md.mesh.numberofvertices,1);
     63        md.inversion.min_parameters    = cuffey(273)*ones(md.mesh.numberofvertices,1);
     64        md.inversion.max_parameters    = cuffey(200)*ones(md.mesh.numberofvertices,1);
    7165
    7266        %Go solve!
    7367        md.verbose=verbose(0);
    74         md=solve(md,StressbalanceSolutionEnum);
    75 
    76         plotmodel(md,'data',md.results.StressbalanceSolution.MaterialsRheologyBbar,'caxis',[ 1.3 1.9]*10^8,'figure',1);
    77         plotmodel(md,'data',md.results.StressbalanceSolution.Vel,'figure',2); 
     68        md=solve(md,StressbalanceSolutionEnum());
     69        plotmodel(md,'axis#all','tight','data',md.results.StressbalanceSolution.MaterialsRheologyBbar,'caxis',[ 1.3 1.9]*10^8,'title','inferred B',...
     70                'data',md.results.StressbalanceSolution.Vel,'title','modeled velocities')
    7871end
Note: See TracChangeset for help on using the changeset viewer.