
	real*8 d
10      write(*,*)' year,month,day,hour(beijing_time)'
        read(*,*,err=100)iy,im,id,hour
        if(iy.lt.100)iy=iy+1900
        call fjd(d,iy,im,id)
        write(*,"('        --- >   ',f11.2)")d+(hour-8.)/24.
        goto 10
100     stop
        end
 
        subroutine fjd(d,iy,im,id)
        real*8 d
        i=(im+9)*0.09
        i=(i+iy-1900)*1.75+0.01
        j=30.56*im
        j=(iy-1950)*367+id+j-i
        d=2433338.5d0+j
        end

