關於時間轉換公式問題

max533 · 发表于 2014-10-8 20:20:21

由於我有2筆數據，時間是一模一樣的，只是表示方式有點不同

1.變數意義如下：
yr(代表西元年分) dayofyear(代表一年365天的第幾天)
Days_Of_year(代表一年有365或366天) hr(代表一天中第幾個小時)

2.轉換公式：

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

`1`	`timeall` `=` `REAL(yr)+(REAL(dayofyear)-0.99)/Days_Of_Year+REAL(hr)/(Days_Of_Year*24)`

3.數據大概如下：

yr    dayofyear Days_Of_Year hr

1994 001          365                24
1994 002          365                00

透過此算式在fortran計算轉換後，輸出的time_all 照理來說該是要一模一樣，但卻不太一樣，
令人感到困惑，是否有哪位高手能幫忙解答呢?!

楚香饭 · 发表于 2014-10-8 21:23:15

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

01Program www_fcode_cn
02  Implicit none
03  write(*,*) timeall( 1994 , 1 , 365 , 24 )
04  write(*,*) timeall( 1994 , 2 , 365 , 0 )
05contains
06  real function timeall( yr , dayofyear , Days_Of_Year , hr )
07    integer :: yr , dayofyear , Days_Of_Year , hr
08    timeall = REAL(yr)+(REAL(dayofyear)-0.99)/Days_Of_Year+REAL(hr)/(Days_Of_Year*24)
09  end function  timeall
10End Program www_fcode_cn

运行后结果为
1994.00281
1994.00281

感觉没有问题

max533 · 发表于 2014-10-8 21:46:57

這是我的程式碼，不過算出來卻有點怪異，是因為我少注意到哪邊嗎?!
他是一段副程式

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

01subroutine convertday (yr,mon,day,hr,timeall,dayofyear)
02   
03  INTEGER :: yr,mon,day,hr,dayofyear
04  REAL*8 :: timeall
05  Integer :: DAY_OF_MONTH(12) , Days_Of_Year
06  DAY_OF_MONTH = (/31,28,31,30,31,30,31,31,30,31,30,31/)  
07  if (MOD(yr,4)==0 .and. MOD(yr,100)/=0) then
08     DAY_OF_MONTH(2) = 29
09     Days_Of_Year = 366
10  else if(MOD(yr,4)==0 .and. MOD(yr,100)==0 .and. MOD(yr,400)==0) then
11     DAY_OF_MONTH(2) = 29
12     Days_Of_Year = 366
13  else !(MOD(yr,4)==0 .and. MOD(yr,100)==0 .and. MOD(yr,400)/=0) then
14     Days_Of_Year = 365    
15  end if
16     
17  if ( mon > 1 ) then
18    dayofyear = day + sum( DAY_OF_MONTH(1:mon-1) )
19  else
20    dayofyear = day
21  end if
22    timeall = REAL(yr)+(REAL(dayofyear)-0.99)/Days_Of_Year+REAL(hr)/(Days_Of_Year*24)
23 end subroutine

楚香饭 · 发表于 2014-10-8 22:14:19

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

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

`1`	`timeall` `=` `yr1.0D0+(dayofyear-0.99D0)/Days_Of_Year` `+` `1.0D0hr/(Days_Of_Year*24.0D0)`

max533 · 发表于 2014-10-9 14:00:08

楚香饭发表于 2014-10-8 22:14
如果你用双精度，需要把常数也改成双精度。

[mw_shl_code=fortran,true]timeall = yr*1.0D0+(dayofyear-0. ...

感謝您的解答，楚香饭大大。

[通用算法] 關於時間轉換公式問題

新人勋章

爱心勋章

水王勋章

元老勋章

热心勋章