#gamit 错误提示:及问题解决
#gamit probles summary by
[email protected]
#遇到问题,请先搜索本博客或者互联网,如果没有答案,请把问题详细的说明(gamit版本、更新包日期、详细fatal信息)发送到
[email protected]。
#我们需要的不是重复性劳动!!!
#update 2013-07-19
>>符号后面是问题处理方法
1.我安装的是GAMIT10.40,GMT4.5版本;使用sh_glred时会有警告"GMT warning:old GMT3 .gmtdefaults file found. May not be fully compatible with GMT 4.",我知道GMT4的默认设置文件是.gmtfaults4,我查看了gamit/com好几个脚本文件都提到了“.gmtdefaults ” ,比如sh_pltvel sh_plot_block等,不知道应该怎么改,这个地方应该如何把 .gmtdefaults 改成 .gmtdefaults4 呢?
>>可能之前因为你安装过gmt3,根据网上的方法,删除所有其它版本的.gmt*文件,在我的电脑上未出现该问题。
2.使用sh_gamti时我加入了参数-pres ELEV ,但是在!文件中出现了“Calling sh_make_sky_gifs with elevation plots
convert: Command not found.”,我查看了gamit/com确实有脚本sh_make_sky_gifs;!文件中还有 pshistogram: Warning: -C for centering is obsolete, please use -F 导致绘制elevation plots 和 phase sky_plots都失败,这个地方应该在哪改呢?或者说GAMIT和GLOBK调用GMT的脚本或程序在哪里?
>>convert 是个程序,是把ps图转为gif图的程序,需要你自己确定系统是否存在,否则自行安装。
3.GLOBK定义参考框架的问题,我看资料感觉至少有两种方法,一个是在globk_comb.cmd中对站点的坐标和速率加约束值;另一个是在glorg_comb.cmd中在stab_site中添加站点,那这两种方法都可用吗?有什么联系吗?
>>首先要对globk和glorg做个小区别,日平差用glorg,而几年的数据合并平差用globk
4.我今天又尝试了一下,选择的数据是wuhn、kunm、bjfs三个站点2000年第231天的数据,运行sh_gamit之后,出现了一个提示:aopts: Subscript out of range。但是结果都出来了,不知道这是为什么。然后我输入sh_gamit_baseline -o o*,系统就会提示昨天的那个问题x min=x max,如果后面加上基线名sh_gamit_baseline -o o* -b wuhn_bjfs,系统就会提示start:Subscript out of range。 date:Subscript out of range。 offset:Subscript out of range。不知道您在运行sh_gamit和sh_gamit_baseline时有没有遇到过这样的问题,求解答。
>>输入数据存在问题,x的最大、最小值相同?是否还有其他信息,另外该命令,需要多天的o文件存在,而不能只有一天数据,否则不可能统计。
5.FATAL :110427:2024:55.0 GRDTAB/grdtab: No site match in otl.list and otl.grid missing
FATAL :110427:2024:55.0 MODEL/model: GAMIT.fatal exists: MODEL not executed
FATAL :110427:2024:55.0 MODEL/model: GAMIT.fatal exists: MODEL not executed
>>出现这个原因,是因为没有连接700多M的那个otl_FES2004.grid文件,你可以在tables中查看一下这个文件是否是链接已坏,如果是的话,到everest.mit.edu,匿名ftp下载(在pub/GRIDS目录下),下载后,做连接。
6.运行sh_gamit后,总是这样的警告
WARNING:090826:1152:15.0 MSTINF2/lib/read_rcvant: Receiver name -------------------- not found in rcvant.dat
WARNING:090826:1152:15.0 MSTINF2/lib/read_rcvant: Receiver name -------------------- not found in rcvant.dat
最后:
FATAL :090825:1202: 1.0 FIXDRV/dcheck: Only one or no existing X-files
Fatal errors occured in MAKEXP, MAKEX, MAKEK or FIXDRV
FATAL :090825:1202: 1.0 MAKEX/set_dcb_flag: PCN-code missing for receiver type in rcvant.dat, cannot set C1/P2 correction flag
FATAL :090825:1202: 1.0 FIXDRV/dcheck: Only one or no existing X-files
查询rcvant.dat与guess_rcvant.dat后,在rcvant.dat里面都能找到,guess_rcvant.dat找不到,修改添加进去后,还是一样的错误。
例如:: Receiver code 5700 not found in rcvant.dat
>>手动制作station.info文件。如果rcvant.dat与guess_rcvant.dat中根本没有这个接收机或天线型号,需要手工修改添加这个型号(修改为近似的仪器型号,或者告诉king,让其添加该仪器型号),修改rcvant.dat和hi.dat文件,特别注意水平偏差和垂直偏差的设置。另外也可能是guess_rcvant.dat中没有该天线型号,但在rcvant.dat中确存在,这就需要你参照rcvant.dat文件中的天线、接收机等型号信息,手动添加到guess_rcvant.dat中。或者参照rcvant.dat中信息手动制作station.info文件,并采用不更新station.info设置。
7.tubi.makex.infor里面的提示是:
No RINEX files found for d062 from ./d062201?.09o
自己观测的测站数据不能参与计算(生成不了X文件),而IGS站确可以。
>>查看o文件等文件名,要小写,而不是大写。
8.我从自己测站中提取,station 信息时失败! 提示:Error reading day-of-year from RINEX filename dxbz.10o
>>RINEX文件,不但内容有一套严格要求,就连文件名也有要求哦。站名+doy+n.yr+类型
9.cannnot get antcod from RINEX entry
>>rinex中天线类型不正确,或者没有,如果自己确切知道天线类型,可以自己手动制作station.info文件,然后在处理中,不更新该文件。或者自己手动在rinex中添加天线类型。
10.FATAL :110430:0413:38.0 GRDTAB/get_otl_grid: Error reading grid file ERROR 3
FATAL :110430:0413:38.0 MODEL/model: GAMIT.fatal exists: MODEL not executed
>>这个问题是otl.grid海潮文件没有配置好。查看链接是否完好,否则重新制作。
11.我有Ubuntu10.04装的gamit10.34,解算数据都可以,但是方向键出现乱码,不能切换到前面使用过的命令,请问是怎么回事?
>>这个问题和gamit一点关系都没有了。出现这种情况,主要是终端缺少插件(安装系统过程中,出现错误,而少安装了东西),如果安装过oracle在linux上,也会遇到这个问题,只能安装插件了。
12.用自己作的station.info文件,怎么设置他才不自动更新呢?
>>在 sites.defaults中设置,不更新..例如下行:
all_sites expt xstinfo
所有测站都不更新。
13.FATAL :110503:1845:18.0 MAKEXP/lib/read_rcvant: Antenna name IAX1202GG not found in rcvant.dat
STOP FATAL Error: Stop from report_stat
MAKEXP fatal ocurred.
>>首先确定gamit提取的天线型号“IAX1202GG”,是你的天线型号,因为有时可能出现提取错误,但这里就是这种情况。LEIAX1202GG是正确的天线,还有就是自己完全手动制作station.info文件时,注意格式等,正确,这样才能被gamit正确识别,并在rcvant.dat中找到。另外,天线型号问题确实很常见,通常做法,在guess_rcvant.dat中查看你的天线是否确实存在,然后再rcvant.dat中也确定存在,如果都存在的话,那就是你的station.info文件存在问题。之后采用sh_upd_stnfo手动生成station.info文件。
14.请问查询天线参数据网址是多少?
>>GPS天线参数提供网站:http://www.ngs.noaa.gov/ANTCAL/#,当然这里基本全是老外的天线。
15.请问为什么igs精密星历文件下载不了呢?
>>这个问题基本不可能出现,因为提供igs精密星历的网站有好几个,不能都同时瘫痪,可在CDDIS(US-MD)、IGS CB(US-CA)、SOPAC(US-CA)、IGN(FR)、KASI (KOREA)机构网站下载。具体网址,这里有:http://igscb.jpl.nasa.gov/components/prods.html
16.gamit的解算中岁差模型可不可一改成IAU2006,或者IAU2000,我解出来的只是IAU76的
>>在sh_sp3fit中可以设置,如下面得命令,其中-i即为惯性参考系统,默认是J2000 IAU76
sh_sp3fit -f
-d -o -i -r -t
17.观测了24小时的数据,横跨两天,想请问一下,用gamit处理,而且是手动处理,怎样进行呢?
>>把精密星历也连在一起就可以,把ORBIT前面那个96改为192。广播星历也类似。
18.遇到了这样的错误 NGSTOT/LIB/NUTRED:JD=2455680。out of range of Nutation Table
>>这是章动表的时间段,不能满足你工程的需要。可以下载最新的章动表(nutabl.IAU00.2011)。(garner.ucsd.edu/pub/gamit/tables)
19.运行gamit,每次输入sh_get_nav命令,都显示command not found,安装好后,输入doy命令,结果是正常的。请问这是为什么呢?
>>这很可能是没有把那个com目录加进环境变量路径中去,gamit最终安装完成后,没有提示那个com目录,只有两个bin目录,所以可以找到kf目录下的doy目录,但不能找到com目录下的sh_get_nav等脚本命令。解决办法是,把com目录加进去(我用的是bash),如下:
export PATH=/home/zhao/program/gamit/gamit/bin:/home/zhao/program/gamit/kf/bin:/home/zhao/program/gamit/com:$PATH
把相应的目录前缀修改下即可。
20.FATAL :110524:0934:33.0 NGSTOT/lib/ut1red: JD= 2455683 out of range of ut1 table, JDT1= 2448622 JDT2= 2455469?
>>从这个错误提示,我们清楚看到,utl表的数据不能满足项目需要,去http://sopac.ucsd.edu/processing/gamit/下载最新的ut1.usno表。
21.如何用gamit处理南方、中海达等国产GPS接收机数据,另外,如果仪器高已经转为相位中心,该使用哪种量高方式?
>>如何处理国产接收机数据,只能采用替代方法了(参见第22条),但精度、误差变换情况不了解。如果仪器高以用转换软件转到相位中心,则在sestbl.表中设置Antenna Model=NONE,不采用相位中心变化,然后再hi.dat文件中设置天线高量测参考点,当然最好采用直高。
22.遇到未识别的天线LEIGS15,在rcvant.dat中可以找到,但是在guess_rcvant.dat中和antmod.dat找不到。我处理时是将观测文件中天线高归化至相位中心,用别的已知天线类型来代替其它未知天线类型,这样处理是否妥当?
>>在处理时,遇到LEIGS15天线未识别,可更新gamit,更新incremental_updates.110531.tar.gz后,在rcvant.dat、antmod.dat中可以找到LEIGS15天线,但guess_rcvant.dat中还是没有,需要自己手动添加下(能和rcvant.dat对应起来即可)。如果不进行天线相位中心改正,可在sestbl.中设置“Antenna Model = NONE”。在别无它法时,采取这种方法当然可以,但gamit的天线改正模型显然精度更高,相对仅进行简单的相位中心改正。
23.在gamit解算中,要让最后的l坐标文件中的坐标 与初始坐标文件中的坐标 相比 变或不变,是不是只需在sittbl.文件中设置就可以了?试过后,但是我发现结果不受我控制,我在sittbl.文件中变换 fix 这一选项的 YYY 或 NNN,但是结果并不是按照我所设置的来的~
>>
24.在电脑上安装了gamit10.35,(服务器操作系统CentOS 5.5 final)配置好安装路径后,在终端输入doy,给出如下的提示:
[zhenghui@node10 2008100]$ doy
**ERROR** Opening help file. Check that enviroment variable HELP_DIR points to
IOSTAT error 1 occurred opening doy.hlp in doy
WARNING:110616:1659:15.0 DOY/doy: IOSTAT error opening file doy.hlp ERROR 1
***TODAY*** IS:
Date 2011/06/16 16:59 hrs, DOY 167 JD 2455729.2078 MJD 55728.7078
GPS Week 1640 Day of week 4, GPS Seconds 406755 Day of Week Thu
Decimal Year 2011.456734
[zhenghui@node10 2008100]$
错误提示是系统HELP_DIR的问题,但是我运行其他的命令不带参数,那么可以正常运行。
我不知道问题出在哪里了,并且在网上也没有找到相关的解决办法。
>>HELP_DIR环境变量是globk的命令需要,采用该环境变量来给出帮助信息,当然,如果不配置该变量,完全可以正常使用gamit。如果添加该环境变量的话,可以在.bashrc或者.cshrc中添加。例如在.cshrc文件中,添加
setenv HELP_DIR /usr/local/gamit/help/
后面的路径当然要修改为你电脑中的gamit目录了。
25.我单步处理当天的GPS数据,到最后执行批处理文件的时候给出如下提示:
[zhenghui@node10 2008100]$ csh bpgga8.bat
STATUS :110616:1703:12.0 YAWTAB/orbits/yawtab: Epoch 500
STATUS :110616:1703:12.0 YAWTAB/orbits/yawtab: Epoch 1000
....省略....
STATUS :110616:1703:13.0 GRDTAB/grdtab: Normal stop in GRDTAB - created upgga8.100
STATUS :110616:1703:13.0 MODEL/open: Site CHAN: Started MODEL version 10.29 2009/7/15 11:15:00 (Linux)
WARNING:110616:1703:13.0 MODEL/open: Site CHAN: Started MODEL version 10.29 2009/7/15 11:15:00 (Linux)
STATUS :110616:1703:13.0 MODEL/open: Site rename File : eq_rename
STATUS :110616:1703:13.0 MODEL/open: Input Observation File : xchan8.100
STATUS :110616:1703:13.0 MODEL/open: Output C-file : /tmp/cchan8.100
STATUS :110616:1703:13.0 MODEL/open: Ephemeris (T-) File : tiggf8.100
model/open ionfil,magfield ierrNONE 0
STATUS :110616:1703:13.0 MODEL/open: Loading/Met (U-) File : upgga8.100
FATAL :110616:1703:13.0 MODEL/lib/read_antex_head: ANTEX version > 1.3
STOP FATAL Error: Stop from report_stat
FATAL :110616:1703:13.0 MODEL/model: GAMIT.fatal exists: MODEL not executed
STOP FATAL Error: Stop from report_stat
FATAL :110616:1703:13.0 MODEL/model: GAMIT.fatal exists: MODEL not executed
fatal出现,不能正常执行,网上说要删除执行目录下的GAMIT*,然后把rtx文件中的COMMENT删除几行了就可以了,但是我不知道运行csh bpgga8.bat命令后调用的是哪一个rtx。
>>是10.35版本问题了。如果是gamit10.4则没有这个问题。但如果你只能用gamit10.35,那么可以将read_antex_head.f中的 1.3改成1.4 重编译gamit
26.在服务器上安装了gamit10.4,我查看了一下网站上下载gamit10.4目录里面的readme,那个文档说明好像说的是Linux下的gfortran(gcc4.2或4.4,4.3可能有bugs),推荐使用g77(gcc<4.0)进行安装。我在笔记本上安装了ubuntu10.04,用apt-get install gfortran 安装的gfortran为4.4.3,我不知道怎么安装上g77,用g77对gamit源码进行编译,而是用默认的gfortran4.4.3进行编译安装的,不知道后面会不会有问题。
>>个人认为,目前linux各版本自带的gfortran可以正常编译gamit,那么就不存在问题。如果出现编译不通过,则需要考虑更改编译器版本,或则更改编译器,毕竟linux下的编译器还是毕竟丰富的。
27.我安装好gamit10.4以后,仿照你的视频例子,进行了批处理,虽然最后在每个单天目录下能生成day.summary,但是终端还是给出如下的提示:
[zhenghui@node10 batchproc]$ ls
brdc coor.apr igs rinex tables
[zhenghui@node10 batchproc]$ sh_gamit -s 2008 275 279 -expt c008 -noftp >sh_gamit.log
ls: No match.
out of memnory
STOP Normal finish of autcln
STOP Normal finish of autcln
aopts: Subscript out of range.
ls: ut1.: No such file or directory
ls: pole.: No such file or directory
STOP Normal finish of autcln
....省略....
aopts: Subscript out of range.
[zhenghui@node10 batchproc]$
我不知道上面的红色部分是否正常,我在超算上(内存比较大)计算的时候没有out of menmory的提示。此外出现的utl.pole.错误也比较奇怪。因为我使用sh_gamit -yr 2008 的时候会自动生成一个tables目录,里面有一些拷贝或者链接文件。虽然终端同是也给出了提示ls: ut1.: No such file or directory ls: pole.: No such file or directory 但是在tables目录下却能够找到lu1.pole.文件,不知道这是为什么,揣测可能是这两个文件的权限问题。
>>出现内存不足情况,是否电脑内存配置太低?可考虑加个内存条了。后面提到出现ls 文件不存在,这很正常,因为这是个确认文件是否存在,然后进行复制的过程。
28.假如我处理的是2008年第100天的数据,当我用执行了sh_sp3fit命令,生成的g文件和t文件名为:gigs08.100 tigs08.100,但是后续的处理需要使用gigsf8.100和tigsf8.1000,如果不将sh_sp3fit生成的文件改名,好像是允许通不过,这点有点奇怪,不知道你碰见过没。
>>出现这中情况,首先你要搞懂sh_sp3fit的每一选项(注意-O)的意思,和之前sh_makexp各选项(注意-orbt)的含义。gamit对命名很敏感。
29.在中国地壳运动观测网络上下载了一天的观测数据,在处理的时候出现如下的问题:
zhenghui@zhenghui-desktop:~/zhenghui/project/chinacrust$ sh_sp3fit -f igs16342.sp3 -t -u
Generating fitted g-file for: 2011 123
....省略....
STATUS :110616:1302:50.0 ARC/arc: Integrating satellite 24 PRN 24
FATAL :110616:1302:50.0 ARC/redsat: Block # = 7for PRN 25 in svnav.dat not valid
STOP FATAL Error: Stop from report_stat
GAMIT.fatal exists, stop in sh_sp3fit
根据错误终端提示,我揣测是没有更新tables里面的svnav.dat,但是我重新更新了以后还是出现上面的提示,不知错在哪里?
>>是的,根据提示,在svnav.dat中找不到最新的25号卫星。请打开该文件,确定有如下这行:
25, 62 7 1555256. P 0.1230 2010 5 28 00 00 0.394 0. 1.407 #launched 2010 5 28 GPS II-F
30.出现错误FATAL :110624:2138:17.0 FIXDRV/bmake: Ocean loading requested no list or grid file,如何解决
>>主要是没有海潮文件,可以关闭海潮模型,或者去网上下载该模型文件。可到everest.mit.edu下载ot1_FES2004(730M),然后做相应的链接。
31.是否碰到过这种问题:接收机钟漂较大的情况下如何解基线;比如5700的接收机,我以前遇到过这种问题,总是在GAMIT解算是最后一步无法进行,提示的出错信息是Edit结束非正常,任何接收机均存在卫星钟的钟漂,只是改正的方式不同而已,有实时改正,也有半秒跳变.在5700设置中可选择施加接收机钟改正,这样就可保证与卫星钟同步.但那次不是我设的,所以没有加改正.后来请教姜卫平老师,他说要修改Edit模块的相关代码.设置里有clock drift,你可能是把它勾上了,就会加入改正,否则会存在钟漂.现在的接收机一般都自动加入钟差改正,所以出现的频率较少,但老机子不排除.接收机钟必须加入改正才能与卫星钟同步,改正方式有两种:一种是实时改正,一种是半秒改正,关于这一问题,可参考UNB的Kim在2004年ION发的一篇文章.
>>
32.FATAL :110714:1513:51.0 FIXDRV/armake: Invalid T-file interval = 6 23 2011 6 23 2011 --need either T-file or one X-file to run FIXDRV
>>Remove bogus test on t-file year > 2010. Copy gamit/fixdrv/armake.f. King 110101/110103,可以自己手动删除,或者下载更新文件。
33.请问如何在GAMIT中设置对流层随机过程模型?它与分段线性函数模型的设置有什么区别?如何设置?
>>你所指的“对流层随机过程模型”是否是“对流层延迟模型”?在sittbl.long 中可以设置对流层延迟模型,默认是SAAS(Saastamoinen模型),并且对干湿可以分别指定(但实际上,gamit中只有Saastamoien模型,别无二选)。与实测气象资料对比显示,采用SAAS模型效果更好,并且目前更关注的是映射函数(GMF或者VMF1)。分段线性函数或者阶跃函数均为参数求解的策略,可在sestbl.中设置。
34.我算的11年的数据,以前用10.31版本能正常解算,现在换成10.40执行最后一个命令csh bspan1.bat提示致命性错误:
FATAL :110924:1724:33.0 SOLVE/elem: S/C coordinates imply hyperbolic trajectory
FATAL GAMIT sh_chksolve: Solve failed to complete normally - check screen or log file
>>
35.我施工的地区没有已知点,所以我通过长时段静态观测,然后用Gamit与周边igs站进行联测解算,经过sh_gamit sh_glred两步,得到的itrf坐标,经转换至x,y,H, 这个坐标当做84坐标来用,有多大差别。
>>itrf坐标与WGS84有区别,采用最新的itrf框架,会与WGS84坐标相差在几个厘米级别,但这种精度在工程中完全可以忽略,因为由于地壳等运动,坐标的变动每年也在厘米级,所以重要的是内部精度要很好的满足。
36.平差过程中,是不是想固定哪个站点(例如6个igs站中的3个或者几个),就在sites.defaults文件该站最后标注 glreps 就可以固定了?
>>定义框架点,是在sites.defaults中用参数glreps表示,另外还有几个文件中有约束,可以修改。
37. sh_gamit -d 2011 297 expt scal -noftp -orbit IGSR Sh_gamit Version 10.4 (2011/1/17),出现错误:No g- or sp3-files available and noftp = Y, cannot continue
>>出现该问题的原因是处理快速星历,shell代码不完善。治标的解决办法:把下载精密星历时生成的g文件(gigsr*.***)复制到g文件夹(gfiles)中,即可解决。另外在下载快速星历时,不要加参数 -makeg no
38.在ubuntu下运行sh_get_rinex -archive sopac -yr 2012 -doy 050 -ndays 1 -sites wuhn,而不能下载数据,并用错误信息“ftp: bind: Address already in use”。
>>该原因主要是ubuntu系统问题,可以采用wget下载模式,而不采用默认的ftp下载。