Layer that invokes individual GFS physics routines. More...
Layer that invokes individual GFS physics routines.
Functions/Subroutines | |
subroutine | gbphys |
Parameter descriptions include intent, name, description, and size. More... | |
subroutine gbphys | ( | ) |
Parameter descriptions include intent, name, description, and size.
gbphys is the driver subroutine to invoke GFS physics (except radiation but radiative heating is applied here) at physics time steps
[in] | ix(1),im(1) | - integer, horiz dimension and num of used pts |
[in] | levs(1) | - integer, vertical layer dimension |
[in] | lsoil(1) | - integer, number of soil layers |
[in] | lsm(1) | - integer, flag for land surface model to use =0 for osu lsm; =1 for noah lsm |
[in] | ntrac(1) | - integer, number of tracers |
[in] | ncld(1) | - integer, number of cloud species |
[in] | ntoz(1) | - integer, ozone location in the tracer array |
[in] | ntcw(1) | - integer, cloud condensate location in the tracer array |
[in] | ntke(1) | - integer, tke location in the tracer array |
[in] | nmtvr(1) | - integer, number of topographic variables such as variance etc used in the GWD parameterization |
[in] | nrcm(1) | - integer, second dimension for the random number array rann |
[in] | ko3(1) | - integer, number of layers for ozone data |
[in] | lonr(1),latr(1) | - integer, number of lon/lat points |
[in] | jcap(1) | - integer, number of spectral wave trancation used only by sascnv shalcnv |
[in] | num_p3d(1) | - integer, number of 3D arrays needed for microphysics |
[in] | num_p2d(1) | - integer, number of 2D arrays needed for microphysics |
[in] | npdf3d(1) | - integer, number of 3d arrays associated with pdf based clouds/microphysics |
[in] | ncnvcld3d(1) | - integer, number of 3d arrays associated with convective cloudiness enhancement |
[in] | kdt(1) | - integer, number of the current time step |
[in] | lat(1) | - integer, latitude index - used for debug prints |
[in] | me(1) | - integer, pe number - used for debug prints |
[in] | pl_coeff(1) | - integer, number coefficients in ozone forcing |
[in] | nlons(im) | - integer, number of total grid points in a latitude circle through a point |
[in] | ncw(2) | - integer, range of droplet number concentrations for Ferrier microphysics |
[in] | flgmin(2) | - real, range of minimum large ice fraction for Ferrier microphys |
[in] | crtrh(3) | - real, critical relative humidity at the surface, PBL top and at the top of the atmosphere |
[in] | cdmbgwd(2) | - real, multiplication factors for cdmb and gwd |
[in] | ccwf(2) | - real, multiplication factor for critical cloud workfunction for RAS |
[in] | dlqf(2) | - real, factor for cloud condensate detrainment from cloud edges (RAS) |
[in] | ctei_rm(2) | - real, critical cloud top entrainment instability criteria (used if mstrat=.true.) |
[in] | clstp(1) | - real, index used by cnvc90 (for convective clouds) legacy stuff - does not affect forecast |
[in] | cgwf(2) | - real, multiplication factor for convective GWD |
[in] | prslrd0(1) | - real, pressure level (Pa) from which Rayleigh Damping is applied |
[in] | ral_ts(1) | - real time scale for Rayleigh damping in days |
[in] | dtp(1) | - real, physics time step in seconds |
[in] | dtf(1) | - real, dynamics time step in seconds |
[in] | fhour(1) | - real, forecast hour |
[in] | solhr(1) | - real, fcst hour at the end of prev time step |
[in] | slag(1) | - real, equation of time ( radian ) |
[in] | sdec(1),cdec(1) | - real, sin and cos of the solar declination angle |
[in] | sinlat(im) | - real, sin of latitude |
[in] | coslat(im) | - real, cos of latitude |
[in] | pgr(im) | - real, surface pressure (Pa) |
[in] | ugrs(ix,levs),vgrs(ix,levs) | - real, u/v component of layer wind |
[in] | tgrs(ix,levs) | - real, layer mean temperature ( k ) |
[in] | qgrs(ix,levs,ntrac) | - real, layer mean tracer concentration |
[in] | vvel(ix,levs) | - real, layer mean vertical velocity (Pa/s) |
[in] | prsi(ix,levs+1) | - real, pressure at layer interfaces |
[in] | prsl(ix,levs) | - real, mean layer presure |
[in] | prsik(ix,levs+1) | - real, Exner function at layer interface |
[in] | prslk(ix,levs) | - real, Exner function at layer |
[in] | phii(ix,levs+1) | - real, interface geopotential (m^2/s^2) |
[in] | phil(ix,levs) | - real, layer geopotential (m^2/s^2) |
[in] | rann(ix,nrcm) | - real, random number array (0-1) |
[in] | prdout(ix,ko3,pl_coeff) | - real, ozone forcing data |
[in] | poz(ko3) | - real, ozone forcing data level pressure (ln(Pa)) |
[in] | dpshc(im) | - real, maximum pressure depth for shallow convection |
[in] | hprime(ix,nmtvr) | - real, orographic std dev |
[in] | xlon(im),xlat(im) | - real, longitude and latitude ( radian ) |
[in] | slope(im) | - real, sfc slope type for lsm |
[in] | shdmin(im1) | - real, min fractional coverage of green veg |
[in] | shdmax(im1) | - real, max fractnl cover of green veg (not used) |
[in] | snoalb(im1) | - real, max snow albedo over land (for deep snow) |
[in] | tg3(im) | - real, deep soil temperature |
[in] | slmsk(im | - real, sea/land/ice mask (=0/1/2) |
[in] | vfrac(im) | - real, vegetation fraction |
[in] | vtype(im) | - real, vegetation type |
[in] | stype(im) | - real, soil type |
[in] | uustar(im) | - real, boundary layer parameter |
[in] | oro(im) | - real, orography |
[in] | oro_uf(im) | - real, unfiltered orography |
[in] | coszen(im) | - real, avg cosz over daytime sw radiation interval |
[in] | sfcdsw(im) | - real, total sky sfc downward sw flux ( w/m**2 ) |
[in] | sfcnsw(im) | - real, total sky sfc netsw flx into ground(w/m**2) |
[in] | sfcdlw(im) | - real, total sky sfc downward lw flux ( w/m**2 ) |
[in] | tsflw(im) | - real, sfc air (layer 1) temp over lw interval (k) |
[in] | sfcemis(im) | - real, sfc lw emissivity ( fraction ) |
[in] | sfalb(im) | - real, mean sfc diffused sw albedo |
[in] | sfcnirbmu(im) | - real, sfc nir-beam sw upward flux (w/m2) |
[in] | sfcnirdfu(im) | - real, sfc nir-diff sw upward flux (w/m2) |
[in] | sfcvisbmu(im) | - real, sfc uv+vis-beam sw upward flux (w/m2) |
[in] | sfcvisdfu(im) | - real, sfc uv+vis-diff sw upward flux (w/m2) |
[in] | sfcnirbmd(im) | - real, sfc nir-beam sw downward flux (w/m2) |
[in] | sfcnirdfd(im) | - real, sfc nir-diff sw downward flux (w/m2) |
[in] | sfcvisbmd(im) | - real, sfc uv+vis-beam sw downward flux (w/m2) |
[in] | sfcvisdfd(im) | - real, sfc uv+vis-diff sw downward flux (w/m2) |
[in] | slimskin_cpl(im) | - real, |
[in] | ulwsfcin_cpl(im) | - real, |
[in] | dusfcin_cpl(im) | - real, |
[in] | dvsfcin_cpl(im) | - real, |
[in] | dtsfcin_cpl(im) | - real, |
[in] | dqsfcin_cpl(im) | - real, |
[in] | swh(ix,levs) | - real, total sky sw heating rates ( k/s ) |
[in] | swhc(ix,levs) | - real, clear sky sw heating rates ( k/s ) |
[in] | hlw(ix,levs) | - real, total sky lw heating rates ( k/s ) |
[in] | hlwc(ix,levs) | - real, clear sky lw heating rates ( k/s ) |
[in] | hlwd(ix,levs) | - real, idea sky lw heating rates ( k/s ) |
[in] | ras(1) | - logical, flag for ras convection scheme |
[in] | cscnv(1) | - logical, flag for Chikira-Sugiyama convection |
[in] | nctp(1) | - integer, number of cloud types in CS scheme |
[in] | do_shoc(1) | - logical, flag for SHOC |
[in] | shocaftcnv(1) | - logical, flag for SHOC |
[in] | ntot3d(1) | - integer, number of total 3d fields for phy_f3d |
[in] | ntot2d(1) | - integer, number of total 2d fields for phy_f2d |
[in] | pre_rad(1) | - logical, flag for testing purpose |
[in] | ldiag3d(1) | - logical, flag for 3d diagnostic fields |
[in] | lgocart(1) | - logical, flag for 3d diagnostic fields for gocart |
[in] | lssav(1) | - logical, flag controls data store and output |
[in] | lssav_cpl(1) | - logical, flag for save data for A/O/I coupling |
[in] | flipv(1) | - logical, flag for vertical direction flip (ras) |
[in] | xkzm_m(1) | - real, background vertical diffusion for momentum |
[in] | xkzm_h(1) | - real, background vertical diffusion for heat, q |
[in] | xkzm_s(1) | - real, sigma threshold for background mom. diffusn |
[in] | psautco(2) | - real, auto conversion coeff from ice to snow |
[in] | prautco(2) | - real, auto conversion coeff from cloud to rain |
[in] | evpco(1) | - real, coeff for evaporation of largescale rain |
[in] | wminco(1) | - real, water and ice minimum threshold for Zhao |
[in] | pdfcld(1) | - logical, flag for pdfcld |
[in] | shcnvcw(1) | - logical, flag for shallow convective cloud |
[in] | sup(1) | - real, supsaturation for ho. nucleation of ice |
[in] | redrag(1) | - logical, flag for reduced drag coeff. over sea |
[in] | hybedmf(1) | - logical, flag for hybrid edmf pbl scheme |
[in] | dspheat(1) | - logical, flag for tke dissipative heating |
[in] | old_monin(1) | - logical, flag for diff monin schemes |
[in] | cnvgwd(1) | - logical, flag for conv gravity wave drag |
[in] | shal_cnv(1) | - logical, flag for calling shallow convection |
[in] | imfshalcnv(1) | - integer, flag for mass-flux shallow conv scheme 1: July 2010 version of mass-flux shallow conv scheme current operational version as of 2016 2: scale- & aerosol-aware mass-flux shallow conv scheme (2017) 0: modified Tiedtke's eddy-diffusion shallow conv scheme -1: no shallow convection used |
[in] | imfdeepcnv(1) | - integer, flag for mass-flux deep conv scheme 1: July 2010 version of SAS conv scheme current operational version as of 2016 2: scale- & aerosol-aware mass-flux deep conv scheme (2017) 0: old SAS Convection scheme before July 2010 |
[in] | cal_pre(1) | - logical, flag controls precip type algorithm |
[in] | mom4ice(1) | - logical, flag controls mom4 sea-ice |
[in] | mstrat(1) | - logical, flag for moorthi approach for stratus |
[in] | trans_trac(1) | - logical, flag for convective transport of tracers |
[in] | nstf_name(1) | - integer array, NSST related flag parameters nstf_name(1) : 0 = NSSTM off 1 = NSSTM on but uncoupled 2 = NSSTM on and coupled nstf_name(2) : 1 = NSSTM spin up on 0 = NSSTM spin up off nstf_name(3) : 1 = NSST analysis on 0 = NSSTM analysis off nstf_name(4) : zsea1 in mm nstf_name(5) : zsea2 in mm |
[in] | moist_adj(1) | - logical, flag for moist convective adjustment |
[in] | thermodyn_id(1) | - integer, valid for GFS only for get_prs/phi |
[in] | sfcpress_id(1) | - integer, valid for GFS only for get_prs/phi |
[in] | gen_coord_hybrid(1) | - logical for pk=ak+bk*ps+ck*theta (Henry) |
[in] | levr(1) | - integer, the number of layers GFS Radiative heating calculted at 1 |
[in] | adjtrc(ntrac) | - real, dynamics adjustments to tracers |
[in] | nnp(1) | - integer, physics substep number |
[in,out] | hice(im) | - real, sea-ice thickness |
[in,out] | fice(im) | - real, sea-ice concentration |
[in,out] | tisfc(im) | - real, sea-ice temperature |
[in,out] | tsea(im) | - real, ground surface temperature ( k ) |
[in,out] | tprcp(im) | - real, total precipitation the following three variables do not affect the forecast |
[in,out] | cv(im) | - real, convective clouds amountt |
[in,out] | cvb(im) | - real, convective clouds base pressure (kPa) |
[in,out] | cvt(im) | - real, convective clouds top pressure (kPa) |
[in,out] | srflag(im) | - real, snow/rain flag for precipitation |
[in,out] | snwdph(im) | - real, actual snow depth (mm) over land/sea ice |
[in,out] | weasd(im) | - real, water equiv of accumulated snow depth (kg/m**2) over land and sea ice |
[in,out] | sncovr(im) | - real, snow cover over land |
[in,out] | zorl(im) | - real, surface roughness |
[in,out] | canopy(im) | - real, canopy water |
[in,out] | ffmm(im) | - real, fm parameter from PBL scheme |
[in,out] | ffhh(im) | - real, fh parameter from PBL scheme |
[in,out] | f10m(im) | - real, fm at 10m |
[in,out] | srunoff(im) | - real, surface water runoff (from lsm) |
[in,out] | evbsa(im) | - real, noah lsm diagnostics |
[in,out] | evcwa(im) | - real, noah lsm diagnostics |
[in,out] | snohfa(im) | - real, noah lsm diagnostics |
[in,out] | transa(im) | - real, noah lsm diagnostics |
[in,out] | sbsnoa(im) | - real, noah lsm diagnostics |
[in,out] | snowca(im) | - real, noah lsm diagnostics |
[in,out] | soilm(im) | - real, soil moisture |
[in,out] | tmpmin(im) | - real, min temperature at 2m height (k) |
[in,out] | tmpmax(im) | - real, max temperature at 2m height (k) |
[in,out] | dusfc(im) | - real, u component of surface stress |
[in,out] | dvsfc(im) | - real, v component of surface stress |
[in,out] | dtsfc(im) | - real, sensible heat flux (w/m2) |
[in,out] | dqsfc(im) | - real, latent heat flux (w/m2) |
[in,out] | totprcp(im) | - real, accumulated total precipitation (kg/m2) |
[in,out] | gflux(im) | - real, groud conductive heat flux |
[in,out] | dlwsfc(im) | - real, time accumulated sfc dn lw flux ( w/m**2 ) |
[in,out] | ulwsfc(im) | - real, time accumulated sfc up lw flux ( w/m**2 ) |
[in,out] | suntim(im) | - real, sunshine duration time (s) |
[in,out] | runoff(im) | - real, total water runoff |
[in,out] | ep(im) | - real, potential evaporation |
[in,out] | cldwrk(im) | - real, cloud workfunction (valid only with sas) |
[in,out] | dugwd(im) | - real, vertically integrated u change by OGWD |
[in,out] | dvgwd(im) | - real, vertically integrated v change by OGWD |
[in,out] | psmean(im) | - real, surface pressure (kPa) |
[in,out] | cnvprcp(im) | - real, accumulated convective precipitation (kg/m2) |
[in,out] | spfhmin(im) | - real, minimum specific humidity |
[in,out] | spfhmax(im) | - real, maximum specific humidity |
[in,out] | rain(im) | - real, total rain at this time step |
[in,out] | rainc(im) | - real, convective rain at this time step |
[in,out] | dt3dt(ix,les,6) | - real, temperature change due to physics |
[in,out] | dq3dt(ix,levs,5+pl_coeff) | - real, moisture change due to physics |
[in,out] | du3dt(ix,levs,4) | - real, u momentum change due to physics |
[in,out] | dv3dt(ix,levs,4) | - real, v momentum change due to physics |
[in,out] | dqdt_v(ix,levs) | - real, total moisture tendency (kg/kg/s) |
[in,out] | cnvqc_v(ix,levs) | - real, total convective conensate (kg/kg) |
[in,out] | acv(im) | - real, array containing accumulated convective clouds |
[in,out] | acvb(im),acvt(im) | - real, arrays used by cnvc90 |
[in,out] | slc(ix,soil) | - real, liquid soil moisture |
[in,out] | smc(ix,soil) | - real, total soil moisture |
[in,out] | stc(ix,soil) | - real, soil temperature |
[in,out] | upd_mf(ix,levs) | - real, convective updraft mass flux |
[in,out] | dwn_mf(ix,levs) | - real, convective downdraft mass flux |
[in,out] | det_mf(ix,levs) | - real, convective detrainment mass flux |
[in,out] | dkh(ix,levs) | - real, vertical diffusion coefficient (gocart), not used |
[in,out] | rnp(ix,levs) | - real, n cloud precipitation rate (gocart), not used |
[in,out] | phy_f3d(ix,levs,num_p3d) | - real, 3d arrays saved for restart |
[in,out] | phy_f2d(ix,num_p2d) | - real, 2d arrays save for restart |
[in,out] | dusfc_cpl(im) | - real, sfc u-momentum flux for A/O/I coupling |
[in,out] | dvsfc_cpl(im) | - real, sfc v-momentum flux for A/O/I coupling |
[in,out] | dtsfc_cpl(im) | - real, sfc sensible heat flux for A/O/I coupling |
[in,out] | dqsfc_cpl(im) | - real, sfc latent heat flux for A/O/I coupling |
[in,out] | dlwsfc_cpl(im) | - real, sfc dnwd lw flux (w/m**2) for A/O/I coupling |
[in,out] | dswsfc_cpl(im) | - real, sfc dnwd sw flux (w/m**2) for A/O/I coupling |
[in,out] | dnirbm_cpl(im) | - real, sfc nir beam dnwd sw rad flux (w/m**2) |
[in,out] | dnirdf_cpl(im) | - real, sfc nir diff dnwd sw rad flux (w/m**2) |
[in,out] | dvisbm_cpl(im) | - real, sfc uv+vis beam dnwd sw rad flux (w/m**2) |
[in,out] | dvisdf_cpl(im) | - real, sfc uv+vis diff dnwd sw rad flux (w/m**2) |
[in,out] | nlwsfc_cpl(im) | - real, net dnwd lw flux (w/m**2) for A/O/I coupling |
[in,out] | nswsfc_cpl(im) | - real, net dnwd sw flux (w/m**2) for A/O/I coupling |
[in,out] | nnirbm_cpl(im) | - real, net nir beam dnwd sw rad flux (w/m**2) |
[in,out] | nnirdf_cpl(im) | - real, net nir diff dnwd sw rad flux (w/m**2) |
[in,out] | nvisbm_cpl(im) | - real, net uv+vis beam dnwd sw rad flux (w/m**2) |
[in,out] | nvisdf_cpl(im) | - real, net uv+vis diff dnwd sw rad flux (w/m**2) |
[in,out] | rain_cpl(im) | - real, total precipitation rain for A/O/I coupling |
[in,out] | snow_cpl(im) | - real, total precipitation snow for A/O/I coupling |
[in,out] | xt(im) | - real, heat content in DTL |
[in,out] | xs(im) | - real, salinity content in DTL |
[in,out] | xu(im) | - real, u-current content in DTL |
[in,out] | xv(im) | - real, v-current content in DTL |
[in,out] | xz(im) | - real, DTL thickness |
[in,out] | zm(im) | - real, MXL thickness |
[in,out] | xtts(im) | - real, d(xt)/d(ts) |
[in,out] | xzts(im) | - real, d(xz)/d(ts) |
[in,out] | d_conv(im) | - real, thickness of Free Convection Layer (FCL) |
[in,out] | ifd(im) | - real, index to start DTM run or not |
[in,out] | dt_cool(im) | - real, Sub-layer cooling amount |
[in,out] | Qrain(im) | - real, sensible heat flux due to rainfall (watts) |
[in,out] | phy_fctd(ix,nctp) | - real, cloud base mass flux for CScnv |
[out] | gt0(ix,levs) | - real, updated temperature |
[out] | gq0(ix,levs,ntrac) | - real, updated tracers |
[out] | gu0(ix,levs) | - real, updated zonal wind |
[out] | gv0(ix,levs) | - real, update meridional wind |
[out] | t2m(im),q2m(im) | - real, 2 meter temperature and humidity |
[out] | u10m(im),v10m(im) | - real, 10 meater u/v wind speed |
[out] | zlvl(im) | - real, layer 1 height (m) |
[out] | psurf(im) | - real, surface pressure (Pa) |
[out] | hpbl(im) | - real, pbl height (m) |
[out] | pwat(im) | - real, precipitable water |
[out] | t1(im) | - real, layer 1 temperature (K) |
[out] | q1(im) | - real, layer 1 specific humidity (kg/kg) |
[out] | u1(im) | - real, layer 1 zonal wind (m/s) |
[out] | v1(im) | - real, layer 1 merdional wind (m/s) |
[out] | chh(im) | - real, thermal exchange coefficient |
[out] | cmm(im) | - real, momentum exchange coefficient |
[out] | dlwsfci(im) | - real, instantaneous sfc dnwd lw flux ( w/m**2 ) |
[out] | ulwsfci(im) | - real, instantaneous sfc upwd lw flux ( w/m**2 ) |
[out] | dswsfci(im) | - real, instantaneous sfc dnwd sw flux ( w/m**2 ) |
[out] | uswsfci(im) | - real, instantaneous sfc upwd sw flux ( w/m**2 ) |
[out] | dusfci(im) | - real, instantaneous u component of surface stress |
[out] | dvsfci(im) | - real, instantaneous v component of surface stress |
[out] | dtsfci(im) | - real, instantaneous sfc sensible heat flux |
[out] | dqsfci(im) | - real, instantaneous sfc latent heat flux |
[out] | gfluxi(im) | - real, instantaneous sfc ground heat flux |
[out] | epi(im) | - real, instantaneous sfc potential evaporation |
[out] | smcwlt2(im) | - real, wilting point (volumetric) |
[out] | smcref2(im) | - real, soil moisture threshold (volumetric) |
[out] | dusfci_cpl(im) | - real, sfc u-momentum flux at time step AOI cpl |
[out] | dvsfci_cpl(im) | - real, sfc v-momentum flux at time step AOI cpl |
[out] | dtsfci_cpl(im) | - real, sfc sensib heat flux at time step AOI cpl |
[out] | dqsfci_cpl(im) | - real, sfc latent heat flux at time step AOI cpl |
[out] | dlwsfci_cpl(im) | - real, sfc dnwd lw flux at time step AOI cpl |
[out] | dswsfci_cpl(im) | - real, sfc dnwd sw flux at time step AOI cpl |
[out] | dnirbmi_cpl(im) | - real, sfc nir beam dnwd sw flx rad at time step |
[out] | dnirdfi_cpl(im) | - real, sfc nir diff dnwd sw flx rad at time step |
[out] | dvisbmi_cpl(im) | - real, sfc uv+vis beam dnwd sw flx at time step |
[out] | dvisdfi_cpl(im) | - real, sfc uv+vis diff dnwd sw flx at time step |
[out] | nlwsfci_cpl(im) | - real, net sfc dnwd lw flux at time step AOI cpl |
[out] | nswsfci_cpl(im) | - real, net sfc dnwd sw flux at time step AOI cpl |
[out] | nnirbmi_cpl(im) | - real, net nir beam dnwd sw flx rad at time step |
[out] | nnirdfi_cpl(im) | - real, net nir diff dnwd sw flx rad at time step |
[out] | nvisbmi_cpl(im) | - real, net uv+vis beam dnwd sw flx at time step |
[out] | nvisdfi_cpl(im) | - real, net uv+vis diff dnwd sw flx at time step |
[out] | ocalnirbm_cpl(im) | - real, ocean alb nir beam (no ice) at time step |
[out] | ocalnirdf_cpl(im) | - real, ocean alb nir diff (no ice) at time step |
[out] | ocalvisbm_cpl(im) | - real, ocean alb vis beam (no ice) at time step |
[out] | ocalvisdf_cpl(im) | - real, ocean alb vis diff (no ice) at time step |
[out] | t2mi_cpl(im) | - real, T2m at time step AOI cpl |
[out] | q2mi_cpl(im) | - real, Q2m at time step AOI cpl |
[out] | u10mi_cpl(im) | - real, U10m at time step AOI cpl |
[out] | v10mi_cpl(im) | - real, V10m at time step AOI cpl |
[out] | tseai_cpl(im) | - real, sfc temp at time step AOI cpl |
[out] | psurfi_cpl(im) | - real, sfc pressure at time step AOI cpl |
[out] | tref(im) | - real, Reference Temperature |
[out] | z_c(im) | - real, Sub-layer cooling thickness |
[out] | c_0(im) | - real, coefficient1 to calculate d(Tz)/d(Ts) |
[out] | c_d(im) | - real, coefficient2 to calculate d(Tz)/d(Ts) |
[out] | w_0(im) | - real, coefficient3 to calculate d(Tz)/d(Ts) |
[out] | w_d(im) | - real, coefficient4 to calculate d(Tz)/d(Ts) |
[out] | rqtk(im) | - real, mass change due to moisture variation |
[out] | dtdtr(ix,levs) | - real, temperature change due to radiative heating per time step (K) |
Definition at line 529 of file gbphys.f.
References physcons::con_cp, physcons::con_g, physcons::con_hfus, physcons::con_hvap, physcons::con_pi, physcons::con_rd, physcons::con_rerth, physcons::con_rv, gscond(), gwdc(), gwdps(), moninedmf(), precpd(), sascnvn(), and shalcnv().