Changeset 18269
- Timestamp:
- 07/18/14 14:49:27 (11 years ago)
- Location:
- issm/trunk-jpl/examples
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/examples/Pig/runme.m
r18267 r18269 188 188 189 189 % Load Model 190 md = loadmodel('./Models/PIG.Control_drag'); 191 md.inversion.iscontrol=0; 192 193 disp(' Extruding mesh') 194 number_of_layers=3; 195 md=extrude(md,number_of_layers,0.9); 190 191 % Disable inversion 196 192 197 193 % Extrude Mesh 198 194 195 % Set Flowequation 196 199 197 % Solve 200 md=solve(md,StressbalanceSolutionEnum);201 198 202 199 % Save Model 203 save ./Models/PIG.ModelHO md;204 200 205 201 end … … 219 215 'ylim#all',[min(mdHO.mesh.y) max(mdHO.mesh.y)]/10^3,... 220 216 'FontSize#all',12,... 221 'data',md SSA.initialization.vel,'title','Observed velocity',...217 'data',mdHO.initialization.vel,'title','Observed velocity',... 222 218 'data',(mdHO.results.StressbalanceSolution.Vel(surf)-mdHO.initialization.vel(surf)),'title','(HO-observed) velocities',... 223 219 'data',mdSSA.results.StressbalanceSolution.Vel,'title','Modeled SSA Velocity',... 224 220 'data',(mdHO.results.StressbalanceSolution.Vel(surf)-mdSSA.results.StressbalanceSolution.Vel),'title','(HO-SSA) velocities',... 225 'data', project2d(mdHO,mdHO.results.StressbalanceSolution.Vel,1),'title','Modeled HO surface Velocities',...221 'data',mdHO.results.StressbalanceSolution.Vel,'title','Modeled HO surface Velocities',... 226 222 'data',(mdHO.results.StressbalanceSolution.Vel(surf)-mdHO.results.StressbalanceSolution.Vel(basal)),'title','(HOsurf-HO base) velocities',... 227 223 'caxis#1',([1.5,4000]),'caxis#3',([1.5,4000]),'caxis#5',([1.5,4000]),... 228 224 'colorbar#all','on','view#all',2,... 229 225 'colorbartitle#all','[m/yr]',... 230 'l og#1', 10,'log#3', 10,'log#5', 10,'gridded#all',1);231 end 232 % }}} 226 'layer#5',1, 'log#1', 10,'log#3', 10,'log#5', 10); 227 end 228 % }}} -
issm/trunk-jpl/examples/UncertaintyQuantification/runme.m
r18267 r18269 1 1 %PIG Uncertainty Quantification Application 2 steps=[ 7];2 steps=[1:3]; 3 3 4 4 if any(steps==1) … … 162 162 163 163 %variables 164 md.qmu.variables.Thickness=normal_uncertain('scaled_Thickness',1,0.05);165 164 md.qmu.variables.DragCoefficient=normal_uncertain(... 166 165 'scaled_FrictionCoefficient',1,0.05); 167 166 md.qmu.variables.rheology_B=normal_uncertain(... 168 167 'scaled_MaterialsRheologyB',1,0.05); 168 md.qmu.variables.Thickness=normal_uncertain('scaled_Thickness',1,0.05); 169 169 170 170 %responses
Note:
See TracChangeset
for help on using the changeset viewer.