1、CSR官网下载NC文件：GRACE/GRACE-FO - Gravity Recovery and Climate Experiment (utexas.edu)

GLDAS数据也是同样的处理方式；

2、matlab读取NC内容

2.1 matlab将NC文件中的gregorian日期转换成年月日ymd

dt = datetime((time)*24*3600,'ConvertFrom','epochtime','Epoch','2002-01-01');

NC文件中时间time表示从2002年1月1日至目前所要读取文件时间为止累计的天数

% 提取给定时间的索引

t0=datetime(2020,06,15,00,00,00);

[~,ind1] = min(abs(datenum(dt)-datenum(t0)));

closest_time = dt(ind1);

2.2 NC读取数据需要翻转90°

rot90(csr_mascon(:,:,idx));

（不知道nc为什么是经纬度翻转了）matlab直接绘图能看到它是反过来的，所以要转置下；

下面是翻转前后对比；

2.3 获取指定经纬度范围数据

% 经纬度范围

lat_range = [20, 30];

lon_range = [100, 110];

% 查找经纬度范围在nc文件中的索引

indices_lat = find(lat >= lat_range(1) & lat <= lat_range(2));

indices_lon = find(lon >= lon_range(1) & lon <= lon_range(2));

% 提取指定范围内的数据

data_subset = data(indices_lon, indices_lat, :);

3、完整代码截图