<p><b>laura@ucar.edu</b> 2012-08-06 09:42:25 -0600 (Mon, 06 Aug 2012)</p><p>added kzh,kzq, and kzm as diagnostics to debug the YSU PBL scheme.<br>
</p><hr noshade><pre><font color="gray">Modified: branches/atmos_physics/src/core_atmos_physics/mpas_atmphys_driver_pbl.F
===================================================================
--- branches/atmos_physics/src/core_atmos_physics/mpas_atmphys_driver_pbl.F        2012-08-06 15:32:18 UTC (rev 2086)
+++ branches/atmos_physics/src/core_atmos_physics/mpas_atmphys_driver_pbl.F        2012-08-06 15:42:25 UTC (rev 2087)
@@ -46,6 +46,10 @@
  if(.not.allocated(rqvblten_p)) allocate(rqvblten_p(ims:ime,kms:kme,jms:jme) )
  if(.not.allocated(rqcblten_p)) allocate(rqcblten_p(ims:ime,kms:kme,jms:jme) )
  if(.not.allocated(rqiblten_p)) allocate(rqiblten_p(ims:ime,kms:kme,jms:jme) )
+!temporary for debugging the YSU PBL scheme:
+ if(.not.allocated(kzh_p)     ) allocate(kzh_p(ims:ime,kms:kme,jms:jme)      )
+ if(.not.allocated(kzm_p)     ) allocate(kzm_p(ims:ime,kms:kme,jms:jme)      )
+ if(.not.allocated(kzq_p)     ) allocate(kzq_p(ims:ime,kms:kme,jms:jme)      )
 
  end subroutine allocate_pbl
 
@@ -77,6 +81,10 @@
  if(allocated(rqvblten_p)) deallocate(rqvblten_p )
  if(allocated(rqcblten_p)) deallocate(rqcblten_p )
  if(allocated(rqiblten_p)) deallocate(rqiblten_p )
+!temporary for debugging the YSU PBL scheme:
+ if(allocated(kzh_p)     ) deallocate(kzh_p      )
+ if(allocated(kzm_p)     ) deallocate(kzm_p      )
+ if(allocated(kzq_p)     ) deallocate(kzq_p      )
 
  end subroutine deallocate_pbl
 
@@ -122,6 +130,10 @@
     rqvblten_p(i,k,j) = 0.
     rqcblten_p(i,k,j) = 0.
     rqiblten_p(i,k,j) = 0.
+!temporary for debugging the YSU PBL scheme:
+    kzh_p(i,k,j)    = 0._RKIND
+    kzm_p(i,k,j)    = 0._RKIND
+    kzq_p(i,k,j)    = 0._RKIND
  enddo
  enddo
  enddo
@@ -155,6 +167,10 @@
     tend_physics % rqvblten % array(k,i) = rqvblten_p(i,k,j)
     tend_physics % rqcblten % array(k,i) = rqcblten_p(i,k,j)
     tend_physics % rqiblten % array(k,i) = rqiblten_p(i,k,j)
+!temporary for debugging the YSU PBL scheme:
+    diag_physics % kzh % array(k,i) = kzh_p(i,k,j)
+    diag_physics % kzm % array(k,i) = kzm_p(i,k,j)
+    diag_physics % kzq % array(k,i) = kzq_p(i,k,j)
  enddo
  enddo
  enddo
@@ -199,6 +215,8 @@
                  qfx      = qfx_p      , gz1oz0   = gz1oz0_p    , wspd     = wspd_p     , &amp;
                  br       = br_p       , dt       = dt_pbl      , kpbl2d   = kpbl_p     , &amp;
                  exch_h   = exch_p     , u10      = u10_p       , v10      = v10_p      , &amp;
+                 kzhout   = kzh_p      , kzmout   = kzm_p       , kzqout   = kzq_p      , &amp;
+                 rho      = rho_p      ,                                                  &amp;
                  ids = ids , ide = ide , jds = jds , jde = jde , kds = kds , kde = kde  , &amp;
                  ims = ims , ime = ime , jms = jms , jme = jme , kms = kms , kme = kme  , &amp;
                  its = its , ite = ite , jts = jts , jte = jte , kts = kts , kte = kte    &amp;
@@ -221,32 +239,13 @@
                  qfx      = qfx_p      , gz1oz0   = gz1oz0_p    , wspd     = wspd_p     , &amp;
                  br       = br_p       , dt       = dt_pbl      , kpbl2d   = kpbl_p     , &amp;
                  exch_h   = exch_p     , u10      = u10_p       , v10      = v10_p      , &amp;
+                 kzhout   = kzh_p      , kzmout   = kzm_p       , kzqout   = kzq_p      , &amp;
+                 rho      = rho_p      ,                                                  &amp;
                  ids = ids , ide = ide , jds = jds , jde = jde , kds = kds , kde = kde  , &amp;
                  ims = ims , ime = ime , jms = jms , jme = jme , kms = kms , kme = kme  , &amp;
                  its = its , ite = ite , jts = jts , jte = jte , kts = kts , kte = kte    &amp;
                 )
 #endif
-!... CALL FROM REVISION 1721:
-!      call ysu ( &amp;
-!                u3d      = u_p        , v3d      = v_p         , th3d     = th_p       , &amp;
-!                t3d      = t_p        , qv3d     = qv_p        , qc3d     = qc_p       , &amp;
-!                qi3d     = qi_p       , p3d      = pres_p      , p3di     = pres2_p    , &amp;
-!                pi3d     = pi_p       , rublten  = rublten_p   , rvblten  = rvblten_p  , &amp;
-!                rthblten = rthblten_p , rqvblten = rqvblten_p  , rqcblten = rqcblten_p , &amp;
-!                rqiblten = rqiblten_p , flag_qi  = f_qi        , cp       = cp         , &amp;
-!                g        = g          , rovcp    = rcp         , rd       = R_d        , &amp;
-!                rovg     = rdg        , ep1      = ep_1        , ep2      = ep_2       , &amp;
-!                karman   = karman     , xlv      = xlv         , rv       = R_v        , &amp;
-!                dz8w     = dz_p       , psfc     = psfc_p      , znt      = znt_p      , &amp;
-!                ust      = ust_p      , hpbl     = hpbl_p      , psim     = psim_p     , &amp;
-!                psih     = psih_p     , xland    = xland_p     , hfx      = hfx_p      , &amp;
-!                qfx      = qfx_p      , gz1oz0   = gz1oz0_p    , wspd     = wspd_p     , &amp;
-!                br       = br_p       , dt       = dt_pbl      , kpbl2d   = kpbl_p     , &amp;
-!                exch_h   = exch_p     , u10      = u10_p       , v10      = v10_p      , &amp;
-!                ids = ids , ide = ide , jds = jds , jde = jde , kds = kds , kde = kde  , &amp;
-!                ims = ims , ime = ime , jms = jms , jme = jme , kms = kms , kme = kme  , &amp;
-!                its = its , ite = ite , jts = jts , jte = jte , kts = kts , kte = kte    &amp;
-!               )
 
      case default
 
@@ -254,7 +253,7 @@
 
 !copy all arrays back to the MPAS grid:
  call pbl_to_MPAS(diag_physics,tend_physics)
- write(0,*) '--- enter subroutine driver_pbl'
+ write(0,*) '--- end subroutine driver_pbl'
 
  end subroutine driver_pbl
 

Modified: branches/atmos_physics/src/core_atmos_physics/mpas_atmphys_vars.F
===================================================================
--- branches/atmos_physics/src/core_atmos_physics/mpas_atmphys_vars.F        2012-08-06 15:32:18 UTC (rev 2086)
+++ branches/atmos_physics/src/core_atmos_physics/mpas_atmphys_vars.F        2012-08-06 15:42:25 UTC (rev 2087)
@@ -205,6 +205,11 @@
     rqcblten_p,       &amp;!
     rqiblten_p         !
 
+ real(kind=RKIND),dimension(:,:,:),allocatable:: &amp;
+    kzh_p,            &amp;!
+    kzm_p,            &amp;!
+    kzq_p              !
+
 !=============================================================================================
 !... variables and arrays related to parameterization of surface layer:
 !=============================================================================================

</font>
</pre>