关于Bernese的海潮负荷BLQ文件

Bernese的海潮负荷文件BLQ生成比较复杂，需要向指定互联网填写站位置信息，然后结果通过邮件形式返回。
主要流程如下：

1. 登录免费的海潮负载提供的网站
   http://holt.oso.chalmers.se/loading
2. 选择海潮模式（例如FES2004）
3. 选择“垂直和水平位移”（vertical and horizontal displacements）
4. 其他默认即可。注意：坐标系原点在固体地球中心。
5. 填写站位置。格式如下：

Name of station_________|   |Longitude (deg)    | Latitude (deg)    | Height (m)       OR
Name of station_________|   X (m)   Y (m)   Z (m)

这里站名的字符数是24个，第26个字符就是坐标位置了。
可以从Bernese的crd坐标文件复制需要获得海潮负荷的站坐标信息行。但是需要把第一位的X坐标向后退2格，也就是如下：

Name of station_________|   |Longitude (deg)    | Latitude (deg)    | Height (m)       OR
Name of station_________|   X (m)   Y (m)   Z (m)

339  ZIMM 14001M004      4331296.99513   567555.96734  4633133.99316    IGS14
下面这个是从crd文件中复制出来的。    
339  ZIMM 14001M004    4331296.99513   567555.96734  4633133.99316    IGS14     

其他不需要修改。一次最多可以填入100个站。

6. 填写邮箱，然后submit提交。等待了，他会每分钟更新提示，这个队列里就是提交的邮箱前缀。

BLQ等待

7. 等待一定时间，在邮件中就发现结果.
   [email protected]（发送邮件地址）
   邮件内容如下：

$$ Ocean loading displacement
$$
$$ Calculated on holt using olfg/olmpp of H.-G. Scherneck
$$
$$ Greens function: mc00egbc
$$ 
$$ COLUMN ORDER:  M2  S2  N2  K2  K1  O1  P1  Q1  MF  MM SSA
$$
$$ ROW ORDER:
$$ AMPLITUDES (m)
$$   RADIAL
$$   TANGENTL    EW
$$   TANGENTL    NS
$$ PHASES (degrees)
$$   RADIAL
$$   TANGENTL    EW
$$   TANGENTL    NS
$$
$$ Displacement is defined positive in Upwards, South and West direction.
$$ The phase lag is relative to Greenwich and lags positive. The 
$$ Gutenberg-Bullen Greens function is used. In the ocean tide model the 
$$ deficit of tidal water mass has been corrected by subtracting a uniform 
$$ layer of water with a certain phase lag globally.
$$
$$ Complete  : No interpolation of ocean model was necessary
$$ _PP       : Ocean model has been interpolated near the station 
$$                         (PP = Post-Processing)
$$ 
$$ Ocean tide model: FES2004
$$
$$ END HEADER
$$
  339  ZIMM 14001M004
$$ Complete FES2004
$$ Computed by OLFG, H.-G. Scherneck, Onsala Space Observatory 2019-Dec-27
$$ 339  ZIMM 14001M004,       RADI TANG  lon/lat:    7.4653   46.8771  956.336
  .00677 .00230 .00146 .00058 .00221 .00091 .00073 .00005 .00048 .00027 .00022
  .00273 .00075 .00061 .00019 .00036 .00025 .00012 .00006 .00004 .00002 .00003
  .00059 .00025 .00010 .00006 .00025 .00004 .00008 .00003 .00001 .00000 .00000
   -72.5  -41.7  -91.6  -52.0  -58.9  -97.2  -59.7 -159.7    7.5    5.6    1.1
    84.5  117.1   62.5  107.6   97.5   21.5   95.4  -40.2 -177.2 -169.8 -177.7
   -25.7    8.4  -28.0   -1.9   41.5  -30.8   40.5  165.6 -114.6 -175.9  -38.7
$$ END TABLE
Errors:
Warnings:

8. 把站内容复制粘贴到Bernese的BLQ文件就可以了。注意，那个序号（339）需要删除，复制粘贴的内容如下：

$$
  ZIMM 14001M004
$$ Complete FES2004
$$ Computed by OLFG, H.-G. Scherneck, Onsala Space Observatory 2019-Dec-27
$$ 339  ZIMM 14001M004,       RADI TANG  lon/lat:    7.4653   46.8771  956.336
  .00677 .00230 .00146 .00058 .00221 .00091 .00073 .00005 .00048 .00027 .00022
  .00273 .00075 .00061 .00019 .00036 .00025 .00012 .00006 .00004 .00002 .00003
  .00059 .00025 .00010 .00006 .00025 .00004 .00008 .00003 .00001 .00000 .00000
   -72.5  -41.7  -91.6  -52.0  -58.9  -97.2  -59.7 -159.7    7.5    5.6    1.1
    84.5  117.1   62.5  107.6   97.5   21.5   95.4  -40.2 -177.2 -169.8 -177.7
   -25.7    8.4  -28.0   -1.9   41.5  -30.8   40.5  165.6 -114.6 -175.9  -38.7
$$ END TABLE