Changeset 10981 for issm/trunk/test/NightlyRun/test1501.m
- Timestamp:
- 11/29/11 19:29:01 (13 years ago)
- Location:
- issm/trunk
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk
- Property svn:mergeinfo changed
/issm/trunk-jpl merged: 10962,10965,10969-10976,10980
- Property svn:mergeinfo changed
-
issm/trunk/test/NightlyRun/test1501.m
r9766 r10981 21 21 for i=1:10 22 22 md=solve(md,PrognosticSolutionEnum); 23 md.surfaceforcings.mass_balance= md.surfaceforcings.mass_balance - ( PatchToVec(md.results.PrognosticSolution.Thickness)-md.geometry.thickness);23 md.surfaceforcings.mass_balance= md.surfaceforcings.mass_balance - ((md.results.PrognosticSolution.Thickness)-md.geometry.thickness); 24 24 end 25 25 … … 47 47 1e-09,1e-09,1e-09,1e-10,1e-10,1e-10,1e-10,1e-10}; 48 48 field_values={... 49 PatchToVec(md.results.TransientSolution(400).Vx),...50 PatchToVec(md.results.TransientSolution(400).Vy),...51 PatchToVec(md.results.TransientSolution(400).Vel),...52 PatchToVec(md.results.TransientSolution(400).Pressure),...53 PatchToVec(md.results.TransientSolution(400).Bed),...54 PatchToVec(md.results.TransientSolution(400).Surface),...55 PatchToVec(md.results.TransientSolution(400).Thickness),...56 PatchToVec(md.results.TransientSolution(400).SurfaceforcingsMassBalance),...57 PatchToVec(md.results.TransientSolution(800).Vx),...58 PatchToVec(md.results.TransientSolution(800).Vy),...59 PatchToVec(md.results.TransientSolution(800).Vel),...60 PatchToVec(md.results.TransientSolution(800).Pressure),...61 PatchToVec(md.results.TransientSolution(800).Bed),...62 PatchToVec(md.results.TransientSolution(800).Surface),...63 PatchToVec(md.results.TransientSolution(800).Thickness),...64 PatchToVec(md.results.TransientSolution(800).SurfaceforcingsMassBalance),...65 PatchToVec(md.results.TransientSolution(1200).Vx),...66 PatchToVec(md.results.TransientSolution(1200).Vy),...67 PatchToVec(md.results.TransientSolution(1200).Vel),...68 PatchToVec(md.results.TransientSolution(1200).Pressure),...69 PatchToVec(md.results.TransientSolution(1200).Bed),...70 PatchToVec(md.results.TransientSolution(1200).Surface),...71 PatchToVec(md.results.TransientSolution(1200).Thickness),...72 PatchToVec(md.results.TransientSolution(1200).SurfaceforcingsMassBalance),...73 PatchToVec(md.results.TransientSolution(1600).Vx),...74 PatchToVec(md.results.TransientSolution(1600).Vy),...75 PatchToVec(md.results.TransientSolution(1600).Vel),...76 PatchToVec(md.results.TransientSolution(1600).Pressure),...77 PatchToVec(md.results.TransientSolution(1600).Bed),...78 PatchToVec(md.results.TransientSolution(1600).Surface),...79 PatchToVec(md.results.TransientSolution(1600).Thickness),...80 PatchToVec(md.results.TransientSolution(1600).SurfaceforcingsMassBalance),...81 PatchToVec(md.results.TransientSolution(2000).Vx),...82 PatchToVec(md.results.TransientSolution(2000).Vy),...83 PatchToVec(md.results.TransientSolution(2000).Vel),...84 PatchToVec(md.results.TransientSolution(2000).Pressure),...85 PatchToVec(md.results.TransientSolution(2000).Bed),...86 PatchToVec(md.results.TransientSolution(2000).Surface),...87 PatchToVec(md.results.TransientSolution(2000).Thickness),...88 PatchToVec(md.results.TransientSolution(2000).SurfaceforcingsMassBalance),...49 (md.results.TransientSolution(400).Vx),... 50 (md.results.TransientSolution(400).Vy),... 51 (md.results.TransientSolution(400).Vel),... 52 (md.results.TransientSolution(400).Pressure),... 53 (md.results.TransientSolution(400).Bed),... 54 (md.results.TransientSolution(400).Surface),... 55 (md.results.TransientSolution(400).Thickness),... 56 (md.results.TransientSolution(400).SurfaceforcingsMassBalance),... 57 (md.results.TransientSolution(800).Vx),... 58 (md.results.TransientSolution(800).Vy),... 59 (md.results.TransientSolution(800).Vel),... 60 (md.results.TransientSolution(800).Pressure),... 61 (md.results.TransientSolution(800).Bed),... 62 (md.results.TransientSolution(800).Surface),... 63 (md.results.TransientSolution(800).Thickness),... 64 (md.results.TransientSolution(800).SurfaceforcingsMassBalance),... 65 (md.results.TransientSolution(1200).Vx),... 66 (md.results.TransientSolution(1200).Vy),... 67 (md.results.TransientSolution(1200).Vel),... 68 (md.results.TransientSolution(1200).Pressure),... 69 (md.results.TransientSolution(1200).Bed),... 70 (md.results.TransientSolution(1200).Surface),... 71 (md.results.TransientSolution(1200).Thickness),... 72 (md.results.TransientSolution(1200).SurfaceforcingsMassBalance),... 73 (md.results.TransientSolution(1600).Vx),... 74 (md.results.TransientSolution(1600).Vy),... 75 (md.results.TransientSolution(1600).Vel),... 76 (md.results.TransientSolution(1600).Pressure),... 77 (md.results.TransientSolution(1600).Bed),... 78 (md.results.TransientSolution(1600).Surface),... 79 (md.results.TransientSolution(1600).Thickness),... 80 (md.results.TransientSolution(1600).SurfaceforcingsMassBalance),... 81 (md.results.TransientSolution(2000).Vx),... 82 (md.results.TransientSolution(2000).Vy),... 83 (md.results.TransientSolution(2000).Vel),... 84 (md.results.TransientSolution(2000).Pressure),... 85 (md.results.TransientSolution(2000).Bed),... 86 (md.results.TransientSolution(2000).Surface),... 87 (md.results.TransientSolution(2000).Thickness),... 88 (md.results.TransientSolution(2000).SurfaceforcingsMassBalance),... 89 89 }; 90 90 … … 106 106 velocity = []; 107 107 for t=starttime:endtime 108 thickness = [thickness PatchToVec(md.results.TransientSolution(t).Thickness)];108 thickness = [thickness (md.results.TransientSolution(t).Thickness)]; 109 109 volume = [volume mean(md.results.TransientSolution(t).Thickness.value,2).*areas]; 110 massbal = [massbal PatchToVec(md.results.TransientSolution(t).SurfaceforcingsMassBalance)];111 velocity = [velocity PatchToVec(md.results.TransientSolution(t).Vel)];110 massbal = [massbal (md.results.TransientSolution(t).SurfaceforcingsMassBalance)]; 111 velocity = [velocity (md.results.TransientSolution(t).Vel)]; 112 112 end 113 113
Note:
See TracChangeset
for help on using the changeset viewer.