diff --git a/ref/src/cu_gf_deep.F90 b/ref/src/cu_gf_deep.F90
index ae11695..be9cf0c 100644
--- a/ref/src/cu_gf_deep.F90
+++ b/ref/src/cu_gf_deep.F90
@@ -494,6 +494,7 @@ subroutine cu_gf_deep_run(        &
          !- zws for shallow convection closure (grant 2001)
          !- height of the pbl
          zws(i) = max(0.,.001-flux_tun(i)*0.41*buo_flux*zo(i,kpbl(i))*g/t(i,kpbl(i)))
+         write(*,'(I5,E20.10)') i, zws(i)
          zws(i) = 1.2*zws(i)**.3333
          zws(i) = zws(i)*rho(i,kpbl(i)) !check if zrho is correct
       enddo