Changeset 21057 for issm/trunk-jpl/examples/Inversion/runme.m
- Timestamp:
- 08/03/16 08:48:38 (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/examples/Inversion/runme.m
r20951 r21057 8 8 md = setflowequation(md,'SSA','all'); 9 9 md.cluster = generic('np',2); 10 md = solve(md, StressbalanceSolutionEnum());10 md = solve(md,'Stressbalance'()); 11 11 plotmodel(md,'axis#all','tight','data',md.materials.rheology_B,'caxis',[ 1.3 1.9]*10^8,'title','"True" B',... 12 12 'data',md.results.StressbalanceSolution.Vel,'title','"observed velocities"') … … 24 24 md.inversion.vel_obs = md.results.StressbalanceSolution.Vel; 25 25 26 md = solve(md, StressbalanceSolutionEnum());26 md = solve(md,'Stressbalance'()); 27 27 plotmodel(md,'axis#all','tight','data',md.materials.rheology_B,'caxis',[ 1.3 1.9]*10^8,'title','B first guess',... 28 28 'data',md.results.StressbalanceSolution.Vel,'title','modeled velocities') … … 45 45 %Go solve! 46 46 md.verbose=verbose(0); 47 md=solve(md, StressbalanceSolutionEnum());47 md=solve(md,'Stressbalance'()); 48 48 plotmodel(md,'axis#all','tight','data',md.results.StressbalanceSolution.MaterialsRheologyBbar,'caxis',[ 1.3 1.9]*10^8,'title','inferred B',... 49 49 'data',md.results.StressbalanceSolution.Vel,'title','modeled velocities') … … 66 66 %Go solve! 67 67 md.verbose=verbose(0); 68 md=solve(md, StressbalanceSolutionEnum());68 md=solve(md,'Stressbalance'()); 69 69 plotmodel(md,'axis#all','tight','data',md.results.StressbalanceSolution.MaterialsRheologyBbar,'caxis',[ 1.3 1.9]*10^8,'title','inferred B',... 70 70 'data',md.results.StressbalanceSolution.Vel,'title','modeled velocities')
Note:
See TracChangeset
for help on using the changeset viewer.