1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
|
!**********************************************************************
! Copyright 1998,1999,2000,2001,2002,2005,2007,2008,2009,2010 *
! Andreas Stohl, Petra Seibert, A. Frank, Gerhard Wotawa, *
! Caroline Forster, Sabine Eckhardt, John Burkhart, Harald Sodemann *
! *
! This file is part of FLEXPART. *
! *
! FLEXPART is free software: you can redistribute it and/or modify *
! it under the terms of the GNU General Public License as published by*
! the Free Software Foundation, either version 3 of the License, or *
! (at your option) any later version. *
! *
! FLEXPART is distributed in the hope that it will be useful, *
! but WITHOUT ANY WARRANTY; without even the implied warranty of *
! MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the *
! GNU General Public License for more details. *
! *
! You should have received a copy of the GNU General Public License *
! along with FLEXPART. If not, see <http://www.gnu.org/licenses/>. *
!**********************************************************************
function juldate(yyyymmdd,hhmiss)
!*****************************************************************************
! *
! Calculates the Julian date *
! *
! AUTHOR: Andreas Stohl (15 October 1993) *
! *
! Variables: *
! dd Day *
! hh Hour *
! hhmiss Hour, minute + second *
! ja,jm,jy help variables *
! juldate Julian Date *
! julday help variable *
! mi Minute *
! mm Month *
! ss Second *
! yyyy Year *
! yyyymmddhh Date and Time *
! *
! Constants: *
! igreg help constant *
! *
!*****************************************************************************
use par_mod, only: dp
implicit none
integer :: yyyymmdd,yyyy,mm,dd,hh,mi,ss,hhmiss
integer :: julday,jy,jm,ja
integer,parameter :: igreg=15+31*(10+12*1582)
real(kind=dp) :: juldate
yyyy=yyyymmdd/10000
mm=(yyyymmdd-10000*yyyy)/100
dd=yyyymmdd-10000*yyyy-100*mm
hh=hhmiss/10000
mi=(hhmiss-10000*hh)/100
ss=hhmiss-10000*hh-100*mi
if (yyyy.eq.0) then
print*, 'there is no year zero.'
stop
end if
if (yyyy.lt.0) yyyy=yyyy+1
if (mm.gt.2) then
jy=yyyy
jm=mm+1
else
jy=yyyy-1
jm=mm+13
endif
julday=int(365.25*jy)+int(30.6001*jm)+dd+1720995
if (dd+31*(mm+12*yyyy).ge.igreg) then
ja=int(0.01*jy)
julday=julday+2-ja+int(0.25*ja)
endif
juldate=real(julday,kind=dp) + real(hh,kind=dp)/24._dp + &
real(mi,kind=dp)/1440._dp + real(ss,kind=dp)/86400._dp
end function juldate
|
