ice_dyn_vp: remove unused local variables

phil-blain · Jul 15, 2020 · 71eba9e · 71eba9e
1 parent 083b611
commit 71eba9e
Showing 1 changed file with 3 additions and 20 deletions.
diff --git a/cicecore/cicedynB/dynamics/ice_dyn_vp.F90 b/cicecore/cicedynB/dynamics/ice_dyn_vp.F90
@@ -754,7 +754,6 @@ subroutine anderson_solver (icellt,   icellu,  &
          res_diff   , & ! difference between current and previous residuals
          fpfunc     , & ! current value of fixed point function
          fpfunc_old , & ! previous value of fixed point function
-         Fvec       , & ! (Fx,Fy) (nonlinear residual) as vector
          tmp            ! temporary vector for BLAS calls
 
       real (kind=dbl_kind), dimension(ntot,im_andacc) :: &
@@ -769,7 +768,7 @@ subroutine anderson_solver (icellt,   icellu,  &
          coeffs       ! coeffs used to combine previous solutions
 
       real (kind=dbl_kind) :: & 
-         tol         , & ! tolerance for fixed point convergence: reltol_andacc * (initial fixed point residual norm)
+         ! tol         , & ! tolerance for fixed point convergence: reltol_andacc * (initial fixed point residual norm)  [unused for now]
          tol_nl      , & ! tolerance for nonlinear convergence: reltol_nonlin * (initial nonlinear residual norm)
          fpres_norm  , & ! norm of current fixed point residual : f(x) = g(x) - x
          prog_norm   , & ! norm of difference between current and previous solution
@@ -2771,9 +2770,7 @@ subroutine fgmres (zetaD,               &
 
       real (kind=dbl_kind), dimension (nx_block,ny_block,max_blocks) :: &
          workspace_x , & ! work vector (x components)
-         workspace_y , & ! work vector (y components)
-         Fx          , & ! residual vector (x components), Fx = bx - Au (N/m^2)
-         Fy              ! residual vector (y components), Fy = by - Av (N/m^2)
+         workspace_y     ! work vector (y components)
 
       real (kind=dbl_kind), dimension (max_blocks) :: &
          norm_squared   ! array to accumulate squared norm of grid function over blocks
@@ -2810,12 +2807,6 @@ subroutine fgmres (zetaD,               &
 
       real (kind=dbl_kind) :: relative_tolerance, r0 !phb DESCRIBE if we keep
 
-      real (kind=dbl_kind) :: &
-         local_dot         ! local value to accumulate dot product computations
-
-      real (kind=dbl_kind), dimension(maxinner) :: &
-         dotprod_local     ! local array to accumulate several dot product computations
-
       character(len=*), parameter :: subname = '(fgmres)'
 
       ! Here we go !
@@ -3172,9 +3163,7 @@ subroutine pgmres (zetaD,                &
 
       real (kind=dbl_kind), dimension (nx_block,ny_block,max_blocks) :: &
          workspace_x , & ! work vector (x components)
-         workspace_y , & ! work vector (y components)
-         Fx          , & ! residual vector (x components), Fx = bx - Au (N/m^2)
-         Fy              ! residual vector (y components), Fy = by - Av (N/m^2)
+         workspace_y     ! work vector (y components)
 
       real (kind=dbl_kind), dimension (max_blocks) :: &
          norm_squared   ! array to accumulate squared norm of grid function over blocks
@@ -3208,12 +3197,6 @@ subroutine pgmres (zetaD,                &
 
       real (kind=dbl_kind) :: relative_tolerance, r0 !phb DESCRIBE if we keep
 
-      real (kind=dbl_kind) :: &
-         local_dot         ! local value to accumulate dot product computations
-
-      real (kind=dbl_kind), dimension(maxinner) :: &
-         dotprod_local     ! local array to accumulate several dot product computations
-
       character(len=*), parameter :: subname = '(pgmres)'
 
       ! Here we go !