Changeset 1215 for LMDZ4/branches/LMDZ4-dev/libf/phylmd/ozonecm_m.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/ozonecm_m.F90 (modified) (4 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