source: LMDZ4/branches/LMDZ4-dev-20091210/libf/dyn3d/juldate.F @ 2370

Last change on this file since 2370 was 524, checked in by lmdzadmin, 20 years ago

Initial revision

  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 710 bytes
Line 
1!
2! $Header$
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 real (a-h,o-z)
10        frac=((os/60.+om)/60.+oh)/24.
11        ojou=dfloat(ijou)+frac
12            year=dfloat(ian)
13            rmon=dfloat(imoi)
14        if (imoi .le. 2) then
15            year=year-1.
16            rmon=rmon+12.
17        endif
18        cf=year+(rmon/100.)+(ojou/10000.)
19        if (cf .ge. 1582.1015) then
20            a=int(year/100)
21            b=2-a+int(a/4)
22        else
23            b=0
24        endif
25        tjd=int(365.25*year)+int(30.6001*(rmon+1))+int(ojou)
26     +   +1720994.5+b
27        tjdsec=(ojou-int(ojou))+(tjd-int(tjd))
28        tjd=int(tjd)+int(tjdsec)
29        tjdsec=tjdsec-int(tjdsec)
30        return
31        end
32
33
Note: See TracBrowser for help on using the repository browser.