Skip to content
Snippets Groups Projects
Commit d514f2d6 authored by Anna Shestakova's avatar Anna Shestakova
Browse files

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

parents 2c6f3887 039d3258
Branches
No related tags found
No related merge requests found
......@@ -82,7 +82,11 @@ use sfx_most, only: &
real,intent(in):: AR2(11),ARDIN(6)
real,intent(out)::ARDOUT(3)
<<<<<<< HEAD
real,parameter::zetalim = 1. !maximum value of z/L for stable SL
=======
real,parameter::zetalim = 2. !maximum value of z/L for stable SL
>>>>>>> 039d3258b422f93c9a44345c5152b600029afbf8
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
......@@ -99,7 +103,10 @@ use sfx_most, only: &
print *, zeta
if(zeta.gt.zetalim) zeta=zetalim
L = HIN/zeta
<<<<<<< HEAD
print *, zeta
=======
>>>>>>> 039d3258b422f93c9a44345c5152b600029afbf8
call get_psi_most(psi_m_hs,psi_h_hs,HIN/L)
call get_psi_most(psi_m,psi_h,HWIND/L)
......@@ -122,7 +129,11 @@ print *, zeta
subroutine pbldia_new_esm(AR2,ARDIN,ARDOUT)
real,intent(in):: AR2(11),ARDIN(6)
real,intent(out)::ARDOUT(3)
<<<<<<< HEAD
real,parameter::zetalim = 1. !maximum value of z/L for stable SL
=======
real,parameter::zetalim = 2. !maximum value of z/L for stable SL
>>>>>>> 039d3258b422f93c9a44345c5152b600029afbf8
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
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment