Changeset 1215 for LMDZ4/branches/LMDZ4-dev/libf/phylmd

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

Location:

LMDZ4/branches/LMDZ4-dev/libf/phylmd

Files:

• ozonecm_m.F90 (modified) (4 diffs)
• phys_output_mod.F90 (modified) (1 diff)
• phys_output_write.h (modified) (1 diff)
• phys_state_var_mod.F90 (modified) (1 diff)
• physiq.F (modified) (2 diffs)
• radiation_AR4.F (modified) (5 diffs)
• radlwsw.F90 (modified) (5 diffs)
• sw_aeroAR4.F90 (modified) (2 diffs)

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

@@ -6 +6 @@
 contains
 
-  SUBROUTINE ozonecm(rjour,rlat,paprs,o3)
+  function ozonecm(rjour,rlat,paprs)
 
     ! The ozone climatology is based on an analytic formula which fits the
@@ -22 +22 @@
 
     USE dimphy, only: klon, klev
+    use assert_m, only: assert
 
     REAL, INTENT (IN) :: rjour
-    REAL, INTENT (IN) :: rlat(klon), paprs(klon,klev+1)
-    REAL o3(klon,klev)
-    ! "o3(j, k)" is the column-density of ozone in cell "(j, k)", that is
+    REAL, INTENT (IN) :: rlat(:) ! (klon)
+    REAL, INTENT (IN) :: paprs(:, :) ! (klon,klev+1)
+
+    REAL ozonecm(klon,klev)
+    ! "ozonecm(j, k)" is the column-density of ozone in cell "(j, k)", that is
     ! between interface "k" and interface "k + 1", in kDU.
 
@@ -47 +50 @@
     !----------------------------------------------------------
 
+    call assert((/size(rlat), size(paprs, 1)/) == klon, "ozonecm klon")
+    call assert(size(paprs, 2) == klev + 1, "ozonecm klev")
+
     pi = 4. * atan(1.)
     DO k = 1, klev
        DO i = 1, klon
           zslat = sin(pi / 180. * rlat(i))
-          zsint = sin(2.*pi*(rjour+15.)/an)
-          zcost = cos(2.*pi*(rjour+15.)/an)
+          zsint = sin(2 * pi * (rjour + 15.) / an)
+          zcost = cos(2 * pi * (rjour + 15.) / an)
           z = 0.0531 + zsint * (-0.001595+0.009443*zslat) &
                + zcost * (-0.001344-0.00346*zslat) &
@@ -79 +85 @@
 
     field(:,klev+1) = 0.
-    forall (k = 1: klev) o3(:,k) = field(:,k) - field(:,k+1)
+    forall (k = 1: klev) ozonecm(:,k) = field(:,k) - field(:,k+1)
 
-  END SUBROUTINE ozonecm
+  END function ozonecm
 
 end module ozonecm_m

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

@@ -907 +907 @@
  CALL histdef3d(iff,o_rnebcon%flag,o_rnebcon%name, "Convective Cloud Fraction", "-")
  CALL histdef3d(iff,o_rhum%flag,o_rhum%name, "Relative humidity", "-")
- CALL histdef3d(iff,o_ozone%flag,o_ozone%name, "Ozone concentration", "ppmv")
+ CALL histdef3d(iff,o_ozone%flag,o_ozone%name, "Ozone mole fraction", "-")
  CALL histdef3d(iff,o_dtphy%flag,o_dtphy%name, "Physics dT", "K/s")
  CALL histdef3d(iff,o_dqphy%flag,o_dqphy%name, "Physics dQ", "(kg/kg)/s")

LMDZ4/branches/LMDZ4-dev/libf/phylmd/phys_output_write.h

@@ -909 +909 @@
 
       IF (o_ozone%flag(iff)<=lev_files(iff)) THEN
-      DO k=1, klev
-       DO i=1, klon
-         zx_tmp_fi3d(i,k)=wo(i,k)*RG/46.6968
-     $                    /(paprs(i,k)-paprs(i,k+1))
-     $                    *(paprs(i,1)/101325.0)
-       ENDDO !i
-      ENDDO !k
-      CALL histwrite_phy(nid_files(iff),o_ozone%name,itau_w,zx_tmp_fi3d)
+         CALL histwrite_phy(nid_files(iff), o_ozone%name, itau_w,
+     $        wo * dobson_u * 1e3 / zmasse / rmo3 * rmd)
       ENDIF
 

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

@@ -202 +202 @@
       REAL,ALLOCATABLE,SAVE :: albsol1(:), albsol2(:)
 !$OMP THREADPRIVATE(albsol1,albsol2)
-      REAL,ALLOCATABLE,SAVE :: wo(:,:)
-!$OMP THREADPRIVATE(wo)
-!
+
+      REAL, ALLOCATABLE, SAVE:: wo(:,:)
+      ! column-density of ozone in a layer, in kilo-Dobsons
+      !$OMP THREADPRIVATE(wo)
+
 ! heat : chauffage solaire
 ! heat0: chauffage solaire ciel clair

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

@@ -1722 +1722 @@
 !           Convert from mole fraction of ozone to column density of ozone in a
 !           cell, in kDU:
-            wo = wo * 48. / 29. * zmasse / dobson_u / 1e3
+            wo = wo * rmo3 / rmd * zmasse / dobson_u / 1e3
 C           (By regridding ozone values for LMDZ only once per day, we
 C           have already neglected the variation of pressure in one
@@ -1728 +1728 @@
 C           "zmasse" changes a little.)
          else
-            CALL ozonecm(real(days_elapsed+1), rlat, paprs, wo)
-
+            wo = ozonecm(real(days_elapsed+1), rlat, paprs)
          end if
       ENDIF

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,

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

@@ -40 +40 @@
   ! t--------input-R- temperature (K)
   ! q--------input-R- vapeur d'eau (en kg/kg)
-  ! wo-------input-R- contenu en ozone (en kg/kg) correction MPL 100505
   ! cldfra---input-R- fraction nuageuse (entre 0 et 1)
   ! cldtaupd---input-R- epaisseur optique des nuages dans le visible (present-day value)
@@ -101 +100 @@
   REAL,    INTENT(in)  :: paprs(KLON,KLEV+1), pplay(KLON,KLEV)
   REAL,    INTENT(in)  :: alb1(KLON), alb2(KLON), tsol(KLON)
-  REAL,    INTENT(in)  :: t(KLON,KLEV), q(KLON,KLEV), wo(KLON,KLEV)
+  REAL,    INTENT(in)  :: t(KLON,KLEV), q(KLON,KLEV)
+
+  REAL, INTENT(in)::wo(KLON,KLEV)
+  ! column-density of ozone in a layer, in kilo-Dobsons
+
   LOGICAL, INTENT(in)  :: ok_ade, ok_aie                                 ! switches whether to use aerosol direct (indirect) effects or not
   REAL,    INTENT(in)  :: cldfra(KLON,KLEV), cldemi(KLON,KLEV), cldtaupd(KLON,KLEV)
@@ -150 +153 @@
   REAL*8 PTL(kdlon,kflev+1), PPMB(kdlon,kflev+1)
   REAL*8 PTAVE(kdlon,kflev)
-  REAL*8 PWV(kdlon,kflev), PQS(kdlon,kflev), POZON(kdlon,kflev)
+  REAL*8 PWV(kdlon,kflev), PQS(kdlon,kflev)
+  real(kind=8) POZON(kdlon,kflev) ! mass fraction of ozone
   REAL*8 PAER(kdlon,kflev,5)
   REAL*8 PCLDLD(kdlon,kflev)
@@ -177 +181 @@
   REAL*8 ztopsw_aero(kdlon,9), ztopsw0_aero(kdlon,9)
   REAL*8 zsolsw_aero(kdlon,9), zsolsw0_aero(kdlon,9)
+  real, parameter:: dobson_u = 2.1415e-05 ! Dobson unit, in kg m-2
+
 
   ! initialisation
@@ -236 +242 @@
         PWV(i,k) = MAX (q(iof+i,k), 1.0e-12)
         PQS(i,k) = PWV(i,k)
-        ! wo:    cm.atm (epaisseur en cm dans la situation standard)
-        ! POZON: kg/kg
-        POZON(i,k) = MAX(wo(iof+i,k),1.0e-12)*RG/46.6968 &
-           /(paprs(iof+i,k)-paprs(iof+i,k+1))&
-           *(paprs(iof+i,1)/101325.0)
+        POZON(i,k) = MAX(wo(iof+i, k), 1e-12) * RG * dobson_u * 1e3 &
+             /(paprs(iof+i, k) - paprs(iof+i, k+1)) &
+             * (paprs(iof+i, 1) / 101325.)
         PCLDLD(i,k) = cldfra(iof+i,k)*cldemi(iof+i,k)
         PCLDLU(i,k) = cldfra(iof+i,k)*cldemi(iof+i,k)

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