matlab读取GRACE的代码

以下是读取GRACE数据的MATLAB代码：

% 设置GRACE数据的路径和文件名
file_path = 'path/to/GRACE/data/';
file_name = 'GRCTellus.JPL.200204_202012.GLO.RL06M.MSCNv02CRI.nc';

% 读取数据
ncid = netcdf.open(fullfile(file_path, file_name), 'NC_NOWRITE');

% 获取变量ID
lon_id = netcdf.inqVarID(ncid, 'lon');
lat_id = netcdf.inqVarID(ncid, 'lat');
time_id = netcdf.inqVarID(ncid, 'time');
data_id = netcdf.inqVarID(ncid, 'lwe_thickness');

% 读取变量
lon = netcdf.getVar(ncid, lon_id);
lat = netcdf.getVar(ncid, lat_id);
time = netcdf.getVar(ncid, time_id);
data = netcdf.getVar(ncid, data_id);

% 关闭文件
netcdf.close(ncid);

% 对时间变量进行转换
time = double(time)/365.25 + datenum('2002-01-01');

% 绘制数据图
figure
imagesc(lon, lat, data(:,:,1)');
axis xy
colorbar
title(datestr(time(1)))