关于hadley center上的nc数据读取

happyhui

mofangbao 发表于 2012-7-30 11:21

                               
登录/注册后可看大图

不会又是set lon 和set lat的问题吧，可能和你数据里面的0 360的表示方法不一致，改成对应的 set x set y试 ...

刚才改了试试，还是0000，抓狂了

做个霸气的木头

happyhui 发表于 2012-7-30 13:35

                               
登录/注册后可看大图

就是用fwrite转化成二进制文件之后，直接写了一个ctl，出图——0！  然后我用fortran将数据读出来，也满屏 ...

如果不是细节问题，这个问题就复杂了，要不你传个小的nc文件给我（xiaomu136@163.com）研究下，或者下载地址给我~

易小凯

happyhui 发表于 2012-7-30 13:37

                               
登录/注册后可看大图

能呀，我一拿到nc就先用sdfopen 打开出了图。非常正常。。。  您可以去Hadley center下一个小数据试一下， ...

你把数据的网址发上来

易小凯

我下了一个hadley的数据，发现有不少缺测值

chaseechoes

我的Hadley的数据中有一段时间的值是一样的。。

happyhui

做个霸气的木头 发表于 2012-7-30 14:10

                               
登录/注册后可看大图

如果不是细节问题，这个问题就复杂了，要不你传个小的nc文件给我（xiaomu136@163.com）研究下，或者下载地 ...

http://www.metoffice.gov.uk/hadobs/hadisst/ 谢谢帮忙！

happyhui

易小凯 发表于 2012-7-30 21:16

                               
登录/注册后可看大图

你把数据的网址发上来

http://www.metoffice.gov.uk/hadobs/hadisst/   谢谢帮忙！

happyhui

易小凯 发表于 2012-7-30 22:32

                               
登录/注册后可看大图

我下了一个hadley的数据，发现有不少缺测值

我们老师说hadley的海温数据不错的。。 你的数据能用fwrite转化么？

happyhui

chaseechoes 发表于 2012-7-31 13:19

                               
登录/注册后可看大图

我的Hadley的数据中有一段时间的值是一样的。。

恩，请问你的能转化成二进制数据吗？

做个霸气的木头

happyhui 发表于 2012-7-31 15:08

                               
登录/注册后可看大图

http://www.metoffice.gov.uk/hadobs/hadisst/ 谢谢帮忙！

我下了HadISST1_SST_update.nc.gz数据，先fwrite成二进制，再fortran读取并转换，很成功啊！因为是海温资料，陆地上都是缺测值，所以前面很多都是缺测值（-1E+30），后面会有温度20多度，一切很正常；
数据只有1个时次，一个层次，所以数据量是360*180
对应的ctl：
dset c:\s01.dat
title sst
undef -1e+30
xdef 360 linear -179.5 1
ydef 180 linear -89.5 1
zdef 1 linear 0 1
tdef 1 linear 000Z01may2012 1mn
vars 1
sst 0 99 monthly sst
endvars
对应的fortran程序：
program ex01
real ss(360,180)
integer i,j
open(100,file='s01.dat',form='binary')
open(200,file='data0.txt')
do j=1,180
        read(100) (ss(i,j),i=1,360)
        write(200,*) (ss(i,j),i=1,360)
        write(*,*) j,"   ",ss(1,j)
enddo
end
你对比看看，是哪里的问题~