GRACE 2 级 GSM 数据处理可以通过运行 Matlab 脚本文件GRACE_Matlab_Toolbox_preprocessing_core.m来执行，其中控制文件作为输入，处理后的网格值或 SH 系数作为输出。在运行程序前需修改控制文件。

控制文件的第一行包含 GSM 文件的总数。

第二行表示高斯滤波的半径。如果在此行中给出"0"，则不会对 SH 系数应用高斯滤波。

第三行指定了所使用的去条带方法，选项为 NONE、SWENSON、CHAMBERS2007、CHAMBERS2012、CHENP3M4、CHENP4M6 和 DUAN。

第四行指定是否使用两个预定义的选项（即GIA_notRemoved和GIA_Removed_Geru）来消除 GIA 效果。

第五行指定了GRACE数据的格式，有两个选项，即ICGEM和GRACE。

第六行指定输出数据的格式和名称。此行中的第一个参数表示输出数据的格式，其形式为 SH 系数 （SH_MAT） 或保存为 Matlab MAT 文件的网格值 （GRID_MAT）。该行中的第二个参数表示 SH 系数的最大程度将另存为输出或用于创建网格化值，该值应小于输入 GSM 数据的最大程度。第三个参数指定输出文件的名称。如果此行中的第一个参数为"GRID_MAT"，则将有第四个参数用于定义具有 0.25、0.5 和 1 选项的网格值的空间分辨率。

第七行指定 C20 文件的完整路径名。如果在此行中给出"NAN"，则原始 C20 值将不会替换为基于 SLR 的值。

第八行和第九行通常设置为"NAN"，但也可用于用具有完整路径名的基于 SLR 的系数替换其他 2 级 SH 系数。

第十行指定 degree-1 文件的完整路径名。

第十一行和第十二行分别指定输入 GSM 文件和输出文件的目录。

从控制文件的第十三行到末尾给出输入GSM文件的名称。

可以使用函数gmt_harmonic进行时间序列的谐波分析。GRACE数据处理过程中的偏差和泄漏可以通过使用gmt_grid2cs和gmt_cs2grid函数来估计。由于GRACE测量误差引起的质量变化估计的不确定性可以基于函数gmt_cs_error来估计。

GRAMAT中的主要Matlab功能：