cu_gf_deep Module Reference

This module contains the Grell_Freitas deep convection scheme. More...

Functions/Subroutines
integer function	my_maxloc1d (a, n)
subroutine	cu_gf_deep_run (itf, ktf, its, ite, kts, kte, dicycle, ichoice, ipr, ccn, ccnclean, dtime, imid, kpbl, dhdt, xland, zo, forcing, t, q, z1, tn, qo, po, psur, us, vs, rho, hfx, qfx, dx, mconv, omeg, csum, cnvwt, zuo, zdo, zdm, edto, edtm, xmb_out, xmbm_in, xmbs_in, pre, outu, outv, outt, outq, outqc, kbcon, ktop, cupclw, frh_out, ierr, ierrc, nchem, fscav, chem3d, wetdpc_deep, do_smoke_transport, rand_mom, rand_vmas, rand_clos, nranflag, do_capsuppress, cap_suppress_j, k22, jmin, kdt, mc_thresh)
	Driver for the deep or congestus GF routine.
subroutine	fct1d3 (ktop, n, dt, z, tracr, massflx, trflx_in, dellac, g)
	Calculates tracer fluxes due to subsidence, only up-stream differencing is currently used but flux corrected transport can be turn on.
subroutine	rain_evap_below_cloudbase (itf, ktf, its, ite, kts, kte, ierr, kbcon, xmb, psur, xland, qo_cup, po_cup, qes_cup, pwavo, edto, pwevo, pre, outt, outq)
	Calculates rain evaporation below cloud base.
subroutine	cup_dd_edt (ierr, us, vs, z, ktop, kbcon, edt, p, pwav, pw, ccn, ccnclean, pwev, edtmax, edtmin, edtc, psum2, psumh, rho, aeroevap, pefc, xland1, itf, ktf, its, ite, kts, kte)
	Calculates strength of downdraft based on windshear and/or aerosol content.
subroutine	cup_dd_moisture (ierrc, zd, hcd, hes_cup, qcd, qes_cup, pwd, q_cup, z_cup, dd_massentr, dd_massdetr, jmin, ierr, gamma_cup, pwev, bu, qrcd, p_cup, q, he, iloop, itf, ktf, its, ite, kts, kte)
	Calcultes moisture properties of downdrafts.
subroutine	cup_env (z, qes, he, hes, t, q, p, z1, psur, ierr, tcrit, itest, itf, ktf, its, ite, kts, kte)
	Calculates environmental moist static energy, saturation moist static energy, heights, and saturation mixing ratio.
subroutine	cup_env_clev (t, qes, q, he, hes, z, p, qes_cup, q_cup, he_cup, hes_cup, z_cup, p_cup, gamma_cup, t_cup, psur, ierr, z1, itf, ktf, its, ite, kts, kte)
	Calculates environmental values on cloud levels.
subroutine	cup_forcing_ens_3d (closure_n, xland, aa0, aa1, xaa0, mbdt, dtime, ierr, ierr2, ierr3, xf_ens, axx, forcing, maxens3, mconv, rand_clos, p_cup, ktop, omeg, zd, zdm, k22, zu, pr_ens, edt, edtm, kbcon, ichoice, imid, ipr, itf, ktf, its, ite, kts, kte, dicycle, tau_ecmwf, aa1_bl, xf_dicycle)
	Calculates an ensemble of closures and the resulting ensemble average to determine cloud base mass-flux.
subroutine	cup_kbcon (ierrc, cap_inc, iloop_in, k22, kbcon, he_cup, hes_cup, hkb, ierr, kbmax, p_cup, cap_max, ztexec, zqexec, jprnt, itf, ktf, its, ite, kts, kte, z_cup, entr_rate, heo, imid)
	Calculates the level of convective cloud base.
subroutine	cup_maximi (array, ks, ke, maxx, ierr, itf, ktf, its, ite, kts, kte)
	Calculates the level at which the maximum value in an array occurs.
subroutine	cup_minimi (array, ks, kend, kt, ierr, itf, ktf, its, ite, kts, kte)
	Calculates the level at which the minimum value in an array occurs.
subroutine	cup_up_aa0 (aa0, z, zu, dby, gamma_cup, t_cup, kbcon, ktop, ierr, itf, ktf, its, ite, kts, kte)
	Calculates the cloud work functions for updrafts.
subroutine	neg_check (name, j, dt, q, outq, outt, outu, outv, outqc, pret, its, ite, kts, kte, itf, ktf, ktop)
	Checks for negative or excessive tendencies and corrects in a mass conversing way by adjusting the cloud base mass-flux.
subroutine	cup_output_ens_3d (xff_mid, xf_ens, ierr, dellat, dellaq, dellaqc, outtem, outq, outqc, dx, zu, pre, pw, xmb, ktop, edt, pwd, name, ierr2, ierr3, p_cup, pr_ens, maxens3, sig, closure_n, xland1, xmbm_in, xmbs_in, rrfs_factor, ichoice, imid, ipr, itf, ktf, its, ite, kts, kte, dicycle, xf_dicycle)
	This subroutine calculates final output fields including physical tendencies, precipitation, and mass-flux.
subroutine	cup_up_moisture (name, ierr, z_cup, qc, qrc, pw, pwav, p_cup, kbcon, ktop, dby, clw_all, xland1, q, gamma_cup, zu, qes_cup, k22, qe_cup, c0, c0t3d, zqexec, ccn, ccnclean, rho, c1d, t, autoconv, up_massentr, up_massdetr, psum, psumh, itest, itf, ktf, its, ite, kts, kte)
	Calculates moisture properties of the updraft.
real function	satvap (temp2)
	Calculates saturation vapor pressure.
subroutine	get_cloud_bc (mzp, array, x_aver, k22, add)
	Calculates the average value of a variable at the updraft originating level.
subroutine	rates_up_pdf (rand_vmas, ipr, name, ktop, ierr, p_cup, entr_rate_2d, hkbo, heo, heso_cup, z_cup, xland, kstabi, k22, kbcon, its, ite, itf, kts, kte, ktf, zuo, kpbl, ktopdby, csum, pmin_lev)
	Driver for the normalized mass-flux routine.
subroutine	get_zu_zd_pdf_fim (kklev, p, rand_vmas, zubeg, ipr, xland, zuh2, draft, ierr, kb, kt, zu, kts, kte, ktf, max_mass, kpbli, csum, pmin_lev)
	Calculates a normalized mass-flux profile for updrafts and downdrafts using the beta function.
subroutine	cup_up_aa1bl (aa0, t, tn, q, qo, dtime, z_cup, zu, dby, gamma_cup, t_cup, kbcon, ktop, ierr, itf, ktf, its, ite, kts, kte)
	Calculates the cloud work function based on boundary layer forcing.
subroutine	get_inversion_layers (ierr, p_cup, t_cup, z_cup, qo_cup, qeso_cup, k_inv_layers, kstart, kend, dtempdz, itf, ktf, its, ite, kts, kte)
	Finds temperature inversions using the first and second derivative of temperature.
subroutine	get_lateral_massflux (itf, ktf, its, ite, kts, kte, ierr, ktop, zo_cup, zuo, cd, entr_rate_2d, up_massentro, up_massdetro, up_massentr, up_massdetr, draft, kbcon, k22, up_massentru, up_massdetru, lambau)
	Calculates mass entranment and detrainment rates.
subroutine	get_partition_liq_ice (ierr, tn, po_cup, p_liq_ice, melting_layer, itf, ktf, its, ite, kts, kte, cumulus)
	Calculates the partition between cloud water and cloud ice.
subroutine	get_melting_profile (ierr, tn_cup, po_cup, p_liq_ice, melting_layer, qrco, pwo, edto, pwdo, melting, itf, ktf, its, ite, kts, kte, cumulus)
	Calculates the melting profile.
subroutine	get_cloud_top (name, ktop, ierr, p_cup, entr_rate_2d, hkbo, heo, heso_cup, z_cup, kstabi, k22, kbcon, its, ite, itf, kts, kte, ktf, zuo, kpbl, klcl, hcot)
	Calculates the cloud top height.

Variables
real(kind=kind_phys), parameter	g =9.81
real(kind=kind_phys), parameter	cp =1004.
real(kind=kind_phys), parameter	xlv =2.5e6
real(kind=kind_phys), parameter	r_v =461.
real(kind=kind_phys), parameter	tcrit =258.
real(kind=kind_phys), parameter	c1 = 0.003
	tuning constant for cloudwater/ice detrainment
integer, parameter	irainevap =1
	parameter to turn on or off evaporation of rainwater as done in sas
real(kind=kind_phys), parameter	frh_thresh = .9
	max allowed fractional coverage (frh_thresh)
real(kind=kind_phys), parameter	rh_thresh = .97
	rh threshold. if fractional coverage ~ frh_thres, do not use cupa any further
real(kind=kind_phys), parameter	betajb =1.2
	tuning constant for j. brown closure (ichoice = 4,5,6)
integer, parameter	use_excess =0
	tuning for shallow and mid convection. ec uses 1.5
real(kind=kind_phys), parameter	fluxtune =1.5
real(kind=kind_phys), parameter	pgcd = 0.1
	flag to turn off or modify mom transport by downdrafts
integer, parameter	autoconv =1
	aerosol awareness, do not use yet!
integer, parameter	aeroevap =1
real(kind=kind_phys), parameter	scav_factor = 0.5
real(kind=kind_phys), parameter	dx_thresh = 6500.
integer, parameter	maxens3 =16
	still 16 ensembles for clousres
logical, parameter	melt_glac = .true.
real(kind=kind_phys), parameter	t_0 = 273.16
	k
real(kind=kind_phys), parameter	t_ice = 250.16
	k
real(kind=kind_phys), parameter	xlf = 0.333e6
	latent heat of freezing (j k-1 kg-1)
real(kind=kind_phys), parameter	qrc_crit = 2.e-4

Detailed Description

This module contains the Grell_Freitas deep convection scheme.

Variable Documentation

aeroevap

integer, parameter cu_gf_deep::aeroevap =1

Definition at line 31 of file cu_gf_deep.F90.

autoconv

integer, parameter cu_gf_deep::autoconv =1

aerosol awareness, do not use yet!

Definition at line 30 of file cu_gf_deep.F90.

betajb

real(kind=kind_phys), parameter cu_gf_deep::betajb =1.2

tuning constant for j. brown closure (ichoice = 4,5,6)

Definition at line 22 of file cu_gf_deep.F90.

c1

real(kind=kind_phys), parameter cu_gf_deep::c1 = 0.003

tuning constant for cloudwater/ice detrainment

Definition at line 14 of file cu_gf_deep.F90.

cp

real(kind=kind_phys), parameter cu_gf_deep::cp =1004.

Definition at line 9 of file cu_gf_deep.F90.

dx_thresh

real(kind=kind_phys), parameter cu_gf_deep::dx_thresh = 6500.

Definition at line 34 of file cu_gf_deep.F90.

fluxtune

real(kind=kind_phys), parameter cu_gf_deep::fluxtune =1.5

Definition at line 25 of file cu_gf_deep.F90.

frh_thresh

real(kind=kind_phys), parameter cu_gf_deep::frh_thresh = .9

max allowed fractional coverage (frh_thresh)

Definition at line 18 of file cu_gf_deep.F90.

g

real(kind=kind_phys), parameter cu_gf_deep::g =9.81

Definition at line 8 of file cu_gf_deep.F90.

irainevap

integer, parameter cu_gf_deep::irainevap =1

parameter to turn on or off evaporation of rainwater as done in sas

Definition at line 16 of file cu_gf_deep.F90.

maxens3

integer, parameter cu_gf_deep::maxens3 =16

still 16 ensembles for clousres

Definition at line 36 of file cu_gf_deep.F90.

melt_glac

logical, parameter cu_gf_deep::melt_glac = .true.

• turn on/off ice phase/melting

Definition at line 39 of file cu_gf_deep.F90.

pgcd

real(kind=kind_phys), parameter cu_gf_deep::pgcd = 0.1

flag to turn off or modify mom transport by downdrafts

Definition at line 27 of file cu_gf_deep.F90.

qrc_crit

real(kind=kind_phys), parameter cu_gf_deep::qrc_crit = 2.e-4

Definition at line 46 of file cu_gf_deep.F90.

r_v

real(kind=kind_phys), parameter cu_gf_deep::r_v =461.

Definition at line 11 of file cu_gf_deep.F90.

rh_thresh

real(kind=kind_phys), parameter cu_gf_deep::rh_thresh = .97

rh threshold. if fractional coverage ~ frh_thres, do not use cupa any further

Definition at line 20 of file cu_gf_deep.F90.

scav_factor

real(kind=kind_phys), parameter cu_gf_deep::scav_factor = 0.5

Definition at line 32 of file cu_gf_deep.F90.

t_0

real(kind=kind_phys), parameter cu_gf_deep::t_0 = 273.16

k

Definition at line 40 of file cu_gf_deep.F90.

t_ice

real(kind=kind_phys), parameter cu_gf_deep::t_ice = 250.16

k

Definition at line 40 of file cu_gf_deep.F90.

tcrit

real(kind=kind_phys), parameter cu_gf_deep::tcrit =258.

Definition at line 12 of file cu_gf_deep.F90.

use_excess

integer, parameter cu_gf_deep::use_excess =0

tuning for shallow and mid convection. ec uses 1.5

Definition at line 24 of file cu_gf_deep.F90.

xlf

real(kind=kind_phys), parameter cu_gf_deep::xlf = 0.333e6

latent heat of freezing (j k-1 kg-1)

Definition at line 40 of file cu_gf_deep.F90.

xlv

real(kind=kind_phys), parameter cu_gf_deep::xlv =2.5e6

Definition at line 10 of file cu_gf_deep.F90.