chaoshu的個人博客分享 http://www.ueservicedoffices.com/u/chaoshu

博文

GAMIT學習中遇到的問題

已有 12770 次閱讀 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的坐標轉換是存在問題的,會出現這種莫名其妙的錯誤,F在還不知道這個問題的原因。

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文件!



http://www.ueservicedoffices.com/blog-858128-902021.html

上一篇:gamit的station.info中的天線高量取方式介紹
下一篇:[轉載]利用gamit結果做基線重復率----預處理命令

0

該博文允許注冊用戶評論 請點擊登錄 評論 (1 個評論)

數據加載中...

Archiver|手機版|科學網 ( 京ICP備14006957 )

GMT+8, 2019-5-26 01:21

Powered by ScienceNet.cn

Copyright © 2007- 中國科學報社

返回頂部
时时彩平台