nc转成dat后，用Fortran读取dat出错

温眸Cutey

第一次发帖，跪求大神帮助！
首先，原始数据是nc数据，我用gs将数据提取成了dat，新的dat和原始nc数据画图能够重合，所以提取的dat数据应该是没问题的。
之后，就是用Fortran处理数据，为了保险起见，我想看一下Fortran提取dat是不是对的，结果还真是不对……


这是gs，用于将nc转成dat：
*gs
'reinit'
'sdfopen e:\0800000.nc'
'set gxout fwrite'
'set fwrite e:\08t.dat'
'set x 1 401'
'set y 1 241'
i=27
while(i<=32)
'set t 'i
k=1
while(k<=20)
'set z 'k
'd t'
k=k+1
endwhile
i=i+1
endwhile
'disable fwrite'



这是Fortran程序，只简单的提取了（2,2,2,2）与grads画的值比较
  program main
  real t(2,2,2,2)
  open(1,file='e:\08t.dat',form='binary',status='old')
  do it=1,2
   do k=1,2
    do j=1,2
     do i=1,2
      read(1)t(i,j,k,it)
      write(*,*)t(i,j,k,it)
     enddo
    enddo
   enddo
  enddo
  close(1)
end


这是Fortran的运行结果：

运行Fortran结果

这是grads的运行结果，只验证了t（1,1,1,1），t（2,1,1,1），t（1,2,1,1），t（2,2,1,1），结果只有前两个数一样，后面都不一样。

grads结果

请大神们赐教啊~~~~我Fortran实在不好……不知道哪里出问题了。

温眸Cutey

木有人么……{:eb302:}{:eb302:}{:eb302:}{:eb302:}

温眸Cutey

发了一天都木有人回我，好忧桑~~~~

lqouc

fortran 写的不对，数据提取是没问题的。

温眸Cutey

lqouc 发表于 2016-8-16 15:20
fortran 写的不对，数据提取是没问题的。

可不可以具体说说呀

lqouc

fortran 里面设置的数组和你想要提取的不一样，想要2222的话必须吧整个文件全部读取，然后print2222才行。你这样输出的是第一时刻第一层的前八个点，也就是1，1，1，1:8。grads里面你设置x 1 8，其他维度全是1输出试试看吧。

温眸Cutey

lqouc 发表于 2016-8-16 19:20
fortran 里面设置的数组和你想要提取的不一样，想要2222的话必须吧整个文件全部读取，然后print2222才行。 ...

哇~~~是真的耶！谢谢@lqouc

Zhangzh

LZ您好，问一下您转化后的dat文件可以用grads打开吗？为啥我的打不开？请您看看，谢谢