爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
12
返回列表 发新帖
楼主: 快乐枫

[求助] Fortran求助呢~~

[复制链接]

新浪微博达人勋

发表于 2012-3-11 00:33:28 | 显示全部楼层
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 成长值: 0
发表于 2012-3-11 08:59:09 | 显示全部楼层
易小凯 发表于 2012-3-11 00:33
复制的时候忘记改了,囧。。。。

密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2012-3-11 23:15:41 | 显示全部楼层
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

评分

参与人数 1金钱 +5 收起 理由
mofangbao + 5

查看全部评分

密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表