4 |
|
|
5 |
contains |
contains |
6 |
|
|
7 |
SUBROUTINE swclr(knu, flag_aer, tauae, pizae, cgae, palbp, pdsig, prayl, & |
SUBROUTINE swclr(knu, flag_aer, palbp, pdsig, prayl, psec, pcgaz, ppizaz, & |
8 |
psec, pcgaz, ppizaz, pray1, pray2, prefz, prj, prk, prmu0, ptauaz, & |
pray1, pray2, prefz, prj, prk, prmu0, ptauaz, ptra1, ptra2) |
|
ptra1, ptra2) |
|
9 |
|
|
10 |
USE raddim, only: kdlon, kflev |
USE raddim, only: kdlon, kflev |
11 |
USE radepsi, only: repsct, zepsec |
USE radepsi, only: repsct, zepsec |
35 |
|
|
36 |
INTEGER knu |
INTEGER knu |
37 |
! -OB |
! -OB |
38 |
DOUBLE PRECISION flag_aer |
logical, intent(in):: flag_aer |
|
DOUBLE PRECISION tauae(kdlon, kflev, 2) |
|
|
DOUBLE PRECISION pizae(kdlon, kflev, 2) |
|
|
DOUBLE PRECISION cgae(kdlon, kflev, 2) |
|
39 |
DOUBLE PRECISION palbp(kdlon, 2) |
DOUBLE PRECISION palbp(kdlon, 2) |
40 |
DOUBLE PRECISION pdsig(kdlon, kflev) |
DOUBLE PRECISION pdsig(kdlon, kflev) |
41 |
DOUBLE PRECISION prayl(kdlon) |
DOUBLE PRECISION prayl(kdlon) |
84 |
|
|
85 |
DO jk = 1, kflev |
DO jk = 1, kflev |
86 |
DO jl = 1, kdlon |
DO jl = 1, kdlon |
87 |
ptauaz(jl, jk) = flag_aer*tauae(jl, jk, knu) |
ptauaz(jl, jk) = 0d0 |
88 |
ppizaz(jl, jk) = flag_aer*pizae(jl, jk, knu) |
ppizaz(jl, jk) = 0d0 |
89 |
pcgaz(jl, jk) = flag_aer*cgae(jl, jk, knu) |
pcgaz(jl, jk) = 0d0 |
90 |
END DO |
END DO |
91 |
|
|
92 |
IF (flag_aer>0) THEN |
IF (flag_aer) THEN |
93 |
! -OB |
! -OB |
94 |
DO jl = 1, kdlon |
DO jl = 1, kdlon |
|
! PCGAZ(JL,JK)=PCGAZ(JL,JK)/PPIZAZ(JL,JK) |
|
|
! PPIZAZ(JL,JK)=PPIZAZ(JL,JK)/PTAUAZ(JL,JK) |
|
95 |
ztray = prayl(jl)*pdsig(jl, jk) |
ztray = prayl(jl)*pdsig(jl, jk) |
96 |
zratio = ztray/(ztray+ptauaz(jl,jk)) |
zratio = ztray/(ztray+ptauaz(jl,jk)) |
97 |
zgar = pcgaz(jl, jk) |
zgar = pcgaz(jl, jk) |
108 |
pcgaz(jl, jk) = 0. |
pcgaz(jl, jk) = 0. |
109 |
ppizaz(jl, jk) = 1. - repsct |
ppizaz(jl, jk) = 1. - repsct |
110 |
END DO |
END DO |
111 |
END IF ! check flag_aer |
END IF |
112 |
END DO |
END DO |
113 |
|
|
114 |
! ------------------------------------------------------------------ |
! ------------------------------------------------------------------ |