diff --git a/diag/pbldia_new_sfx.f90 b/diag/pbldia_new_sfx.f90
index 3ac1d887f040723ffd9dcf0a2a449082594e9d86..3ccf023a9ba3e6470ffa5ad3391339540e1fa69f 100755
--- a/diag/pbldia_new_sfx.f90
+++ b/diag/pbldia_new_sfx.f90
@@ -82,7 +82,7 @@ use sfx_most, only: &
 
       real,intent(in):: AR2(11),ARDIN(6)
       real,intent(out)::ARDOUT(3)
-      real,parameter::zetalim = 2. !maximum value of z/L for stable SL
+      real,parameter::zetalim = 1. !maximum value of z/L for stable SL
       real psi_m,psi_h,psi_m_hs,psi_h_hs
       real hwind,htemp,ustar,tstar,qstar,&
 & ufwind,uftemp,dteta,dq,wind,L,HIN,zeta
@@ -96,8 +96,10 @@ use sfx_most, only: &
       TSTAR = AR2(9) * ARDIN(2)
       QSTAR = AR2(9) * ARDIN(3)
 
+print *, zeta
      if(zeta.gt.zetalim) zeta=zetalim
       L = HIN/zeta
+print *, zeta
 
      call get_psi_most(psi_m_hs,psi_h_hs,HIN/L)
      call get_psi_most(psi_m,psi_h,HWIND/L)
@@ -120,7 +122,7 @@ use sfx_most, only: &
 subroutine pbldia_new_esm(AR2,ARDIN,ARDOUT)
       real,intent(in):: AR2(11),ARDIN(6)
       real,intent(out)::ARDOUT(3)
-      real,parameter::zetalim = 2. !maximum value of z/L for stable SL
+      real,parameter::zetalim = 1. !maximum value of z/L for stable SL
       real psi_m,psi_h,psi_m_hs,psi_h_hs
       real hwind,htemp,ustar,tstar,qstar,&
 & ufwind,uftemp,dteta,dq,wind,L,HIN,zeta