Changeset 1955 for LMDZ5/trunk/libf/phylmd/physiq.F90

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

File:

• LMDZ5/trunk/libf/phylmd/physiq.F90 (modified) (2 diffs)

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
   !