Changeset 2551

Timestamp:

Jun 9, 2016, 11:30:18 AM (8 years ago)

Author:

jyg

Message:

Some improvements to commits 2548 and 2549

Location:

LMDZ5/trunk/libf/phylmd

Files:

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

LMDZ5/trunk/libf/phylmd/phys_cal_mod.F90

@@ -76 +76 @@
     year_len = ioget_year_len(year_cur)
 
+    print *,'phys_cal_update: jH_cur, hour ',jH_cur, hour
+
   END SUBROUTINE phys_cal_update
 

LMDZ5/trunk/libf/phylmd/phys_output_mod.F90

@@ -56 +56 @@
     include "clesphys.h"
     include "thermcell.h"
+    include "YOMCST.h"
 
     ! ug Nouveaux arguments n\'ecessaires au histwrite_mod:
@@ -263 +264 @@
      ELSE
          CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)
-         CALL ymds2ju(annee_ref, 1, day_ini, start_time*86400, zjulian_start)
+         CALL ymds2ju(annee_ref, 1, day_ini, start_time*rday, zjulian_start)
      END IF
 

LMDZ5/trunk/libf/phylmd/phys_output_write_mod.F90

@@ -331 +331 @@
 
     ! On calcul le nouveau tau:
-!!    itau_w = itau_phy + itap + start_time * day_step_phy      !! essais jyg
-    itau_w = itau_phy + itap                                    !! essais jyg
+    itau_w = itau_phy + itap 
     ! On le donne à iophy pour que les histwrite y aient accès:
     CALL set_itau_iophy(itau_w)

LMDZ5/trunk/libf/phylmd/time_phylmdz_mod.F90

@@ -83 +83 @@
   USE print_control_mod, ONLY: lunout
   IMPLICIT NONE
+  INCLUDE 'YOMCST.h'
     REAL,INTENT(IN) :: pdtphys_
     REAL            :: julian_date
@@ -94 +95 @@
     
     ! Update elapsed time since begining of run:
-    current_time=current_time+pdtphys/86400
+    current_time=current_time+pdtphys/rday
+       print *,'update_time: start_time, current_time ', start_time, current_time
 
     ! Compute corresponding Julian date and update calendar
-    CALL ymds2ju(annee_ref,1,day_ini,(start_time+current_time)*86400,julian_date)
+    CALL ymds2ju(annee_ref,1,day_ini,(start_time+current_time)*rday,julian_date)
     CALL phys_cal_update(julian_date)