Changeset 485


Ignore:
Timestamp:
05/19/09 11:28:00 (16 years ago)
Author:
seroussi
Message:

fixed problem with loadresults for stokes

File:
1 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk/src/m/solutions/cielo/loadresults.m

    r305 r485  
    2222
    2323else
    24 
    25         if ismacayealpattyn | ishutter,
    2624       
     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
    2731                md.vx=u_g(1:3:end)*md.yts;
    2832                md.vy=u_g(2:3:end)*md.yts;
    2933                md.vz=u_g(3:3:end)*md.yts;
    3034                md.pressure=p_g;
    31         else
    32                 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)
    3635        end
    3736        md.vel=sqrt(md.vx.^2+md.vy.^2+md.vz.^2);
Note: See TracChangeset for help on using the changeset viewer.