code changed to pass RTs

physics/satmedmfvdifq.F

@@ -56,7 +56,7 @@ end subroutine satmedmfvdifq_finalize
 !! -# A mass-flux approach is also used to represent the stratocumulus-top-induced turbulence
 !! (mfscuq.f).
 !! \section detail_satmedmfvidfq GFS satmedmfvdifq Detailed Algorithm 
-!! @{
+!> @{
       subroutine satmedmfvdifq_run(ix,im,km,ntrac,ntcw,ntiw,ntke,       &
      &     grav,rd,cp,rv,hvap,hfus,fv,eps,epsm1,                        &
      &     dv,du,tdt,rtg,u1,v1,t1,q1,swh,hlw,xmu,garea,                 &
@@ -196,7 +196,7 @@ subroutine satmedmfvdifq_run(ix,im,km,ntrac,ntcw,ntiw,ntke,       &
      &                     rlmn,    rlmn1,  rlmx,   elmx,
      &                     ttend,   utend,  vtend,  qtend,
      &                     zfac,    zfmin,  vk,     spdk2,
-     &                     tkmin,   tkminx, xkzinv, xkgdx,
+     &                     tkmin,   xkzinv, xkgdx,
      &                     zlup,    zldn,   bsum,
      &                     tem,     tem1,   tem2,
      &                     ptem,    ptem0,  ptem1,  ptem2
@@ -215,11 +215,11 @@ subroutine satmedmfvdifq_run(ix,im,km,ntrac,ntcw,ntiw,ntke,       &
       parameter(prmin=0.25,prmax=4.0)
       parameter(pr0=1.0,prtke=1.0,prscu=0.67)
       parameter(f0=1.e-4,crbmin=0.15,crbmax=0.35)
-      parameter(tkmin=1.e-9,tkminx=0.2,dspmax=10.0)
+      parameter(tkmin=1.e-9,dspmax=10.0)
       parameter(qmin=1.e-8,qlmin=1.e-12,zfmin=1.e-8)
       parameter(aphi5=5.,aphi16=16.)
       parameter(elmfac=1.0,elefac=1.0,cql=100.)
-      parameter(dw2min=1.e-4,dkmax=1000.,xkgdx=5000.)
+      parameter(dw2min=1.e-4,dkmax=1000.,xkgdx=25000.)
       parameter(qlcr=3.5e-5,zstblmax=2500.,xkzinv=0.1)
       parameter(h1=0.33333333)
       parameter(ck0=0.4,ck1=0.15,ch0=0.4,ch1=0.15)
@@ -335,20 +335,20 @@ subroutine satmedmfvdifq_run(ix,im,km,ntrac,ntcw,ntiw,ntke,       &
           xkzo(i,k)  = 0.0
           xkzmo(i,k) = 0.0
           if (k < kinver(i)) then
-!                               minimum turbulent mixing length
+!                                  vertical background diffusivity
+            ptem      = prsi(i,k+1) * tx1(i)
+            tem1      = 1.0 - ptem
+            tem2      = tem1 * tem1 * 10.0
+            tem2      = min(1.0, exp(-tem2))
+            xkzo(i,k) = xkzm_hx(i) * tem2
+!
             ptem      = prsl(i,k) * tx1(i)
             tem1      = 1.0 - ptem
             tem2      = tem1 * tem1 * 2.5
             tem2      = min(1.0, exp(-tem2))
             rlmnz(i,k)= rlmn * tem2
             rlmnz(i,k)= max(rlmnz(i,k), rlmn1)
-!                               vertical background diffusivity
-            ptem      = prsi(i,k+1) * tx1(i)
-            tem1      = 1.0 - ptem
-            tem2      = tem1 * tem1 * 10.0
-            tem2      = min(1.0, exp(-tem2))
-            xkzo(i,k) = xkzm_hx(i) * tem2
-!                               vertical background diffusivity for momentum
+!                                 vertical background diffusivity for momentum
             if (ptem >= xkzm_s) then
               xkzmo(i,k) = xkzm_mx(i)
               kx1(i)     = k + 1
@@ -542,7 +542,7 @@ subroutine satmedmfvdifq_run(ix,im,km,ntrac,ntcw,ntiw,ntke,       &
       enddo
       enddo
 !
-! Find pbl height based on bulk richardson number (mrf pbl scheme)
+! find pbl height based on bulk richardson number (mrf pbl scheme)
 !   and also for diagnostic purpose
 !
       do i=1,im
@@ -808,23 +808,23 @@ subroutine satmedmfvdifq_run(ix,im,km,ntrac,ntcw,ntiw,ntke,       &
 !
 !  background diffusivity decreasing with increasing surface layer stability
 !
-!     do i = 1, im
-!       if(.not.sfcflg(i)) then
-!         tem = (1. + 5. * rbsoil(i))**2.
-!!        tem = (1. + 5. * zol(i))**2.
-!         frik(i) = 0.1 + 0.9 / tem
-!       endif
-!     enddo
-!
-!     do k = 1,km1
-!       do i=1,im
-!         xkzo(i,k) = frik(i) * xkzo(i,k)
-!         xkzmo(i,k)= frik(i) * xkzmo(i,k)
-!       enddo
-!     enddo
+      do i = 1, im
+        if(.not.sfcflg(i)) then
+          tem = (1. + 5. * rbsoil(i))**2.
+!         tem = (1. + 5. * zol(i))**2.
+          frik(i) = 0.1 + 0.9 / tem
+        endif
+      enddo
+!
+      do k = 1,km1
+        do i=1,im
+          xkzo(i,k) = frik(i) * xkzo(i,k)
+          xkzmo(i,k)= frik(i) * xkzmo(i,k)
+        enddo
+      enddo
 !
 !>  ## The background vertical diffusivities in the inversion layers are limited 
-!!     to be less than or equal to xkzinv
+!    to be less than or equal to xkzminv
 !
       do k = 1,km1
         do i=1,im
@@ -1013,14 +1013,13 @@ subroutine satmedmfvdifq_run(ix,im,km,ntrac,ntcw,ntiw,ntke,       &
         do i = 1, im
           if(k == 1) then
             tem = ckz(i,1)
-            tem1 = 0.5 * xkzmo(i,1)
+            tem1 = xkzmo(i,1)
           else
             tem = 0.5 * (ckz(i,k-1) + ckz(i,k))
             tem1 = 0.5 * (xkzmo(i,k-1) + xkzmo(i,k))
           endif
           ptem = tem1 / (tem * elm(i,k))
           tkmnz(i,k) = ptem * ptem
-          tkmnz(i,k) = min(tkmnz(i,k), tkminx)
           tkmnz(i,k) = max(tkmnz(i,k), tkmin)
         enddo
       enddo