Changeset 2491 for LMDZ5/trunk/libf/phylmd

Timestamp:

Apr 5, 2016, 4:29:27 PM (9 years ago)

Author:

Ehouarn Millour

Message:

Fix calendar issue in outputs when in 1D.
EM+MPL+LF

Location:

LMDZ5/trunk/libf/phylmd

Files:

• phys_cal_mod.F90 (modified) (1 diff)
• phys_output_mod.F90 (modified) (2 diffs)

LMDZ5/trunk/libf/phylmd/phys_cal_mod.F90

@@ -15 +15 @@
   INTEGER,SAVE :: year_len      ! number of days in the current year
 !$OMP THREADPRIVATE(year_len)
-  REAL,SAVE    :: hour
+  REAL,SAVE    :: hour         ! seconds elapsed (in the current day) since midnight
 !$OMP THREADPRIVATE(hour)
   REAL,SAVE    :: jD_1jan

LMDZ5/trunk/libf/phylmd/phys_output_mod.F90

@@ -39 +39 @@
     USE ioipsl
     USE phys_cal_mod, only : hour, calend
+    USE time_phylmdz_mod   : pdtphys
     USE mod_phys_lmdz_para
     USE aero_mod, only : naero_spc,name_aero
@@ -255 +256 @@
      idayref = day_ref
      IF (klon_glo==1) THEN
-         CALL ymds2ju(annee_ref, 1, idayref, hour, zjulian)     !jyg
+       ! current_time (used to compute hour) is updated at the begining of
+       ! the physics; to set the correct outputs "initial time" we thus 
+       ! have to use (hour-dtphys).
+         CALL ymds2ju(annee_ref, 1, idayref, hour-pdtphys, zjulian)
+         print *,'phys_output_mod: annee,iday,hour,zjulian=',annee_ref,idayref, hour, zjulian
      ELSE
          CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)