<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 , &
br = br_p , dt = dt_pbl , kpbl2d = kpbl_p , &
exch_h = exch_p , u10 = u10_p , v10 = v10_p , &
+ kzhout = kzh_p , kzmout = kzm_p , kzqout = kzq_p , &
+ rho = rho_p , &
ids = ids , ide = ide , jds = jds , jde = jde , kds = kds , kde = kde , &
ims = ims , ime = ime , jms = jms , jme = jme , kms = kms , kme = kme , &
its = its , ite = ite , jts = jts , jte = jte , kts = kts , kte = kte &
@@ -221,32 +239,13 @@
qfx = qfx_p , gz1oz0 = gz1oz0_p , wspd = wspd_p , &
br = br_p , dt = dt_pbl , kpbl2d = kpbl_p , &
exch_h = exch_p , u10 = u10_p , v10 = v10_p , &
+ kzhout = kzh_p , kzmout = kzm_p , kzqout = kzq_p , &
+ rho = rho_p , &
ids = ids , ide = ide , jds = jds , jde = jde , kds = kds , kde = kde , &
ims = ims , ime = ime , jms = jms , jme = jme , kms = kms , kme = kme , &
its = its , ite = ite , jts = jts , jte = jte , kts = kts , kte = kte &
)
#endif
-!... CALL FROM REVISION 1721:
-! call ysu ( &
-! u3d = u_p , v3d = v_p , th3d = th_p , &
-! t3d = t_p , qv3d = qv_p , qc3d = qc_p , &
-! qi3d = qi_p , p3d = pres_p , p3di = pres2_p , &
-! pi3d = pi_p , rublten = rublten_p , rvblten = rvblten_p , &
-! rthblten = rthblten_p , rqvblten = rqvblten_p , rqcblten = rqcblten_p , &
-! rqiblten = rqiblten_p , flag_qi = f_qi , cp = cp , &
-! g = g , rovcp = rcp , rd = R_d , &
-! rovg = rdg , ep1 = ep_1 , ep2 = ep_2 , &
-! karman = karman , xlv = xlv , rv = R_v , &
-! dz8w = dz_p , psfc = psfc_p , znt = znt_p , &
-! ust = ust_p , hpbl = hpbl_p , psim = psim_p , &
-! psih = psih_p , xland = xland_p , hfx = hfx_p , &
-! qfx = qfx_p , gz1oz0 = gz1oz0_p , wspd = wspd_p , &
-! br = br_p , dt = dt_pbl , kpbl2d = kpbl_p , &
-! exch_h = exch_p , u10 = u10_p , v10 = v10_p , &
-! ids = ids , ide = ide , jds = jds , jde = jde , kds = kds , kde = kde , &
-! ims = ims , ime = ime , jms = jms , jme = jme , kms = kms , kme = kme , &
-! its = its , ite = ite , jts = jts , jte = jte , kts = kts , kte = kte &
-! )
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, &!
rqiblten_p !
+ real(kind=RKIND),dimension(:,:,:),allocatable:: &
+ kzh_p, &!
+ kzm_p, &!
+ kzq_p !
+
!=============================================================================================
!... variables and arrays related to parameterization of surface layer:
!=============================================================================================
</font>
</pre>