This subroutine computes grid-scale condensation and evaporation of cloud condensate. More...
This subroutine computes grid-scale condensation and evaporation of cloud condensate.
There are two sources of condensation, one from large-scale processes and the other from convective processes. Both of them produce either cloud water or cloud ice, depending on the cloud substance at and above the grid point at current and previous time steps, and on the temperature. Evaporation of cloud is allowed at points where the relative humidity is lower than the critical value required for condensation.
subroutine | gscond (im, ix, km, dt, dtf, prsl, ps, q, cwm, t , tp, qp, psp, tp1, qp1, psp1, u, lprnt, ipr) |
subroutine gscond | ( | integer | im, |
integer | ix, | ||
integer | km, | ||
real (kind=kind_phys) | dt, | ||
real (kind=kind_phys) | dtf, | ||
real (kind=kind_phys), dimension(ix,km) | prsl, | ||
real (kind=kind_phys), dimension(im) | ps, | ||
real (kind=kind_phys), dimension(ix,km) | q, | ||
real (kind=kind_phys), dimension(ix,km) | cwm, | ||
real (kind=kind_phys), dimension(ix,km) | t, | ||
real (kind=kind_phys), dimension(ix,km) | tp, | ||
real (kind=kind_phys), dimension(ix,km) | qp, | ||
real (kind=kind_phys), dimension(im) | psp, | ||
real (kind=kind_phys), dimension(ix,km) | tp1, | ||
real (kind=kind_phys), dimension(ix,km) | qp1, | ||
real (kind=kind_phys), dimension(im) | psp1, | ||
real (kind=kind_phys), dimension(im,km) | u, | ||
logical | lprnt, | ||
integer | ipr | ||
) |
[in] | ix | horizontal dimension |
[in] | im | horizontal number of used pts |
[in] | km | vertical layer dimension |
[in] | dt | physics time step in seconds |
[in] | dtf | dynamics time step in seconds |
[in] | prsl | pressure values for model layers |
[in] | ps | surface pressure (Pa) |
[in,out] | q | model layer specific humidity (gm/gm) |
[in,out] | cwm | model layer cloud condensate |
[in,out] | t | model layer mean temperature (K) |
[in,out] | tp | model layer mean temperature (K) saved for restart |
[in,out] | qp | model layer specific humidity (gm/gm) saved for restart |
[in,out] | psp | surface pressure (Pa) saved for restart |
[in,out] | tp1 | updated model layer mean temperature (K) saved for restart |
[in,out] | qp1 | updated model layer specific humidity (gm/gm) saved for restart |
[in,out] | psp1 | updated surface pressure (Pa) saved for restart |
[in] | u | the critical value of relative humidity for large-scale condensation |
[in] | lprnt | logical print flag |
[in] | ipr | check print point for debugging |
\[ A_{t}=(t-tp)/dt \]
\[ A_{q}=(q-qp)/dt \]
\[ A_{p}=(prsl-\frac{prsl}{ps} \times psp)/dt \]
\[ b=1-\left ( \frac{f_{s}-f}{f_{s}-u} \right )^{1/2} \]
for \(f>u\); and \(b=0\) for \(f<u\). where \(f_{s}=1.0\) is the relative humidity in a cloud region and \(u\) ,which is an input parameter accounts for the effects of subgrid-scale variations in moisture on large-scale condensation. Since both temperature and moisture may vary at scales smaller than the model grid scale, it is possible for condensation to occur before the grid-average relative humidity reaches 100%. Therefore \(u\) needs to be less than 1.0 to account for the subgrid-scale variation of temperature and moisture fields and allow subgrid-scale condensation.\[ q_{0}=uq_{s} \]
\[ E_{c}=\frac{q_{0}-q}{dt} \]
\[ E_{c}=\frac{q_{s}}{dt}(u-f) \]
where \(dt\) is the time step for precipitation calculation in the model. It is a simplified version of a higher-order cloud evaporation algorithm (Rutledge and Hobbs 1983 [rutledge_and_hobbs_1983]). In the case where all clouds will evaporate before \(u\) is reached, the following equation is used:\[ E_{c}=\frac{cwm}{dt} \]
\[ C_{g}=\frac{M-q_{s}f_{t}}{1+(f\epsilon L^{2}q_{s}/RC_{p}T^{2})}+E_{c} \]
where\[ M=A_{q}-\frac{f\epsilon Lq_{s}}{RT^{2}}A_{t}+\frac{fq_{s}}{p}A_{p} \]
To close the system, an equation for the relative humidity tendency \(f_{t}\) was derived by Sundqvist et al. (1989) [51] using the hypothesis that the quantity \(M+E_{c}\) is divided into one part, \(bM\),which condenses in the already cloudy portion of a grid square, and another part, \((1-b)M+E_{c}\),which is used to increase the relative humidity of the cloud-free portion and the cloudiness in the square. The equation is written as\[ f_{t}=\frac{2(1-b)(f_{s}-u)[(1-b)M+E_{c}]}{2q_{s}(1-b)(f_{s}-u)+cwm/b} \]
\[ cwm=cwm+(C_{g}-E_{c})\times dt \]
\[ q=q-(C_{g}-E_{c})\times dt \]
\[ t=t+\frac{L}{C_{p}}(C_{g}-E_{c})\times dt \]
Definition at line 87 of file gscond.f.
References physcons::con_cp, physcons::con_g, physcons::con_hfus, physcons::con_hvap, physcons::con_psat, physcons::con_rd, physcons::con_rv, and physcons::con_ttp.
Referenced by gbphys().