看到網(wǎng)上很多的同學(xué)和粉絲都提到了馮偉老師GRACE程序包的數(shù)據(jù)讀取問題。這里我將詳細(xì)介紹具體的代碼調(diào)試過程。第一個(gè)專欄是如何讀取自帶的RL05數(shù)據(jù)，目前數(shù)據(jù)已經(jīng)發(fā)布到RL06了，因此我將在這里介紹本人改進(jìn)的可以讀取RL06數(shù)據(jù)的一套代碼。其實(shí)大多數(shù)函數(shù)都沒有改變，只是修改了讀取GSM數(shù)據(jù)、一階項(xiàng)、二階項(xiàng)的函數(shù)。下載地址：

我們下載馮偉老師的工具包，我在這里把我調(diào)好的可以讀取RL05數(shù)據(jù)的程序放在以下的鏈接，感興趣的可以下載調(diào)試，后續(xù)會(huì)上傳可以讀取最新的RL06數(shù)據(jù)的程序包。?

1. 前期準(zhǔn)備：

接著放入對(duì)應(yīng)的目錄下，C:\GRACE_Matlab_Toolbox\GRACE_data\RL06

?2.?控制文件的修改

相比于RL05數(shù)據(jù)的控制文件，讀取ICGEM網(wǎng)站的CSR RL06數(shù)據(jù)需要對(duì)應(yīng)修改如下：

注意：控制文件中需要讀取目錄下所有GSM的文件名，手動(dòng)修改是不可能的，可以采用以下的方法：在文件對(duì)應(yīng)得目錄下，新建txt文件，里面寫入

并將后綴名改為.bat，運(yùn)行即可得到所有文件得txt文件，復(fù)制粘貼即可。

需要注意的是，我沒有采用具有GUI界面的那個(gè)程序。下面切入正題：直接運(yùn)行即可。

最主要得函數(shù)修改如下：

（1）讀取GRACE GSM數(shù)據(jù)

（2）讀取一階項(xiàng)

（3）讀取二階項(xiàng)，C20沒有改變，只是修改了C21S21C22S22的函數(shù)。

希望對(duì)大家有幫助！有問題，歡迎加入交流群