请选择 进入手机版 | 继续访问电脑版
爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

新浪微博登陆

只需一步, 快速开始

QQ登录

只需一步,快速开始

搜索
查看: 377|回复: 8

[求助] Fortran合并grd文件

[复制链接] |关注本帖

新浪微博达人勋

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

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

x
我有676个grd文件,现在想把这些文件合并到一个文件中,之前试过用动态数组没弄出来,按照下面的程序会报错,想求教各位大神该如何改啊?
implicit none
integer i,j,k
integer,parameter:: n=676
integer,parameter::xn=72,yn=128
real,dimension(xn,yn)::prec
character*20 ff
open(800,file='E:\grddata\prec.grd',access='append',form='binary')
do k=1,n
   write(ff,*)k
open(k,file='E:\grddata\'//trim(adjustl(ff))//'.grd',form='binary')
  do i=1,xn
      do j=1,yn
     read(k,*,end=900)prec
  write(800,*)prec
900 continue
enddo
enddo
close(k)
enddo
end
报错提示:

QQ图片20180707231305.png
密码修改失败请联系qq:937062711

新浪微博达人勋

发表于 2018-7-8 08:36:14 | 显示全部楼层 |取消关注该作者的回复
这些资料怎么排列的,是grads可以识别的吗?grads可以用的话,直接用grads好了,不要这么麻烦
密码修改失败请联系qq:937062711

新浪微博达人勋

 楼主| 发表于 2018-7-8 08:42:48 | 显示全部楼层 |取消关注该作者的回复
river 发表于 2018-7-8 08:36
这些资料怎么排列的,是grads可以识别的吗?grads可以用的话,直接用grads好了,不要这么麻烦

river大神,这些grd文件grads是可以读的,每个文件只是时间不一样,其他的都是一样的,我需要提取每个文件里面特定格点的数据,因此我想把所有文件合并后再读取,我贴个ctl:
dset E:\grddata\2.grd
title rainfall
undef  -9999.0
options yrev
xdef 128 linear 72.25  0.5
ydef 72 linear  18.25 0.5
zdef 1 linear 0 1
tdef 1 linear Feb1961 1mo
vars 1
prec 0 99  rainfall
endvars
密码修改失败请联系qq:937062711

新浪微博达人勋

发表于 2018-7-8 08:57:48 | 显示全部楼层 |取消关注该作者的回复
原来数据是二维的,676个文件合并以后应该存一个三维的数组吧:
每读一个文件it的数据prec(ix,iy)就复制给数组hebing(ix,iy,it),最后得到一个xn*yn*n的数组,然后再操作就好了
密码修改失败请联系qq:937062711

新浪微博达人勋

发表于 2018-7-8 09:53:40 来自手机 | 显示全部楼层 |取消关注该作者的回复
把ctl改一下,就可以用grads读取
密码修改失败请联系qq:937062711

新浪微博达人勋

 楼主| 发表于 2018-7-8 10:05:00 | 显示全部楼层 |取消关注该作者的回复
qxsdjt 发表于 2018-7-8 09:53
把ctl改一下,就可以用grads读取

请问是用options template吗?怎么改啊
密码修改失败请联系qq:937062711

新浪微博达人勋

 楼主| 发表于 2018-7-8 10:28:59 | 显示全部楼层 |取消关注该作者的回复
nuisttony 发表于 2018-7-8 08:57
原来数据是二维的,676个文件合并以后应该存一个三维的数组吧:
每读一个文件it的数据prec(ix,iy)就复制给 ...

谢谢,已经解决了,二进制文件应该是无格式读取read(文件号)。。。。
密码修改失败请联系qq:937062711

新浪微博达人勋

发表于 2018-7-9 10:20:50 | 显示全部楼层 |取消关注该作者的回复
好东西,谢谢分享
密码修改失败请联系qq:937062711

新浪微博达人勋

发表于 2018-7-9 10:20:56 | 显示全部楼层 |取消关注该作者的回复
好东西,谢谢分享
密码修改失败请联系qq:937062711
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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