爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 6388|回复: 12

[源代码] 刚学习fortran,最近拿到.grd文件,被强制要求用fortran求年平均值

[复制链接]

新浪微博达人勋

发表于 2015-12-11 14:31:41 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
今天边看书边自己写了段fortran的程序,检查下来没有问题,只是用ifortran运行的时候被提示错误,如图

错误提示

错误提示
。同时附上程序文件:
program main
    implicit none
integer ::x,y,z,n,t,i,j
integer,parameter :: neibian=119603200
integer,parameter :: waibian=119603200
real :: b(neibian)
real :: av(waibian)
real :: a
  open(1,file='j:\modeling\ma_mbc.grd',status='old',form='binary')
read(1,*) (b(i),i=1,119603200)
i=1
j=1
do x=1,128,1
     do y=1,64,1
         do z=1,40,1
             av(j)=a/365
             a=0
             j=j+1
             do t=1,365,1
                 a=a+b(i)
                 i=i+1
             end do
         end do
     end do
end do

open(2,file='j:\ma_mbc.txt',status='new')
write(2,*) 'data=',av(j)

    stop
    end

附件是grd程序的ctl文件 ma_mbc.ctl (2.72 KB, 下载次数: 4)
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-12-11 15:41:25 | 显示全部楼层
read里面*去掉
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-12-11 15:56:09 | 显示全部楼层

是这样的么?
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-12-11 16:14:56 | 显示全部楼层

回帖奖励 +19 金钱


应该是啊,别吓我嘛,图上报错是这么说的啦,其他的代码我根本没看。。。。。。
话说你这学期还教这个么?
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2015-12-11 19:34:47 | 显示全部楼层

求大神帮忙!!!谢谢
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2015-12-11 19:35:32 | 显示全部楼层

read里面你确定*能去掉么?我去试试
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-12-12 10:12:08 | 显示全部楼层
lqouc 发表于 2015-12-11 16:14
应该是啊,别吓我嘛,图上报错是这么说的啦,其他的代码我根本没看。。。。。。
话说你这学期还教这个么 ...

应该不是这个,这个用法是对的,问题估计出在后面的循环和数据存的方式不符合,我根本没有下载那个ctl,所以具体也不知道该怎么改。教啊,我每学期都有,呵呵哒。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-12-12 11:05:10 | 显示全部楼层
晓风残月 发表于 2015-12-12 10:12
应该不是这个,这个用法是对的,问题估计出在后面的循环和数据存的方式不符合,我根本没有下载那个ctl, ...

唉,姐咱私下聊吧。
相信我,不管还有没有别的错,那个地方肯定不对。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2015-12-12 14:50:17 | 显示全部楼层
lqouc 发表于 2015-12-12 11:05
唉,姐咱私下聊吧。
相信我,不管还有没有别的错,那个地方肯定不对。

那关键是改成什么呢?rec=1?
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2015-12-12 14:51:44 | 显示全部楼层
晓风残月 发表于 2015-12-12 10:12
应该不是这个,这个用法是对的,问题估计出在后面的循环和数据存的方式不符合,我根本没有下载那个ctl, ...

麻烦你看一下那个ctl文件吧,谢谢谢谢
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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