请教MODIS的MYD08-D3月资料的云水路径读取问题

ljj

用MATLAB读取MODIS的MYD08-D3月资料的云水路径变量后，得到中国区域夏季云水路径的空间分布不合适，并且量值偏大，不知道哪里出错了。望前辈们多多指教，谢谢！在https://ladsweb.modaps.eosdis.nasa.gov/search/order/4这个网址下载的MODIS月资料
以下为MATLAB的代码：
filepath='D:\MODIS\MYD08-M3\';
filename = dir(strcat(filepath,'*.hdf'));
for i=7:186
  clw = hdfread([filepath,filename(i).name], 'Cloud_Water_Path_Liquid_Mean_Mean');
  ciw = hdfread([filepath,filename(i).name], 'Cloud_Water_Path_Ice_Mean_Mean');%读取文件
clw1(:,:,i-6)=clw(:,:);
ciw1(:,:,i-6)=ciw(:,:);
end
v1=double(clw1);
v1(v1==-9999)=nan;
  v2=double(ciw1);
v2(v2==-9999)=nan;
cwc=v2+v1;
[mp,np,tp] = size(cwc);
cwy = reshape (cwc,mp,np,12,tp/12);
cwyy= squeeze(nanmean(cwy(:,:,6:8,:),3));%去掉维数为1的行/列
[lat,lon,year] = size(cwyy);
cwcsum= squeeze(nanmean(cwyy(:,:,:),3));%去掉维数为1的行/列
lon=-180:179;
lat=90:-1:-89;
lonmin=find(lon==70)
lonmax=find(lon==135)
latmin=find(lat==55)
latmax=find(lat==15)
vchina=cwcsum(latmin:latmax,lonmin:lonmax);
nx=max(max(vchina))
ny=min(min(vchina))
lonc=70:135;
latc=55:-1:15;
figure
x1=70
x2=135
y1=55
y2=15
[c,h]=contourf(lonc,latc,vchina)
set (h,'linestyle','none');
  hold on;
  ChinaX(x1,x2,y1,y2) ;
  colorbar

ljj

不好意思，是MYD08-M3，不是MYD08-D3