diff --git a/diag/pbldia_new_sfx.f90 b/diag/pbldia_new_sfx.f90
index bbcc8cdaa498c4b03e3c4cc29579de2d419831dd..83f1b4000fbd79e3485cad36144ff8dd10b5f854 100755
--- a/diag/pbldia_new_sfx.f90
+++ b/diag/pbldia_new_sfx.f90
@@ -1,6 +1,5 @@
 module pbldia_new_sfx
 
-use sfx_phys_const
 
 implicit none 
 real,parameter,private::undef=-999.0,R=287.,appa=0.287
@@ -54,7 +53,7 @@ use sfx_sheba, only: &
       real,intent(in):: AR2(11),ARDIN(6)!ARDIN(10)
       real,intent(out)::ARDOUT(3) !ARDOUT(4)
 !      real,parameter::Llim = 50.0 !minimum value of L for stable SL (25)
-      real,parameter::zetalim = 2. !minimum value of L for stable SL (25)
+      real,parameter::zetalim = 2. !maximum value of z/L for stable SL
       real,parameter::kappa = 0.4 
       real psi_m,psi_h,psi_m_hs,psi_h_hs
       real hwind,htemp,ustar,tstar,qstar,&
@@ -100,7 +99,7 @@ use sfx_most, only: &
       real,intent(in):: AR2(11),ARDIN(6)!ARDIN(10)
       real,intent(out)::ARDOUT(3) !ARDOUT(4)
 !      real,parameter::Llim = 50.0 !minimum value of L for stable SL (25)
-      real,parameter::zetalim = 2. !minimum value of L for stable SL (25)
+      real,parameter::zetalim = 2. !maximum value of z/L for stable SL
       real,parameter::kappa = 0.4 
       real psi_m,psi_h,psi_m_hs,psi_h_hs
       real hwind,htemp,ustar,tstar,qstar,&
@@ -147,7 +146,7 @@ subroutine pbldia_new_esm(AR2,ARDIN,ARDOUT)
       real,intent(in):: AR2(11),ARDIN(6)!ARDIN(10)
       real,intent(out)::ARDOUT(3) !ARDOUT(4)
 !      real,parameter::Llim = 50.0 !minimum value of L for stable SL (25)
-      real,parameter::zetalim = 2. !minimum value of L for stable SL (25)
+      real,parameter::zetalim = 2. !maximum value of z/L for stable SL
       real,parameter::kappa = 0.4 
       real psi_m,psi_h,psi_m_hs,psi_h_hs
       real hwind,htemp,ustar,tstar,qstar,&