同一份数据，输出grd和txt不一样？

dzxconan

river 发表于 2013-4-20 15:39

                               
登录/注册后可看大图

换个文件是什么意思，难道不是上面那个程序生成的那个二进制文件？

就是分别读上面生成的txt和grd滴~~

dzxconan

river 发表于 2013-4-20 15:39

                               
登录/注册后可看大图

换个文件是什么意思，难道不是上面那个程序生成的那个二进制文件？

就是分别读上面生成的txt和grd滴~~

river

dzxconan 发表于 2013-4-20 19:42

                               
登录/注册后可看大图

就是分别读上面生成的txt和grd滴~~

那就很奇怪了。加上close（12）  close（13)试试

易小凯

dzxconan 发表于 2013-4-20 11:58

                               
登录/注册后可看大图

parameter (n=160,m=55)
      real ts(n,m)
        open(1,file='ts.grd',form='binary')

grd是binary文件
读的时候，直接read（1）

蔚蓝的天

应该还是程序有问题

dzxconan

易小凯 发表于 2013-4-19 23:07

                               
登录/注册后可看大图

你grd文件是怎么读的呢？

parameter (n=160,m=55)
      real ts(n,m)
        open(1,file='ts.grd',form='binary')
        read(1)((ts(i,j),i=1,n),j=1,m)
        print*,ts
        close(1)

dzxconan

蔚蓝的天 发表于 2013-4-21 11:15

                               
登录/注册后可看大图

应该还是程序有问题

program t
parameter(n=160,m=55,tt=12)
integer it(n,m,tt),ave(m,n)
CHARACTER aa*5,mon(tt)*2
data mon/'01','02','03','04','05','06','07','08','09','10','11','12'/
data aa/'t1601'/
do k=1,tt
aa(4:5)=mon(k)
open(11,file=aa)
read(11,*)((it(i,j,k),i=1,n),j=1,m)
close(11)
enddo
do i=1,m
do j=1,n
do k=6,8
ave(i,j)=ave(i,j)+it(j,i,k)/3
enddo
enddo
enddo
  open(12,file='d:\160\ts.grd',form='binary')
write(12)((ave(i,j),i=1,m),j=1,n)
  open(13,file='d:\160\ts.txt')
write(13,*)((ave(i,j),i=1,m),j=1,n)
end
是的 我现在运行这么程序就变成file not found,unit 12,file d:\160\ts.grd

dzxconan

易小凯 发表于 2013-4-21 00:43

                               
登录/注册后可看大图

grd是binary文件
读的时候，直接read（1）

program t
parameter(n=160,m=55,tt=12)
integer it(n,m,tt),ave(m,n)
CHARACTER aa*5,mon(tt)*2
data mon/'01','02','03','04','05','06','07','08','09','10','11','12'/
data aa/'t1601'/
do k=1,tt
aa(4:5)=mon(k)
open(11,file=aa)
read(11,*)((it(i,j,k),i=1,n),j=1,m)
close(11)
enddo
do i=1,m
do j=1,n
do k=6,8
ave(i,j)=ave(i,j)+it(j,i,k)/3
enddo
enddo
enddo
  open(12,file='d:\160\ts.grd',form='binary')
write(12)((ave(i,j),i=1,m),j=1,n)
  open(13,file='d:\160\ts.txt')
write(13,*)((ave(i,j),i=1,m),j=1,n)
end
我现在运行这么程序就变成file not found,unit 12,file d:\160\ts.grd

dzxconan

river 发表于 2013-4-20 20:13

                               
登录/注册后可看大图

那就很奇怪了。加上close（12）  close（13)试试

program t
parameter(n=160,m=55,tt=12)
integer it(n,m,tt),ave(m,n)
CHARACTER aa*5,mon(tt)*2
data mon/'01','02','03','04','05','06','07','08','09','10','11','12'/
data aa/'t1601'/
do k=1,tt
aa(4:5)=mon(k)
open(11,file=aa)
read(11,*)((it(i,j,k),i=1,n),j=1,m)
close(11)
enddo
do i=1,m
do j=1,n
do k=6,8
ave(i,j)=ave(i,j)+it(j,i,k)/3
enddo
enddo
enddo
  open(12,file='d:\160\ts.grd',form='binary')
write(12)((ave(i,j),i=1,m),j=1,n)
  open(13,file='d:\160\ts.txt')
write(13,*)((ave(i,j),i=1,m),j=1,n)
end
好奇怪 我现在运行这么程序就变成file not found,unit 12,file d:\160\ts.grd

river

dzxconan 发表于 2013-4-21 19:43

                               
登录/注册后可看大图

program t
parameter(n=160,m=55,tt=12)
integer it(n,m,tt),ave(m,n)

把编译、运行过程中形成的那些东西都删掉试试