eof后的特征场txt转化为grd时出错

nuistzhou

egvt.txt (1.29 MB, 下载次数: 12)

2012-4-29 13:18 上传

点击文件名下载附件
特征场
由于原场是39*32*54（39是nx，32是ny，54是年数），分辨率为（1.875*1.904）
按照lijianping老师的eof程序分解时，是按照二维数组读入的原场。
因此上述txt里横向是54个模态，纵向是39*32（1248）个格点
我需要把输出成txt的特征场用grads画图，于是用f90转换，但是一直显示end of file during read 。。。
我怀疑可能的问题是：
   虽然我没有选择“自动换行”，但是每个站点的54个模态都是写了两行，即前一行是前51个模态，后3个模态就放到下一行了，不知道这样的格式会不会造成f90读取的错误。。。
现贴上f90程序

program ceshi
implicit none
integer,parameter::nx=39,ny=32,nt=54
integer i,j,k
real m(nx,ny,nt)
open(1,file='e:\bs\eof\egvt.txt')
do j=1,ny
  do i=1,nx
      do k=1,nt
            read(1,*) m(i,j,k)
          enddo
  enddo
enddo

open(2,file='e:\bs\eof\egvt.grd',form='binary')
do k=1,nt
   do j=1,ny
      do i=1,nx
             write(2) m(i,j,k)
          enddo
        enddo
enddo
end

Nafisa6

do j=1,ny
  do i=1,nx
      do k=1,nt
            read(1,*) m(i,j,k)
          enddo
  enddo
enddo

改成这样试试
do j=1,ny
  do i=1,nx
      
            read(1,*) (m(i,j,k),k=1,nt)
         
  enddo
enddo

nuistzhou

Nafisa6 发表于 2012-4-29 14:04

                               
登录/注册后可看大图

do j=1,ny
  do i=1,nx
      do k=1,nt

感谢！！！
可以读进去，并且输出了一个grd。
但是写了个ctl后，画出的图全是重复一个样式。。。
这是我的ctl：
dset e:\bs\eof\egvt.grd
undef 100000000
xdef 39 linear  69.375  1.875
ydef 32 levels 0.9524 2.8571 4.7618 6.6666 8.5713 10.476 12.3808 14.2855 16.1902 18.095 19.9997 21.9044 23.8092 25.7139 27.6186 29.5234 31.4281 33.3328 35.2375 37.1422 39.047 40.9517 42.8564 44.7611 46.6658 48.5705 50.4752 52.3799 54.2846 56.1893 58.094 59.9986
zdef  1 linear 0 1
tdef  54 linear 00dec1957 1yr
vars 1
seof 0 99
endvars
。第一次画eof图，完全不会。求教

Nafisa6

看了一下用你这个数据画出的图，再看看你计算EOF的程序吧，不一定是ctl的问题

nuistzhou

Nafisa6 发表于 2012-4-29 15:15

                               
登录/注册后可看大图

看了一下用你这个数据画出的图，再看看你计算EOF的程序吧，不一定是ctl的问题

eof程序我直接用的lijianping老师的，应该不会有问题。