source: LMDZ5/branches/testing/libf/misc/juldate.F @ 5455

Last change on this file since 5455 was 2298, checked in by Laurent Fairhead, 10 years ago

Merged trunk changes -r2237:2291 into testing branch

  • 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
Line 
1!
2! $Id: juldate.F 2298 2015-06-14 19:13:32Z aborella $
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
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       
16        frac=((os/60.+om)/60.+oh)/24.
17        ojou=dble(ijou)+frac
18            year=dble(ian)
19            rmon=dble(imoi)
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.