Changeset 1215 for LMDZ4/branches/LMDZ4-dev/libf/phylmd/sw_aeroAR4.F90

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/sw_aeroAR4.F90 (modified) (2 diffs)

LMDZ4/branches/LMDZ4-dev/libf/phylmd/sw_aeroAR4.F90

@@ -89 +89 @@
   !* LOCAL VARIABLES:
   !
-  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)
@@ -203 +207 @@
         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