關於時間轉換公式問題

楚香饭

[Fortran] 纯文本查看 复制代码

Program www_fcode_cn
  Implicit none
  write(*,*) timeall( 1994 , 1 , 365 , 24 )
  write(*,*) timeall( 1994 , 2 , 365 , 0 )
contains
  real function timeall( yr , dayofyear , Days_Of_Year , hr )
    integer :: yr , dayofyear , Days_Of_Year , hr
    timeall = REAL(yr)+(REAL(dayofyear)-0.99)/Days_Of_Year+REAL(hr)/(Days_Of_Year*24)
  end function  timeall
End Program www_fcode_cn  

运行后结果为
   1994.00281   
   1994.00281   

感觉没有问题

楚香饭

如果你用双精度，需要把常数也改成双精度。

[Fortran] 纯文本查看 复制代码

timeall = yr*1.0D0+(dayofyear-0.99D0)/Days_Of_Year + 1.0D0*hr/(Days_Of_Year*24.0D0)