gamit中用到的GRID文件

在使用 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/

發表評論
所有評論
還沒有人評論,想成為第一個評論的人麼? 請在上方評論欄輸入並且點擊發布.
相關文章