Fortran求助呢~~

易小凯

言深深 发表于 2012-3-10 20:09

                               
登录/注册后可看大图

哈哈，我发小呢的一个bug，rn的二月是28天？

复制的时候忘记改了，囧。。。。

言深深

易小凯 发表于 2012-3-11 00:33

                               
登录/注册后可看大图

复制的时候忘记改了，囧。。。。

哈

maplesnow82

How about this one? I think it is more straightforward.
     function dayofy(year,month,day)
     ! give the year, month and day=> the number of days in this year
      integer year,month,day,leapyr
      integer daylis(0:13,0:1)                ! (month,year_common:leap)

      data days/0,31,59,90,120,151,181,212,243,273,304,334,365,365
     $           ,0,31,60,91,121,152,182,213,244,274,305,335,366,366/

      dayofy = days(month-1,leapyr(year)) + day
      return
      end
     
!-----This is the function for leap or common year justification.     
      function leapyr(iyr)

      logical y4 , y100 , y400

      y4 = mod(iyr,4) .eq. 0
      y100 = mod(iyr,100) .eq.0
      y400 = mod(iyr,400) .eq.0
      if ((y4 .and. .not.y100) .or. y400) then
         leapyr = 1
      else
         leapyr = 0
      endif
      return
      end