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

Last change on this file since 2152 was 2098, checked in by lguez, 10 years ago

Replaced 360 in calbeta_clim by length of current year according to
chosen calendar. Length of current year is given by
ioget_year_len. Since we already need this for ozone, moved the call
to ioget_year_len from physiq to phys_cal_mod and created variable
year_len of module phys_cal_mod.

Control the output from minmaxqfi.

Non-ASCII characters in comments are not always rendered properly and
they risk being lost. See revision 1740.
  • 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: 1.3 KB
Line 
1! $Id:$
2MODULE phys_cal_mod
3! This module contains information on the calendar at the actual time step
4
5  SAVE
6
7  INTEGER :: year_cur      ! current year
8  INTEGER :: mth_cur       ! current month
9  INTEGER :: day_cur       ! current day
10  INTEGER :: days_elapsed  ! number of whole days since start of the simulation
11  INTEGER :: mth_len       ! number of days in the current month
12  INTEGER year_len ! number of days in the current year
13  REAL    :: hour
14  REAL    :: jD_1jan
15  REAL    :: jH_1jan
16  REAL    :: xjour
17
18
19CONTAINS
20 
21  SUBROUTINE phys_cal_update(jD_cur, jH_cur)
22    ! This subroutine updates the module saved variables.
23
24    USE IOIPSL, only: ju2ymds, ymds2ju, ioget_mon_len, ioget_year_len
25   
26    REAL, INTENT(IN) :: jD_cur ! jour courant a l'appel de la physique (jour julien)
27    REAL, INTENT(IN) :: jH_cur ! heure courante a l'appel de la physique (jour julien)
28   
29    CALL ju2ymds(jD_cur+jH_cur, year_cur, mth_cur, day_cur, hour)
30    CALL ymds2ju(year_cur, 1, 1, 0., jD_1jan)
31   
32    jH_1jan = jD_1jan - int (jD_1jan)
33    jD_1jan = int (jD_1jan)
34    xjour = jD_cur - jD_1jan
35    days_elapsed = jD_cur - jD_1jan
36
37    ! Get lenght of acutual month
38    mth_len = ioget_mon_len(year_cur,mth_cur)
39
40    year_len = ioget_year_len(year_cur)
41
42  END SUBROUTINE phys_cal_update
43
44END MODULE phys_cal_mod
Note: See TracBrowser for help on using the repository browser.