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