source: LMDZ5/trunk/libf/dyn3d/juldate.F @ 1930

Last change on this file since 1930 was 1930, checked in by lguez, 11 years ago

abort, dfloat and pause are not in the Fortran standard. Replaced
abort by abort_gcm and dfloat by dble. Note: I modified dyn3dpar files
that were identical to dyn3d modified files.

  • 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: 700 bytes
Line 
1!
2! $Id: juldate.F 1930 2014-01-17 16:45:09Z lguez $
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=dble(ijou)+frac
12            year=dble(ian)
13            rmon=dble(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.