source: LMDZ5/trunk/libf/phylmd/phys_cal_mod.F90 @ 5362

Last change on this file since 5362 was 2802, checked in by oboucher, 8 years ago

Changed comment for days_elapsed (upon recommendation of Laurent F)

  • Property copyright set to
    Name of program: LMDZ
    Creation date: 1984
    Version: LMDZ5
    License: CeCILL version 2
    Holder: Laboratoire de m\'et\'eorologie dynamique, CNRS, UMR 8539
    See the license file in the root directory
File size: 2.6 KB
RevLine 
[1227]1! $Id:$
2MODULE phys_cal_mod
[2344]3! This module contains information on the calendar at the current time step
[1227]4
[2344]5  INTEGER,SAVE :: year_cur      ! current year
6!$OMP THREADPRIVATE(year_cur)
7  INTEGER,SAVE :: mth_cur       ! current month
8!$OMP THREADPRIVATE(mth_cur)
9  INTEGER,SAVE :: day_cur       ! current day
10!$OMP THREADPRIVATE(day_cur)
[2802]11  INTEGER,SAVE :: days_elapsed  ! number of whole days since start of the current year
[2344]12!$OMP THREADPRIVATE(days_elapsed)
13  INTEGER,SAVE :: mth_len       ! number of days in the current month
14!$OMP THREADPRIVATE(mth_len)
15  INTEGER,SAVE :: year_len      ! number of days in the current year
16!$OMP THREADPRIVATE(year_len)
[2491]17  REAL,SAVE    :: hour         ! seconds elapsed (in the current day) since midnight
[2344]18!$OMP THREADPRIVATE(hour)
19  REAL,SAVE    :: jD_1jan
20!$OMP THREADPRIVATE(jD_1jan)
21  REAL,SAVE    :: jH_1jan
22!$OMP THREADPRIVATE(jH_1jan)
23  REAL,SAVE    :: xjour
24!$OMP THREADPRIVATE(xjour)
25  REAL,SAVE    :: jD_cur  ! jour courant a l'appel de la physique (jour julien)
26!$OMP THREADPRIVATE(jD_cur)
27  REAL,SAVE    :: jH_cur  ! heure courante a l'appel de la physique (jour julien)
28!$OMP THREADPRIVATE(jH_cur)
29  REAL,SAVE    :: jD_ref  ! jour du demarage de la simulation (jour julien)
30!$OMP THREADPRIVATE(jD_ref)
[2358]31 CHARACTER (len=10) :: calend ! type of calendar to use
32                              ! (can be earth_360d, earth_365d or earth_366d)
33!$OMP THREADPRIVATE(calend)
[1227]34
35CONTAINS
36 
[2344]37  SUBROUTINE phys_cal_init(annee_ref,day_ref)
[2802]38
39    USE IOIPSL, ONLY:  ymds2ju
40    USE ioipsl_getin_p_mod, ONLY: getin_p
41
42    IMPLICIT NONE
[2344]43    INTEGER,INTENT(IN) :: annee_ref
44    INTEGER,INTENT(IN) :: day_ref
45
[2358]46    ! Find out which type of calendar we are using
47    calend = 'earth_360d' ! default
48    CALL getin_p("calend",calend)
49     
[2344]50    CALL ymds2ju(annee_ref, 1, day_ref, 0., jD_ref)
51    jD_ref=INT(jD_ref)
52 
53  END SUBROUTINE  phys_cal_init
54
[2422]55  SUBROUTINE phys_cal_update(julian_date)
[1227]56    ! This subroutine updates the module saved variables.
57
[2098]58    USE IOIPSL, only: ju2ymds, ymds2ju, ioget_mon_len, ioget_year_len
[2422]59    IMPLICIT NONE
60    REAL, INTENT(IN) :: julian_date
61
62    jD_cur=INT(julian_date)
63    jH_cur=julian_date-jD_cur
[1227]64   
65    CALL ju2ymds(jD_cur+jH_cur, year_cur, mth_cur, day_cur, hour)
66    CALL ymds2ju(year_cur, 1, 1, 0., jD_1jan)
67   
68    jH_1jan = jD_1jan - int (jD_1jan)
69    jD_1jan = int (jD_1jan)
70    xjour = jD_cur - jD_1jan
71    days_elapsed = jD_cur - jD_1jan
72
[2422]73    ! Get lenght of current month
[1227]74    mth_len = ioget_mon_len(year_cur,mth_cur)
75
[2422]76    ! Get length of current year
[2098]77    year_len = ioget_year_len(year_cur)
78
[1227]79  END SUBROUTINE phys_cal_update
80
81END MODULE phys_cal_mod
Note: See TracBrowser for help on using the repository browser.