Changeset 1955 for LMDZ5/trunk/libf/phylmd

Timestamp:

Jan 31, 2014, 9:17:49 AM (10 years ago)

Author:

fhourdin

Message:

Options for idealized ozone (read_climoz=0) :
Imposing ozone be constant in case were solarlong0> -999
(defaut value is -999.999 when seasonal cycle is compiuted) and
to be symetric with respect to equator when read_climoz=-1
FH

Location:

LMDZ5/trunk/libf/phylmd

Files:

• ozonecm_m.F90 (modified) (5 diffs)
• physiq.F90 (modified) (2 diffs)

LMDZ5/trunk/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
 

LMDZ5/trunk/libf/phylmd/physiq.F90

@@ -870 +870 @@
   integer iostat
 
+  REAL zzz
+
   !======================================================================
   ! Gestion calendrier : mise a jour du module phys_cal_mod
@@ -1604 +1606 @@
         co3i = ro3i
      end if
-  elseif (MOD(itap-1,lmt_pas) == 0) THEN
+  ELSEIF (MOD(itap-1,lmt_pas) == 0) THEN
      !        Once per day, update ozone from Royer:
-     wo(:, :, 1) = ozonecm(rlat, paprs, rjour=real(days_elapsed+1))
+
+     IF (solarlong0<-999.) then
+     ! Generic case with evolvoing season
+        zzz=real(days_elapsed+1)
+     ELSE IF (abs(solarlong0-1000.)<1.e-4) then
+     ! Particular case with annual mean insolation
+        zzz=real(90) ! could be revisited
+        IF (read_climoz/=-1) THEN
+           abort_message ='read_climoz=-1 is recommended when solarlong0=1000.'
+           CALL abort_gcm (modname,abort_message,1)
+        ENDIF
+     ELSE
+     ! Case where the season is imposed with solarlong0
+        zzz=real(90) ! could be revisited
+     ENDIF
+     wo(:,:,1)=ozonecm(rlat, paprs,read_climoz,rjour=zzz)
   ENDIF
   !