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

Timestamp:

Aug 21, 2015, 9:23:13 AM (9 years ago)

Author:

Ehouarn Millour

Message:

Physics/dynamics separation: get rid of all the 'include "temps.h"' in the physics; variables in module time_phylmdz_mod must be used instead. Also added JD_cur, JH_cur and JD_ref in module phys_cal_mod, in preparation for having physics handle its calendar internally.
EM

File:

• LMDZ5/trunk/libf/phylmd/phys_cal_mod.F90 (modified) (1 diff)

LMDZ5/trunk/libf/phylmd/phys_cal_mod.F90

@@ -1 +1 @@
 ! $Id:$
 MODULE phys_cal_mod
-! This module contains information on the calendar at the actual time step
+! This module contains information on the calendar at the current time step
 
-  SAVE
-
-  INTEGER :: year_cur      ! current year
-  INTEGER :: mth_cur       ! current month
-  INTEGER :: day_cur       ! current day
-  INTEGER :: days_elapsed  ! number of whole days since start of the simulation 
-  INTEGER :: mth_len       ! number of days in the current month
-  INTEGER year_len ! number of days in the current year
-  REAL    :: hour
-  REAL    :: jD_1jan
-  REAL    :: jH_1jan
-  REAL    :: xjour
+  INTEGER,SAVE :: year_cur      ! current year
+!$OMP THREADPRIVATE(year_cur)
+  INTEGER,SAVE :: mth_cur       ! current month
+!$OMP THREADPRIVATE(mth_cur)
+  INTEGER,SAVE :: day_cur       ! current day
+!$OMP THREADPRIVATE(day_cur)
+  INTEGER,SAVE :: days_elapsed  ! number of whole days since start of the simulation 
+!$OMP THREADPRIVATE(days_elapsed)
+  INTEGER,SAVE :: mth_len       ! number of days in the current month
+!$OMP THREADPRIVATE(mth_len)
+  INTEGER,SAVE :: year_len      ! number of days in the current year
+!$OMP THREADPRIVATE(year_len)
+  REAL,SAVE    :: hour
+!$OMP THREADPRIVATE(hour)
+  REAL,SAVE    :: jD_1jan
+!$OMP THREADPRIVATE(jD_1jan)
+  REAL,SAVE    :: jH_1jan
+!$OMP THREADPRIVATE(jH_1jan)
+  REAL,SAVE    :: xjour
+!$OMP THREADPRIVATE(xjour)
+  REAL,SAVE    :: jD_cur  ! jour courant a l'appel de la physique (jour julien)
+!$OMP THREADPRIVATE(jD_cur)
+  REAL,SAVE    :: jH_cur  ! heure courante a l'appel de la physique (jour julien)
+!$OMP THREADPRIVATE(jH_cur)
+  REAL,SAVE    :: jD_ref  ! jour du demarage de la simulation (jour julien)
+!$OMP THREADPRIVATE(jD_ref)
 
 
 CONTAINS
   
+  SUBROUTINE phys_cal_init(annee_ref,day_ref)
+  USE IOIPSL, ONLY:  ymds2ju
+  IMPLICIT NONE
+    INTEGER,INTENT(IN) :: annee_ref
+    INTEGER,INTENT(IN) :: day_ref
+
+    CALL ymds2ju(annee_ref, 1, day_ref, 0., jD_ref)
+    jD_ref=INT(jD_ref)
+  
+  END SUBROUTINE  phys_cal_init
+
   SUBROUTINE phys_cal_update(jD_cur, jH_cur)
     ! This subroutine updates the module saved variables.