将txt格式的数据读成dat数据

luom

program main
implicit none

integer, parameter :: ylat=72,xlon=144,day=2,z=1,k=20736

real  :: var(xlon,ylat,z,day)
integer t,irec,l,i,j,m
irec=0
open(1,file='dong2.dat',status='new',access='direct',FORM='unformatted',recl=144*72)
   
open (4,file='d:\dong1.txt',status='old')
do t=1,day
do l= 1,1
irec=irec+1
do j=1,72
do i=1,144

   read(4,100) var(i,j,l,t)
   write(1,rec=irec)var(i,j,l,t)
print*,var(i,j,l,t)
enddo
enddo
enddo
enddo
close(1)
close(4)
100 format(1x,f4.3)
end
运行结果 得到了一个dat数据大小为 81k  字节为82944也符合情况，但是我写的ctl文件怎么打不开这个dat数据请各位指点指点
dset d:\dong2.dat
title monthly jiangshui
undef -999.
xdef 144 linear 1.25 2.5
ydef 72 -88.75 2.5
zdef 1 linear 1000 1
tdef 2 linear 00Z01jan2005 1mo
vars 1
pre       0   99    rainfall
endvars

artmunich

楼主可以在终端上输出少量数据测试一下dat是否写准确了

mofangbao

ctl的ydef仔细检查下，这么粗心的问题自己应该先检查的！！

luom

mofangbao 发表于 2011-12-7 19:21

                               
登录/注册后可看大图

ctl的ydef仔细检查下，这么粗心的问题自己应该先检查的！！

ydef没有错误，只有72个

luom

artmunich 发表于 2011-12-7 18:35

                               
登录/注册后可看大图

楼主可以在终端上输出少量数据测试一下dat是否写准确了

是从得到的dat数据 中输出数据么，怎么检验？？

mofangbao

luom 发表于 2011-12-7 21:24

                               
登录/注册后可看大图

ydef没有错误，只有72个

没发现你ydef少了写 linear么

artmunich

luom 发表于 2011-12-7 21:25

                               
登录/注册后可看大图

是从得到的dat数据 中输出数据么，怎么检验？？

应该是你的fortran程序写数据有问题吧。

luom

artmunich 发表于 2011-12-7 23:07

                               
登录/注册后可看大图

应该是你的fortran程序写数据有问题吧。

我的fortran程序哪里有错 请你指出来吧， 诚心请教

luom

artmunich 发表于 2011-12-7 23:07

                               
登录/注册后可看大图

应该是你的fortran程序写数据有问题吧。

你能给我举一个正确的列子么，

artmunich

你的recl是144*72，这样写可能更好一些吧
irec=1
do t=1,day
do l=1,1
write(1,rec=irec)(var(i,j,l,t),i=1,144),j=1,72)
irec=irec+1
enddo
enddo
你觉得呢