Skip to content

+Add wind-stress acceleration diagnostics #1437

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
adcroft merged 2 commits into from
Jul 27, 2021
Merged

+Add wind-stress acceleration diagnostics #1437

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
2 commits
Select commit Hold shift + click to select a range
d7dece9
+Add wind-stress acceleration diagnostics
Hallberg-NOAA Jul 2, 2021
2137fce
Merge branch 'dev/gfdl' into stress_accel_diag
Hallberg-NOAA Jul 26, 2021
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
4 changes: 4 additions & 0 deletions src/core/MOM_variables.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -169,6 +169,10 @@ module MOM_variables
PFv => NULL(), & !< Meridional acceleration due to pressure forces [L T-2 ~> m s-2]
du_dt_visc => NULL(), &!< Zonal acceleration due to vertical viscosity [L T-2 ~> m s-2]
dv_dt_visc => NULL(), &!< Meridional acceleration due to vertical viscosity [L T-2 ~> m s-2]
du_dt_str => NULL(), & !< Zonal acceleration due to the surface stress (included
!! in du_dt_visc) [L T-2 ~> m s-2]
dv_dt_str => NULL(), & !< Meridional acceleration due to the surface stress (included
!! in dv_dt_visc) [L T-2 ~> m s-2]
du_dt_dia => NULL(), & !< Zonal acceleration due to diapycnal mixing [L T-2 ~> m s-2]
dv_dt_dia => NULL(), & !< Meridional acceleration due to diapycnal mixing [L T-2 ~> m s-2]
u_accel_bt => NULL(), &!< Pointer to the zonal barotropic-solver acceleration [L T-2 ~> m s-2]
Expand Down
54 changes: 43 additions & 11 deletions src/diagnostics/MOM_diagnostics.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -107,6 +107,7 @@ module MOM_diagnostics
!! of this spurious Coriolis source.
KE_adv => NULL(), & !< KE source from along-layer advection [H L2 T-3 ~> m3 s-3]
KE_visc => NULL(), & !< KE source from vertical viscosity [H L2 T-3 ~> m3 s-3]
KE_stress => NULL(), & !< KE source from surface stress (included in KE_visc) [H L2 T-3 ~> m3 s-3]
KE_horvisc => NULL(), & !< KE source from horizontal viscosity [H L2 T-3 ~> m3 s-3]
KE_dia => NULL() !< KE source from diapycnal diffusion [H L2 T-3 ~> m3 s-3]

Expand All @@ -121,8 +122,8 @@ module MOM_diagnostics
integer :: id_col_ht = -1, id_dh_dt = -1
integer :: id_KE = -1, id_dKEdt = -1
integer :: id_PE_to_KE = -1, id_KE_BT = -1
integer :: id_KE_Coradv = -1
integer :: id_KE_adv = -1, id_KE_visc = -1
integer :: id_KE_Coradv = -1, id_KE_adv = -1
integer :: id_KE_visc = -1, id_KE_stress = -1
integer :: id_KE_horvisc = -1, id_KE_dia = -1
integer :: id_uh_Rlay = -1, id_vh_Rlay = -1
integer :: id_uhGM_Rlay = -1, id_vhGM_Rlay = -1
Expand Down Expand Up @@ -1060,7 +1061,7 @@ subroutine calculate_energy_diagnostics(u, v, h, uh, vh, ADp, CDp, G, GV, US, CS
KE_h(i,j) = CS%KE(i,j,k) * CS%dh_dt(i,j,k)
enddo ; enddo
if (.not.G%symmetric) &
call do_group_pass(CS%pass_KE_uv, G%domain)
call do_group_pass(CS%pass_KE_uv, G%domain)
do j=js,je ; do i=is,ie
CS%dKE_dt(i,j,k) = KE_h(i,j) + 0.5 * G%IareaT(i,j) &
* (KE_u(I,j) + KE_u(I-1,j) + KE_v(i,J) + KE_v(i,J-1))
Expand All @@ -1078,7 +1079,7 @@ subroutine calculate_energy_diagnostics(u, v, h, uh, vh, ADp, CDp, G, GV, US, CS
KE_v(i,J) = vh(i,J,k) * G%dyCv(i,J) * ADp%PFv(i,J,k)
enddo ; enddo
if (.not.G%symmetric) &
call do_group_pass(CS%pass_KE_uv, G%domain)
call do_group_pass(CS%pass_KE_uv, G%domain)
do j=js,je ; do i=is,ie
CS%PE_to_KE(i,j,k) = 0.5 * G%IareaT(i,j) &
* (KE_u(I,j) + KE_u(I-1,j) + KE_v(i,J) + KE_v(i,J-1))
Expand All @@ -1096,7 +1097,7 @@ subroutine calculate_energy_diagnostics(u, v, h, uh, vh, ADp, CDp, G, GV, US, CS
KE_v(i,J) = vh(i,J,k) * G%dyCv(i,J) * ADp%v_accel_bt(i,J,k)
enddo ; enddo
if (.not.G%symmetric) &
call do_group_pass(CS%pass_KE_uv, G%domain)
call do_group_pass(CS%pass_KE_uv, G%domain)
do j=js,je ; do i=is,ie
CS%KE_BT(i,j,k) = 0.5 * G%IareaT(i,j) &
* (KE_u(I,j) + KE_u(I-1,j) + KE_v(i,J) + KE_v(i,J-1))
Expand All @@ -1118,7 +1119,7 @@ subroutine calculate_energy_diagnostics(u, v, h, uh, vh, ADp, CDp, G, GV, US, CS
* (uh(I,j,k) - uh(I-1,j,k) + vh(i,J,k) - vh(i,J-1,k))
enddo ; enddo
if (.not.G%symmetric) &
call do_group_pass(CS%pass_KE_uv, G%domain)
call do_group_pass(CS%pass_KE_uv, G%domain)
do j=js,je ; do i=is,ie
CS%KE_CorAdv(i,j,k) = KE_h(i,j) + 0.5 * G%IareaT(i,j) &
* (KE_u(I,j) + KE_u(I-1,j) + KE_v(i,J) + KE_v(i,J-1))
Expand Down Expand Up @@ -1146,7 +1147,7 @@ subroutine calculate_energy_diagnostics(u, v, h, uh, vh, ADp, CDp, G, GV, US, CS
* (uh(I,j,k) - uh(I-1,j,k) + vh(i,J,k) - vh(i,J-1,k))
enddo ; enddo
if (.not.G%symmetric) &
call do_group_pass(CS%pass_KE_uv, G%domain)
call do_group_pass(CS%pass_KE_uv, G%domain)
do j=js,je ; do i=is,ie
CS%KE_adv(i,j,k) = KE_h(i,j) + 0.5 * G%IareaT(i,j) &
* (KE_u(I,j) + KE_u(I-1,j) + KE_v(i,J) + KE_v(i,J-1))
Expand All @@ -1164,7 +1165,7 @@ subroutine calculate_energy_diagnostics(u, v, h, uh, vh, ADp, CDp, G, GV, US, CS
KE_v(i,J) = vh(i,J,k) * G%dyCv(i,J) * ADp%dv_dt_visc(i,J,k)
enddo ; enddo
if (.not.G%symmetric) &
call do_group_pass(CS%pass_KE_uv, G%domain)
call do_group_pass(CS%pass_KE_uv, G%domain)
do j=js,je ; do i=is,ie
CS%KE_visc(i,j,k) = 0.5 * G%IareaT(i,j) &
* (KE_u(I,j) + KE_u(I-1,j) + KE_v(i,J) + KE_v(i,J-1))
Expand All @@ -1173,6 +1174,24 @@ subroutine calculate_energy_diagnostics(u, v, h, uh, vh, ADp, CDp, G, GV, US, CS
if (CS%id_KE_visc > 0) call post_data(CS%id_KE_visc, CS%KE_visc, CS%diag)
endif

if (associated(CS%KE_stress)) then
do k=1,nz
do j=js,je ; do I=Isq,Ieq
KE_u(I,j) = uh(I,j,k) * G%dxCu(I,j) * ADp%du_dt_str(I,j,k)
enddo ; enddo
do J=Jsq,Jeq ; do i=is,ie
KE_v(i,J) = vh(i,J,k) * G%dyCv(i,J) * ADp%dv_dt_str(i,J,k)
enddo ; enddo
if (.not.G%symmetric) &
call do_group_pass(CS%pass_KE_uv, G%domain)
do j=js,je ; do i=is,ie
CS%KE_stress(i,j,k) = 0.5 * G%IareaT(i,j) * &
((KE_u(I,j) + KE_u(I-1,j)) + (KE_v(i,J) + KE_v(i,J-1)))
enddo ; enddo
enddo
if (CS%id_KE_stress > 0) call post_data(CS%id_KE_stress, CS%KE_stress, CS%diag)
endif

if (associated(CS%KE_horvisc)) then
do k=1,nz
do j=js,je ; do I=Isq,Ieq
Expand All @@ -1182,7 +1201,7 @@ subroutine calculate_energy_diagnostics(u, v, h, uh, vh, ADp, CDp, G, GV, US, CS
KE_v(i,J) = vh(i,J,k) * G%dyCv(i,J) * ADp%diffv(i,J,k)
enddo ; enddo
if (.not.G%symmetric) &
call do_group_pass(CS%pass_KE_uv, G%domain)
call do_group_pass(CS%pass_KE_uv, G%domain)
do j=js,je ; do i=is,ie
CS%KE_horvisc(i,j,k) = 0.5 * G%IareaT(i,j) &
* (KE_u(I,j) + KE_u(I-1,j) + KE_v(i,J) + KE_v(i,J-1))
Expand All @@ -1203,7 +1222,7 @@ subroutine calculate_energy_diagnostics(u, v, h, uh, vh, ADp, CDp, G, GV, US, CS
KE_h(i,j) = CS%KE(i,j,k) * (CDp%diapyc_vel(i,j,k) - CDp%diapyc_vel(i,j,k+1))
enddo ; enddo
if (.not.G%symmetric) &
call do_group_pass(CS%pass_KE_uv, G%domain)
call do_group_pass(CS%pass_KE_uv, G%domain)
do j=js,je ; do i=is,ie
CS%KE_dia(i,j,k) = KE_h(i,j) + 0.5 * G%IareaT(i,j) &
* (KE_u(I,j) + KE_u(I-1,j) + KE_v(i,J) + KE_v(i,J-1))
Expand Down Expand Up @@ -1894,6 +1913,11 @@ subroutine MOM_diagnostics_init(MIS, ADp, CDp, Time, G, GV, US, param_file, diag
'm3 s-3', conversion=GV%H_to_m*(US%L_T_to_m_s**2)*US%s_to_T)
if (CS%id_KE_visc>0) call safe_alloc_ptr(CS%KE_visc,isd,ied,jsd,jed,nz)

CS%id_KE_stress = register_diag_field('ocean_model', 'KE_stress', diag%axesTL, Time, &
'Kinetic Energy Source from Surface Stresses or Body Wind Stress', &
'm3 s-3', conversion=GV%H_to_m*(US%L_T_to_m_s**2)*US%s_to_T)
if (CS%id_KE_stress>0) call safe_alloc_ptr(CS%KE_stress,isd,ied,jsd,jed,nz)

CS%id_KE_horvisc = register_diag_field('ocean_model', 'KE_horvisc', diag%axesTL, Time, &
'Kinetic Energy Source from Horizontal Viscosity', &
'm3 s-3', conversion=GV%H_to_m*(US%L_T_to_m_s**2)*US%s_to_T)
Expand Down Expand Up @@ -2294,7 +2318,7 @@ subroutine set_dependent_diagnostics(MIS, ADp, CDp, G, GV, CS)

if (associated(CS%dKE_dt) .or. associated(CS%PE_to_KE) .or. &
associated(CS%KE_BT) .or. associated(CS%KE_CorAdv) .or. &
associated(CS%KE_adv) .or. associated(CS%KE_visc) .or. &
associated(CS%KE_adv) .or. associated(CS%KE_visc) .or. associated(CS%KE_stress) .or. &
associated(CS%KE_horvisc) .or. associated(CS%KE_dia)) &
call safe_alloc_ptr(CS%KE,isd,ied,jsd,jed,nz)

Expand Down Expand Up @@ -2323,6 +2347,11 @@ subroutine set_dependent_diagnostics(MIS, ADp, CDp, G, GV, CS)
call safe_alloc_ptr(ADp%dv_dt_visc,isd,ied,JsdB,JedB,nz)
endif

if (associated(CS%KE_stress)) then
call safe_alloc_ptr(ADp%du_dt_str,IsdB,IedB,jsd,jed,nz)
call safe_alloc_ptr(ADp%dv_dt_str,isd,ied,JsdB,JedB,nz)
endif

if (associated(CS%KE_dia)) then
call safe_alloc_ptr(ADp%du_dt_dia,IsdB,IedB,jsd,jed,nz)
call safe_alloc_ptr(ADp%dv_dt_dia,isd,ied,JsdB,JedB,nz)
Expand Down Expand Up @@ -2353,6 +2382,7 @@ subroutine MOM_diagnostics_end(CS, ADp, CDp)
if (associated(CS%KE_Coradv)) deallocate(CS%KE_Coradv)
if (associated(CS%KE_adv)) deallocate(CS%KE_adv)
if (associated(CS%KE_visc)) deallocate(CS%KE_visc)
if (associated(CS%KE_stress)) deallocate(CS%KE_stress)
if (associated(CS%KE_horvisc)) deallocate(CS%KE_horvisc)
if (associated(CS%KE_dia)) deallocate(CS%KE_dia)
if (associated(CS%dv_dt)) deallocate(CS%dv_dt)
Expand All @@ -2368,6 +2398,8 @@ subroutine MOM_diagnostics_end(CS, ADp, CDp)
if (associated(ADp%gradKEu)) deallocate(ADp%gradKEu)
if (associated(ADp%du_dt_visc)) deallocate(ADp%du_dt_visc)
if (associated(ADp%dv_dt_visc)) deallocate(ADp%dv_dt_visc)
if (associated(ADp%du_dt_str)) deallocate(ADp%du_dt_str)
if (associated(ADp%dv_dt_str)) deallocate(ADp%dv_dt_str)
if (associated(ADp%du_dt_dia)) deallocate(ADp%du_dt_dia)
if (associated(ADp%dv_dt_dia)) deallocate(ADp%dv_dt_dia)
if (associated(ADp%du_other)) deallocate(ADp%du_other)
Expand Down
Loading