65 |
DOUBLE PRECISION zscat(kdlon) |
DOUBLE PRECISION zscat(kdlon) |
66 |
DOUBLE PRECISION ztr(kdlon, 2, kflev+1) |
DOUBLE PRECISION ztr(kdlon, 2, kflev+1) |
67 |
|
|
68 |
INTEGER jl, jk, ja, jkl, jklp1, jaj, jkm1, in |
INTEGER jl, jk, ja, jkl, jklp1, jaj, jkm1 |
69 |
DOUBLE PRECISION ztray, zgar, zratio, zff, zfacoa, zcorae |
DOUBLE PRECISION ztray, zgar, zratio, zff, zfacoa, zcorae |
70 |
DOUBLE PRECISION zmue, zgap, zww, zto, zden, zmu1, zden1 |
DOUBLE PRECISION zmue, zgap, zww, zto, zden, zmu1, zden1 |
71 |
DOUBLE PRECISION zbmu0, zbmu1, zre11 |
DOUBLE PRECISION zbmu0, zbmu1, zre11 |
72 |
|
|
|
! * Prescribed Data for Aerosols: |
|
|
|
|
|
DOUBLE PRECISION taua(2, 5), rpiza(2, 5), rcga(2, 5) |
|
|
SAVE taua, rpiza, rcga |
|
|
DATA ((taua(in,ja),ja=1,5), in=1, 2)/.730719, .912819, .725059, .745405, & |
|
|
.682188, .730719, .912819, .725059, .745405, .682188/ |
|
|
DATA ((rpiza(in,ja),ja=1,5), in=1, 2)/.872212, .982545, .623143, .944887, & |
|
|
.997975, .872212, .982545, .623143, .944887, .997975/ |
|
|
DATA ((rcga(in,ja),ja=1,5), in=1, 2)/.647596, .739002, .580845, .662657, & |
|
|
.624246, .647596, .739002, .580845, .662657, .624246/ |
|
73 |
! ------------------------------------------------------------------ |
! ------------------------------------------------------------------ |
74 |
|
|
75 |
! * 1. OPTICAL PARAMETERS FOR AEROSOLS AND RAYLEIGH |
! * 1. OPTICAL PARAMETERS FOR AEROSOLS AND RAYLEIGH |
86 |
END DO |
END DO |
87 |
|
|
88 |
DO jk = 1, kflev |
DO jk = 1, kflev |
|
! -OB |
|
|
! DO 104 JL = 1, KDLON |
|
|
! PCGAZ(JL,JK) = 0. |
|
|
! PPIZAZ(JL,JK) = 0. |
|
|
! PTAUAZ(JL,JK) = 0. |
|
|
! 104 CONTINUE |
|
|
! -OB |
|
|
! DO 106 JAE=1,5 |
|
|
! DO 105 JL = 1, KDLON |
|
|
! PTAUAZ(JL,JK)=PTAUAZ(JL,JK) |
|
|
! S +PAER(JL,JK,JAE)*TAUA(KNU,JAE) |
|
|
! PPIZAZ(JL,JK)=PPIZAZ(JL,JK)+PAER(JL,JK,JAE) |
|
|
! S * TAUA(KNU,JAE)*RPIZA(KNU,JAE) |
|
|
! PCGAZ(JL,JK) = PCGAZ(JL,JK) +PAER(JL,JK,JAE) |
|
|
! S * TAUA(KNU,JAE)*RPIZA(KNU,JAE)*RCGA(KNU,JAE) |
|
|
! 105 CONTINUE |
|
|
! 106 CONTINUE |
|
|
! -OB |
|
89 |
DO jl = 1, kdlon |
DO jl = 1, kdlon |
90 |
ptauaz(jl, jk) = flag_aer*tauae(jl, jk, knu) |
ptauaz(jl, jk) = flag_aer*tauae(jl, jk, knu) |
91 |
ppizaz(jl, jk) = flag_aer*pizae(jl, jk, knu) |
ppizaz(jl, jk) = flag_aer*pizae(jl, jk, knu) |
114 |
ppizaz(jl, jk) = 1. - repsct |
ppizaz(jl, jk) = 1. - repsct |
115 |
END DO |
END DO |
116 |
END IF ! check flag_aer |
END IF ! check flag_aer |
|
! 107 CONTINUE |
|
|
! PRINT 9107,JK,((PAER(JL,JK,JAE),JAE=1,5) |
|
|
! $ ,PTAUAZ(JL,JK),PPIZAZ(JL,JK),PCGAZ(JL,JK),JL=1,KDLON) |
|
|
! 9107 FORMAT(1X,'SWCLR_107',I3,8E12.5) |
|
|
|
|
117 |
END DO |
END DO |
118 |
|
|
119 |
! ------------------------------------------------------------------ |
! ------------------------------------------------------------------ |