Changeset 485
- Timestamp:
- 05/19/09 11:28:00 (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk/src/m/solutions/cielo/loadresults.m
r305 r485 22 22 23 23 else 24 25 if ismacayealpattyn | ishutter,26 24 25 if isstokes, 26 md.vx=u_g(1:4:end)*md.yts; 27 md.vy=u_g(2:4:end)*md.yts; 28 md.vz=u_g(3:4:end)*md.yts; 29 md.pressure=p_g; 30 else 27 31 md.vx=u_g(1:3:end)*md.yts; 28 32 md.vy=u_g(2:3:end)*md.yts; 29 33 md.vz=u_g(3:3:end)*md.yts; 30 34 md.pressure=p_g; 31 else32 md.vx=u_g(1:4:end)*md.yts;33 md.vy=u_g(2:4:end)*md.yts;34 md.vz=u_g(3:4:end)*md.yts;35 md.pressure=u_g(4:4:end)36 35 end 37 36 md.vel=sqrt(md.vx.^2+md.vy.^2+md.vz.^2);
Note:
See TracChangeset
for help on using the changeset viewer.