Add new PBL diagnostic variables dkt and dku

physics/moninedmf.f

@@ -64,7 +64,7 @@ subroutine hedmf_run (im,km,ntrac,ntcw,dv,du,tau,rtg,             &
      &   prsi,del,prsl,prslk,phii,phil,delt,dspheat,                    &
      &   dusfc,dvsfc,dtsfc,dqsfc,hpbl,hgamt,hgamq,dkt,                  &
      &   kinver,xkzm_m,xkzm_h,xkzm_s,lprnt,ipr,                         &
-     &   xkzminv,moninq_fac,hurr_pbl,islimsk,var_ric,                   &
+     &   xkzminv,moninq_fac,hurr_pbl,islimsk,dkudiagnostic,var_ric,     &
      &   coef_ric_l,coef_ric_s,lssav,ldiag3d,qdiag3d,ntoz,              &
      &   du3dt_PBL,dv3dt_PBL,dt3dt_PBL,dq3dt_PBL,do3dt_PBL,             &
      &   flag_for_pbl_generic_tend,errmsg,errflg)
@@ -154,7 +154,7 @@ subroutine hedmf_run (im,km,ntrac,ntcw,dv,du,tau,rtg,             &
      &                     ti(im,km-1),  shr2(im,km-1),                 &
      &                     al(im,km-1),  ad(im,km),                     &
      &                     au(im,km-1),  a1(im,km),                     &
-     &                     a2(im,km*ntrac)
+     &                     a2(im,km*ntrac), dkudiagnostic(im,km-1)
 !
       real(kind=kind_phys) tcko(im,km),  qcko(im,km,ntrac),             &
      &                     ucko(im,km),  vcko(im,km),  xmf(im,km)
@@ -1402,6 +1402,13 @@ subroutine hedmf_run (im,km,ntrac,ntcw,dv,du,tau,rtg,             &
 !
         enddo
       enddo
+
+      do k = 1,km1
+        do i=1,im
+           dkudiagnostic(i,k) = dku(i,k)
+        enddo
+      enddo
+
 !
 !     solve tridiagonal problem for momentum
 !

physics/moninedmf.meta

@@ -513,6 +513,15 @@
   type = integer
   intent = in
   optional = F
+[dkudiagnostic]
+  standard_name = atmosphere_momentum_diffusivity
+  long_name = diffusivity for momentum
+  units = m2 s-1
+  dimensions = (horizontal_loop_extent,vertical_dimension_minus_one)
+  type = real
+  kind = kind_phys
+  intent = out
+  optional = F
 [var_ric]
   standard_name = flag_variable_bulk_richardson_number
   long_name = flag for calculating variable bulk richardson number for hurricane PBL
@@ -538,6 +547,8 @@
   dimensions = ()
   type = real
   kind = kind_phys
+  intent = in
+  optional = F
 [lssav]
   standard_name = flag_diagnostics
   long_name = logical flag for storing diagnostics

physics/satmedmfvdifq.F

@@ -65,6 +65,7 @@ subroutine satmedmfvdifq_run(im,km,ntrac,ntcw,ntiw,ntke,          &
      &     prsi,del,prsl,prslk,phii,phil,delt,                          &
      &     dspheat,dusfc,dvsfc,dtsfc,dqsfc,hpbl,                        &
      &     kinver,xkzm_m,xkzm_h,xkzm_s,dspfac,bl_upfr,bl_dnfr,          &
+     &     dkudiagnostic,                                               &
      &     ntoz,du3dt,dv3dt,dt3dt,dq3dt,do3dt,gen_tend,ldiag3d,qdiag3d, &
      &     errmsg,errflg)
 !
@@ -131,7 +132,8 @@ subroutine satmedmfvdifq_run(im,km,ntrac,ntcw,ntiw,ntke,          &
      &                     slx(im,km),  svx(im,km),   qtx(im,km),
      &                     tvx(im,km),  pix(im,km),   radx(im,km-1),
      &                     dku(im,km-1),dkt(im,km-1), dkq(im,km-1),
-     &                     cku(im,km-1),ckt(im,km-1)
+     &                     cku(im,km-1),ckt(im,km-1),
+     &                     dkudiagnostic(im,km-1)
 !
       real(kind=kind_phys) plyr(im,km), rhly(im,km),  cfly(im,km),
      &                     qstl(im,km)
@@ -1538,6 +1540,12 @@ subroutine satmedmfvdifq_run(im,km,ntrac,ntcw,ntiw,ntke,          &
 !
         enddo
       enddo
+
+      do k = 1,km1
+        do i=1,im
+           dkudiagnostic(i,k) = dku(i,k)
+        enddo
+      enddo
 c
 !> - Call tridi2() to solve tridiagonal problem for momentum
 c

physics/satmedmfvdifq.meta

@@ -593,6 +593,15 @@
   kind = kind_phys
   intent = in
   optional = F
+[dkudiagnostic]
+  standard_name = atmosphere_momentum_diffusivity
+  long_name = diffusivity for momentum
+  units = m2 s-1
+  dimensions = (horizontal_dimension,vertical_dimension_minus_one)
+  type = real
+  kind = kind_phys
+  intent = out
+  optional = F
 [ntoz]
   standard_name = index_for_ozone
   long_name = tracer index for ozone mixing ratio