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

博文

解决GAMIT不识别最新的O-file文件中的观测值代码

已有 3608 次阅读 2019-10-28 23:33 |个人分类:GAMIT/GLOBK|系统分类:科研笔记

一、基线结算出现如下错误:

FATAL  :191028:1849:57.0 MODEL/setup: GPS frequencies on x-/c-file file not L1 and L2

打开生成的x文件查看发现,只有L1载波和C1伪距值,L2均为0,如下:

x00019.240 (~-project-shenda-liuh-240) - VIM_007.png

说明一下,本次结算采用RINEX3.04格式的o文件,观测数据如下:

00012400.19o (~-project-shenda-liuh-240) - VIM_008.png

GPS的观测值有6个,其中拥有C1C,C2Y,L1C和L2Y,然而GAMIT只读取了单频的数据;


二、寻找问题根源

单步运行:

makex demo.makex.batch

发现出现如下错误:

WARNING:191028:2050:29.0 MAKEX/settyp: Single-frequency observations
WARNING:191028:2050:29.0 MAKEX/settyp: No L2 pseudorange: codeless tracking


根据上面提示,可以看到是在~/gamit/makex/settyp.f出现错误;

打开这个文件,查找发现是由于iobtypx这个数组中2和4元素都为0,导致错误发生;

image.png

继续探索,iobtypx是在哪里赋值的?

发现iobtypx是settyp函数的一个参数,由外部传入,因此需要找寻在哪里调用了settyp;

      Subroutine settyp ( debug,rxver,nobtyp,rxobtyp,iobtypx

     .                  , nwave2,nsat,rcvrsw,swver,ndat,dattyp,lambda )

查找调用settyp,发现在makex中调用了此函数;

image.png

打开makex文件发现iobtypx赋值是在sel_obtyp函数中完成的!

image.png

所以找到并打开~/gamit/lib/sel_obtyp.f函数文件

image.png

发现G中并不包含C2Y和L2Y观测,至此,终于找到问题所在!


三、解决办法,添加我们需要的观测值,然后重新编译;

image.png



https://blog.sciencenet.cn/blog-858128-1203824.html

上一篇:[转载]RINEX3.04观测值代码
下一篇:FATAL:SV antenna offsets for SVN ** not found in antmod.dat
收藏 IP: 119.79.228.*| 热度|

0

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

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

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

GMT+8, 2024-12-22 21:18

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部