diff --git a/diag/pbldia_new_sfx.f90 b/diag/pbldia_new_sfx.f90
index 3ccf023a9ba3e6470ffa5ad3391339540e1fa69f..4a645f13f1b8d1848fde66e36f6b1653600eddc9 100755
--- a/diag/pbldia_new_sfx.f90
+++ b/diag/pbldia_new_sfx.f90
@@ -96,10 +96,8 @@ 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)