/[lmdze]/trunk/phylmd/Radlwsw/sw.f
ViewVC logotype

Diff of /trunk/phylmd/Radlwsw/sw.f

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

revision 117 by guez, Wed Mar 5 14:57:53 2014 UTC revision 118 by guez, Thu Dec 18 17:30:24 2014 UTC
# Line 31  contains Line 31  contains
31      ! 95-01-01 J.-J. Morcrette direct/diffuse albedo      ! 95-01-01 J.-J. Morcrette direct/diffuse albedo
32      ! 03-11-27 J. Quaas Introduce aerosol forcings (based on Boucher)      ! 03-11-27 J. Quaas Introduce aerosol forcings (based on Boucher)
33    
     USE clesphys, ONLY: bug_ozone  
34      USE raddim, ONLY: kdlon, kflev      USE raddim, ONLY: kdlon, kflev
35      USE suphec_m, ONLY: rcpd, rday, rg, md, rmo3      USE suphec_m, ONLY: rcpd, rday, rg
36    
37      ! ARGUMENTS:      ! ARGUMENTS:
38    
# Line 127  contains Line 126  contains
126      DOUBLE PRECISION, save:: ZFSDNAI(KDLON, KFLEV+1)      DOUBLE PRECISION, save:: ZFSDNAI(KDLON, KFLEV+1)
127    
128      logical:: initialized = .false.      logical:: initialized = .false.
129        REAL, PARAMETER :: dobson_u = 2.1415E-05 ! Dobson unit, in kg m-2
130    
131      !-------------------------------------------------------------------      !-------------------------------------------------------------------
132    
# Line 150  contains Line 150  contains
150         DO JK = 1 , KFLEV         DO JK = 1 , KFLEV
151            DO JL = 1, KDLON            DO JL = 1, KDLON
152               ZCLDSW0(JL, JK) = 0.0               ZCLDSW0(JL, JK) = 0.0
153               IF (bug_ozone) then               ZOZ(JL, JK) = POZON(JL, JK) / (dobson_u * 1E3 * rg) * PDP(JL, JK)
                 ZOZ(JL, JK) = POZON(JL, JK)*46.6968/RG &  
                      *PDP(JL, JK)*(101325.0/PPSOL(JL))  
              ELSE  
                 ! Correction MPL 100505  
                 ZOZ(JL, JK) = POZON(JL, JK)*MD/RMO3*46.6968/RG*PDP(JL, JK)  
              ENDIF  
154            ENDDO            ENDDO
155         ENDDO         ENDDO
156    
Legend:
Removed from v.117  
changed lines
  Added in v.118
  ViewVC Help
Powered by ViewVC 1.1.21