Changeset 134

Timestamp:

04/29/09 15:12:25 (16 years ago)

Author:

seroussi

Message:

added velocity in inputs if exists

File:

• issm/trunk/src/m/solutions/cielo/diagnostic.m (modified) (1 diff)

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

@@ -22 +22 @@
         disp(sprintf('\n%s',['computing horizontal velocities...']));
 
-        u_g=diagnostic_core_nonlinear(m_dh,{});
+        %plug existing velocity in inputs
+        if ~isnan(md.vx) & ~isnan(md.vy),
+                indx=m_dh.part(1:2:end);
+                indy=m_dh.part(2:2:end);
+                velocity=zeros(m_dh.nodesets.gsize,1);
+                velocity(indx)=md.vx/md.yts;
+                velocity(indy)=md.vy/md.yts;
+                inputs=struct('velocity',velocity);
+        else
+                inputs={};
+        end
+
+        u_g=diagnostic_core_nonlinear(m_dh,inputs);
 
         if strcmpi(md.type,'3d'),