/[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

trunk/phylmd/Radlwsw/swclr.f revision 105 by guez, Thu Sep 4 10:40:24 2014 UTC trunk/Sources/phylmd/Radlwsw/swclr.f revision 157 by guez, Mon Jul 20 16:01:49 2015 UTC
# Line 65  SUBROUTINE swclr(knu, paer, flag_aer, ta Line 65  SUBROUTINE swclr(knu, paer, flag_aer, ta
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
# Line 96  SUBROUTINE swclr(knu, paer, flag_aer, ta Line 86  SUBROUTINE swclr(knu, paer, flag_aer, ta
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)
# Line 142  SUBROUTINE swclr(knu, paer, flag_aer, ta Line 114  SUBROUTINE swclr(knu, paer, flag_aer, ta
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    ! ------------------------------------------------------------------    ! ------------------------------------------------------------------
Legend:
Removed from v.105  
changed lines
  Added in v.157
  ViewVC Help
Powered by ViewVC 1.1.21