@@ -27,7 +27,8 @@ subroutine diagnoseKdWork(G, GV, US, N2, Kd, Bdif_flx, dz, Bdif_flx_dz)
27
27
type (verticalGrid_type), intent (in ) :: GV ! < ocean vertical grid structure
28
28
type (unit_scale_type), intent (in ) :: US ! < A dimensional unit scaling type
29
29
real , dimension (SZI_(G),SZJ_(G),SZK_(GV)+ 1 ), &
30
- intent (in ) :: N2, & ! < Buoyancy frequency [T-2 ~> s-2]
30
+ intent (in ) :: N2 ! < Buoyancy frequency [T-2 ~> s-2]
31
+ real , dimension (SZI_(G),SZJ_(G),SZK_(GV)+ 1 ), &
31
32
Kd ! < Diffusivity [H2 T-1 ~> m2 s-1]
32
33
real , dimension (SZI_(G),SZJ_(G),SZK_(GV)+ 1 ), &
33
34
intent (out ) :: Bdif_flx ! < Buoyancy flux [H2 T-3 ~> m2 s-3]
@@ -41,7 +42,7 @@ subroutine diagnoseKdWork(G, GV, US, N2, Kd, Bdif_flx, dz, Bdif_flx_dz)
41
42
! $OMP parallel do default(shared)
42
43
do j= G% jsc,G% jec ; do i= G% isc,G% iec
43
44
Bdif_flx(i,j,1 ) = 0.0
44
- Bdif_flx(i,j,GV% ke+1 ) = 0.0
45
+ Bdif_flx(i,j,GV% ke+1 ) = 0.0
45
46
do K= 2 ,GV% ke
46
47
Bdif_flx(i,j,K) = - N2(i,j,K) * Kd(i,j,K)
47
48
enddo
@@ -57,7 +58,7 @@ subroutine diagnoseKdWork(G, GV, US, N2, Kd, Bdif_flx, dz, Bdif_flx_dz)
57
58
enddo
58
59
enddo ; enddo
59
60
endif
60
-
61
+
61
62
end subroutine diagnoseKdWork
62
63
63
64
! > \namespace mom_diagnose_kdwork
0 commit comments