Changeset 1231

Timestamp:

Aug 27, 2009, 2:50:15 PM (15 years ago)

Author:

lguez

Message:

These are changes related to ozone in radiative transfer computations:

-- Moved "max(ozonecm, 1e-12)" from procedure "radlwsw" to procedure

"ozonecm".

-- Removed ratio "p_surface / p_reference" in the computation of

"pozon" in "radlwsw".

-- Removed ratio "p_reference / p_surface" in the computation of

"zoz" in "SW_LMDAR4" and "SW_AEROAR4".

-- Removed ratio "mass of air / mass of ozone" in the computation of

"zoz" in "LW_LMDAR4".

Location:

LMDZ4/branches/LMDZ4-dev/libf

Files:

• dyn3dpar/calfis_p.F (modified) (1 diff)
• phylmd/ozonecm_m.F90 (modified) (1 diff)
• phylmd/radiation_AR4.F (modified) (2 diffs)
• phylmd/radlwsw.F90 (modified) (1 diff)
• phylmd/sw_aeroAR4.F90 (modified) (1 diff)

LMDZ4/branches/LMDZ4-dev/libf/dyn3dpar/calfis_p.F

@@ -210 +210 @@
       SAVE firstcal,debut
 c$OMP THREADPRIVATE(firstcal,debut)
-      REAL :: jD_cur, jH_cur
+      REAL, intent(in):: jD_cur, jH_cur
       
       REAL,SAVE,dimension(1:iim,1:llm):: du_send,du_recv,dv_send,dv_recv

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

@@ -86 +86 @@
     field(:,klev+1) = 0.
     forall (k = 1: klev) ozonecm(:,k) = field(:,k) - field(:,k+1)
+    ozonecm = max(ozonecm, 1e-12)
 
   END function ozonecm

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

@@ -183 +183 @@
          ZCLDSW0(JL,JK) = 0.0
          ZOZ(JL,JK) = POZON(JL,JK) / dobson_u / 1e3 / RG * PDP(JL,JK)
-     $        * (101325. / PPSOL(JL))
       ENDDO
       ENDDO
@@ -2560 +2559 @@
       DO k = 1, KFLEV
          DO i = 1, KDLON
-            ZOZ(i,k) = POZON(i,k)*PDP(i,k) * RMD/RMO3
+            ZOZ(i,k) = POZON(i,k)*PDP(i,k)
          ENDDO
       ENDDO

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

@@ -241 +241 @@
         PWV(i,k) = MAX (q(iof+i,k), 1.0e-12)
         PQS(i,k) = PWV(i,k)
-        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.)
+        POZON(i,k) = wo(iof+i, k) * RG * dobson_u * 1e3 &
+             / (paprs(iof+i, k) - paprs(iof+i, k+1))
         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

@@ -207 +207 @@
         DO JL = 1, KDLON
            ZCLDSW0(JL,JK) = 0.0
-           ZOZ(JL,JK) = POZON(JL,JK) / dobson_u / 1e3 / RG * PDP(JL,JK) &
-                * (101325. / PPSOL(JL))
+           ZOZ(JL,JK) = POZON(JL,JK) / dobson_u / 1e3 / RG * PDP(JL,JK)
         ENDDO
      ENDDO