Changeset 2819 for LMDZ5/trunk/libf/phylmd/physiq_mod.F90

Timestamp:

Mar 8, 2017, 3:34:12 PM (7 years ago)

Author:

dcugnet

Message:

• Modification to allow previous convention climoz files (360 records for all calendars) to be handled again:
  • read_climoz=1 / 2 (tro3 only / tro3 and tro3_daylight) => 360 records file
  • read_climoz=3 / 4 (same) => daily (calendar compliant) file or monthly (14 records) file
• Few fixes in the ozone forcing reading routine

File:

• LMDZ5/trunk/libf/phylmd/physiq_mod.F90 (modified) (3 diffs)

LMDZ5/trunk/libf/phylmd/physiq_mod.F90

@@ -1021 +1021 @@
                                   = ["tro3         ","tro3_daylight"]
     ! vars_climoz(1:read_climoz): variables names in climoz file.
+    ! vars_climoz(1:read_climoz-2) if read_climoz>2 (temporary)
+    REAL :: ro3i ! 0<=ro3i<=360 ; required time index in NetCDF file for
+                 ! the ozone fields, old method.
 
     include "YOMCST.h"
@@ -1972 +1975 @@
 
           wo(:,:,1)=ozonecm(latitude_deg, paprs,read_climoz,rjour=zzz)
+       ELSE IF(ANY([1,2]==read_climoz)) THEN
+          ro3i = REAL(INT((days_elapsed+jh_cur-jh_1jan)/year_len*360.))
+          IF(INT(ro3i) == 360) ro3i = 359.
+          ! (This should never occur, except perhaps because of roundup
+          ! error. See documentation.)
+          WRITE(*,*)' >> Using 360 records ozone files (old method).'
+          CALL regr_pr_time_av(ncid_climoz, vars_climoz(1:read_climoz),      &
+                               ro3i+1., press_edg_climoz, paprs, wo)
+          FORALL (l = 1: read_climoz) wo(:, :, l) = wo(:, :, l) * rmo3 / rmd &
+               * zmasse / dobson_u / 1e3
        ELSE
           IF(adjust_tropopause) THEN
@@ -1977 +1990 @@
             CALL dyn_tropopause(t_seri, ztsol, paprs, pplay, presnivs, rot,  &
                  pr_tropopause)
-            CALL regr_pr_time_av(ncid_climoz, vars_climoz(1:read_climoz),    &
+            CALL regr_pr_time_av(ncid_climoz, vars_climoz(1:read_climoz-2),  &
                  days_elapsed+jh_cur-jh_1jan, press_edg_climoz, paprs, wo,   &
                  time_climoz,  latitude_deg,  press_cen_climoz, pr_tropopause)
           ELSE
             WRITE(*,*)' >> Interpolating O3 field directly on gcm levels.'
-            CALL regr_pr_time_av(ncid_climoz, vars_climoz(1:read_climoz),    &
+            CALL regr_pr_time_av(ncid_climoz, vars_climoz(1:read_climoz-2),  &
                  days_elapsed+jh_cur-jh_1jan, press_edg_climoz, paprs, wo,   &
                  time_climoz)