diff --git a/diag/pbldia_new_sfx.f90 b/diag/pbldia_new_sfx.f90
index bbcc8cdaa498c4b03e3c4cc29579de2d419831dd..fddff40b5d96a9b91eea69d0e71b278db51c57c8 100755
--- a/diag/pbldia_new_sfx.f90
+++ b/diag/pbldia_new_sfx.f90
@@ -54,7 +54,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 +100,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 +147,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,&