WaterGAP Global Hydrology Model (WGHM)数据的下载、读取和全球格网图

一、数据的下载

下载链接：Müller Schmied, H et al. (2022): The global water resources and use model WaterGAP v2.2d - Alternative model output driven with gswp3-w5e5年）

数据详细解释文档：GMD - The global water resources and use model WaterGAP v2.2d: model description and evaluation

二、数据读取

下载对应水文分量的”nc4”文件，以TWS为例，利用如下matlab代码读取

clearvars
addpath('E:\RSE\WGHM\data');
file='watergap_22d_WFDEI-GPCC_histsoc_tws_monthly_1901_2016.nc4';
lon=ncread(file,'lon');lat=ncread(file,'lat');
tws=rot90(ncread(file,'tws'));
index_str=1;index_end=12;num_file=index_end-index_str+1;
TWS=tws(:,:,index_str:index_end);%截取1901年的数据

for kk=1:num_file
tep(1,1)=TWS(1,1,1);temp(:,:)=TWS(:,:,kk);
for jj=1:size(lon,1)
    for ii=1:size(lat,1)
        if temp(ii,jj)==tep
            temp(ii,jj)=0;%将数据中的缺省值替换为0
        end
    end
end
TWS_new(kk,:,:)=temp(:,:);
end
TWS_mean=mean(TWS_new);%去平均
for i=1:num_file
    TWS_new(i,:,:)=TWS_new(i,:,:)-TWS_mean;
end

三、绘制全球格网图

利用冯伟老师的“GRACE_Matlab_Toolbox-master”工具包绘制全球格网图

for j=1:num_file
    gws_1(:,:)=TWS_new(j,:,:);
    gws_1Q=gws_1(:,1:360);gws_1H=gws_1(:,361:720);
    gws_2=[gws_1H,gws_1Q]; %gldas_2(find(isnan(gldas_2)==1))=0;%找到缺失值NaN赋值为0
    TWS(:,:,j)=gws_2(:,:);
end

for kk=1:num_file
temp(:,:)=TWS(:,:,kk);
WGHM(:,:,kk)=flipud(temp);
end
figure(1);
temp(:,:)=WGHM(:,:,1);
gmt_grid2map((temp)*0.1,-30,30,0,0,'cm',['WGHM TWS 1901.01']);