在使用 GAMIT/GLOBK 程序進行數據處理時,你肯定已經注意到了其 tables 文件夾中存在許多 .grid 文件。如果你檢查過 GAMIT/GLOBK 安裝目錄的 tables 文件夾,你會發現更多這樣的文件。它們多數體積都很大,一般爲數百 MB。這樣的文件稱爲格網模型文件,一般在 MIT 的 FTP 服務器上提供下載。
本文將試圖逐個地介紹其中的格網模型文件,可以算是介紹其格網模型文件的一份文檔。
atl.grid
全球的大氣潮格網模型文件,由澳大利亞國家大學的 Paul Tregoning 創建。該文件包含在 GAMIT/GLOBK 程序的安裝包中且不需要更新。也就是說,不管你處理哪一年的觀測數據,都可以引入這個文件。
要在數據處理中應用大氣潮格網模型,你首先需要保證在進行數據處理的工程文件夾的 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
atml.grid
全球的無潮汐大氣負載參數格網模型文件,由澳大利亞國家大學的 Paul Tregoning 根據盧森堡大學的 Tonie van Dam 提供的數據創建。該模型一般每年對應一個文件,供 GAMIT/GLOBK 使用的文件命名方式一般爲:atmdisp_cm.。如用於 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
gpt.grid
全球氣壓和溫度模型文件,該模型由 Boehm 和 Schuh 提出。在估計對流層參數時,GAMIT/GLOBK 從該模型中內插獲取觀測所在地區的氣壓和溫度。該模型一般不會更新,並且已經被包含在 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/GLOBK 就將從該模型中獲取氣壓和溫度參數。但內插出的精度畢竟不如實測數據,因此如果你有氣象觀測數據的話,最好還是將其引入。所以我推薦如下的設置:
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/GLOBK 將首先檢查有沒有對應的 RINEX 格式氣象數據,如果沒有的話再檢查有沒有對應的 U-文件,最後纔會嘗試使用全球氣壓與溫度模型提供的數據。
map.grid
全球的大氣映射函數模型文件,該文件由澳大利亞國家大學的 Paul Tregoning 根據維也納科技大學的 Johannes Boehm 提供的數據創建。GAMIT/GLOBK 目前支持 GMF、NMFW 和 VMF1 三種映射函數,但只有 VMF1 需要引入模型文件。該模型逐年更新,每年生成一個文件。你可以根據需要在 MIT 的 FTP 服務器上下載,其命名方式爲:vmf1grd.。
要在數據處理時應用 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
met.grid
雖然各版本的 GAMIT/GLOBK 一直保留了這個模型文件鏈接,但該模型並沒有什麼用。有 sestbl. 中對該項配置的說明爲證:
Use met.list = N ; not yet supported
Use met.grid = N ; not yet supported
otl.grid
全球的海潮模型文件,GAMIT/GLOBK 現在支持多種海潮模型。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;
補充
文行至此,已經將 GAMIT/GLOBK 程序中所有的 .grid 文件都已做過簡要的介紹。這裏以一個小技巧作爲最後的補充:或許你還發現,在 tables 文件夾內還存在一些與 .grid 文件同名的 .list 文件(除了 gpt.grid)。這是爲了方便處理 IGS 跟蹤站的觀測數據,.list 文件是從模型文件中提取出的所有 IGS 站點對應的參數值。如果你要處理的全都是 IGS 站的觀測數據,完全可以使用這些 .list 文件,只需在對應的配置項中將使用 .list 文件的選項設置爲 “Y”,並將使用 .grid 文件的選項設置爲 “N”。
轉自GNSS helper:http://www.gnss.help/2016/10/02/grids-in-gamit/