CCPP SciDoc  v6.0.0
Common Community Physics Package Developed at DTC
GFDL In-Core Fast Saturation Adjustment Module

The subroutine 'fv_sat_adj' implements the fast processes in the GFDL Cloud MP. It is part of the GFDL Cloud MP. More...

Functions/Subroutines

subroutine, public fv_sat_adj::fv_sat_adj_run (mdt, zvir, is, ie, isd, ied, kmp, km, kmdelz, js, je, jsd, jed, ng, hydrostatic, fast_mp_consv, te0_2d, te0, ngas, qvi, qv, ql, qi, qr, qs, qg, hs, peln, delz, delp, pt, pkz, q_con, akap, cappa, area, dtdt, out_dt, last_step, do_qa, qa, nthreads, errmsg, errflg)
 
real(kind=kind_dyn) function fv_sat_adj::wqs1 (ta, den)
 the function 'wqs1' computes the saturated specific humidity for table ii. More...
 
real(kind=kind_dyn) function fv_sat_adj::iqs1 (ta, den)
 the function 'wqs1' computes the saturated specific humidity for table iii More...
 
real(kind=kind_dyn) function fv_sat_adj::wqs2 (ta, den, dqdt)
 The function 'wqs2'computes the gradient of saturated specific humidity for table ii. More...
 
subroutine fv_sat_adj::wqs2_vect (is, ie, ta, den, wqsat, dqdt)
 The function wqs2_vect computes the gradient of saturated specific humidity for table ii. It is the same as "wqs2", but written as vector function. More...
 
real(kind=kind_dyn) function fv_sat_adj::iqs2 (ta, den, dqdt)
 The function 'iqs2' computes the gradient of saturated specific humidity for table iii. More...
 
subroutine fv_sat_adj::qs_table (n)
 saturation water vapor pressure table i More...
 
subroutine fv_sat_adj::qs_tablew (n)
 saturation water vapor pressure table ii. More...
 
subroutine fv_sat_adj::qs_table2 (n)
 saturation water vapor pressure table iii. More...
 
subroutine fv_sat_adj::fv_sat_adj_work (mdt, zvir, is, ie, js, je, ng, hydrostatic, consv_te, te0, ifdef MULTI_GASES
 This subroutine includes the entity of the fast saturation adjustment processes. More...
 

Detailed Description

Author
Shian-Jiann Lin, Linjiong Zhou

The subroutine 'fv_sat_adj' performs the fast processes in the GFDL microphysics.

This is designed for single-moment 6-class cloud microphysics schemes. It handles the heat release due to in situ phase changes.

Argument Table

fv_sat_adj_run argument table
local_namestandard_namelong_nameunitstypedimensionskindintent
mdttime_step_for_remapping_for_fast_physicsremapping time step for fast physicssreal()kind_dynin
zvirratio_of_vapor_to_dry_air_gas_constants_minus_one_default_kindzvir=rv/rd-1.0nonereal()kind_dynin
isstarting_x_direction_indexstarting X direction indexcountinteger()in
ieending_x_direction_indexending X direction indexcountinteger()in
isdstarting_x_direction_index_domainstarting X direction index for domaincountinteger()in
iedending_x_direction_index_domainending X direction index for domaincountinteger()in
kmptop_layer_index_for_fast_physicstop layer index for GFDL mpindexinteger()in
kmvertical_dimension_for_fast_physicsnumber of vertical levelscountinteger()in
kmdelzvertical_dimension_for_thickness_at_lagrangian_surfacevertical dimension for thickness at Lagrangian surfacecountinteger()in
jsstarting_y_direction_indexstarting Y direction indexcountinteger()in
jeending_y_direction_indexending Y direction indexcountinteger()in
jsdstarting_y_direction_index_domainstarting X direction index for domaincountinteger()in
jedending_y_direction_index_domainending X direction index for domaincountinteger()in
ngnumber_of_ghost_zonesnumber of ghost zones defined in fv_mpcountinteger()in
hydrostaticflag_for_hydrostatic_solver_for_fast_physicsflag for use the hydrostatic or nonhydrostatic solverflaglogical()in
fast_mp_consvflag_for_fast_microphysics_energy_conservationflag for fast microphysics energy conservationflaglogical()in
te0_2datmosphere_energy_content_in_columnatmosphere total energy in columnsJ m-2real(starting_x_direction_index:ending_x_direction_index, starting_y_direction_index:ending_y_direction_index)kind_dyninout
te0atmosphere_energy_content_at_lagrangian_surfaceatmosphere total energy at Lagrangian surfaceJ m-2real(starting_x_direction_index_domain:ending_x_direction_index_domain, starting_y_direction_index_domain:ending_y_direction_index_domain, 1:vertical_dimension_for_fast_physics)kind_dynout
ngasnumber_of_gases_for_multi_gases_physicsnumber of gases for multi gases physicscountinteger()in
qvigas_tracers_for_multi_gas_physics_at_lagrangian_surfacegas tracers for multi gas physics at Lagrangian surfacekg kg-1real(starting_x_direction_index_domain:ending_x_direction_index_domain, starting_y_direction_index_domain:ending_y_direction_index_domain, 1:vertical_dimension_for_fast_physics, 1:number_of_gases_for_multi_gases_physics)kind_dyninout
qvwater_vapor_specific_humidity_at_lagrangian_surfacewater vapor specific humidity updated by fast physics at Lagrangian surfacekg kg-1real(starting_x_direction_index_domain:ending_x_direction_index_domain, starting_y_direction_index_domain:ending_y_direction_index_domain, 1:vertical_dimension_for_fast_physics)kind_dyninout
qlcloud_liquid_water_specific_humidity_at_lagrangian_surfacecloud liquid water specific humidity updated by fast physics at Lagrangian surfacekg kg-1real(starting_x_direction_index_domain:ending_x_direction_index_domain, starting_y_direction_index_domain:ending_y_direction_index_domain, 1:vertical_dimension_for_fast_physics)kind_dyninout
qicloud_ice_specific_humidity_at_lagrangian_surfacecloud ice specific humidity updated by fast physics at Lagrangian surfacekg kg-1real(starting_x_direction_index_domain:ending_x_direction_index_domain, starting_y_direction_index_domain:ending_y_direction_index_domain, 1:vertical_dimension_for_fast_physics)kind_dyninout
qrcloud_rain_specific_humidity_at_lagrangian_surfacecloud rain specific humidity updated by fast physics at Lagrangian surfacekg kg-1real(starting_x_direction_index_domain:ending_x_direction_index_domain, starting_y_direction_index_domain:ending_y_direction_index_domain, 1:vertical_dimension_for_fast_physics)kind_dyninout
qscloud_snow_specific_humidity_at_lagrangian_surfacecloud snow specific humidity updated by fast physics at Lagrangian surfacekg kg-1real(starting_x_direction_index_domain:ending_x_direction_index_domain, starting_y_direction_index_domain:ending_y_direction_index_domain, 1:vertical_dimension_for_fast_physics)kind_dyninout
qgcloud_graupel_specific_humidity_at_lagrangian_surfacecloud graupel specific humidity updated by fast physics at Lagrangian surfacekg kg-1real(starting_x_direction_index_domain:ending_x_direction_index_domain, starting_y_direction_index_domain:ending_y_direction_index_domain, 1:vertical_dimension_for_fast_physics)kind_dyninout
hssurface_geopotential_at_lagrangian_surfacesurface geopotential at Lagrangian surfacem2 s-2real(starting_x_direction_index_domain:ending_x_direction_index_domain, starting_y_direction_index_domain:ending_y_direction_index_domain)kind_dynin
pelnlog_pressure_at_lagrangian_surfacelogarithm of pressure at Lagrangian surfacePareal(starting_x_direction_index:ending_x_direction_index, 1:vertical_dimension_for_fast_physics_plus_one, starting_y_direction_index:ending_y_direction_index)kind_dynin
delzthickness_at_lagrangian_surfacethickness at Lagrangian_surfacemreal(starting_x_direction_index:ending_x_direction_index, starting_y_direction_index:ending_y_direction_index, 1:vertical_dimension_for_thickness_at_Lagrangian_surface)kind_dynin
delppressure_thickness_at_lagrangian_surfacepressure thickness at Lagrangian surfacePareal(starting_x_direction_index_domain:ending_x_direction_index_domain, starting_y_direction_index_domain:ending_y_direction_index_domain, 1:vertical_dimension_for_fast_physics)kind_dynin
ptvirtual_temperature_at_lagrangian_surfacevirtual temperature at Lagrangian surfaceKreal(starting_x_direction_index_domain:ending_x_direction_index_domain, starting_y_direction_index_domain:ending_y_direction_index_domain, 1:vertical_dimension_for_fast_physics)kind_dyninout
pkzfinite_volume_mean_edge_pressure_raised_to_the_power_of_kappafinite-volume mean edge pressure in Pa raised to the power of kappa1real(starting_x_direction_index:ending_x_direction_index, starting_y_direction_index:ending_y_direction_index, 1:vertical_dimension_for_fast_physics)kind_dyninout
q_concloud_condensed_water_specific_humidity_at_lagrangian_surfacecloud condensed water specific humidity updated by fast physics at Lagrangian surfacekg kg-1real(starting_x_direction_index_domain:ending_x_direction_index_domain, starting_y_direction_index_domain:ending_y_direction_index_domain, 1:vertical_dimension_for_condensed_water_at_Lagrangian_surface)kind_dyninout
akapkappa_dry_for_fast_physicsmodified kappa for dry air, fast physicsnonereal()kind_dynin
cappacappa_moist_gas_constant_at_lagrangian_surfacecappa(i,j,k) = rdgas / ( rdgas + cvm(i)/(1.+r_vir*q(i,j,k,sphum)) )nonereal(starting_x_direction_index_domain:ending_x_direction_index_domain, starting_y_direction_index_domain:ending_y_direction_index_domain, 1:vertical_dimension_for_cappa_at_Lagrangian_surface)kind_dyninout
areacell_area_for_fast_physicsarea of the grid cell for fast physicsm2real(starting_x_direction_index_domain:ending_x_direction_index_domain, starting_y_direction_index_domain:ending_y_direction_index_domain)kind_gridin
dtdttendency_of_air_temperature_at_lagrangian_surfaceair temperature tendency due to fast physics at Lagrangian surfaceK s-1real(starting_x_direction_index:ending_x_direction_index, starting_y_direction_index:ending_y_direction_index, 1:vertical_dimension_for_fast_physics)kind_dyninout
out_dtflag_for_tendency_of_air_temperature_at_lagrangian_surfaceflag for calculating tendency of air temperature due to fast physicsflaglogical()in
last_stepflag_for_the_last_step_of_k_split_remappingflag for the last step of k-split remappingflaglogical()in
do_qaflag_for_inline_cloud_fraction_calculationflag for the inline cloud fraction calculationflaglogical()in
qacloud_fraction_at_lagrangian_surfacecloud fraction at Lagrangian surfacenonereal(starting_x_direction_index_domain:ending_x_direction_index_domain, starting_y_direction_index_domain:ending_y_direction_index_domain, 1:vertical_dimension_for_fast_physics)kind_dynout
nthreadsomp_threads_for_fast_physicsnumber of OpenMP threads available for fast physics schemescountinteger()in
errmsgccpp_error_messageerror message for error handling in CCPPnonecharacter()len=*out
errflgccpp_error_codeerror code for error handling in CCPP1integer()out