用grads将nc转成grd出现问题 求解决啊

潇悦

max020729 发表于 2012-6-25 22:41

                               
登录/注册后可看大图

原来你也在啊  也是一个小月亮啊  我今天发帖子才升级成为月亮的

在这等着呢，主要是太郁闷了，完全看不出来哪里错了啊

max020729

潇悦 发表于 2012-6-25 22:44

                               
登录/注册后可看大图

在这等着呢，主要是太郁闷了，完全看不出来哪里错了啊

坐等大神回复

言深深

额，我用你的gs

：

'reinit'

'sdfopen D:\ypc2012\j\nc\precip.mon.mean.nc'
'set fwrite D:\ypc2012\j\nc\jiangshui.grd'
'set gxout fwrite'

'set lon 1.25 358.75'
'set lat -88.75 88.75'
'set z 1'
'set t 1 369'
'd precip'
'disable fwrite'
;

搞出来了，但是grd文件用fortran读取的时候发现部分出现无法显示的问题，不知道为什么。读取结果用的程序 f.f90 (0 Bytes, 下载次数: 0)

2012-6-25 23:09 上传

点击文件名下载附件 ,注意，此处指输出第一个时间点上的个点数据，需要多时次输出自行在write(2`````)行修改
!************************************
parameter(m=72,n=144,o=369)
dimension  x(m,n,o)

open(1,file='jiangshui.grd',access='direct',form='binary',recl=n*m*o*4)
read(1,rec=1) (((x(i,j,k),j=1,n),i=1,m),k=1,o)
close(1)     

open(2,file='output.txt')
do i=1,m
print*,i
do j=1,n
write(2,"(f16.6)")x(i,j,1)
enddo
enddo
close(2)

end

!************************************
结果是这样的 output.txt (0 Bytes, 下载次数: 16)

2012-6-25 23:09 上传

点击文件名下载附件

潇悦

言深深 发表于 2012-6-25 23:04

                               
登录/注册后可看大图

额，我用你的gs

：

请问这个gs不对的话应该怎么写呀？

言深深

潇悦 发表于 2012-6-25 23:10

                               
登录/注册后可看大图

请问这个gs不对的话应该怎么写呀？

会不会是你的原来的数据有问题呢？或者你原来数据有缺测？我不知道有缺测数据的gs应该咋写呢。不过这个gs应该是可以的啊，总之我能读出来一部分。

max020729

言深深 发表于 2012-6-25 23:12

                               
登录/注册后可看大图

会不会是你的原来的数据有问题呢？或者你原来数据有缺测？我不知道有缺测数据的gs应该咋写呢。不过这个gs ...

不了解啊  之前直接用nc文件画出的图还是能看的  但是转成grd就不能看了  我不想1个1个时间的分开  因为之后还要用来求相关什么的  肿么办啊

言深深

max020729 发表于 2012-6-25 23:16

                               
登录/注册后可看大图

不了解啊  之前直接用nc文件画出的图还是能看的  但是转成grd就不能看了  我不想1个1个时间的分开  因为之 ...

额，不需要分开的啊，直接写成txt文件，我写出来一个时间是为了方便看对不对的。write那一句改成
write(2,"(<o>f16.6)")(x(i,j,k),k=1,o)

就可以了啊

max020729

言深深 发表于 2012-6-25 23:20

                               
登录/注册后可看大图

额，不需要分开的啊，直接写成txt文件，我写出来一个时间是为了方便看对不对的。write那一句改成
write( ...

看txt发现好多格点都没有的数据的啊  这是不是无法用grd文件画出图的原因呢  可是为什么nc的可以画出来呢？  

言深深

max020729 发表于 2012-6-25 23:29

                               
登录/注册后可看大图

看txt发现好多格点都没有的数据的啊  这是不是无法用grd文件画出图的原因呢  可是为什么nc的可以画出来呢 ...

nc画图的时候有没有缺值之类的呢？感觉有些奇怪呢

做个霸气的木头

直觉告诉我，深哥生成数据的中的﹡是缺测值较大造成的！建议查下nc文件的缺测值，对应修改ctl的缺测值