From 9af882c63615925e75a6b812ab4ab227a113f454 Mon Sep 17 00:00:00 2001 From: JosephMouallem <75327266+JosephMouallem@users.noreply.github.com> Date: Fri, 13 Oct 2023 07:51:14 -0400 Subject: [PATCH] fix couple issues with diss_est (#295) --- model/sw_core.F90 | 5 ++++- 1 file changed, 4 insertions(+), 1 deletion(-) diff --git a/model/sw_core.F90 b/model/sw_core.F90 index b1c21facd..058463dff 100644 --- a/model/sw_core.F90 +++ b/model/sw_core.F90 @@ -1439,7 +1439,7 @@ subroutine d_sw(delpc, delp, ptc, pt, u, v, w, uc,vc, & endif - if ( d_con > 1.e-5 ) then + if ( d_con > 1.e-5 .or. flagstruct%do_diss_est) then do j=js,je+1 do i=is,ie ub(i,j) = vort(i,j) - vort(i+1,j) @@ -1493,6 +1493,9 @@ subroutine d_sw(delpc, delp, ptc, pt, u, v, w, uc,vc, & if ( damp_v>1.E-5 ) then damp4 = (damp_v*gridstruct%da_min_c)**(nord_v+1) call del6_vt_flux(nord_v, npx, npy, damp4, wk, vort, ut, vt, gridstruct, bd) + elseif (flagstruct%do_diss_est) then + ut=0. + vt=0. endif if ( d_con > 1.e-5 .or. flagstruct%do_diss_est ) then