5 |
contains |
contains |
6 |
|
|
7 |
SUBROUTINE SW(PSCT, PRMU0, PFRAC, PPMB, PDP, PPSOL, PALBD, PALBP, PTAVE, & |
SUBROUTINE SW(PSCT, PRMU0, PFRAC, PPMB, PDP, PPSOL, PALBD, PALBP, PTAVE, & |
8 |
PWV, PQS, POZON, PCLDSW, PTAU, POMEGA, PCG, PHEAT, PHEAT0, PALBPLA, & |
PWV, PQS, POZON, PCLDSW, PTAU, POMEGA, PCG, PHEAT, PHEAT0, PTOPSW, & |
9 |
PTOPSW, PSOLSW, PTOPSW0, PSOLSW0, ZFSUP, ZFSDN, ZFSUP0, ZFSDN0, & |
PSOLSW, PTOPSW0, PSOLSW0, ZFSUP, ZFSDN, ZFSUP0, ZFSDN0, PTOPSWAD, & |
10 |
PTOPSWAD, PSOLSWAD, ok_ade) |
PSOLSWAD, ok_ade) |
11 |
|
|
12 |
! Purpose. |
! Purpose. |
13 |
! This routine computes the shortwave radiation fluxes in two |
! This routine computes the shortwave radiation fluxes in two |
38 |
|
|
39 |
! ARGUMENTS: |
! ARGUMENTS: |
40 |
|
|
41 |
DOUBLE PRECISION PSCT ! constante solaire (valeur conseillee: 1370) |
DOUBLE PRECISION PSCT ! constante solaire |
42 |
DOUBLE PRECISION PRMU0(KDLON) ! COSINE OF ZENITHAL ANGLE |
DOUBLE PRECISION PRMU0(KDLON) ! COSINE OF ZENITHAL ANGLE |
43 |
DOUBLE PRECISION PFRAC(KDLON) ! fraction de la journee |
DOUBLE PRECISION PFRAC(KDLON) ! fraction de la journee |
44 |
DOUBLE PRECISION PPMB(KDLON, KFLEV+1) ! HALF-LEVEL PRESSURE (MB) |
DOUBLE PRECISION PPMB(KDLON, KFLEV+1) ! HALF-LEVEL PRESSURE (MB) |
56 |
DOUBLE PRECISION PCG(KDLON, 2, KFLEV) ! ASYMETRY FACTOR |
DOUBLE PRECISION PCG(KDLON, 2, KFLEV) ! ASYMETRY FACTOR |
57 |
DOUBLE PRECISION PHEAT(KDLON, KFLEV) ! SHORTWAVE HEATING (K/DAY) |
DOUBLE PRECISION PHEAT(KDLON, KFLEV) ! SHORTWAVE HEATING (K/DAY) |
58 |
DOUBLE PRECISION PHEAT0(KDLON, KFLEV)! SHORTWAVE HEATING (K/DAY) clear-sky |
DOUBLE PRECISION PHEAT0(KDLON, KFLEV)! SHORTWAVE HEATING (K/DAY) clear-sky |
|
DOUBLE PRECISION PALBPLA(KDLON) ! PLANETARY ALBEDO |
|
59 |
DOUBLE PRECISION PTOPSW(KDLON) ! SHORTWAVE FLUX AT T.O.A. |
DOUBLE PRECISION PTOPSW(KDLON) ! SHORTWAVE FLUX AT T.O.A. |
60 |
DOUBLE PRECISION PSOLSW(KDLON) ! SHORTWAVE FLUX AT SURFACE |
DOUBLE PRECISION PSOLSW(KDLON) ! SHORTWAVE FLUX AT SURFACE |
61 |
DOUBLE PRECISION PTOPSW0(KDLON) ! SHORTWAVE FLUX AT T.O.A. (CLEAR-SKY) |
DOUBLE PRECISION PTOPSW0(KDLON) ! SHORTWAVE FLUX AT T.O.A. (CLEAR-SKY) |
200 |
ENDDO |
ENDDO |
201 |
ENDDO |
ENDDO |
202 |
DO i = 1, KDLON |
DO i = 1, KDLON |
|
PALBPLA(i) = ZFSUP(i, KFLEV+1)/(ZFSDN(i, KFLEV+1)+1.0e-20) |
|
|
|
|
203 |
PSOLSW(i) = ZFSDN(i, 1) - ZFSUP(i, 1) |
PSOLSW(i) = ZFSDN(i, 1) - ZFSUP(i, 1) |
204 |
PTOPSW(i) = ZFSDN(i, KFLEV+1) - ZFSUP(i, KFLEV+1) |
PTOPSW(i) = ZFSDN(i, KFLEV+1) - ZFSUP(i, KFLEV+1) |
205 |
|
|