59 |
/ (p(:, :, l) * unpl2k + p(:, :, l+1) * (beta(:, :, l+1) - unpl2k)) |
/ (p(:, :, l) * unpl2k + p(:, :, l+1) * (beta(:, :, l+1) - unpl2k)) |
60 |
ENDDO |
ENDDO |
61 |
|
|
62 |
pk(:, :, 1) = p(:, :, 1) * pks & |
pk(:, :, 1) = ps * pks & |
63 |
/ (p(:, :, 1) * (1. + kappa) & |
/ (ps * (1. + kappa) + 0.5 * (beta(:, :, 2) - unpl2k) * p(:, :, 2)) |
|
+ 0.5 * (beta(:, :, 2) - unpl2k) * p(:, :, 2)) |
|
64 |
DO l = 2, llm |
DO l = 2, llm |
65 |
pk(:, :, l) = beta(:, :, l) * pk(:, :, l-1) |
pk(:, :, l) = beta(:, :, l) * pk(:, :, l - 1) |
66 |
ENDDO |
ENDDO |
67 |
|
|
68 |
if (present(pkf)) then |
if (present(pkf)) then |