最近出一批降水数据处理专题,主要提供代码及处理步骤,数据说明请自己查阅。
另外,最近qq问问题的太多,故暂时停止qq免费答疑。有需求,请付费。
APHRODITE 官网:http://aphrodite.st.hirosaki-u.ac.jp/download/
http://aphrodite.st.hirosaki-u.ac.jp/faq_en.html#Q01 FAQ 这些问题先读一遍
目前已经发展到第二代,具体的说明请自己查阅reademe。
打小横线的代表暂时不公开。
重点看一下EOD mix 和 unmixd的区别。
EOD 是 end of the day 的简称。
了解下世界统一时间UTC如何表示 http://blog.sina.com.cn/s/blog_6575a3c20100jdqe.html (引自新浪博客)
数据处理部分:
ncdisp 在matlab 中查看nc文件的属性
可以看出,降水数据维度为3维,其中时间的表示方法为自1998年第一天之后的分钟数。另外,务必仔细检查降水数据的单位。
接下来开始处理:(复制到matlab即可运行)
clc
clear all
trackin='I:\Forcing\APHRODITE\APHRO_1901\'; %数据输入路径
trackout = 'I:\Forcing\APHRODITE\APHRO_1901\test\' %输出路径
%----经纬度范围
lat_max = 80;
lat_min = 70;
lon_max = 100;
lon_min = 80;
%----年份范围
y_min=1998;
y_max=2015;
vname='precip'; %变量名
offset=0; %变量要加上的值
% offset=86400 % to mm
a=dir([trackin,'*.nc']);
init=0; %初始化标识
for i=1:length(a)
str=a(i).name;
fn1=[trackin,str];
disp(['