chaoshu的个人博客分享 http://blog.sciencenet.cn/u/chaoshu

博文

GAMIT学习中遇到的问题

已有 37912 次阅读 2015-7-1 16:57 |个人分类:GAMIT/GLOBK|系统分类:科研笔记|关键词:学者

1、如果安装之后doy能用其他的命令(如sh_get_nav,makexp等)出现command not found 很有可能是环境变量没有配置好,gamit的程序命令都是在/com /gamit/bin /kf/bin 三个目录下的,环境变量一定要包含这三个目录,不然部分命令用不了。

2、在用在gamit分布解算的时候,用sh_rx2apr生成测站近似坐标的时候,出现如下错误:

Running svdiff...

At line 1336 of file svdiff.f

Fortran runtime error : Bad real number in item 1 of list input

生成的.apr文件中坐标都是0。

待解决!


第2个问题,没有找到很好的解决办法,检查了很多地方,都不知道究竟如何解决,感觉应该是安装过程的问题,我的配置是:

ubuntu12.04+GAMIT10.5+bash

在网上看过很多教程,都说可以在bash下安装gamit,我就是在bash下安装的,确实安装成功,但是总是会有各种小问题,解决不了啊!所以还是重装了,改成csh安装后,这个问题没有出现了。


3、有一次在解算一个站点W263时,过程没有大错误,很顺利的完成,但最后却没有这个站的结果。

进行分布解算的时候才发现原来在进行makex的运算时发生了错误。

具体是这样的:makex会生成k-file和x-file,而先生成k-file时,本来应该文件名为kW2634.193,却生成了

kw2634.193,这就出现了一个致命问题,get_rxfiles:在rinex目录下查找不到w263对应的o文件(即大写字母变成了小写字母,W2631930.14o不被认识了),后面的该站对应的xW2634.193文件也就没有生成,后面的批量解算也就自然没有该站的结果了。

解答:请教过一些同仁,了解到gamit 所有的测站文件、扩展名都应该是小写。

大写字母change小写字母: sh_casefold -dir d -files ./rinex/*


4、最近尝试了一下最新的版本gamit10.6,这个版本较以前有比较大的变化,说明一下几个问题,

1)安装的基本方法和之前的博文介绍差不多,但有几个注意地方,

修改Makefile.config里面的X11路径,去掉后面的/X11

详细参考

http://wenku.baidu.com/link?url=qrB8HuJ3O_X9L3MH3bsOjjQXVhgIN2eph73Zq0K_MIiFq6KEHQuc9DzFBPM4YAq1hi498eG3-ODt4cMUbP7Lw-pvxOSDskjpZBwuf3FrEKe

2)安装完之后,运行sh_get_orbits,只能下载igs的sp3文件,得不到g-file。原因是sh_sp3fit出现问题,根据错误提示,新版本的gamit的tables里面的svnav.dat文件版本太低。

很奇怪,我明明昨天才更新了tables里的所有文件,不可能是表文件太久的问题。但我还是又一次从这个地方ftp://garner.ucsd.edu/archive/garner/gamit/tables/

又更新了一次tables,结果还是出现同样的问题。

没办法,去官网查看最新的使用说明,看到Documentation下面又几个tables文件,就打开来看一下,发现

svnav.dat  Version  2.0 ,这才明白原来gamit10.6开始使用svnav.dat等文件的新版本,

旧的格式已经不能再用了,果断把这几个文件下载,替换了tables里面的文件,再运行,问题解决。


因此,这里提醒一下,以后这几个文件不能在从上面的ftp更新了,要直接从官网更新。


5、10.6新版本进行简单的实例运算,出现fixdrv如下错误

Running fixdrv
STATUS :150722:2144:30.0 FIXDRV/fixdrv: Started v.10.44 of 2015/6/17 21:10:00 (Linux)
STATUS :150722:2144:30.0 FIXDRV/fixdrv: New Clock-polynomial (I-) file being written--see fixdrv.out
WARNING:150722:2144:30.0 FIXDRV/lib/rsesfo: Session number is zero; setting to 1
WARNING:150722:2144:30.0 FIXDRV//lib/rstnfo: Station.info entry BJNM 2015 137  0  0  0  2015 137 23 59 30 ends early for session but may be ok for station
WARNING:150722:2144:30.0 FIXDRV//lib/rstnfo: Station.info entry CHAN 2015 137  0  0  0  2015 137 23 59 30 ends early for session but may be ok for station
WARNING:150722:2144:30.0 FIXDRV//lib/rstnfo: Station.info entry WUHN 2015 137  0  0  0  2015 137 23 59 30 ends early for session but may be ok for station
STATUS :150722:2144:30.0 FIXDRV/bmake: Setting numzen =  13 from zenint =  2.0 hr
FATAL  :150722:2144:30.0 FIXDRV/bmake: Ocean loading requested no list or grid file
Fatal errors occured in MAKEXP, MAKEX, MAKEK or FIXDRV

FATAL  :150722:2144:30.0 FIXDRV/bmake: Ocean loading requested no list or grid file

从上面的错误可以看出,是由于海潮改正模型otl.grid文件没有更新。

一般解算实例都是最简单的设置,能运行通过就行了。所以基本什么设置都没改。

ftp://garner.ucsd.edu/archive/garner/gamit/tables/

上面的海潮改正模型文件很大,有600多M,所以我也就在这个实例中没有更新,以往的老板本都是这样处理的,但是这次其他的所有过程都很正常的结束,可就是在最后的批处理fixdrv出现上面的错误。

查找了各种原因,后来发现,在新版本中sestbl.文件中

Use otl.grid = Y                 ; Ocean tidal loading grid file, GAMIT-format converted from OSO

设置成了默认为 Y,而旧版本中这一项默认是N;

所有,要不将Y改成N,要不就更新到最新的海潮改正模型otl_FES2004.grid。建议还是更新模型最好。

从这里也可以看出来,新版本在细节上和旧的版本有不少的改变,在使用时,一定要多加注意。

6、10.6运行fixdrv dyzfc5.211出现如下错误:

STATUS :150826:2217:25.0 FIXDRV/fixdrv: Started v.10.44 of 2015/6/17 21:10:00 (Linux)
STATUS :150826:2217:25.0 FIXDRV/fixdrv: Old I-file used, print file fixdrv.out not written
WARNING:150826:2217:25.0 FIXDRV/lib/rsesfo: Session number is zero; setting to 1
WARNING:150826:2217:25.0 FIXDRV//lib/rstnfo: Station.info entry HYDE 2015 211  0  0  0  2015 211 23 59 30 ends early for session but may be ok for station
WARNING:150826:2217:25.0 FIXDRV//lib/rstnfo: Station.info entry LHAZ 2015 211  0  0 30  2015 211 23 59 30 starts late for session but may be ok for station
WARNING:150826:2217:25.0 FIXDRV//lib/rstnfo: Station.info entry URUM 2015 211  0  0  0  2015 211 23 59 30 ends early for session but may be ok for station
WARNING:150826:2217:25.0 FIXDRV//lib/rstnfo: Station.info entry YZF1 2015 211  8 38  0  2015 211  9 30 30 starts late for session but may be ok for station
WARNING:150826:2217:25.0 FIXDRV/lib/check_y2k: Unexpected 2-digit year
WARNING:150826:2217:25.0 FIXDRV/lib/check_y2k: Unexpected 2-digit year
WARNING:150826:2217:25.0 FIXDRV/lib/check_y2k: Unexpected 2-digit year
FATAL  :150826:2217:25.0 FIXDRV/lib/julday: Unreasonable month:        0
STOP FATAL Error: Stop from report_stat

解決:上面同樣的數據,用gamit10.5處理沒有出現這樣的錯誤,初步感覺可能是最新版本的一個degud。又重新用gamit10.6解算了另一天的數據,沒有出現上面的問題。所以現在只能這麼認爲,這就是gamit10.6的degud了。

王东振-湖北局提出的解决办法:

可能有的人装的10.6没有这种问题。我在执行sh_sp3fit的时候输出的tigsf二进制文件比正常的要小很多,通过查看代码发现10.6在执行sh_sp3fit之前运行了sh_bcfit这个命令来得到tigsf,所以在执行makej之前执行一下sh_bcfit

sh_sp3fit得到的T文件不行,所以用sh_bcfit来得到T文件就行了

出现这种情况不是sh_bcfit可以解决的,而是由于没有运行两次sh_check_sess导致的,也可能是运行makex expt.makex.batch中最好写成makex expt.makex.batch 1 0

10.6   T文件的生成是在最后一步生成的,10.4是在sh_sp3fit这步生成的

把sh_sp3fit这条命令里面的-t删掉试一下

7、csh bzang5.bat 出現Geodetic height unreasonable: check p- and l-files錯誤

STATUS :151027:1104:35.0 MODEL/open: Site QHGE: Started MODEL version 10.45 2013/4/19 12:50:00 (Linux)
WARNING:151027:1104:35.0 MODEL/open: Site QHGE: Started MODEL version 10.45 2013/4/19 12:50:00 (Linux)
STATUS :151027:1104:35.0 MODEL/open:  Site rename File         : eq_rename
STATUS :151027:1104:35.0 MODEL/open:  Input Observation File   : xqhge5.226
STATUS :151027:1104:35.0 MODEL/open:  Output C-file            : /tmp/cqhge5.226.2796801000
STATUS :151027:1104:35.0 MODEL/open:  Ephemeris (T-) File      : tigsf5.226
STATUS :151027:1104:35.0 MODEL/open:  Loading/Met (U-) File : uzang5.226
STATUS :151027:1104:36.0 MODEL/setup: Yaw modelling is implemented
FATAL  :151027:1104:36.0 MODEL/setup: Geodetic height unreasonable: check p- and l-files


根據上面說明,很可能就是l-file出現問題,也就是先驗座標有問題,檢查發現:

l-fiel文件中

QHGE QHGE GPS    N37 40 51.62936 E 91 13 34.55720 6249761.4582 Ref. Epoch 2015.6200 QHGE_GPS

這裏看不出太大問題,所以在向上查找,生成這個文件的是Epoch 2015.6164: From file lfile.apr

打開lfile.apr文件

QHGE_GPS   -105852.89     4945092.12     3820258.81   0.0000 0.0000 0.0000 2015.62

打開o-file,qhge2260.15o

 -429152.7055  5141120.9335  3743045.1321                  APPROX POSITION XYZ

發現兩者差別太大了,這就是問題所在啊!

爲了說明問題:

sh_rx2apr -site qhge2260.15o

生成: QHGE_GPS   -429152.705500     5141120.933500     3743045.132100   0.0000 0.0000 0.0000 2015.62

sh_rx2apr -site qhge2260.15o -nav brdc2260.15n

生成:  QHGE_GPS   -105852.89     4945092.12     3820258.81   0.0000 0.0000 0.0000 2015.62

???爲什麼會出現這麼詭異的問題呢???


8、sh_gamit錯誤

Failure in sh_preproc. -- sh_gamit terminated

在批量自動解算時,遇到這個錯誤,查看官方說明發現問題原因:

sh_preproc: Fixed bug when no sites are to be excluded.  Herring 150805.

很明顯,是sites.defaults設置有問題,本來數據都準備好了,不需要在設置這個啦,我通常是全部註釋掉,但是現在10.5版本漏洞,要不換新版本,要不就採用默認設置。


9、gamit出现如下错误

Warning: ridiculously long PATH truncated

分布解算,发现sh_makexp就出现错误了,估计应该是brdc导航文件出现问题,使用sh_get_nav,

sh_get_nav -archive sopac -yr 2015 -doy 001 -ndays 1

Information extracted from ftp_info
########################################
ftpsite garner.ucsd.edu
ftplogin anonymous chaoshu@cc.
ftpdir /pub/rinex/YYYY/DDD
ftpcmd ftp -inv garner.ucsd.edu
wgetsite ftp://garner.ucsd.edu
wlogin anonymous chaoshu@cc.
########################################
Getting auto or brdc files for 2015 001 for 1 days from garner.ucsd.edu
Warning: ridiculously long PATH truncated
Warning: ridiculously long PATH truncated
Word too long.

如此,基本可以判断就是这个问题了。

解决办法:关闭此终端,重新打开一个,再次输入上面的命令,如果没有出现这个问题,那就解决了!这个问题的原因暂时不知。

10、POS8 POS8_GPS    N** **      NaN E  0  0  0.00000       0.0000

sh_rx2apr -site pos81310.15o 出现如上错误,无法获得近似坐标。

原因是

    2.11           OBSERVATION DATA    M (MIXED)           RINEX VERSION / TYPE
teqc  2015Jun23     chaoshu             20160115 03:02:44UTCPGM / RUN BY / DATE
Linux 2.4.21-27.ELsmp|Opteron|gcc|Linux x86_64|=+           COMMENT
teqc  2015Jun23     chaoshu             20160115 03:02:43UTCCOMMENT
teqc  2015Jun23     chaoshu             20160115 03:02:41UTCCOMMENT
teqc  2015Jun23     chaoshu             20160115 03:02:40UTCCOMMENT
teqc  2015Jun23     chaoshu             20160115 03:02:39UTCCOMMENT
teqc  2015Jun23     chaoshu             20160115 03:02:38UTCCOMMENT
teqc  2015Jun23     chaoshu             20160115 03:02:37UTCCOMMENT
teqc  2015Jun23     chaoshu             20160115 03:02:35UTCCOMMENT
teqc  2015Jun23     chaoshu             20160115 03:02:34UTCCOMMENT
teqc  2015Jun23     chaoshu             20160115 03:02:33UTCCOMMENT
teqc  2015Jun23     chaoshu             20160115 03:02:32UTCCOMMENT
teqc  2015Jun23     chaoshu             20160115 03:02:31UTCCOMMENT
teqc  2015Jun23     chaoshu             20160115 03:02:30UTCCOMMENT
teqc  2015Jun23     chaoshu             20160115 03:02:29UTCCOMMENT
teqc  2015Jun23     chaoshu             20160115 03:02:29UTCCOMMENT
teqc  2015Jun23     chaoshu             20160115 03:02:28UTCCOMMENT
teqc  2015Jun23     chaoshu             20160115 03:02:27UTCCOMMENT
teqc  2015Jun23     chaoshu             20160115 03:02:27UTCCOMMENT
teqc  2015Jun23     chaoshu             20160115 03:02:26UTCCOMMENT
teqc  2015Jun23     chaoshu             20160115 03:02:26UTCCOMMENT
teqc  2015Jun23     chaoshu             20160115 03:02:25UTCCOMMENT
teqc  2015Jun23     chaoshu             20160115 03:02:25UTCCOMMENT
teqc  2015Jun23     chaoshu             20160115 03:02:25UTCCOMMENT
teqc  2015Jun23     chaoshu             20160115 03:02:24UTCCOMMENT
teqc  2015Jun23     chaoshu             20160115 03:02:24UTCCOMMENT

  0.076      (antenna height)                              COMMENT
+28.55021313 (latitude)                                    COMMENT
+86.65622741 (longitude)                                   COMMENT
4382.949      (elevation)                                   COMMENT
BIT 2 OF LLI FLAGS DATA COLLECTED UNDER A/S CONDITION       COMMENT
pos8 (COGO code)                                            COMMENT
POS8                                                        MARKER NAME
POS8                                                        MARKER NUMBER
-Unknown-           -Unknown-                               OBSERVER / AGENCY
5440R49418          TRIMBLE NETR9       4.81                REC # / TYPE / VERS
21353166            TRM59800.00     SCIS                    ANT # / TYPE
  327253.3626  5601133.3713  3032303.0357                  APPROX POSITION XYZ
       0.0760        0.0000        0.0000                  ANTENNA: DELTA H/E/N
    1     1                                                WAVELENGTH FACT L1/2
    7    L1    L2    C1    P1    P2    S1    S2            # / TYPES OF OBSERV
   30.0000                                                 INTERVAL
   16                                                      LEAP SECONDS
SNR is mapped to RINEX snr flag value [0-9]                COMMENT
 L1 & L2: min(max(int(snr_dBHz/6), 0), 9)                  COMMENT
Forced Modulo Decimation to 30 seconds                      COMMENT


头文件中包含了一些不规则的内容,上面红色标记的。删除这些内容之后就可以了!


11、globk使用sh_exglk提取.org文件的itrf速度场文件,使用cvframe转行成欧亚速度场

cvframe fjxm.vel fjxm_oy.vel ITRF08 EURA

12、BCSUM/multiread:IOSTAT error decoding file AKE THE STABILIZATION MORE ROBUST OR MORE PRECISE ERROR 5010

IOSTAT error 5010 occurred decoding

AKE THE STABILIZATION MORE ROBUST OR MORE PRECISE

运行sh_globk_scatter时遇到如上错误。

sh_globk_scatter -f globk_rep.org
注意到,globk_rep.org,是glred命令生成的,而glorg_rep.cmd中有如下内容

* Controls for removing sites from the stabilization #   Vary these to make the stabilization more robust or more precise                                     stab_it 4 0.8 3.0

x stab_it 4 0.5 4.0  

感觉这和错误信息提示的内容有关联。

13、在解算几年的多个站点的数据时,发现存在某条基线在某一单日解q-file中存在错误:

Baseline vector (m ): AIRA        (Site 1) to DZZZ        (Site 8)

X   891556.61074 Y(E)  1056082.07664 Z  -718562.60680  L  1557727.43529

N  -709189.66984 E    -1250005.95802 U  -164385.39672  L  1557727.43529

且,数值比正常值偏差很大,所以造成没有计算精度,估计还是lfile.文件的概略坐标有问题。

在数据量一大之后,如何对某个异常的观测文件进行检测,并查找出来,这很麻烦呢

14、GAMIT自带的sh_rx2apr存在bug,

o文件头文件中的XYZ坐标(460603.7128,5635817.8132,2943450.2693),GAMIT自带的转换工具是27°29′55″,可是我用第三方转换软件是27°39′24″

用rtklib检验,后面的坐标是正确的,也就是说GAMIT的坐标转换是存在问题的,会出现这种莫名其妙的错误。现在还不知道这个问题的原因。

GAMIT/GLOBK群中睿控-彭小青,解释说:通过笑脸转过去的是WGS84的经纬度,gamit的是球坐标,所以就不是椭球的,当然就不一样了

13、在进行GLOBK进行合并多期数据时,发生如下错误信息:

WARNING:160907:1901:38.0 GLOBK/lib/xyz_to_geod: Convergence failure in computing geodetic coordinates


14、解算TRACK,出现如下错误:

STOP TRACK DISASTER: Too many observables

IOSTAT error   5010 occurred decoding

WARNING:161116:1332:59.0 TRACK/rxhead: IOSTAT error decoding file  S2    S5          # / TYPES OF OBSERV ERROR  5010

WARNING:161116:1332:59.0 TRACK/rxhead: IOSTAT error decoding file    30.0000 INTERVAL ERROR  5010

从上面的截图可以看出,rinex的观测文件中存在大量的空的数据;想办法把这些冗余的数据剔除。

选择观测值类型。

例如:从test.08o文件中提取单频观测值L1(不区分大小写),可运行以下命令:

teqc -O.obs L1 test.08o >test.L1.08o

例如:要提取L1,L2,P1,P2(不区分大小写,+可省略),可运行以下命令:

teqc -O.obs L1+L2+P1+P2 test.08o >test.L1L2P1P2.08o

15、TRACK解算遇到如下错误

**DISASTER** No overlapping data at site ml15

16、MAKEX/makex: # epochs ( 5760) gt maxepc { 2880

这是由于安装的时候MAXEPC 2880,应该重新编译,将2880修改成5760

17、Error opening file for COM ocean loading correction (Name otlcmc.dat)

解决办法,在gamit安装目录下,更新软件。

chmod +x install_updates

./install_updates

18、FATAL  :171020:0915:58.0 MAKEX/lib/rrxhed: # obs types =33 > maxobt = 27 (Name ./cksv1270.17o)

查找到maxbot的定义,如下:

/opt/gamit10.61/gamit/includes/makex.h:34:      PARAMETER   (MAXOBT=27)

将MAXOBT=27修改更大的数字,然后重新编译。

19、gamit10.61处理北斗数据,单步解算,前面处理都正常结束,最后一步出现如下错误

chaoshu@cc:~/project/test_gamit1061/056$ csh bchen7.bat

STATUS :171020:1154:49.0 ARC/aversn: Started ARC, Version 9.72 of 2013/3/26 10:30 (Linux)         Library ver. 10.96 of 2013/4/24 14:30:00:00 (Linux)
FATAL  :171020:1154:49.0 ARC/filopn: Error, old style svnav.dat, need a new svnav.dat file (Name svnav.dat)
STOP FATAL Error: Stop from report_stat
FATAL  :171020:1154:49.0 YAWTAB/orbits/yawtab: GAMIT.fatal exists: YAWTAB not executed
STOP FATAL Error: Stop from report_stat
FATAL  :171020:1154:49.0 GRDTAB/grdtab: GAMIT.fatal exists: GRDTAB not executed
STOP FATAL Error: Stop from report_stat
FATAL  :171020:1154:49.0 MODEL/model: GAMIT.fatal exists: MODEL not executed
STOP FATAL Error: Stop from report_stat
FATAL  :171020:1154:49.0 MODEL/model: GAMIT.fatal exists: MODEL not executed
STOP FATAL Error: Stop from report_stat
FATAL  :171020:1154:49.0 MODEL/model: GAMIT.fatal exists: MODEL not executed
STOP FATAL Error: Stop from report_stat
FATAL  :171020:1154:49.0 MODEL/model: GAMIT.fatal exists: MODEL not executed
STOP FATAL Error: Stop from report_stat
FATAL  :171020:1154:49.0 MODEL/model: GAMIT.fatal exists: MODEL not executed
STOP FATAL Error: Stop from report_stat
FATAL  :171020:1154:49.0 MODEL/model: GAMIT.fatal exists: MODEL not executed
STOP FATAL Error: Stop from report_stat
AUTCLN is running--see autcln.out for messages
STOP FATAL Error: Stop from report_stat

原因是之前的版本为10.5,后面安装10.61,添加环境变量时,只修改了.bashrc,忘了修改.cshrc;

所以,在最后一步 csh bchen7.bat 调用csh时候出现了问题。

20、sh_check_sess步骤时候出现如下错误:

makexp: ../../../src/libgfortran/io/read.c:1137: read_f: Assertion `exponent < 10000' failed.

Program received signal SIGABRT: Process abort signal.

Backtrace for this error:
#0  0xB76C5133
#1  0xB76C57D0
#2  0xB77CEC8B
#3  0xB77CECB0
#4  0xB7534606
#5  0xB7537A32
#6  0xB752D756
#7  0xB752D806
#8  0xB776F9E8
#9  0xB777299D
#10  0xB77709F2
#11  0x805158F in reade_
#12  0x804FFCF in checke_
#13  0x804BF34 in MAIN__ at makexp.f:?
makej: ../../../src/libgfortran/io/read.c:1137: read_f: Assertion `exponent < 10000' failed.

Program received signal SIGABRT: Process abort signal.

Backtrace for this error:
#0  0xB76B8133
#1  0xB76B87D0
#2  0xB77C1C8B
#3  0xB77C1CB0
#4  0xB7509606
#5  0xB750CA32
#6  0xB7502756
#7  0xB7502806
#8  0xB77629E8
#9  0xB776599D
#10  0xB77639F2
#11  0x806B87F in reade_
#12  0x804C000 in j_from_nav_
#13  0x804A29B in MAIN__ at makej.f:?
STOP FATAL Error: Stop from report_stat
ls: No match.
STOP FATAL Error: Stop from report_stat

sh_check_sess: Removing any PRN's from session.info that are missing from: gigsf7.265

missing session.info file QUIT

上面提示没有session.info文件,这个文件是在第一步sh_makexp中生成的,因此进行单步运行查找原因.

为了查看原因,单步运行,sh_makexp出现如下问题

STATUS :171107:1055:37.0 MAKEXP/makexp: Started  MAXEXP Ver. 9.85 2016/2/18 16:00:00 (Linux) Library Ver. 11.24 of 2017/6/7 10:00 (Linux)
WARNING:171107:1055:37.0 MAKEXP/lib/reade: Bad year on nav-file: G 6 2533141641 23 59 (Name brdc2640.17n)
makexp: ../../../src/libgfortran/io/read.c:1137: read_f: Assertion `exponent < 10000' failed.

Program received signal SIGABRT: Process abort signal.

Backtrace for this error:
#0  0xB763F133
#1  0xB763F7D0
#2  0xB7748C8B
#3  0xB7748CB0
#4  0xB74AE606
#5  0xB74B1A32
#6  0xB74A7756
#7  0xB74A7806
#8  0xB76E99E8
#9  0xB76EC99D
#10  0xB76EA9F2
#11  0x80515A7 in reade_
#12  0x804FFCF in checke_
#13  0x804BF34 in MAIN__ at makexp.f:?
Aborted (core dumped)

有上面基本可以断定是导航文件出现错误导致.

打开当天的brdc文件,发现最后几行出现错乱,删除掉乱行,从新计算,没有出现错误;问题解决

21、gamit解算出现如下情况,正常结束,但是最后的nrms超过了0.5,达到0.929,如何解决问题呢?

ls: No match.
STOP Normal finish of autcln
STOP Normal finish of autcln
STOP Normal finish of autcln
aopts: Subscript out of range.

Prefit nrms:  0.28385E+01    Postfit nrms: 0.92961E+00

查看结果文件q-file,发现Baseline vector前面有如下问题:

Adjustments larger than l-file tolerance: 0.30                                                                              
GCR LFTOL AIRA GEOC LONG     1.077                                                                                          
GCR LFTOL ARTU GEOC LAT      0.682                                                                                          
GCR LFTOL ARTU GEOC LONG     0.905                                                                                          
GCR LFTOL BAKO GEOC LONG     1.361                                                                                          
GCR LFTOL BJFS GEOC LONG     1.194                                                                                          
GCR LFTOL DAEJ GEOC LONG     1.123                                                                                          
GCR LFTOL DAEJ RADIUS        0.355                                                                                          
GCR LFTOL GUAM GEOC LONG     1.017                                                                                          
GCR LFTOL GUAM RADIUS        0.553                                                                                          
GCR LFTOL IISC GEOC LONG     1.223                                                                                                                                                                                  
Adjustments larger than twice the a priori constraint:                                                                      
GCR APTOL AIRA GEOC LAT     -0.247                                                                                          
GCR APTOL AIRA GEOC LONG     1.077                                                                                          
GCR APTOL AIRA RADIUS        0.298                                                                                          
GCR APTOL ARTU GEOC LAT      0.682                                                                                          
GCR APTOL ARTU GEOC LONG     0.905                                                                                          
GCR APTOL ARTU RADIUS       -0.184                                                                                          
GCR APTOL BAKO GEOC LAT      0.149                                                                                          
GCR APTOL BAKO GEOC LONG     1.361                                                                                          
GCR APTOL BAKO RADIUS        0.160

从上面看,基本可以确定应该是近似坐标L文件出现问题!但是查看了gamit官方文档,发现如下:

The Adjustments larger than summary will list sites for which the adjustment was greater than the tolerance you set in the sestbl. (nominally 30 cm) and those for which the adjustment was greater than the a priori constraint. Both of these tests are useful primarily in the preliminary solution (“p” Q-file), which we have not shown. For the first category, the L-file will be updated so that the final solution (shown here) will have adjustments within a linear range. For the second category, the a priori constraint will be relaxed for the final solution so that the solution is not strained (high nrms). These messages are also written into the sh_gamit_[ddd].summary file emailed to the user.

把sittbl.中的IGS站约束更改

0.025 0.025 0.050

改为

0.050 0.050 0.050

从新解算后如下:

Prefit nrms:  0.21040E+01   Postfit nrms: 0.52267E+00

说明修改sittbl.起作用,继续调整COORD.CONSTR,直到符合解算要求;

如果还是不行,删除所有数据,重新来过!

22、gamit基线解算正常结束,q-file中基线也是正常的,但是h-file里面,某一个站点Obs为 0

最后globk的org文件里面也没有这个测站的结果。(待解决

23、出现如下错误:

FATAL  :180106:2238:39.0 MODEL/setup: GPS frequencies on x-/c-file file not L1 and L2

原因:检查o-files发现其中为单频接收机,头文件中有如下语句:

1     0      WAVELENGTH FACT L1/2

gamit在10.6.10.5以后的版本不在支持L1_RECEIVE

SOLVE/read_bfl: L1_RECEIVER option no longer supported: use L1_ONLY

应该在sestbl.里面修改解算策略选择:

Choice of Observable = L1_ONLY

Quick-pre observable = L1_ONLY

并且必须在,必须在 autcln 命令文件中设置将注释取消 L1_only,并同时将注释上noL1only;

* Exclude L1-only (or bad RINEX files) to avoid problems: comment out if you want to process L1 data
x noL1only

* For autcln versions 3.30 or higher, the default is to exclude L1-only data; to include it use

L1only

但是这样设置以后,解算正常结束,并且能出结果,然后会出现如下信息:

ls: No match.

STOP Normal finish of autcln

STOP AUTCLN: No ref site and satellite found

aopts: Subscript out of range.

24、MAKEX/lib/readj: Satellite channel  1 (PRN= 0) not in J-file

更新tables,然后删除gfiles目录下的g文件,重新由sh_sp3fit生成g文件;

另外需要注意导航文件中一定也要包含相应的卫星


25、SOLVE/read_biases: Zero WL biases read from N-file, check processing sequence

相关解决方案坏apriori文件。
本教程中遇到的错误:在这种情况下,还请检查
表格文件夹中的符号链接regional.apr指向正确的位置。

更多信息:
编辑2000 / tables / process.defaults以从itrf05.apr更改为
regional.apr和'mailto'到您自己的电子邮件地址以便接收
sh_gamit摘要文件。


26、gamit处理2000周以后出现错误;

At line 304 of file timcon.f

Fortran runtime error: Expected INTEGER for item 7 in formatted transfer, got REAL

(a,i3,i6,4i4,f6.2,i4,f10.2)

解决办法:

利用gamit处理2000周之后的数据,会出现一个错误,是关于timcon.f的错误,这是一个时间转换的。编软件的时候把文件里的最大处理时间设置为了2000,需要把2000改成5000,然后重新编译一下即可。


27、在csh bbeid8.bat 处理出现如下错误;

Zero WL biases read from N-file, check processing sequence

这是因为导航文件的问题,从新下载文件;(最好下载单一系统的问题,混合系统的文件,很可能会存在问题)。


28、FATAL  :181102:1046:55.0 MODEL/get_dcb2: DCB value for SVN  54 not found on dcb.dat

更新dcb.dat.gps,并重新链接此文件到项目中的tables下dcb.dat文件!



https://m.sciencenet.cn/blog-858128-902021.html

上一篇:gamit的station.info中的天线高量取方式介绍
下一篇:[转载]利用gamit结果做基线重复率----预处理命令

0

该博文允许注册用户评论 请点击登录 评论 (1 个评论)

数据加载中...
扫一扫,分享此博文

Archiver|手机版|科学网 ( 京ICP备07017567号-12 )

GMT+8, 2024-5-4 09:45

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部