[2014年 9月 12日 金曜日 09:41:27 JST]
[~/to_manda_sensei/Tools/MGDSST.AlongTrack]
[kunoki@aofd165]
$ ifort lst2utc.f90 -o lst2utc
$ lst2utc
LST: 2012 5 15 23
UTC: 2012 5 15 14
LST: 2012 5 15 5
UTC: 2012 5 14 20
LST: 2012 5 1 5
UTC: 2012 4 30 20
LST: 2012 3 1 5
UTC: 2012 2 29 20
LST: 2011 3 1 5
UTC: 2011 2 28 20
program check_lst2utc
integer lyr,lmo,ldy, lhr, uyr, umo, udy, uhr, tdh
tdh=+9
lyr=2012; lmo=5; ldy=15; lhr=23
call lst2utc(lyr,lmo,ldy, lhr, uyr, umo, udy, uhr, tdh)
print *,"LST: ",lyr, lmo, ldy,lhr
print *,"UTC: ",uyr, umo, udy,uhr
print *
lyr=2012; lmo=5; ldy=15; lhr=5
call lst2utc(lyr,lmo,ldy, lhr, uyr, umo, udy, uhr, tdh)
print *,"LST: ",lyr, lmo, ldy,lhr
print *,"UTC: ",uyr, umo, udy,uhr
print *
lyr=2012; lmo=5; ldy=1; lhr=5
call lst2utc(lyr,lmo,ldy, lhr, uyr, umo, udy, uhr, tdh)
print *,"LST: ",lyr, lmo, ldy,lhr
print *,"UTC: ",uyr, umo, udy,uhr
print *
lyr=2012; lmo=3; ldy=1; lhr=5
call lst2utc(lyr,lmo,ldy, lhr, uyr, umo, udy, uhr, tdh)
print *,"LST: ",lyr, lmo, ldy,lhr
print *,"UTC: ",uyr, umo, udy,uhr
print *
lyr=2011; lmo=3; ldy=1; lhr=5
call lst2utc(lyr,lmo,ldy, lhr, uyr, umo, udy, uhr, tdh)
print *,"LST: ",lyr, lmo, ldy,lhr
print *,"UTC: ",uyr, umo, udy,uhr
print *
end program check_lst2utc
subroutine lst2utc(lyr,lmo,ldy, lhr, uyr, umo, udy, uhr, tdh)
! Description:
!
! Author: am
!
! Host: aofd165.fish.nagasaki-u.ac.jp
! Directory: /work2/kunoki/to_manda_sensei/Tools/MGDSST.AlongTrack
!
! Revision history:
! This file is created by /usr/local/mybin/nff.sh at 09:19 on 09-12-2014.
! Reference
! 1行で書けるうるう年判別法
! http://d.hatena.ne.jp/Kappuccino/20081025/1224906495
!
implicit none
integer,intent(in):: lyr,lmo,ldy, lhr
integer,intent(inout)::uyr,umo,udy, uhr
integer,intent(in):: tdh
integer y
integer,dimension(12)::last_day
data last_day/31, 28, 31, 30, 31, 30, 31, 31, 30, 31, 30, 31/
uyr=lyr
umo=lmo
udy=ldy
uhr=lhr
uhr=lhr-tdh
if(uhr<=0)then
uhr=uhr+24
udy=udy-1
endif
if(uhr==24)then
uhr=0
udy=udy+1
endif
if(udy<=0)then
umo=umo-1
udy=last_day(umo)
if(umo == 2)then
y=uyr
udy = 28 + (1 / (y - y / 4 * 4 + 1)) * &
& (1 - 1 / (y - y / 100 * 100 + 1)) &
& + (1 / (y - y / 400 * 400 + 1));
endif
endif
if(umo<=0)then
umo=1
uyr=uyr-1
endif
end subroutine lst2utc