Index | Rundeck: P1SoM40 | Created: Fri May 10 03:30:18 EDT 2024 |
k_gcm | Module: | File: ATURB_E1.f |
Author : Ye Cheng/G. Hartke
Version:
Declaration:
subroutine k_gcm(tvflx, qflx, ustar, wstar, dbl, ze, lscale, e, qturb, an2, as2, dtdz, dqdz, dudz, dvdz, kh, km, ke, wt, wq, w2, uw, vw, wt_nl, wq_nl, n)
- tvflx :
real*8, intent(in)
- virtual potential temperature flux at surface
- qflx :
real*8, intent(in)
- moisture flux at surface
- ustar :
real*8, intent(in)
- friction velocity
- wstar :
real*8, intent(in)
- Deardorff convective velocity scale
- dbl :
real*8, intent(in)
- the height of the PBL (real*8, in meters)
- kh :
real*8, dimension(n), intent(out)
- turbulent diffusivity for scalars (heat, moisture,...)
- ze :
real*8, dimension(n+1), intent(in)
- height (in meters) of layer edges
- lscale :
real*8, dimension(n), intent(in)
- turbulent length scale
- e :
real*8, dimension(n), intent(in)
- turbulent kinetic energy
- qturb :
real*8, dimension(n), intent(in)
- sqrt(2*e)
- an2 :
real*8, dimension(n), intent(in)
- g_alpha*dTdz
- as2 :
real*8, dimension(n), intent(in)
- (dUdz)**2+(dVdz)**2
- dtdz :
real*8, dimension(n), intent(in)
- dt/dz
- dqdz :
real*8, dimension(n), intent(in)
- dq/dz
- dudz :
real*8, dimension(n), intent(in)
- du/dz
- dvdz :
real*8, dimension(n), intent(in)
- dv/dz
- km :
real*8, dimension(n), intent(out)
- turbulent diffusivity for momentun
- ke :
real*8, dimension(n), intent(out)
- turbulent diffusivity for e
- w2 :
real*8, dimension(n), intent(out)
- vertical component of 2*e
- vw :
real*8, dimension(n), intent(out)
- turbulent fluxes
- uw :
real*8, dimension(n), intent(out)
- turbulent fluxes
- wq :
real*8, dimension(n), intent(out)
- turbulent fluxes
- wt :
real*8, dimension(n), intent(out)
- turbulent fluxes
- wt_nl :
real*8, dimension(n), intent(out)
- non-local part of heat flux wt
- wq_nl :
real*8, dimension(n), intent(out)
- non-local part of moisture flux wq
- wc_nl :
- non-local part of tracer flux
- n :
integer, intent(in)
- number of layers
- teeny : used from constant
- by3 : used from constant
- kappa : used from socpbl
- prt : used from socpbl
- ghmin : used from socpbl
- ghmax : used from socpbl
- d1 : used from socpbl
- d2 : used from socpbl
- d3 : used from socpbl
- d4 : used from socpbl
- d5 : used from socpbl
- s0 : used from socpbl
- s1 : used from socpbl
- s2 : used from socpbl
- s4 : used from socpbl
- s5 : used from socpbl
- s6 : used from socpbl
- s7 : used from socpbl
- s8 : used from socpbl
- b1 : used from socpbl
- g5 : used from socpbl
- kmmin :
real*8, parameter
- Initial Value
= 1.5d-5
- khmin :
real*8, parameter
- Initial Value
= 2.5d-5
- kmax :
real*8, parameter
- Initial Value
= 600.d0
- j :
integer
- loop variable
- tmp0 :
real*8
- tmp :
real*8
- tau :
real*8
- gm :
real*8
- gh :
real*8
- gmmax :
real*8
- byden :
real*8
- sm :
real*8
- sh :
real*8
- ustar2 :
real*8
- wstar3 :
real*8
- zzi :
real*8
- tau_pt :
real*8
- w2j :
real*8