引自:http://www.gnss.help/2016/10/02/grids-in-gamit/
原作者:山科研究生jiangyingming@live.com
0、前言
在使用 GAMIT/GLOBK 程序進行數據處理時,你肯定已經注意到了其 tables 文件夾中存在許多 *.grid 文件。如果你檢查過 GAMIT/GLOBK 安裝目錄的 tables 文件夾,你會發現更多這樣的文件。它們多數體積都很大,一般為數百 MB。這樣的文件稱為格網模型文件,一般在 MIT 的 FTP 服務器上提供下載。
本文將試圖逐個地介紹其中的格網模型文件,可以算是介紹其格網模型文件的一份文檔。
1、atl.grid
全球的大氣潮格網模型文件,由澳大利亞國家大學的 Paul Tregoning 創建。
該文件包含在 GAMIT/GLOBK 程序的安裝包中且不需要更新。也就是說,不管你處理哪一年的觀測數據,都可以引入這個文件。
要在 GAMIT 數據處理中應用大氣潮格網模型,你首先需要保證在進行數據處理的工程文件夾的 tables 目錄中,其下的 atl.grid 已經正確鏈接到 GAMIT/GLOBK 安裝目錄中的 tables/atl.grid 文件。然后在 sestbl. 中修改設置為:
Use atl.list = N ; Atmospheric tides, list file, not yet available
Use atl.grid = Y ; Atmospheric tides, grid file
2、atml.grid
全球的無潮汐大氣負載參數格網模型文件,該文件由澳大利亞國家大學的 Paul Tregoning 根據盧森堡大學的 Tonie van Dam 提供的數據創建。
該模型一般每年對應一個文件,GAMIT 使用的該文件的命名方式一般為:atmdisp_cm.<year>。如用於 2015 年的文件名為 atmdisp_cm.2015。
要在數據處理中應用無潮汐大氣負載參數格網模型,首先從 MIT 的 FTP 服務器下載適用於當年的模型文件放入 GAMIT/GLOBK 安裝目錄的 tables 文件夾內,然后確保在進行數據處理的工程文件夾的 tables 目錄中,其下的 atml.grid 正確鏈接到當年的模型文件。最后在 sestbl. 中設置:
Use atml.list = N ; Atmospheric (non-tidal) loading list file from LU
Use atml.grid = Y ; Atmospheric (non-tidal) loading grid file from LU, converted to GAMIT format
3、gpt.grid
全球氣壓和溫度模型文件,該模型由 Boehm 和 Schuh 提出。在估計對流層參數時,GAMIT 從該模型中內插獲取觀測所在地區的氣壓和溫度。
該模型一般不會更新,並且已經被包含在 GAMIT/GLOBK 程序的安裝包中。
要在數據處理中使用全球氣壓與溫度模型,首先確保你進行數據處理的工程文件夾的 tables 目錄中,其下的 gpt.grid 已經正確鏈接到 GAMIT/GLOBK 安裝目錄中的 tables/gpt.grid 文件。然后在 sestbl. 中設置:
Met obs source = GPT 50 ; hierarchical list with humidity value at the end; e.g. RNX UFL GPT 50 ; default GTP 50
這樣,在進行數據處理時,GAMIT 就將從該模型中獲取氣壓和溫度參數。但內插出的精度畢竟不如實測數據,因此如果你有氣象觀測數據的話,最好還是將其引入。所以我推薦如下的設置:
Met obs source = RNX UFL GPT 50 ; hierarchical list with humidity value at the end; e.g. RNX UFL GPT 50 ; default GTP 50
如此設置后,需要測站所在地區的氣壓和溫度數據時,GAMIT 程序將首先檢查有沒有對應的 RINEX 格式氣象數據,如果沒有的話再檢查有沒有對應的 U-文件,最后才會嘗試使用全球氣壓與溫度模型提供的數據。
4、map.grid
全球的大氣映射函數模型文件,該文件由澳大利亞國家大學的 Paul Tregoning 根據維也納科技大學的 Johannes Boehm 提供的數據創建。
GAMIT 程序目前支持 GMF、NMFW 和 VMF1 三種映射函數,但只有 VMF1 需要引入模型文件。該模型逐年更新,每年生成一個文件。你可以根據需要在 MIT 的 FTP 服務器上下載,其命名方式為:vmf1grd.<year>。
要在數據處理時應用 VMF1 映射函數模型,需要首先下載適用於當年的模型文件放入 GAMIT/GLOBK 安裝目錄的 tables 文件夾內。然后確保在進行數據處理的工程文件夾的 tables 目錄中,其下的 map.grid 正確鏈接到當年的模型文件。最后在 sestbl. 中設置:
DMap = VMF1 ; GMF(default)/VMF1/NMFH; GMF now invokes GPT2 if gpt.grid is available (default)
WMap = VMF1 ; GMF(default)/VMF1/NMFW; GMF now invokes GPT2 if gpt.grid is available (default)
Use map.list = N ; VMF1 list file with mapping functions, ZHD, ZWD, P, Pw, T, Ht
Use map.grid = Y ; VMF1 grid file with mapping functions and ZHD
5、met.grid
雖然各版本的 GAMIT 一直保留了這個模型文件鏈接,但該模型並沒有什么用。有 sestbl. 中對該項配置的說明為證:
Use met.list = N ; not yet supported
Use met.grid = N ; not yet supported
6、otl.grid
全球的海潮模型文件,GAMIT 現在支持多種海潮模型。MIT 的 FTP 服務器上提供的海潮模型文件有:
-
otl_CSR4.grid
-
otl_FES2004.grid
-
otl_FES99.grid
-
otl_GOT00.grid
-
otl_NAO99b.grid
要在 GAMIT 數據處理中應用全球海潮模型文件,首先下載對應的海潮模型文件放入 GAMIT/GLOBK 安裝目錄的 tables 文件夾內。然后確保在進行數據處理的工程文件夾的 tables 目錄中,其下的 otl.grid 正確鏈接到對應的模型文件。最后在 sestbl. 中設置:
Tides applied = 31 ; Binary coded: 1 earth 2 freq-dep 4 pole 8 ocean 16 remove mean for pole tide
; 32 atmosphere ; default = 31
Use otl.list = N ; Ocean tidal loading list file from OSO
Use otl.grid = Y ; Ocean tidal loading grid file, GAMIT-format converted from OSO
這里的 “Tides applied”,以二進制的形式表示解算中應用的選項。如果你知道 UNIX/Linux 系統中的權限位表示方法,這一項應該很容易理解。否則,可以簡易地以數字的對應關系計算:
-
1 :earth;
-
2 :freq-dep;
-
4 :pole;
-
8 :ocean;
-
16:remove mean for pole tide;
-
32:atmosphere
對於你要在解算中應用的項,只需將其對應的數字的和作為配置即可。例如這里的 31:
31 = 1 + 2 + 4 + 8 + 16
因此數據處理中應用到的選項為:
-
earth;
-
freq-dep;
-
pole;
-
ocean;
-
remove mean for pole tide;
7、補充
文行至此,已經將 GAMIT/GLOBK 程序中所有的 .grid 文件都已做過簡要的介紹。這里以一個小技巧作為最后的補充。
或許你還發現,在 tables 文件夾內還存在一些與 .grid 同名的 .list 文件(除了 gpt.grid)。這些 .list 文件是從模型文件中提取出的所有 IGS 站點對應的參數值。如果你要處理的全都是 IGS 站的觀測數據,完全可以使用這些 .list 文件,只需在對應的配置項中將使用 .list 文件的選項設置為 “Y”,並將使用 .grid 文件的選項設置為 “N”。