source: trunk/LMDZ.COMMON/libf/misc/juldate.F @ 3553

Last change on this file since 3553 was 1441, checked in by emillour, 10 years ago

Updates in common dynamics (seq and ) to keep up with updates
in LMDZ5 (up to LMDZ5 trunk, rev 2250):

  • compilation:
  • added test in grid/dimension/makdim to check that # of longitudes is a multiple of 8
  • dyn3d_common:

Bug correction concerning zoom (cf LMDZ5 rev 2218)

  • coefpoly.F becomes coefpoly_m.F90 (in misc)
  • fxhyp.F => fxhyp_m.F90 , fyhyp.F => fyhyp_m.F90
  • new routines for zoom: invert_zoom_x_m.F90 and principal_cshift_m.F90
  • inigeom.F adapted to new zoom definition routines
  • fluxstokenc.F : got rid of calls to initial0()
  • dyn3d:
  • advtrac.F90 : got rid of calls to initial0()
  • conf_gcm.F90 : cosmetic changes and change in default dzoomx,dzoomy values
  • guide_mod.F90 : followed updates from Earth Model
  • gcm.F is now gcm.F90
  • dyn3dpar:
  • advtrac_p.F90, covcont_p.F90, mod_hallo.F90 : cosmetic changes
  • conf_gcm.F90 : cosmetic and changed in default dzoomx,dzoomy values
  • parallel_lmdz.F90 : updates to keep up with Earth model
  • misc:
  • arth.F90 becomes arth_m.F90
  • wxios.F90 updated wrt Earth model changes
  • nrtype.F90 and coefpoly_m.F90 added
  • ran1.F, sort.F, minmax.F, minmax2.F, juldate.F moved over from dyn3d_common

EM
File size: 910 bytes
Line 
1!
2! $Id: juldate.F 2197 2015-02-09 07:13:05Z emillour $
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 NONE
10        INTEGER,INTENT(IN) :: ian,imoi,ijou,oh,om,os
11        REAL,INTENT(OUT) :: tjd,tjdsec
12       
13        REAL frac,year,rmon,cf,a,b
14        INTEGER ojou
15       
16        frac=((os/60.+om)/60.+oh)/24.
17        ojou=dble(ijou)+frac
18            year=dble(ian)
19            rmon=dble(imoi)
20        if (imoi .le. 2) then
21            year=year-1.
22            rmon=rmon+12.
23        endif
24        cf=year+(rmon/100.)+(ojou/10000.)
25        if (cf .ge. 1582.1015) then
26            a=int(year/100)
27            b=2-a+int(a/4)
28        else
29            b=0
30        endif
31        tjd=int(365.25*year)+int(30.6001*(rmon+1))+int(ojou)
32     +   +1720994.5+b
33        tjdsec=(ojou-int(ojou))+(tjd-int(tjd))
34        tjd=int(tjd)+int(tjdsec)
35        tjdsec=tjdsec-int(tjdsec)
36        return
37        end
38
39
Note: See TracBrowser for help on using the repository browser.