Index: LMDZ5/trunk/libf/phylmd/phys_cal_mod.F90
===================================================================
--- LMDZ5/trunk/libf/phylmd/phys_cal_mod.F90	(revision 2490)
+++ LMDZ5/trunk/libf/phylmd/phys_cal_mod.F90	(revision 2491)
@@ -15,5 +15,5 @@
   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
Index: LMDZ5/trunk/libf/phylmd/phys_output_mod.F90
===================================================================
--- LMDZ5/trunk/libf/phylmd/phys_output_mod.F90	(revision 2490)
+++ LMDZ5/trunk/libf/phylmd/phys_output_mod.F90	(revision 2491)
@@ -39,4 +39,5 @@
     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,5 +256,9 @@
      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)