source: LMDZ6/branches/Amaury_dev/libf/obsolete/juldate.f90 @ 5208

! Amaury (07/2024): utilisé nulle part, code de <2004
! Je suggère de bouger dans obsolète

! $Id: juldate.f90 5116 2024-07-24 12:54:37Z fairhead $

  SUBROUTINE juldate(ian,imoi,ijou,oh,om,os,tjd,tjdsec)
  ! Sous-routine de changement de date:
  ! gregorien>>>date julienne
  ! En entree:an,mois,jour,heure,min.,sec.
  ! En sortie:tjd
    IMPLICIT NONE
    INTEGER,INTENT(IN) :: ian,imoi,ijou,oh,om,os
    REAL,INTENT(OUT) :: tjd,tjdsec

    REAL :: frac,year,rmon,cf,a,b
    INTEGER :: ojou

    frac=((os/60.+om)/60.+oh)/24.
    ojou=dble(ijou)+frac
        year=dble(ian)
        rmon=dble(imoi)
    if (imoi <= 2) THEN
        year=year-1.
        rmon=rmon+12.
    endif
    cf=year+(rmon/100.)+(ojou/10000.)
    if (cf >= 1582.1015) THEN
        a=int(year/100)
        b=2-a+int(a/4)
    else
        b=0
    endif
    tjd=int(365.25*year)+int(30.6001*(rmon+1))+int(ojou) &
          +1720994.5+b
    tjdsec=(ojou-int(ojou))+(tjd-int(tjd))
    tjd=int(tjd)+int(tjdsec)
    tjdsec=tjdsec-int(tjdsec)

END SUBROUTINE juldate