Ignore:
Timestamp:
02/07/21 00:19:47 (4 years ago)
Author:
jdquinn
Message:

CHG: Cleanup

File:
1 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk-jpl/examples/Inversion/runme.m

    r25915 r25986  
    11steps=[1];
    22
    3 if any(steps==1) %{{{
     3if any(steps==1)
    44        %Generate observations
    55        md = model;
     
    1313                'data',md.results.StressbalanceSolution.Vel,'title','"observed velocities"')
    1414        save model1 md
    15 end %}}}
     15end
    1616
    17 if any(steps==2) %{{{
     17if any(steps==2)
    1818        %Modify rheology, now constant
    1919        loadmodel('model1.mat');
     
    3030                'data',md.results.StressbalanceSolution.Vel,'title','modeled velocities')
    3131        save model2 md
    32 end %}}}
     32end
    3333
    34 if any(steps==3) %{{{
     34if any(steps==3)
    3535        %invert for ice rigidity
    3636        loadmodel('model2.mat');
     
    5151        plotmodel(md,'axis#all','tight','data',md.results.StressbalanceSolution.MaterialsRheologyBbar,'caxis',[ 1.3 1.9]*10^8,'title','inferred B',...
    5252                'data',md.results.StressbalanceSolution.Vel,'title','modeled velocities')
    53 end %}}}
     53end
    5454
    55 if any(steps==4) %{{{
     55if any(steps==4)
    5656        %invert for ice rigidity
    5757        loadmodel('model2.mat');
     
    7373        plotmodel(md,'axis#all','tight','data',md.results.StressbalanceSolution.MaterialsRheologyBbar,'caxis',[ 1.3 1.9]*10^8,'title','inferred B',...
    7474                'data',md.results.StressbalanceSolution.Vel,'title','modeled velocities')
    75 end %}}}
     75end
Note: See TracChangeset for help on using the changeset viewer.