Changeset 1999 for LMDZ5/branches/testing/libf/phylmd/ozonecm_m.F90

Timestamp:

Mar 20, 2014, 10:57:19 AM (10 years ago)

Author:

Laurent Fairhead

Message:

Merged trunk changes r1920:1997 into testing branch

Location:

LMDZ5/branches/testing

Files:

• . (modified) (1 prop)
• libf/phylmd/ozonecm_m.F90 (modified) (5 diffs)

LMDZ5/branches/testing/libf/phylmd/ozonecm_m.F90

@@ -6 +6 @@
 contains
 
-  function ozonecm(rlat, paprs, rjour)
+  function ozonecm(rlat, paprs,read_climoz, rjour)
 
     ! The ozone climatology is based on an analytic formula which fits the
@@ -27 +27 @@
     REAL, INTENT (IN) :: paprs(:, :) ! (klon,klev+1)
     REAL, INTENT (IN) :: rjour
+    INTEGER, INTENT (IN) :: read_climoz
 
     REAL ozonecm(klon,klev)
@@ -45 +46 @@
     REAL, parameter:: an = 360., zo3q3 = 4E-8
     REAL, parameter:: dobson_unit = 2.1415E-5 ! in kg m-2
-    REAL gms, zslat, zsint, zcost, z, ppm, qpm, a
+    REAL gms, zslat,zslat2, zsint, zcost, z, ppm, qpm, a
     REAL asec, bsec, aprim, zo3a3
 
@@ -57 +58 @@
        DO i = 1, klon
           zslat = sin(pi / 180. * rlat(i))
+          zslat2=zslat*zslat
+          IF (read_climoz==-1) zslat=0. ! Imposing hemispheric symetry
           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) &
-               + zslat**2 * (.056222 + zslat**2 &
+               + zslat2 * (.056222 + zslat2 &
                * (-.037609+.012248*zsint+.00521*zcost+.008890*zslat))
           zo3a3 = zo3q3/ps/2.
           z = z - zo3q3*ps
           gms = z
-          ppm = 800. - (500.*zslat+150.*zcost)*zslat
-          qpm = 1.74E-5 - (7.5E-6*zslat+1.7E-6*zcost)*zslat
-          bsec = 2650. + 5000.*zslat**2
+          ppm = 800. - 500.*zslat2 - 150.*zcost*zslat
+          qpm = 1.74E-5 - 7.5E-6*zslat2 - 1.7E-6*zcost*zslat
+          bsec = 2650. + 5000.*zslat2
           a = 4.0*(bsec)**(3./2.)*(ppm)**(3./2.)*(1.0+(bsec/ps)**(3./2.))
           a = a/(bsec**(3./2.)+ppm**(3./2.))**2
@@ -88 +91 @@
     ozonecm = max(ozonecm, 1e-12)
 
+    print*,'ozonecm Version2'
+
   END function ozonecm