Context Navigation

juldate.F @ 3603

Last change on this file since 3603 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

File size: 706 bytes

Line
1	!
2	! $Id: juldate.F 1299 2010-01-20 14:27:21Z musat $
3	!
4	subroutine juldate(ian,imoi,ijou,oh,om,os,tjd,tjdsec)
5	c Sous-routine de changement de date:
6	c gregorien>>>date julienne
7	c En entree:an,mois,jour,heure,min.,sec.
8	c 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.25year)+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.