subroutine gwdps::gwdps_run | ( | integer, intent(in) | im, |
integer, intent(in) | km, | ||
real(kind=kind_phys), dimension(:,:), intent(inout) | a, | ||
real(kind=kind_phys), dimension(:,:), intent(inout) | b, | ||
real(kind=kind_phys), dimension(:,:), intent(inout) | c, | ||
real(kind=kind_phys), dimension(:,:), intent(in) | u1, | ||
real(kind=kind_phys), dimension(:,:), intent(in) | v1, | ||
real(kind=kind_phys), dimension(:,:), intent(in) | t1, | ||
real(kind=kind_phys), dimension(:,:), intent(in) | q1, | ||
integer, dimension(:), intent(in) | kpbl, | ||
real(kind=kind_phys), dimension(:,:), intent(in) | prsi, | ||
real(kind=kind_phys), dimension(:,:), intent(in) | del, | ||
real(kind=kind_phys), dimension(:,:), intent(in) | prsl, | ||
real(kind=kind_phys), dimension(:,:), intent(in) | prslk, | ||
real(kind=kind_phys), dimension(:,:), intent(in) | phii, | ||
real(kind=kind_phys), dimension(:,:), intent(in) | phil, | ||
real(kind=kind_phys), intent(in) | deltim, | ||
integer, intent(in) | kdt, | ||
real(kind=kind_phys), dimension(:), intent(in) | hprime, | ||
real(kind=kind_phys), dimension(:), intent(in) | oc, | ||
real(kind=kind_phys), dimension(:,:), intent(in) | oa4, | ||
real(kind=kind_phys), dimension(:,:), intent(in) | clx4, | ||
real(kind=kind_phys), dimension(:), intent(in) | theta, | ||
real(kind=kind_phys), dimension(:), intent(in) | sigma, | ||
real(kind=kind_phys), dimension(:), intent(in) | gamma, | ||
real(kind=kind_phys), dimension(:), intent(inout) | elvmax, | ||
real(kind=kind_phys), dimension(:), intent(inout) | dusfc, | ||
real(kind=kind_phys), dimension(:), intent(inout) | dvsfc, | ||
real(kind=kind_phys), dimension(:,:), intent(inout), optional | dtaux2d_ms, | ||
real(kind=kind_phys), dimension(:,:), intent(inout), optional | dtauy2d_ms, | ||
real(kind=kind_phys), dimension(:,:), intent(inout), optional | dtaux2d_bl, | ||
real(kind=kind_phys), dimension(:,:), intent(inout), optional | dtauy2d_bl, | ||
real(kind=kind_phys), dimension(:), intent(inout), optional | dusfc_ms, | ||
real(kind=kind_phys), dimension(:), intent(inout), optional | dvsfc_ms, | ||
real(kind=kind_phys), dimension(:), intent(inout), optional | dusfc_bl, | ||
real(kind=kind_phys), dimension(:), intent(inout), optional | dvsfc_bl, | ||
real(kind=kind_phys), intent(in) | g, | ||
real(kind=kind_phys), intent(in) | cp, | ||
real(kind=kind_phys), intent(in) | rd, | ||
real(kind=kind_phys), intent(in) | rv, | ||
integer, intent(in) | imx, | ||
integer, intent(in) | nmtvr, | ||
real(kind=kind_phys), dimension(:), intent(in) | cdmbgwd, | ||
integer, intent(in) | me, | ||
logical, intent(in) | lprnt, | ||
integer, intent(in) | ipr, | ||
real(kind=kind_phys), dimension(:), intent(inout) | rdxzb, | ||
logical, intent(in) | ldiag_ugwp, | ||
character(len=*), intent(out) | errmsg, | ||
integer, intent(out) | errflg ) |
— Subgrid Mountain Blocking Section
UDS=\max(\sqrt{U1^2+V1^2},minwnd)
\frac{U^{2}(h_{d})}{2}=\int_{h_{d}}^{H} N^{2}(z)(H-z)dz
ZLEN=\sqrt{[\frac{h_{d}-z}{z+h'}]}
R=\frac{\cos^{2}\psi+\gamma\sin^{2}\psi}{\gamma\cos^{2}\psi+\sin^{2}\psi}
D_{b}(z)=-C_{d}\max(2-\frac{1}{R},0)\rho\frac{\sigma}{2h'}ZLEN\max(\cos\psi,\gamma\sin\psi)\frac{UDS}{2}
E(OA,F_{r_{0}})=(OA+2)^{\delta}
\delta=C_{E}F_{r_{0}}/F_{r_{c}}
m'(OA,CLX)=C_{m}\triangle x(1+CLX)^{OA+1}
G'(OC,F_{r_{0}})=\frac{F_{r_{0}}^2}{F_{r_{0}}^2+a^{2}}
a^{2}=C_{G}OC^{-1}
\tau_0=E\frac{m'}{\triangle x}\frac{\rho_{0}U_0^3}{N_{0}}G'
R_{scor}=\min \left[\frac{\tau_i}{\tau_{i+1}},1\right]
\tau=\frac{m'}{\triangle x}\rho NUh_d^2
h_{d_i}^2=\frac{\triangle x}{m'}\frac{\tau_{i+1}}{\rho_{i}N_{i}U_{i}}
Ri_{m}=\frac{Ri(1-Fr_{d})}{(1+\sqrt{Ri}\cdot Fr_{d})^{2}}
h_{c}=\frac{U}{N}\left\{2(2+\frac{1}{\sqrt{Ri}})^{1/2}-(2+\frac{1}{\sqrt{Ri}})\right\}
Definition at line 194 of file gwdps.f.
Referenced by cires_ugwp::cires_ugwp_run().