source: LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3d/juldate.F @ 5225

Last change on this file since 5225 was 1299, checked in by Laurent Fairhead, 15 years ago

Nettoyage general pour se rapprocher des normes et éviter des erreurs a la
compilation:

  • tous les FLOAT() sont remplacés par des REAL()
  • tous les STOP dans phylmd sont remplacés par des appels à abort_gcm
  • le common control défini dans le fichier control.h est remplacé par le module control_mod pour éviter des messages sur l'alignement des variables dans les déclarations
  • des $Header$ remplacés par des $Id$ pour svn

Quelques remplacements à faire ont pu m'échapper


General cleanup of the code to try and adhere to norms and to prevent some
compilation errors:

  • all FLOAT() instructions have been replaced by REAL() instructions
  • all STOP instructions in phylmd have been replaced by calls to abort_gcm
  • the common block control defined in the control.h file has been replaced by the control_mod to prevent compilation warnings on the alignement of declared variables
  • $Header$ replaced by $Id$ for svn

Some changes which should have been made might have escaped me

  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 706 bytes
Line 
1!
2! $Id: juldate.F 1299 2010-01-20 14:27:21Z abarral $
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.