nc文件转grd最后一个时次数据问题

海阔天又空

存储哈德莱海温1870年1月到2016年12月，共1764个时次，其他时次值都正常，只有最后一个时次1764，值不全仅有一部分在lat-90 54有值，如果画全球分布图提示如下           后两张图是 存完grd后画的有值的部分 和 原始nc数据画的全球分布图

海阔天又空

gs文件
'reinit'
'sdfopen E:\ziliao\HadISST_sst.nc'

'set fwrite E:\ziliao\sst_1870_2016.grd'
'set gxout fwrite'


'set x 1 360'
'set y 1 180'

i=1
while(i<1765)
'set t ' i
'd sst'
i=i+1
endwhile

'disable fwrite'

ctl文件
DSET E:\ziliao\sst_1870_2016.grd
TITLE test
UNDEF -9.99e+08
XDEF   360 LINEAR -179.5 1
YDEF   180 LINEAR -89.5 1
ZDEF   1   LINEAR 0 1
TDEF  1764   LINEAR 11:59Z16JAN1870 1mo
VARS   1
sst 0 99 1monsst
ENDVARS

海阔天又空

最后一个时次有值的地方 值没有问题

river

分辨率是1*1的？为什么ctl 经纬度是带0.5的？原始的经度不是从0开始吗？

海阔天又空

river 发表于 2017-11-13 13:45
分辨率是1*1的？为什么ctl 经纬度是带0.5的？原始的经度不是从0开始吗？

分辨率是1*1的  这时原始信息
ga-> sdfopen E:\ziliao\HadISST_sst.nc
Scanning self-describing file:  E:\ziliao\HadISST_sst.nc
SDF file E:\ziliao\HadISST_sst.nc is open as file 1
LON set to 0 360
LAT set to -89.5 89.5
LEV set to 0 0
Time values set: 1870:1:16:11 1870:1:16:11
E set to 1 1

ga-> q ctlinfo
dset E:\ziliao\HadISST_sst.nc
title Monthly version of HadISST sea surface temperature component
undef -1e+30
dtype netcdf
xdef 360 linear -179.5 1
ydef 180 linear -89.5 1
zdef 1 linear 0 1
tdef 1766 linear 11:59Z16JAN1870 1mo
vars 1
sst=>sst  0  t,y,x  sst
endvars

海阔天又空

river 发表于 2017-11-13 13:45
分辨率是1*1的？为什么ctl 经纬度是带0.5的？原始的经度不是从0开始吗？

其他时次的值和原始nc数据一样只有最后一个时次有问题    拜托大神帮我看看 谢谢啦

river

海阔天又空 发表于 2017-11-13 15:16
其他时次的值和原始nc数据一样只有最后一个时次有问题    拜托大神帮我看看 谢谢啦

你是取了前1764个时次？最后两个时次为啥不要，原始的时次是1766。你试一下把所有的都转成二进制的看看

海阔天又空

river 发表于 2017-11-13 23:11
你是取了前1764个时次？最后两个时次为啥不要，原始的时次是1766。你试一下把所有的都转成二进制的看看

想只截取到2016年12月 全部提取了1764时次有值了 但最后一个值1766也是同样问题  在lat -90 54有值 再取高点纬度会报错 很费解 提取季节平均有同样的问题

pb880991

请问楼主原始数据那里弄的？最近也需要

海阔天又空

pb880991 发表于 2017-11-23 19:17
请问楼主原始数据那里弄的？最近也需要

同学那里考的