Changeset 113

Timestamp:

04/28/09 16:39:07 (16 years ago)

Author:

seroussi

Message:

Better initialization of velocity (to be fixed in diagfnostic.m tomorrow)

File:

• issm/trunk/src/m/solutions/cielo/diagnostic_core_nonlinear.m (modified) (3 diffs)

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

@@ -2 +2 @@
 %INPUT function [ru_g varargout]=cielodiagnostic_core_nonlinear(m,inputs)
         
-        %check inputs
-        if isempty(inputs),
-                clear inputs;
-        end
-
 %   first off! We are going to modify the loads dataset. We need to shield the loads from those changes once we return;
         loads=m.loads;
@@ -16 +11 @@
         count=1;
         converged=0;
-        soln(count).u_g=[];
+        [velocity_param velocity_is_present]=recover_input(inputs,'velocity');
+        if velocity_is_present
+                soln(count).u_g=velocity_param;
+        else    
+                soln(count).u_g=[];
+        end
         soln(count).u_f=[];
 
@@ -30 +30 @@
                 if count>1,
                         inputs.old_velocity=full(soln(count-1).u_g);
+                else
+                        inputs.old_velocity=full(soln(count).u_g);
                 end