关于hadley center上的nc数据读取

happyhui

各位师兄师姐，之前我发过这样一个类似帖子，但是无果而终。。。 大家的建议都试过，都不太管用。然而此数据要用，不得不再麻烦大家帮我想想对策。
    不知大家有没有用过hadley center的海温数据？
                                                                                                                                                                                                                                                                                                                            我最近在用hadley center的海温数据做分析，从网站上下载了nc格式的文件，以往的做法都是将nc文件用set gxout fwrie来转化成二进制，从而进行fortran计算。 我对这个语句算是非常熟悉，每次都可以成功转化，但是这次却无论如何都无法转化成二进制数据，输出的数据全是0.
     我想请问，nc格式的文件有很多种么？ 为何只有hadley center的nc数据无法成功转化成二进制数据呢？
    希望大家能给予帮助！
    在此附上我的gs文件：
'reinit'
'sdfopen f:\sam\areacc\hadleycenter\HadISST.nc'
'set fwrite f:\sam\areacc\hadleycenter\sst.dat'
'set gxout fwrite'
i=1
while(i<=1706)
'set lon -179.5 179.5'
'set lat -89.5 89.5'
'set z 1'
'set t 'i
'd sst'
i=i+1
endwhile
'disable fwrite'
'reinit'

做个霸气的木头

语句看上去没有问题，nc的ctl截个图传上来看看

易小凯

只有标准的netCDF文件，grads的sdfopen才能够打开的

happyhui

做个霸气的木头 发表于 2012-7-29 15:14

                               
登录/注册后可看大图

语句看上去没有问题，nc的ctl截个图传上来看看

nc 的metadata信息很传统啊。。 截图不太会传
也是 xdef 360 linear -179.5 1
       ydef 180 linear -89.5 1
       zdef 1 linear 0 1
       tdef 1706 linear.   blablah  跟以前做的nc没什么不同

happyhui

易小凯 发表于 2012-7-29 16:55

                               
登录/注册后可看大图

只有标准的netCDF文件，grads的sdfopen才能够打开的

直接用sdfopen 可以打开这个nc, 而且用nc直接可以出图。
问题在于，我想把他转化成二进制文件用fortran 进一步处理，但是用fwrite怎么着也写入不了二进制文件，全是0值。。。

做个霸气的木头

happyhui 发表于 2012-7-29 18:38

                               
登录/注册后可看大图

nc 的metadata信息很传统啊。。 截图不太会传
也是 xdef 360 linear -179.5 1
       ydef 180 linear  ...

格点也木有错啊！你看到的0值是fortran转化成十进制后看到的，还是你写了ctl画的图，如果是前者，就可能是你fortran程序的问题了

易小凯

happyhui 发表于 2012-7-29 18:40

                               
登录/注册后可看大图

直接用sdfopen 可以打开这个nc, 而且用nc直接可以出图。
问题在于，我想把他转化成二进制文件用fortran  ...

那你sdfopen这个nc文件，然后看能否出图

mofangbao

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

happyhui

做个霸气的木头 发表于 2012-7-29 19:15

                               
登录/注册后可看大图

格点也木有错啊！你看到的0值是fortran转化成十进制后看到的，还是你写了ctl画的图，如果是前者，就可能是 ...

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

happyhui

易小凯 发表于 2012-7-29 21:46

                               
登录/注册后可看大图

那你sdfopen这个nc文件，然后看能否出图

能呀，我一拿到nc就先用sdfopen 打开出了图。非常正常。。。  您可以去Hadley center下一个小数据试一下，我觉得可能时数据格式跟ncep上的不大一样所以不好读取。