fix declaration bug

FESOM · Jul 12, 2023 · ae5bf22 · ae5bf22
1 parent cdb4c4e
commit ae5bf22
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 3 deletions.
diff --git a/src/oce_ale_ssh_splitexpl_subcycl.F90 b/src/oce_ale_ssh_splitexpl_subcycl.F90
@@ -1351,8 +1351,8 @@ subroutine compute_vert_vel_transpv(dynamics, partit, mesh)
     type(t_dyn)   , intent(inout), target :: dynamics
     type(t_partit), intent(inout), target :: partit
     type(t_mesh)  , intent(in)   , target :: mesh
-    integer                               :: node, elem, nz, nzmin, nzmax, ednodes(2), edelem(2) 
-    real(kind=WP)                         :: hh_inv 
+    integer                               :: node, elem, nz, ed, nzmin, nzmax, ednodes(2), edelem(2) 
+    real(kind=WP)                         :: hh_inv, deltaX1, deltaX2, deltaY1, deltaY2 
     real(kind=WP)                         :: e1c1(mesh%nl-1), e1c2(mesh%nl-1)
     real(kind=WP)                         :: e2c1(mesh%nl-1), e2c2(mesh%nl-1)
 

diff --git a/src/oce_ale_vel_rhs.F90 b/src/oce_ale_vel_rhs.F90
@@ -173,7 +173,7 @@ subroutine compute_vel_rhs(ice, dynamics, partit, mesh)
             end do
         else
             UVBT_4AB(1:2, elem) = 0.0_WP
-            do nz=nzmin,nzmax-1
+            do nz=nzmin, nzmax-1
                 ! add pressure gradient terms
                 UV_rhs(  1, nz, elem) = UV_rhs(1, nz, elem) + (Fx-pgf_x(nz, elem))*elem_area(elem)*helem(nz,elem)
                 UV_rhs(  2, nz, elem) = UV_rhs(2, nz, elem) + (Fy-pgf_y(nz, elem))*elem_area(elem)*helem(nz,elem)