From 71eba9ecec186f02307b83f93fbae363a3563e7f Mon Sep 17 00:00:00 2001 From: Philippe Blain Date: Wed, 15 Jul 2020 17:19:27 -0400 Subject: [PATCH] ice_dyn_vp: remove unused local variables --- cicecore/cicedynB/dynamics/ice_dyn_vp.F90 | 23 +++-------------------- 1 file changed, 3 insertions(+), 20 deletions(-) diff --git a/cicecore/cicedynB/dynamics/ice_dyn_vp.F90 b/cicecore/cicedynB/dynamics/ice_dyn_vp.F90 index 76163b867..88fdac3d7 100644 --- 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 !