/[lmdze]/trunk/Sources/phylmd/Radlwsw/swclr.f
ViewVC logotype

Diff of /trunk/Sources/phylmd/Radlwsw/swclr.f

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 216 by guez, Wed Dec 7 16:44:53 2016 UTC revision 217 by guez, Thu Mar 30 14:25:18 2017 UTC
# Line 4  module swclr_m Line 4  module swclr_m
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
# Line 36  contains Line 35  contains
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)
# Line 88  contains Line 84  contains
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)
# Line 114  contains Line 108  contains
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      ! ------------------------------------------------------------------      ! ------------------------------------------------------------------
Legend:
Removed from v.216  
changed lines
  Added in v.217
  ViewVC Help
Powered by ViewVC 1.1.21