CCPP SciDoc v7.0.0  v7.0.0
Common Community Physics Package Developed at DTC
 
Loading...
Searching...
No Matches

◆ zhaocarr_gscond_run()

subroutine, public zhaocarr_gscond::zhaocarr_gscond_run ( integer, intent(in) im,
integer, intent(in) km,
real(kind=kind_phys), intent(in) dt,
real(kind=kind_phys), intent(in) dtf,
real(kind=kind_phys), dimension(:,:), intent(in) prsl,
real(kind=kind_phys), dimension(:), intent(in) ps,
real(kind=kind_phys), dimension(:,:), intent(inout) q,
real(kind=kind_phys), dimension(:,:), intent(in) clw1,
real(kind=kind_phys), dimension(:,:), intent(in) clw2,
real(kind=kind_phys), dimension(:,:), intent(out) cwm,
real(kind=kind_phys), dimension(:,:), intent(inout) t,
real(kind=kind_phys), dimension(:,:), intent(inout), optional tp,
real(kind=kind_phys), dimension(:,:), intent(inout), optional qp,
real(kind=kind_phys), dimension(:), intent(inout), optional psp,
real(kind=kind_phys), intent(in) psat,
real(kind=kind_phys), intent(in) hvap,
real(kind=kind_phys), intent(in) grav,
real(kind=kind_phys), intent(in) hfus,
real(kind=kind_phys), intent(in) ttp,
real(kind=kind_phys), intent(in) rd,
real(kind=kind_phys), intent(in) cp,
real(kind=kind_phys), intent(in) eps,
real(kind=kind_phys), intent(in) epsm1,
real(kind=kind_phys), intent(in) rv,
real(kind=kind_phys), dimension(:,:), intent(inout), optional tp1,
real(kind=kind_phys), dimension(:,:), intent(inout), optional qp1,
real(kind=kind_phys), dimension(:), intent(inout), optional psp1,
real(kind=kind_phys), dimension(:,:), intent(in) u,
logical, intent(in) lprnt,
integer, intent(in) ipr,
character(len=*), intent(out) errmsg,
integer, intent(out) errflg )

Argument Table

general_gscond GFS gscond Scheme General Algorithm

  1. Calculate ice-water identification number \(IW\) in order to make a distinction between cloud water and cloud ice (table2 of Zhao and Carr (1997) [223]).
  2. Calculate the changes in \(t\), \(q\) and \(p\) due to all the processes except microphysics.
  3. Calculate cloud evaporation rate ( \(E_c\), eq. 19 of Zhao and Carr (1997)[223]).
  4. Calculate cloud condensation rate ( \(C_g\), eq.8 of Zhao and Carr (1997)[223]).
  5. Update \(t\), \(q\), \(cwm\) due to cloud evaporation and condensation processes.

GFS gscond Scheme Detailed Algorithm

  1. Begining of grid-scale condensation/evaporation loop (start of k-loop, i-loop)
  1. Compute ice-water identification number IW.
    The distinction between cloud water and cloud ice is made by the cloud identification number IW, which is zero for cloud water and unity for cloud ice (Table 2 in Zhao and Carr (1997) [223]):
    • All clouds are defined to consist of liquid water below the freezing level ( \(T\geq 0^oC\)) and of ice particles above the \(T=-15^oC\) level.
    • In the temperature region between \(-15^oC\) and \(0^oC\), clouds may be composed of liquid water or ice. If there are cloud ice particles above this point at the previous or current time step, or if the cloud at this point at the previous time step consists of ice particles, then the cloud substance at this point is considered to be ice particles because of the cloud seeding effect and the memory of its content. Otherwise, all clouds in this region are considered to contain supercooled cloud water.
  1. Condensation and evaporation of cloud
  • Compute the changes in t, q and p ( \(A_{t}\), \(A_{q}\) and \(A_{p}\)) caused by all the processes except grid-scale condensation and evaporation.

    \[ A_{t}=(t-tp)/dt \]

    \[ A_{q}=(q-qp)/dt \]

    \[ A_{p}=(prsl-\frac{prsl}{ps} \times psp)/dt \]

  • Calculate the saturation specific humidity \(q_{s}\) and the relative humidity \(f\) using IW.
  • According to Sundqvist et al. (1989) [198], estimate cloud fraction \(b\) at a grid point from relative humidity \(f\) using the equation

    \[ b=1-\left ( \frac{f_{s}-f}{f_{s}-u} \right )^{1/2} \]

    for \(f>u\)
  1. End of the condensation/evaporation loop (end of i-loop,k-loop).
  1. Store \(t\), \(q\), \(ps\) for next time step.

Definition at line 73 of file zhaocarr_gscond.f.