Changeset 702


Ignore:
Timestamp:
06/02/09 11:32:55 (15 years ago)
Author:
seroussi
Message:

wrong computation of vel

File:
1 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk/src/m/solutions/ice/Loadresults.m

    r674 r702  
    3131        md.results.diagnostic.vx=u_g(indx)*md.yts;
    3232        md.results.diagnostic.vy=u_g(indy)*md.yts;
    33         md.results.diagnostic.vel=sqrt(md.vx.^2+md.vy.^2);
     33        md.results.diagnostic.vel=sqrt(u_g(indx).^2+u_g(indy).^2)*md.yts;
    3434        md.results.diagnostic.pressure=u_g(indp);
    3535
     
    5151        md.results.diagnostic.vy=u_g(indy)*md.yts;
    5252        md.results.diagnostic.vz=u_g(indz)*md.yts;
    53         md.results.diagnostic.vel=sqrt(md.vx.^2+md.vy.^2+md.vz.^2);
     53        md.results.diagnostic.vel=sqrt(u_g(indx).^2+u_g(indy).^2+u_g(indz).^2)*md.yts;
    5454        md.results.diagnostic.pressure=u_g(indp);
    5555end
Note: See TracChangeset for help on using the changeset viewer.