Bernese提供了生成該大氣負(fù)荷文件的工具，位于：
Menu >Service >Coordinate tools >Extract atmospheric tidal loading coefficients
需要注意的是掸宛，輸入文件包括2個(gè)：
第一個(gè)：s1_s2_def_ce.dat，已經(jīng)在 "${X}/GEN" 目錄中了。
第二個(gè)：站信息表刑桑，也就是需要指明那些站計(jì)算大氣負(fù)荷失息，當(dāng)然這個(gè)程序會(huì)根據(jù)站名自動(dòng)去尋找該站的經(jīng)緯度譬淳。
下圖是S1和S2大氣潮汐的徑向表面位移
https://geophy.uni.lu/atmosphere/tide-loading-calculator/
以下翻譯自官方文檔：
================================================== ===============================================
Bernese GNSS軟件，版本5.2盹兢，最新版本：2012年10月3日
安裝S1 / S2大氣潮汐荷載網(wǎng)格文件
S1 / S2位移網(wǎng)格的來源
在全球地球物理流體中心（GGFC）的網(wǎng)站上
http://geophy.uni.lu/ggfc-atmosphere/tide-loading-calculator.html
大氣潮汐負(fù)荷位移網(wǎng)格
-CE：固體地球的中心
（在這種情況下邻梆，也可以用作“ CF：圖形中心”）
-CM：重心
提供。
根據(jù)Ray和Ponte模型[2003]給出了更正
加入IERS 2010大會(huì)绎秒。
文件“s1_s2_def_ce.dat”與Bernese GNSS一起分發(fā)
“${X}/GEN”目錄中的軟件浦妄。
重要的提示：
如果您還打算針對(duì)大氣非潮汐載荷進(jìn)行校正，則可以
使用一致的潮汐荷載貢獻(xiàn)必不可少见芹。格式
預(yù)期網(wǎng)格文件的內(nèi)容在本文檔的最后一部分中進(jìn)行了描述
文件剂娄。
置換網(wǎng)格的實(shí)現(xiàn)
S1/S2負(fù)載系數(shù)存儲(chǔ)在伯爾尼戰(zhàn)役。
使用以下程序GRDS1S2
菜單>服務(wù)>坐標(biāo)工具>提取大氣潮汐負(fù)荷系數(shù)
為給定的一組站點(diǎn)坐標(biāo)生成“.ATL”文件
伯爾尼格式的坐標(biāo)文件玄呛。
程序GRDS1S2.f是程序grdintrp.f的派生自
上述網(wǎng)站阅懦。它已適應(yīng)環(huán)境
的GNSS GNSS軟件可以按預(yù)期的方式讀寫文件-
墊子（由盧森堡大學(xué)的Tonie van Dam提供）。
使用替代的大氣潮汐網(wǎng)格
在文件的開頭徘铝，可以給出數(shù)據(jù)的簡短描述
復(fù)制到生成的伯爾尼語格式的“ .ATL”文件中耳胎。這些
該行的開頭用“ $$”表示注釋。
數(shù)據(jù)部分由14列帶實(shí)數(shù)的列組成（無列
寬度是預(yù)定義的惕它，但必須至少用一個(gè)空格分隔
字符）：
-1/2：經(jīng)度怕午，緯度
-3/4：徑向上S1的余弦和正弦項(xiàng)
-5/6：徑向上S2的余弦和正弦項(xiàng)
-7/8：北向S1的余弦和正弦項(xiàng)
-9/10：北方向S2的余弦和正弦項(xiàng)
-11/12：S1向東的余弦和正弦項(xiàng)
-13/14：S2向東的余弦和正弦項(xiàng)
如果您引入的分辨率高于1度的網(wǎng)格
必須增加MXLON和MXLAT定義的數(shù)組
$FG/GRDS1S2.f（或％FG：\GRDS1S2.f）。
==================================================
英文文檔：
Bernese GNSS Software, Version 5.2 Last mod.: 03-Oct-2012
Installation of the S1/S2 atmospheric tidal loading grid files
Source of the S1/S2 displacement grids
On the web site of the Global Geophysical Fluid Center (GGFC)
http://geophy.uni.lu/ggfc-atmosphere/tide-loading-calculator.html
the atmospheric tidal loading displacement grids refering either to
- CE: center of the solid Earth
(can be used in this context also as "CF: center of figure") - CM: center of the mass
are provided.
The corrections are given for the Ray and Ponte model [2003] according
to the IERS 2010 conventions.
The file "s1_s2_def_ce.dat" is distributed with the Bernese GNSS
Software in the "${X}/GEN" directory.
IMPORTANT NOTE:
If you plan to correct also for atmospheric non-tidal loading, it is
essential to use a consistent tidal loading contribution. The format
of the expected grid file is described in the last section of this
file.
Implementation of the displacement grids
S1/S2 loading coefficients are stored in the STA directory of a
Bernese campaign.
Use the program GRDS1S2 from
Menu >Service >Coordinate tools >Extract atmospheric tidal loading
coefficients
to generate the ".ATL" file for a set of station coordinates given in
a Bernese formatted coordinate file.
The program GRDS1S2.f is a derivative of the program grdintrp.f from
the above mentioned web-site. It has been adapted for the environment
of Bernese GNSS Software to read and write files in the expected for-
mats (courtesy to Tonie van Dam, University of Luxembourg).
Using alternative atmospheric tidal loading grids
At beginning of the file a short description of the data may be given
that is copied into the resulting Bernese formatted ".ATL" file. These
comments are indicated by "$$ " at the beginning of the line.
The data part consists of 14 columns with real numbers (no column
width is predefined but they must be seperated by at least one blank
character):
- 1/ 2: longitude, latitude
- 3/ 4: cosine and sine terms for S1 in the radial direction
- 5/ 6: cosine and sine terms for S2 in the radial direction
- 7/ 8: cosine and sine terms for S1 in the north direction
- 9/10: cosine and sine terms for S2 in the north direction
- 11/12: cosine and sine terms for S1 in the east direction
- 13/14: cosine and sine terms for S2 in the east direction
If you introduce a grid with a higher resolution than 1-degree you
have to increase the arrays defined by MXLON and MXLAT in
$FG/GRDS1S2.f (or %FG:\GRDS1S2.f).
======================================================================
