Changeset 25986 for issm/trunk-jpl/examples/Inversion/runme.m
- Timestamp:
- 02/07/21 00:19:47 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/examples/Inversion/runme.m
r25915 r25986 1 1 steps=[1]; 2 2 3 if any(steps==1) %{{{3 if any(steps==1) 4 4 %Generate observations 5 5 md = model; … … 13 13 'data',md.results.StressbalanceSolution.Vel,'title','"observed velocities"') 14 14 save model1 md 15 end %}}}15 end 16 16 17 if any(steps==2) %{{{17 if any(steps==2) 18 18 %Modify rheology, now constant 19 19 loadmodel('model1.mat'); … … 30 30 'data',md.results.StressbalanceSolution.Vel,'title','modeled velocities') 31 31 save model2 md 32 end %}}}32 end 33 33 34 if any(steps==3) %{{{34 if any(steps==3) 35 35 %invert for ice rigidity 36 36 loadmodel('model2.mat'); … … 51 51 plotmodel(md,'axis#all','tight','data',md.results.StressbalanceSolution.MaterialsRheologyBbar,'caxis',[ 1.3 1.9]*10^8,'title','inferred B',... 52 52 'data',md.results.StressbalanceSolution.Vel,'title','modeled velocities') 53 end %}}}53 end 54 54 55 if any(steps==4) %{{{55 if any(steps==4) 56 56 %invert for ice rigidity 57 57 loadmodel('model2.mat'); … … 73 73 plotmodel(md,'axis#all','tight','data',md.results.StressbalanceSolution.MaterialsRheologyBbar,'caxis',[ 1.3 1.9]*10^8,'title','inferred B',... 74 74 'data',md.results.StressbalanceSolution.Vel,'title','modeled velocities') 75 end %}}}75 end
Note:
See TracChangeset
for help on using the changeset viewer.