Merge branch 'diag_new' of http://tesla.parallel.ru/INMCM60_pbl/INMCM60_sfx into diag_new

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