本有志于学术，在去年这个时候尝试过，静下心来努力学习、认真功课，可今年离该条路越来越远，但这里自己还是想保存一些有关自己学术的东西，也希望对别人有些帮助。

   因为今天写论文的时候在网上查资料又走到了这里，看了两篇别人的博客，发现有时候别人的一言两语便胜过自己的万分辛苦的成果，想着自己是不是也可以帮助别人呢。

   自己硕士的方向是基于GRACE还有其他水文数据分析全国各流域的水资源的时空变化，所以牵涉到与水有关的数据，自己基本上都接触过、整理过。

   GRACE数据因为自己不是搞地球物理的，也尝试过从球谐函数滤波反演，自己也有一些成果了，可再看到nasa的产品数据之后，还是选择了用产品数据；GLDAS数据，为许多科研中心下载的grib文件，尝试了许多种提取方式之后（linux下grib程序）选择在windows下用matlab反演，trmm则是参考师兄的结果在envi下提取。

因为论文还没结束，等毕业了会把一些成果共享出来，以便需要的参考。

   现在只贴一些代码和几张图片：

year = [2002*ones(12,1);2003*ones(12,1);2004*ones(12,1);2005*ones(12,1);2006*ones(12,1);2007*ones(12,1);2008*ones(12,1);2009*ones(12,1);2010*ones(12,1);2011*ones(12,1);2012*ones(12,1);2013*ones(12,1)];
month = repmat((1:12)',12,1);
day = ones(144,1);
date = datenum(year, month, day);

figure('color','white'),
%
% GRACE
%
subplot(3,1,1)
titleString = areaTitle;

plot(date,M_tmp(1:144),'Color','b')
%ylim([-10 10]);
datetick('x','yy');
title(strcat(titleString));
xlabel('year')
ylabel('TWS(m)')

set(gcf,'name',strcat(titleString));
%
% TRMM
%
subplot(3,1,2)
titleString = '';

plot(date,S_t,'Color','r')
%ylim([-10 10]);
datetick('x','yy');
title(strcat(titleString));
xlabel('year')
ylabel('Rainfall(m)')

set(gcf,'name',strcat(titleString));
%
% GLDAS
%
subplot(3,1,3)
titleString = '';

plot(date,S_w,'Color','g')
%ylim([-10 10]);
datetick('x','yy');
title(strcat(titleString));
xlabel('year')
ylabel('SM(m)')

set(gcf,'name',strcat(titleString));

print(gcf,'-dtiff','-r300',strcat('./fig/',areaName,'-all.tif'))  

因为时间原因，后续吧