Changeset 2422 for LMDZ5/trunk/libf/phylmd/phys_cal_mod.F90

Timestamp:

Jan 6, 2016, 12:37:41 PM (9 years ago)

Author:

Ehouarn Millour

Message:

Small modification in the way time and calendar are handled: Now all the time keeping is done in the physics and only the timestep is transfered from the dynamics to the physics. Due to changes in computations and roundoffs this will change reference bench results.
The implementation of this change in phymar is left as future work.
EM

File:

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

LMDZ5/trunk/libf/phylmd/phys_cal_mod.F90

@@ -52 +52 @@
   END SUBROUTINE  phys_cal_init
 
-  SUBROUTINE phys_cal_update(jD_cur, jH_cur)
+  SUBROUTINE phys_cal_update(julian_date)
     ! This subroutine updates the module saved variables.
 
     USE IOIPSL, only: ju2ymds, ymds2ju, ioget_mon_len, ioget_year_len
-    
-    REAL, INTENT(IN) :: jD_cur ! jour courant a l'appel de la physique (jour julien)
-    REAL, INTENT(IN) :: jH_cur ! heure courante a l'appel de la physique (jour julien)
+    IMPLICIT NONE
+    REAL, INTENT(IN) :: julian_date
+
+    jD_cur=INT(julian_date)
+    jH_cur=julian_date-jD_cur
     
     CALL ju2ymds(jD_cur+jH_cur, year_cur, mth_cur, day_cur, hour)
@@ -68 +70 @@
     days_elapsed = jD_cur - jD_1jan
 
-    ! Get lenght of acutual month
+    ! Get lenght of current month
     mth_len = ioget_mon_len(year_cur,mth_cur)
 
+    ! Get length of current year
     year_len = ioget_year_len(year_cur)