Changeset 129

Timestamp:

04/29/09 14:50:09 (16 years ago)

Author:

Eric.Larour

Message:

Fixed problem for vx

File:

• issm/trunk/src/m/solutions/cielo/diagnostic.m (modified) (2 diffs)

issm/trunk/src/m/solutions/cielo/diagnostic.m

@@ -1 +1 @@
 function md=diagnostic(md)
+%DIAGNOSTIC - diagnostic solution sequence
+%
+%   Usage:
+%      md=diagnostic(md)
 
         %timing
@@ -38 +42 @@
                 md.pressure=md.rho_ice*md.g*(md.surface-md.z);
         else
-                
                 %load results onto model
                 md.vx=u_g(1:2:end)*md.yts;
                 md.vy=u_g(2:2:end)*md.yts;
-                md.vz=zeros(m_dv.nodesets.gsize,1);
+                md.vz=zeros(md.numberofgrids,1);
                 md.vel=sqrt(md.vx.^2+md.vy.^2+md.vz.^2);