Changeset 1215 for LMDZ4/branches/LMDZ4-dev/libf/phylmd/radiation_AR4.F

Timestamp:

Jul 23, 2009, 5:52:59 PM (15 years ago)

Author:

lguez

Message:

-- Made "ozonecm" a function instead of a subroutine. Used assumed shape
arguments in "ozonecm".

-- Corrected long name and computation of NetCDF variable "ozone" in
the files "hist*".

-- Corrected comments for ozone variables.

-- In the case "read_climoz", used variables "rmd" and "rmo3" from
"YOMCST.h" instead of writing approximate values.

-- Replaced "real*..." declarations (not conforming to Fortran standard)
by "real(kind=...)" declarations.

-- Replaced value "1./46.6968" in ozone computations by the equivalent
(but clearer) "dobson_u * 1e3" (relative difference ~ 1e-5).

File:

• LMDZ4/branches/LMDZ4-dev/libf/phylmd/radiation_AR4.F (modified) (5 diffs)

LMDZ4/branches/LMDZ4-dev/libf/phylmd/radiation_AR4.F

@@ -89 +89 @@
 C* LOCAL VARIABLES:
 C
-      REAL*8 ZOZ(KDLON,KFLEV)
+      real, parameter:: dobson_u = 2.1415e-05 ! Dobson unit, in kg m-2
+
+      REAL(kind=8) ZOZ(KDLON,KFLEV)
+!     column-density of ozone in layer, in kilo-Dobsons
+
       REAL*8 ZAKI(KDLON,2)     
       REAL*8 ZCLD(KDLON,KFLEV)
@@ -178 +182 @@
       DO JL = 1, KDLON
          ZCLDSW0(JL,JK) = 0.0
-         ZOZ(JL,JK) = POZON(JL,JK)*46.6968/RG
-     .               *PDP(JL,JK)*(101325.0/PPSOL(JL))
+         ZOZ(JL,JK) = POZON(JL,JK) / dobson_u / 1e3 / RG * PDP(JL,JK)
+     $        * (101325. / PPSOL(JL))
       ENDDO
       ENDDO
@@ -2477 +2481 @@
       REAL*8 PPMB(KDLON,KFLEV+1)  ! HALF LEVEL PRESSURE (mb)
       REAL*8 PPSOL(KDLON)         ! SURFACE PRESSURE (Pa)
-      REAL*8 POZON(KDLON,KFLEV)   ! O3 CONCENTRATION (kg/kg)
+      REAL(kind=8) POZON(KDLON,KFLEV)   ! O3 mass fraction
       REAL*8 PTL(KDLON,KFLEV+1)   ! HALF LEVEL TEMPERATURE (K)
       REAL*8 PAER(KDLON,KFLEV,5)  ! OPTICAL THICKNESS OF THE AEROSOLS
@@ -2503 +2507 @@
 C-------------------------------------------------------------------------
       REAL*8 ZABCU(KDLON,NUA,3*KFLEV+1)
-      REAL*8 ZOZ(KDLON,KFLEV)
-c
+
+      REAL(kind=8) ZOZ(KDLON,KFLEV)
+!     equivalent pressure of ozone in a layer, in Pa
+
 cym      REAL*8 ZFLUX(KDLON,2,KFLEV+1) ! RADIATIVE FLUXES (1:up; 2:down)
 cym      REAL*8 ZFLUC(KDLON,2,KFLEV+1) ! CLEAR-SKY RADIATIVE FLUXES
@@ -2550 +2556 @@
 C
       IF (MOD(itaplw0,lw0pas).EQ.0) THEN
-      DO k = 1, KFLEV  ! convertir ozone de kg/kg en pa/pa
-      DO i = 1, KDLON
-c convertir ozone de kg/kg en pa (modif MPL 100505)
-         ZOZ(i,k) = POZON(i,k)*PDP(i,k) * RMD/RMO3
-c        print *,'LW: ZOZ*10**6=',ZOZ(i,k)*1000000.
-      ENDDO
+c     Compute equivalent pressure of ozone from mass fraction:
+      DO k = 1, KFLEV
+         DO i = 1, KDLON
+            ZOZ(i,k) = POZON(i,k)*PDP(i,k) * RMD/RMO3
+         ENDDO
       ENDDO
 cIM ctes ds clesphys.h   CALL LWU(RCO2,RCH4, RN2O, RCFC11, RCFC12,