Skip to content

Commit

Permalink
ice_dyn_vp: matvec: remove unneeded local variable
Browse files Browse the repository at this point in the history
  • Loading branch information
phil-blain committed Jul 16, 2020
1 parent 5ff03ab commit eadcfa6
Showing 1 changed file with 2 additions and 6 deletions.
8 changes: 2 additions & 6 deletions cicecore/cicedynB/dynamics/ice_dyn_vp.F90
Original file line number Diff line number Diff line change
Expand Up @@ -1565,7 +1565,6 @@ subroutine matvec (nx_block, ny_block, &
str

real (kind=dbl_kind) :: &
utp, vtp , & ! utp = uvel, vtp = vvel
ccaimp,ccb , & ! intermediate variables
strintx, strinty ! divergence of the internal stress tensor

Expand Down Expand Up @@ -1754,9 +1753,6 @@ subroutine matvec (nx_block, ny_block, &
i = indxui(ij)
j = indxuj(ij)

utp = uvel(i,j)
vtp = vvel(i,j)

ccaimp = umassdti(i,j) + vrel(i,j) * cosw + Cb(i,j) ! kg/m^2 s

ccb = fm(i,j) + sign(c1,fm(i,j)) * vrel(i,j) * sinw ! kg/m^2 s
Expand All @@ -1767,8 +1763,8 @@ subroutine matvec (nx_block, ny_block, &
strinty = uarear(i,j)* &
(str(i,j,5) + str(i,j+1,6) + str(i+1,j,7) + str(i+1,j+1,8))

Au(i,j) = ccaimp*utp - ccb*vtp - strintx
Av(i,j) = ccaimp*vtp + ccb*utp - strinty
Au(i,j) = ccaimp*uvel(i,j) - ccb*vvel(i,j) - strintx
Av(i,j) = ccaimp*vvel(i,j) + ccb*uvel(i,j) - strinty
enddo ! ij - icellu

end subroutine matvec
Expand Down

0 comments on commit eadcfa6

Please sign in to comment.