ICGEM计算的2190阶次重力异常计算时分辩率设置1分和5分，计算结果除了分辩率不同数值有差异吗？

答：有

实验数据：
计算网站：https://icgem.gfz-potsdam.de/home

实验方法：
1.将5分的数据使用GM的grdsample命令重采样为1分
2.计算第一步的结果与网站直接计算的一分的的结果有何差异。

#!/bin/bash
#Program:
#      测试由icgem网站计算的1分和5分分辨率的EGM08重力异常模型有什么差异
#Author:
#   ---,----
#History:
#       2025/06/26
#--------------------------------------------

#绘制5分的数据
gmt begin GravityAnomailyEGM08_5 png E600
gmt set COLOR_BACKGROUND=darkblue
gmt set COLOR_FOREGROUND=darkred
    gmt basemap -JM15c -R135/145/20/30 -Bafg -BWSen+t"(a) EGM08-5 min"
    #gmt plot simulated_true_ssh_karin_SWOT.txt   -W0.01p,black  -Sc0.05c  -Gblack
    # gmt coast -W1/0.5p,black -A1000 -Dh -Ggrey
    gmt makecpt -T-200/200 -Cjet
    gmt xyz2grd gravity_anomaly_EGM2008_5.gdf -h35 -i0,1,3 -R135/145/20/30 -I5m -GGAEGM08_5.nc
    gmt grdimage GAEGM08_5.nc -C
    gmt colorbar -C -DjBC+w13c/0.5c+o0c/-1.5c+m+e -Bxaf  -Byaf+l"mGal"
    gmt end show

#绘制1分的数据
gmt begin GravityAnomailyEGM08_1 png E600
gmt set COLOR_BACKGROUND=darkblue
gmt set COLOR_FOREGROUND=darkred
    gmt basemap -JM15c -R135/145/20/30 -Bafg -BWSen+t"(b) EGM08-1 min"
    #gmt plot simulated_true_ssh_karin_SWOT.txt   -W0.01p,black  -Sc0.05c  -Gblack
    # gmt coast -W1/0.5p,black -A1000 -Dh -Ggrey
    gmt makecpt -T-200/200 -Cjet
    gmt xyz2grd gravity_anomaly_EGM2008_1.gdf -h35 -i0,1,3 -R135/145/20/30 -I1m -GGAEGM08_1.nc
    gmt grdimage GAEGM08_1.nc -C
    gmt colorbar -C -DjBC+w13c/0.5c+o0c/-1.5c+m+e -Bxaf  -Byaf+l"mGal"
    gmt end show
#绘制差异
gmt begin difference png E600
gmt set COLOR_BACKGROUND=darkblue
gmt set COLOR_FOREGROUND=darkred
    gmt basemap -JM15c -R135/145/20/30 -Bafg -BWSen+t"(c) EGM08-difference"
    gmt grdsample GAEGM08_5.nc -R135/145/20/30 -I1m -GGAEGM08_5_1.nc
    gmt grdmath GAEGM08_1.nc GAEGM08_5_1.nc SUB = dif.nc
    gmt makecpt -T-2/2 -Cjet
    gmt grdimage dif.nc -C
    gmt colorbar -C -DjBC+w13c/0.5c+o0c/-1.5c+m+e -Bxaf  -Byaf+l"mGal
    gmt end show

结论：

我不知道ICGEM如何计算重力异常的，本以为209-阶次的分辨率不高，那么五分和1分应该差异不大，但是对比结果来看，插值引入了误差。因此即使5分计算速度更快，但不能拿来插值后使用，在网站计算时还是用多大的分辨率就把参数修改为多大吧。