Index

Rundeck: P1SoM40

Created: Fri May 10 03:30:18 EDT 2024

pbl

Module: PBL_DRV

File: PBL_DRV.f

Summary: PBL contains code common for all surface. It calculates pbl profiles, for each surface type, to find the values of the various PBL variables at the surface, and accumulates diagnostics and output. It is called from within the subroutine SURFCE (for itype=1, 2 and 3, i.e., surface types ocean, seaice and landice respectively, in SURFACE.f), and from within the subroutine earth (for itype=4, i.e., surface type land, in GHY_DRV.f). Dynamic equations for the mean turbulent variables are integrated over npbl(=8) sublayers between the surface (sublayer 1) and the middle of the first GCM layer (sublayer npbl), using tridiagonal method.
Author : Greg. Hartke/Ye Cheng
Version:

---

Declaration:

subroutine pbl(i, j, ihc, itype, ptype, pbl_args, atm)

Variables:

ddms : used from clouds_com

downdraft mass flux in kg/(m^2 s), (i,j)

tdn1 : used from clouds_com

downdraft temperature in K, (i,j)

qdn1 : used from clouds_com

downdraft humidity in kg/kg, (i,j)

rgas : used from constant

grav : used from constant

omega2 : used from constant

deltx : used from constant

teeny : used from constant

lhe : used from constant

lhs : used from constant

planet_name : used from constant

sinlat2d : used from geom

pk : used from atm_com

dpdx_by_rho : used from atm_com

dpdy_by_rho : used from atm_com

dpdx_by_rho_0 : used from atm_com

dpdy_by_rho_0 : used from atm_com

ddm1 : used from clouds_com

ddml : used from clouds_com

npbl : used from socpbl

Initial Value = >n

zgs : used from socpbl

advanc : used from socpbl

j : integer, intent(in)

grid point

i : integer, intent(in)

grid point

height :

class (or 1)

itype : integer, intent(in)

surface type

ihc : integer, intent(in)

ptype : real*8, intent(in)

percent surface type

pbl_args : type(t_pbl_args)

ts : real*8

tmp : real*8

qsat : real*8

s1byg1 : real*8, parameter

Initial Value = .57735d0

ztop : real*8

coriol : real*8

rhosrf : real*8

qtop : real*8

utop : real*8

vtop : real*8

ufluxs : real*8

vfluxs : real*8

tfluxs : real*8

qfluxs : real*8

psitop : real*8

psisrf : real*8

vocean :

ocean/ice velocities for use in drag calulation

uocean :

ocean/ice velocities for use in drag calulation

evap_max :

maximal evaporation from unsaturated soil

fr_sat :

fraction of saturated soil

zs1 : real*8

= height of the first model layer (m)

tgv :

= virtual potential temperature of the ground (K)

tkv :

= virtual potential temperature of first model layer (K)

ws :

= magn. of GCM surf wind - ocean curr + buoyancy + gust (m/s)

psi : real*8

= angular diff. btw geostrophic and surface winds (rads)

wg :

= magnitude of the geostrophic wind (m/s)

hemi : real*8

= 1 for northern hemisphere, -1 for southern hemisphere

tg :

= bulk ground temperature (K)

elhx :

= latent heat for saturation humidity (J/kg)

dskin :

= skin-bulk SST difference (C)

qsol :

= solar heating (W/m2)

k : integer

pole :

= .TRUE. if at the north or south pole, .FALSE. otherwise

us :

= x component of surface wind, positive eastward (m/s)

vs :

= y component of surface wind, positive northward (m/s)

wsgcm :

= magnitude of the GCM surface wind - ocean currents (m/s)

wspdf :

= mean surface wind calculated from PDF of wind speed (m/s)

ws :

= magn. of GCM surf wind - ocean curr + buoyancy + gust (m/s)

tsv :

= virtual potential temperature of the surface (K)

qs :

= surface value of the specific moisture

dbl : real*8

= boundary layer height (m)

dbls : real*8

= stable boundary layer height (m)

ldbl : integer

= layer immediately above dbl

ldbls : integer

= layer immediately above dbls

kms : real*8

= momentum transport coefficient at ZGS (m**2/s)

khs :

= heat transport coefficient at ZGS (m**2/s)

khq :

= moist transport coefficient at ZGS (m**2/s)

ustar : real*8

= friction speed (square root of momentum flux) (m/s)

cm : real*8

= drag coefficient (dimensionless surface momentum flux)

ch : real*8

= Stanton number (dimensionless surface heat flux)

cq : real*8

= Dalton number (dimensionless surface moisture flux)

z0m : real*8

= roughness length for momentum,

z0h : real*8

= roughness length for temperature (m)

z0q : real*8

= roughness length for water vapor (m)

ug : real*8

= eastward component of the geostrophic wind (m/s)

vg : real*8

= northward component of the geostrophic wind (m/s)

mdf : real*8

= downdraft mass flux (m/s)

wint :

= integrated surface wind speed over sgs wind distribution

dtdt_gcm : real*8

temp. tendency from processes other than turbulence (K/s)

kqs : real*8

w2_1 : real*8

dpdxr : real*8

dpdyr : real*8

dpdxr0 : real*8

dpdyr0 : real*8

mdn : real*8

upbl : real*8, dimension(npbl)

vpbl : real*8, dimension(npbl)

tpbl : real*8, dimension(npbl)

qpbl : real*8, dimension(npbl)

epbl : real*8, dimension(npbl-1)

lmonin : real*8