fix couple issues with diss_est

NOAA-GFDL · Oct 6, 2023 · 9f507e2 · 9f507e2
1 parent 324e3da
commit 9f507e2
Showing 1 changed file with 4 additions and 1 deletion.
diff --git 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