CERA-20C，日平均数据的月平均相关问题

大气小白

第一次用欧洲中心的这个资料，在处理降水数据的时候遇到了一些问题，希望家园里用过这个资料的朋友能给予一些帮助。所下载的数据为forecast，surface，monthly means of daily means的总降水量（total precipitation）
用ncdump查询文件描述信息给出的描述信息如下图

缺测值为-32767 单位是米

下图是简单绘制出了六月份 30°S-60°N 60°-180°E的降水空间分布

降水图

可以看出降水量从-30000开始，降水越多的地方，图上数值越趋近于0
用过这套资料的朋友，能不能介绍下如何将这个降水量的单位转化为mm/day，在欧洲中心的网站上找半天也没有找到相应的介绍。
资料来源地址http://apps.ecmwf.int/datasets/data/cera20c-edmo/levtype=sfc/type=fc/

大气小白

问题已解决

Black_kiki

最近刚刚下载了CERA-20C的数据，monthly means of daily means 选了六个参数。对MATLAB和气象数据都不是很熟悉，能不能看一下你写的这个文件读取程序啊？

大气小白

Black_kiki 发表于 2017-11-25 21:38
最近刚刚下载了CERA-20C的数据，monthly means of daily means 选了六个参数。对MATLAB和气象数据都不是很 ...

你下载的是什么格式的

Black_kiki

大气小白 发表于 2017-12-19 16:31
你下载的是什么格式的

我下在的.nc格式的，最近一直在处理这些数据

Black_kiki

大气小白 发表于 2017-12-19 16:31
你下载的是什么格式的

不知道为什么，我的mean sea level pressure 数据 *scale_factor+add_offset 后 数值全是32676，还没找出什么原因

Black_kiki

Black_kiki 发表于 2017-12-19 18:56
不知道为什么，我的mean sea level pressure 数据 *scale_factor+add_offset 后 数值全是32676，还没找出 ...

这个问题已经解决了，msldata 数据类型是int16
我原先写的是msldata = double(msldata*scale_factor+add_offset);
现在改成了msldata = double(msldata)*scale_factor+add_offset;

还想弱弱的问一下；
对于missingdata ，应该是下面哪一种处理方法
msldata = double(msldata)*scale_factor+add_offset;
msldata(msldata ==-32676) = NaN;

还是-32676*scale_factor+add_offset的值设置为NaN

大气小白

Black_kiki 发表于 2017-12-19 19:20
这个问题已经解决了，msldata 数据类型是int16
我原先写的是msldata = double(msldata*scale_factor+add ...

我下载的时候把10个ensmbles下在一个.nc文件里面了 所以默认没有提供描述文件，我自己写了一个CTL，里面加上unpack scale_factor add_offset 这一句就可以了


另外吐槽一句，现在另有几个新的数据要下，突然发现欧洲中心从可以一次下110年改成一次下55年。到现在一次只能下1年，python脚本调半天调不对，也找不到对于这个数据的例子，已经快把我搞疯了

Black_kiki

大气小白 发表于 2017-12-21 22:47
我下载的时候把10个ensmbles下在一个.nc文件里面了 所以默认没有提供描述文件，我自己写了一个CTL，里面 ...

你的意思是你一次下了10个集成模型出来的数据吗？
我当时不知道select number 什么意思，就写了邮件问他们。
他们说，这10个集成模型差不了太多，任选一个就行。

你是用python处理的？我是matlab.

大气小白

Black_kiki 发表于 2017-12-25 08:45
你的意思是你一次下了10个集成模型出来的数据吗？
我当时不知道select number 什么意思，就写了邮件问他 ...

10个number是10个不同初值的预报结果，一般会用多个初值的结果做一个ensmble，我是简单的平均了一下，没有分配权重。
另外我先用grads处理了一下，python是爬数据用的。