source: LMDZ5/trunk/libf/misc/juldate.F @ 3603

Last change on this file since 3603 was 2239, checked in by Ehouarn Millour, 10 years ago

Reorganizing physics/dynamics interface:

  • what is related to dynamics-physics interface is now in a seperate directory: dynlmdz_phy* for physics in phy*
  • 1d model and related dependencies (including a couple from "dynamics", set up as symbolic links) is now in subdirectory "dyn1d" of phy*.
  • "bibio" directory is now "misc" and should only contain autonomous utilities.
  • "cosp" is now a subdirectory of phylmd.

EM

  • 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
  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 863 bytes
RevLine 
[524]1!
[1403]2! $Id: juldate.F 2239 2015-03-23 07:27:30Z musat $
[524]3!
4        subroutine juldate(ian,imoi,ijou,oh,om,os,tjd,tjdsec)
5c       Sous-routine de changement de date:
6c       gregorien>>>date julienne
7c       En entree:an,mois,jour,heure,min.,sec.
8c       En sortie:tjd
[2197]9        IMPLICIT NONE
10        INTEGER,INTENT(IN) :: ian,imoi,ijou,oh,om,os
11        REAL,INTENT(OUT) :: tjd,tjdsec
12       
13        REAL frac,year,rmon,cf,a,b
14        INTEGER ojou
15       
[524]16        frac=((os/60.+om)/60.+oh)/24.
[1930]17        ojou=dble(ijou)+frac
18            year=dble(ian)
19            rmon=dble(imoi)
[524]20        if (imoi .le. 2) then
21            year=year-1.
22            rmon=rmon+12.
23        endif
24        cf=year+(rmon/100.)+(ojou/10000.)
25        if (cf .ge. 1582.1015) then
26            a=int(year/100)
27            b=2-a+int(a/4)
28        else
29            b=0
30        endif
31        tjd=int(365.25*year)+int(30.6001*(rmon+1))+int(ojou)
32     +   +1720994.5+b
33        tjdsec=(ojou-int(ojou))+(tjd-int(tjd))
34        tjd=int(tjd)+int(tjdsec)
35        tjdsec=tjdsec-int(tjdsec)
36        return
37        end
38
39
Note: See TracBrowser for help on using the repository browser.