From 99c68e3599f7c62cbea916c79342a3e62f9e87d2 Mon Sep 17 00:00:00 2001
From: Anna Shestakova <shestakova.aa.92@gmail.com>
Date: Sun, 27 Oct 2024 23:00:37 +0300
Subject: [PATCH] diag module for interpolation in surface layer
---
diag/pbldia_new_sfx.f90 | 7 +++----
1 file changed, 3 insertions(+), 4 deletions(-)
diff --git a/diag/pbldia_new_sfx.f90 b/diag/pbldia_new_sfx.f90
index bbcc8cd..83f1b40 100755
--- a/diag/pbldia_new_sfx.f90
+++ b/diag/pbldia_new_sfx.f90
@@ -1,6 +1,5 @@
module pbldia_new_sfx
-use sfx_phys_const
implicit none
real,parameter,private::undef=-999.0,R=287.,appa=0.287
@@ -54,7 +53,7 @@ use sfx_sheba, only: &
real,intent(in):: AR2(11),ARDIN(6)!ARDIN(10)
real,intent(out)::ARDOUT(3) !ARDOUT(4)
! real,parameter::Llim = 50.0 !minimum value of L for stable SL (25)
- real,parameter::zetalim = 2. !minimum value of L for stable SL (25)
+ real,parameter::zetalim = 2. !maximum value of z/L for stable SL
real,parameter::kappa = 0.4
real psi_m,psi_h,psi_m_hs,psi_h_hs
real hwind,htemp,ustar,tstar,qstar,&
@@ -100,7 +99,7 @@ use sfx_most, only: &
real,intent(in):: AR2(11),ARDIN(6)!ARDIN(10)
real,intent(out)::ARDOUT(3) !ARDOUT(4)
! real,parameter::Llim = 50.0 !minimum value of L for stable SL (25)
- real,parameter::zetalim = 2. !minimum value of L for stable SL (25)
+ real,parameter::zetalim = 2. !maximum value of z/L for stable SL
real,parameter::kappa = 0.4
real psi_m,psi_h,psi_m_hs,psi_h_hs
real hwind,htemp,ustar,tstar,qstar,&
@@ -147,7 +146,7 @@ subroutine pbldia_new_esm(AR2,ARDIN,ARDOUT)
real,intent(in):: AR2(11),ARDIN(6)!ARDIN(10)
real,intent(out)::ARDOUT(3) !ARDOUT(4)
! real,parameter::Llim = 50.0 !minimum value of L for stable SL (25)
- real,parameter::zetalim = 2. !minimum value of L for stable SL (25)
+ real,parameter::zetalim = 2. !maximum value of z/L for stable SL
real,parameter::kappa = 0.4
real psi_m,psi_h,psi_m_hs,psi_h_hs
real hwind,htemp,ustar,tstar,qstar,&
--
GitLab