fixup! ice_dyn_vp: initialize residual vector to zero

phil-blain · Aug 12, 2022 · be571c5 · be571c5
1 parent a1f6364
commit be571c5
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/cicecore/cicedynB/dynamics/ice_dyn_vp.F90 b/cicecore/cicedynB/dynamics/ice_dyn_vp.F90
@@ -816,6 +816,8 @@ subroutine anderson_solver (icellt  , icellu , &
 
       ! Initialization
       res_num = 0
+      Fx = c0
+      Fy = c0
 
       !$OMP PARALLEL DO PRIVATE(iblk)
       do iblk = 1, nblocks
@@ -874,8 +876,6 @@ subroutine anderson_solver (icellt  , icellu , &
                             vrel     (:,:,iblk))
 
             ! Compute nonlinear residual norm (PDE residual)
-            Fx = c0
-            Fy = c0
             call matvec (nx_block             , ny_block           , &
                          icellu       (iblk)  , icellt       (iblk), &
                          indxui     (:,iblk)  , indxuj     (:,iblk), &