GAMIT软件操作手册

GLORG: GLOBK coordinate frame realization programGLORG:GLOBK坐标框架实现程序GLORG Ver 5.13: Origin resolution program for the GLOBK.GLORG Ver5.13：GLOBK的原始解决程序。

Runstring:% glorg <output> <options> <command_file> <com_file> <OPTION>where <output> is the name of the output file (may be 6 for outputto current window.<output>输出文件的名字（输出到当前窗口为6位）<options> is the bit mapped option selection or it can bespecified by codes separated by : or = (no spaces)(Code feature added at Ver. 4.0)<options>是二进制选项，可以用代码表示，用：或=（不是空格）分开（4.0版本增加代码功能）The options are:CODE BIT Decimal MeaningCORR 0 1 Output correlation matrix输出相关矩阵BLEN 1 2 Output baseline lengths and components输出基线长度和元件BRAT 2 4 Output baseline lengths and components ratesof change.输出基线长度和元件的变化率CMDS 3 8 Write a summary of the markov file to theoutput file.将马尔科夫文件的总结写入输出文件VSUM 4 16 Write the short version of the velocity fieldinformation (one line per site)写出速度场信息的短版本（每个站点一行）5-9 32-512 NO LONGER USED (see POS_ORG and RATE_ORGbelow).不再使用（参见下面POS_ORG 和RATE_ORG）RAFX 10 1024 Fix the Right ascension origin of the system.固定系统的赤经起点MOUT 11 2048 Only output baselines if both sites are Markov.(Used to limit output in large back solutions)如果站点都是马尔科夫，只输出基线（用来限制大型的返回解决方案）COVA 12 4096 Output full precision covariance matrix.输出所有精确协方差矩阵PSUM 13 8192 Output position adjustments in summary form以摘要形式输出位置调整GDLF 14 16384 Output the GDL file used in the solution输出解决方案中的GDL文件DBUG 15 Output matrix details when there are negativevariances and negative chi**2 increments如果有负方差和负chi**2增量，输出矩阵细节ERAS 16 Erase the output file before writing solution写解决方案前改正输出文件NOPR 17 Do not output the file (either crt, prt or orgdendending on opt set).不输出文件（根据选项设置，决定crt，prt或org）SDET 18 Output details of the stabilization calculationsin glorg输出稳定计算的细节在glorgRNRP 19 Report the statistics of the differences in thepositions and velocities of renamed sites. Alsogenerates equate lines than can be filtered withshell script sh_exeqs. The equates lines arewritten to file <org root>.eqs报告重命名站点的位置和速度差异统计。

Gamit软件安装说明SUSE Linux 8.0使用说明一、安装组件选择：1，在categories中选择：●Advanced-Devel●Gnome System●KDE Desktop Environment●Linux Development Tools●Network/Server●Help Support Documentation2，选择Single Package，在其中利用Search选择需要安装的组件，OK3，利用Time Zone改时间/时区二、Login Shell修改安装完毕，利用Yast2控制项（用户及密码）修改为CSH三、使Windows下的磁盘分区（包括NTFS文件系统）可用：1，修改etc/fstab.文件，在其中增加: mount –a (链接指定的磁盘分区) 或2，在root的/sbin/下，运行fdisk /dev/hda 查看分区情况（p 查看，q退出）运行mkdir Dell_D 建立链接目录运行mount /dev/hda5 Dell_D 完成链接。

四、修改启动系统选项编辑lilo.conf（或Redhat的grub.conf）文件，更改启动顺序，并运行lilo（或grub）。

五、Gamit 等软件路径设置1，在一般用户的首级目录(例如：home/cxc/)下，拷贝/etc的csh.cshrc文件为.cshrc：cp /etc/csh.cshrc .cshrc2，编辑.cshrc文件，在其最后一行增加如下路径：Setenv PATH“/home/cxc/gs110/com:/home/cxc/gs110/gamit/bin: /home/cxc/gs110/kf/bin: $PATH”（$PATH表示以前设置的路径）3，运行 source .cshrc 以激活以上设置六、设置vi在/etc/vimrc中设置vi，将其中打开“syntax on”选项(去掉22行号)。

GAMIT/GLOBK软件使用手册一软解介绍GAMIT软件最初由美国麻省理工学院研制, 后与美国SCRIPPS海洋研究所共同开发改进。

该软件是世界上最优秀的GSP定位和定轨软件之一, 采用精密星历和高精度起算点时, 其解算长基线的相对精度能达到10-9量级, 解算短基线的精度能优于1mm, 特点是运算速度快、版木更新周期短以及在精度许可范围内自动化处理程度高等, 因此应用相当广泛。

GAMIT软件由许多不同功能的模块组成, 这些模块可以独立地运行。

按其功能可分成两个部分: 数据准备和数据处理。

此外, 该软件还带有功能强大的shell程序。

目前，比较著名的GPS数据处理软件主要有美国麻省理工学院（MIT）和海洋研究所（SIO）联合研制的GAMIT/GLOBK软件、瑞士伯尔尼大学研制的BERNESE软件、美国喷气推进实验室（JPL）研制的GIPSY软件等。

GAMIT/GLOBK和BERNESE软件采用相位双差数据作为基本解算数据，GIPSY软件采用非差相位数据作为基本解算数据，在精度方面，三个软件没有明显的差异，都可得到厘米级的点位坐标精度。

相比较而言，GIPSY软件为美国军方研制的软件，国内只能得到它的执行程序，在国内，它的用户并不多，BERNESE软件需要购买，它的用户稍微多一点，GAMIT/GLOBK软件接近于自由软件，在国内拥有大量用户。

GLOBK软件核心思想是卡尔曼滤波（卡尔曼滤波理论是一种对动态系统进行数据处理的有效方法, 它利用观测向量来估计随时间不断变化的状态向量），其主要目的是综合处理多元测量数据。

GLOBK的主要输人是经GAMIT处理后的h-file和近似坐标, 当然,它亦己成功地应用于综合处理其它的GPS软件（如Bernese和GIPSY）产生的数据以及其它大地测量和SLR 观测数据。

GLOBK的主要输出有测站坐标的时间序列、测站平均坐标、测站速度和多时段轨道参数，GLOBK可以有效地检验不同约束条件下的影响, 因为单时段分析使用了非常宽松的约束条件，所以在GLOBK中就可以对任一参数强化约束。

GAMITGLOBK软件操作GAMITGLOBK是一种用于精密测量地壳运动和形变的软件工具。

它以其高度准确的结果和灵活可靠的功能而受到许多地质学家和测量师的青睐。

本文将详细介绍GAMITGLOBK的操作步骤和功能。

首先，为了使用GAMITGLOBK，我们需要在电脑上安装该软件。

在安装完成后，我们可以通过命令提示符或终端来打开GAMITGLOBK。

在打开GAMITGLOBK后，我们将进入软件的主界面。

在主界面中，我们可以看到一系列的选项和命令。

首先，我们可以使用“load”命令导入所需的测量数据文件。

GAMITGLOBK支持多种数据格式，包括RINEX、SP3和轨道文件。

我们可以使用命令“load rinex filename rinexdir”来导入RINEX文件。

同样，我们可以使用命令“load sp3 filename”来导入SP3文件。

一旦我们导入了所需的数据文件，我们可以使用GAMITGLOBK来进行数据处理和分析。

其中一个常用的功能是进行基线数据处理。

为了进行基线数据处理，我们可以使用命令“base bx by bz ant”来指定基线的坐标和天线类型。

然后，我们可以使用命令“solve”来执行基线数据处理。

结果将包括各个基线的长度和倾角。

此外，GAMITGLOBK还提供了一些高级功能和选项。

例如，我们可以使用“process”命令来指定数据处理的参数和选项。

我们可以使用“model”命令来选择地球模型。

我们还可以使用“trop”命令来计算和校正对流层延迟误差。

在数据处理和分析完成后，我们可以使用GAMITGLOBK来生成报告和结果。

我们可以使用“report”命令来生成结果的报告文件。

报告文件将包括数据处理的摘要、结果的图表和图像。

我们也可以使用“save”命令来保存数据处理的结果和文件。

综上所述，GAMITGLOBK是一种功能强大且易于使用的地壳运动和形变测量软件工具。

通过了解GAMITGLOBK的操作步骤和功能，我们可以更好地利用该软件进行数据处理和分析，并获得准确可靠的测量结果。

武汉大学测绘学院GAMIT/GLOBK数据处理报告[键入文档副标题]李文文20122021400092012/12/13GAMIT/GLOBK 是一套高精度数据处理软件，主要用于分析研究地壳变形、高精度GPS测量数据处理等领域。

它由美国麻省理工学院( MIT) 和斯克里普斯海洋研究所(SIO) 联合开发，并得到美国哈佛大学和美国国家科学基金会的支持，是目前世界上应用最为广泛的高精度GPS数据处理软件之一。

GAMIT/GLOBK基于UNIX（Linux）系统开发和运行。

本文中所有数据处理工作均是基于Ubuntu9.0与csh SHELL环境下完成的。

一数据预备为了学习使用GAMIT处理GPS数据，本文选择2012.07.01（DOY 183）天如下共15个全球IGS跟踪站建立全球观测网。

由于该网最初是用于评定北斗电离层模型的改正精度，故而在选站上更加偏重中国及周边地区。

在完成跟踪网选择后需要下载相应的导航电文和精密星历数据。

这些数据亦可以通过GAMIT中的sh_get_rinex, sh_get_navs, sh_get_orbits脚本根据指定的站点名称和时间直接从CDDIS，SOPAC等服务器上下载。

这里需要注意的是，由于这些脚本均是基于csh（或tcsh）解释器，故而在bash环境中无法正确执行。

总结准备数据的相关信息如下：二建立工程根据GAMIT软件处理要求，需要建立相关目录。

一个GAMIT工程主要包括如下几个工程目录：DOY: processing data, final solutions, etc.rinex: observation file in RINEX o format.igs: precise orbit file from IGS in sp3 orbit Project Namebrdc: broadcast file in RINEX N formattables: table files linked to ~/ gg/tablesOther directories created during processing最初建立工程只需要在主工程目录下建立相应的DOY, igs, igs, brdc四个目录，并在相应的目录存放数据。

GAMIT数据处理图⽂流程G A M I T10.5数据处理图⽂流程冉启顺⽬录⼀、前⾔⼆、数据准备1.前期准备⾸先想好要处理什么时段、什么地⽅的数据本⽂档以处理bjfs,shao,lhaz,urum四个站点2013年2⽉1⽇的数据为例本⽂档使⽤软件版本GAMIT10.5处理环境：win10上的虚拟机ubuntu15.102.⼯程⽬录设置在任意位置建⽴⼯程⽂件夹（全⽂以test为例），并在其⽬录下建⽴如下的⽂件夹：3.数据下载在终端中打开test，查询2013年2⽉1⽇的年积⽇，和GPS周等信息代码：doy 2013 02 01[1]在test的⽬录下再建⼀个⽂件夹以032命名的⽂件夹，此时，test⽬录下有：[2]根据步骤3的信息，到相关⽹站下载对应的⼴播星历，精密星历，观测值⽂件分别放在⽂件夹brdc,igs,rinex⾥⾯。

当然还涉及到下载的数据解压，.d⽂件到.o⽂件的转换等操作，其中⼴播星历可以不改名字（即可以不⽤改成brdc0320.13n）4.更新tables表，更新准则，下⾯的表格为准：还有个dcb.dat需要每次处理都要更新i.常见⽂件说明ii.更新的准则是：“更新频次”中的红字部分每次处理都更新⼀下，同时结合本次处理与上次处理的年份是否相同来进⾏年表的选择。

原本有些⽂件是每周更新⼀次，但是作为新⼿为了避免出错还是每次更新吧有些时候，有些⽂件找不到，或者⾥⾯没有，则可以忽略。

如svs_exclude，gdetic.dat常常找不到我在更新时，这四个⽂件在服务器上没找到下图是我更新的⽂件。

原本还应该有2013年的各种年表，但是我之前处理过⼀次2013年的数据，并且我的gamit安装路径/opt/gamit/GAMIT10.5/tables已经有了2013年的这些年表。

iii.将上述更新的数据复制到gamit安装⽬录下的tables⽂件夹内，并替换原有⽂件。

我的gamit安装路径是/opt/gamit，即我安装路径下的tables 的路径是：/opt/gamit/GAMIT10.5/tablesiv.然后将更新后的tables⽂件夹，即/opt/gamit/GAMIT10.5/tables拷贝（或者设置连接，我习惯拷贝）到test⽬录下。

gamit 操作步骤摘要：一、引言二、GAMIT 软件介绍三、GAMIT 操作步骤1.准备工作2.数据导入3.参数设置4.模型计算5.结果输出与分析四、GAMIT 软件在实际应用中的优势五、结论正文：【引言】GAMIT（Geodetic Analysis and Modeling Integrated Tool）是一款综合性的地球物理建模分析工具，广泛应用于大地测量、地震学、地壳形变等领域。

本文将详细介绍GAMIT 软件的操作步骤，以帮助用户更好地使用该软件进行科研工作。

【GAMIT 软件介绍】GAMIT 软件由美国麻省理工学院（MIT）地球、大气与行星科学系的研究人员开发，是一款功能强大的地球物理建模分析软件。

GAMIT 软件主要应用于GPS 数据分析、地球物理建模、地壳形变分析等领域，支持多种数据格式，具有丰富的功能和灵活的参数设置。

【GAMIT 操作步骤】1.准备工作在使用GAMIT 软件前，需要确保电脑上安装了GAMIT 软件，并正确配置环境变量。

此外，需要准备相应的数据文件，如观测数据、测站信息、基线文件等。

2.数据导入将所需数据文件导入GAMIT 软件中，包括观测数据（如GPS 观测数据）、测站信息文件（如站点坐标、高程等）和基线文件（如基线向量等）。

3.参数设置GAMIT 软件具有丰富的参数设置，包括地球模型、大气模型、噪声模型等。

用户需要根据实际需求和数据特点选择合适的参数设置。

此外，还可以根据需要设置其他参数，如迭代次数、约束条件等。

4.模型计算在完成参数设置后，启动GAMIT 软件进行模型计算。

软件将根据设定的参数对数据进行处理，生成地球物理模型。

5.结果输出与分析计算完成后，GAMIT 软件会生成一系列结果文件，包括地球物理模型、残差图、相关系数等。

用户可以根据需要对这些结果进行分析和解读，以获取有关地球物理过程的信息。

【GAMIT 软件在实际应用中的优势】GAMIT 软件在实际应用中具有以下优势：（1）功能强大，支持多种数据格式和多种地球物理模型；（2）参数设置灵活，可根据用户需求进行调整；（3）计算速度快，适用于大规模数据处理；（4）结果准确，能够提供可靠的地球物理模型。

解算软件GAMIT处理流程文档.txt小时候觉得父亲不简单，后来觉得自己不简单，再后来觉得自己孩子不简单。

越是想知道自己是不是忘记的时候，反而记得越清楚。

gamit 处理流程流程概览：1、规划工程2、更新table文件夹3、创建工程目录，连接table4、编辑sites.defaults5、手动建立其余必须目录6、下载、准备数据7、制作准备文件8、sh_gamit批处理解算9、具体介绍：1、工程规划：处理中国大陆IGS站8个（wuhn、shao、chan、bjfs、kunm、guao、urum、lhaz、），台湾IGS站2个(tnml、twtf)数据跨度：2008年10月1日~2008年10月5日，2009年10月1日~2009年10月5日框架站（6个）：wuhn、chan、bjfs、kunm、lhaz、twtf2、更新~/gg/tables文件夹(需要一定权限)可到/processing/gamit/ 更新所需表文件，或者到ftp:///archive/garner/gamit/tables/ 下载（用户名：ananymous 密码：邮箱）通常来自IERS的地球自转参数表ut1. 、pole. 需要每日或周更新章动表nutabl. 太阳星历表soltab. 月亮星历表luntab. 跳秒表leap.sec 需要每年更新接收机/天线列表rcvant.dat 在有新仪器使用时需要更新卫星列表svnav.dat 在有新卫星使用时更新码偏文件(P1-C1、P1-P2)dcb.dat 需要月更新如果想方便的话，每次解算时全部更新3、建立工程目录，例如china2008,同时运行sh_setup 命令链接talbes文件夹可以查看ut1. 文件，是否在你处理时间内，处理的是2008年数据，显然在范围内。

4、编辑工程china2008文件夹下的talbes中的sites.defaults文件，把需要处理的测站加进去，并规划参考固定站(globk需要)。

GAMIT软件Track模块使用帮助TrackRT : Realtime GPS kinematic processing programTrackRTB : Realtime GPS kinematic processing program for use with Pre2.5 BNC TrackRTr : Rinex file emulator of trackRT for post-processing evaluationRunstring:% trackRT <options>% trackRTB <options>where options are-m <machine> Host name of system with real-time data port(Normally BNC would be running on this system)(Default local host)-p <port> Port number on host supplying the data stream(Mandatory input).-f <command file> trackRT command file (similar to track command file) -r <ref code> 4-char code for reference site.-d <list of 4-character codes> list of sites to be processed. If notspecified all sites in data stream will be processed(there is a limit on maximum number of sites that can beprocessed depending on max_site in track. The referencecode should NOT be given here.)-n <root name> Specifies a root name to which output will be directed(.out is appended for stdout output). Use ? in thepos_root, sum_file and csv_file names the ? will bereplaced by this string.Command file name must be given.% trackRTr <options> Same list of option although -m, -p have no effect -d <list of rinex files> list of rinex files (site names are generated fromthe first 4-characters of rinex file names)VERSIONS:---------1.13 -- Fixed time-equate iin SaveObsA so that a differnce of morethan 0.1 seconds in needed for a new epoch (trackRT only)Added reference relative humidity for GPT. (120103)1.12 -- Handle bad PRN from BNC where G23 == 203.1.11 -- Version to handle missing satellites in SP3 files (110505)1.10 -- Version to handler BNC2.5 and greater that is now outputingascii rather than binary data (trackRTB handles BNC2.0) (110501) 1.00 -- Initial version 100301 (2010 March 1)INTRODUCTION:-------------TrackRT/TrackRTB and trackRTr have many commands that allow it be controlled both when started and during runs. Below the detailed commands are given and here we give the minimum commands needed.SP3_DIR : (unless sp3 files are in the current directory)SITE_POS: Site positions must be given for every site to be processed.SITE_STATS: To set the apriori sigma and process noise for the sites.NOTE: The reference site position should be set to zerofor the apriori sigma and process noise. Typical site_statswould be where cit1 here is the reference site.site_statsall 0.1 0.1 0.1 0.025 0.025 0.025cit1 0.0 0.0 0.0 0.0 0.0 0.0ATM_STATS: By default atmospheric delays are not estimated, and normally these should be. Normally the reference site is set to zerowhen site separations are less the 500-1000km. Typical atm_statscommand would beatm_statsall 0.20 0.00010 0.000 ! Unit m/sqrt(sec) -> 0.0001 = 0.03 m/sqrt(day) cit1 0.00 0.00000 0.000If mixed antenna and receiver types are usedANTE_OFF: To specify the antenna heights and types and receiver types ANTMOD_FILE: Must be given to get the antenna phase center modelsRCV_TYPE: Needed for mixed receiver types (entry can be specified inANTE_OFF also and this command would not be needed).DCB_FILE: Up to date, data code bias (DCB) file (part of GAMIT ftp area).IncludingUPDATE_FILE is useful so that trackRT can be controlled on the fly.TrackRTB should be used for older versions of BNC. If no files are createdby trackRT or trackRTB (list of files should be listed shortly after theprogram is started) then try using the other version. The program are identical except trackRTB is expecting a binary data stream.COMMANDS--------Input/Output commands---------------------SP3_DIR <Directory> <center><Directory> -- Directory where sp3 files are stored<center> -- Center for orbits (default igs, igr and igu also tested)POS_ROOT <root> <duration>Set the root part of the name for the output files.<root> Root part of name. When ? included in the <root>, the ? is replaced with the -n string.<duration> Duration of data in each file. The designations of d, h, or m may be used to specify the units of days, hours, minutes.Default is d.The position file names take the form: (see commands below)<root>.<outtype>.<site>.<start time>.<datatype>The resolution of the <start time> which is modulo the output interval depends on the output duration. For output durations greater than or equal to 1-day, the time is given as YYYYMMDD. For intervals shorter than 1-day, it is YYYYMMDD:HHMN. Minium output interval is 1 minute.SUM_FILE <root>Sets the root part of the summary file name. Using ? in the namewith be replaced by the -n string. (Default if command is not givenis trackRT or the -n string when -n used). File names are time taggedaccording to the pos_root output interval.CSV_ROOT <root>Set the root part of the name for comma separated values (CSV) outputfile. These files are used for AmCharts web plots.DCB_FILE <file name>Set the name of the data-code-bias (DCB) file. This file is part of theGAMIT tables directory and should be updated regularly. It is used to remove biases in the Melbourne-Wubbena widelanes. The receiver type can be specified with the RCV_TYPE or ANTE_OFF commands.ANTMOD_FILE <file name>Sets the name of a standard IGS antex file with phase center models forthe GPS ground antennas. Antenna types at specific sites are given withthe ANTE_OFF command. This command my be used multiple times for site specific model with new models replacing previously read ones.UPDATE_FILE <file name>Allows new trackRT commands to be issued during a run. Once the file is read it needs to be deleted before the trackRT will re-read it. File isonly read if it exists. NOTE: File should be removed before trackRT isrun or else it will be read when the command file is read (ie., it will overwrite the commands in the command file.OUT_TYPE <NEU+GEOD+XYZ+DHU>Specifies types of output coordinates. All types can specified in a single string with no spaces. The types areNEU -- North, East, Up differences from the reference site or from the coordinates given in the REF_NEU command.GEOD -- Geodetic latitude, longitude and height (in the GEOD format, the total atmosheric delay is given, while in the other formats theadjustment to the apriori delay is given).XYZ -- Cartersian XYZ coordinatesDHU -- Delta horizontal and Up coordinates from the apriori coordinates of each site (default output type)OUT_SIG_LIMIT <sigma (m)>Sets the maximum sigma of a position estimate for it to be output.If pseudorange data types are used, the default value of 1 m needsto increased to 10-100 meters.Analysis commands-----------------DATA_NOISE <L1 (m)> <L2 (m)> <P1 (m)> <P2 (m)> <Elev Weight> [PRN] Allows specificiation of the noise in the L1 phase, L2 phase,P1 range and P2 range, and the weight given to elevation angledepepence (at ver 1.20); variance is scale by (1+(W/sin(el))^2)where W is the <Elev Weight>. These values affect the sigmas printed for the position determinationsOptional: PRN may be added and noise assigned to that PRN (ifnon-PRN form is used, this will replace all PRN specific valuesso use the non-PRN first followed by specific PRN valuesDATA_TYPE <type>Specifies type of data to process. Types supported are L1 L2 LC P1 P1 PC which can be combined in a contiguous string. Files names include thisstring at the end.Vers 1.0: Only LC, LCPC and PC have been tested. Specifically L1 only data have not been fully implemented yet.USE_GPTGMF <Relative humidity (0-1)>Set the used of the GPT temperature pressure model and the GPT dry and wetmapping functions. Optional relative humity added vers 1.27; default 0.00. Default is the older MTT atmospheric model.Site and Antenna informationSite information is entered by first giving the command and then the information by site on the subsequent linesSITE_POSSite <X (m)> <Y (m)> <Z (m)> <Vx (m)> <Vy (m)> <Vz (m)><Epoch (yrs)>Site is the four character name of the site (more characters can be includedbut only the first 4 are checked). Site names that do not appeat in the listof sites to be processed are ignored). The remainder of the line contains positon and velocity and the epoch in deciminal years to which the position refers.ANTE_OFFSite <ARP dN (m)> <ARP dE (m)> <ARP dU (m)> <Antenna Name><Receiver Code>Specifies the type of antenna and its position of antena reference point (ARP) at each site. The antenna name including radome should be specified with theofficial IGS name for a standard ANTEX file or with a unique name that appearsin the ANTEX file for site specific calibrations. (Note: There is oneadditional character in the long antenna and this extra character before the radome name must be removed. The antenna name and radome can becopied directly from the rinex file if present. One more ANTEX files must be specified with the ANTMOD_FILE command for the antenna names to be useful. The receiver type DCB code can be optionally specified here as well (SeeRCV_TYPEcommand).RCV_TYPESite <Receiver code N/P/C>Specifies the type of data-code-bias (DCB) correction needed for the receiver. Code specifies the type of L1 and L2 ranges being measures. The choices are P -- Pcode, C -- C/A and N C/A with cross corelation for L2 range. The codes can be found in gamit/tables/rcvant.dat. These codes can also be given in the ante_off command. An up-to-date DCB_FILE command must be used to specify theDCB biases. The files are available from the MIT ftp site and update once per month.SITE_STATSSite <Apriori Sigmas in XYZ (m)> <RW noises in XYZ (m/sqrt(sec))> Gives statistics to assign to the kinematic station positions. The<Apriori Sigma in XYZ> are the three sigmas in XYZ for the initial postion and<RW noise in XYZ> are the three sigmas for the random walk in position. Units of the ransom walk are m/ssqrt(s).ALL can be used for the station name and the same statistics will be appliedto all kinematic sites (NOTE: the fixed site do not change position).e.g.site_statsall 1.0 1.0 1.0 0.010 0.010 0.010base 0.0 0.0 0.0 0.000 0.000 0.000(1 meters apriori sigmas and changes 10 mm/sqrt(s)=600 mm/sqrt(hr)).NOTE: One station should always set to zero sigmas and random walk otherwise allstation positions will be estimated and noramlly drift by large amounts.The fixed station does NOT need to be the reference site.ATM_STATSSite <Apriori Zenith delay sigma> <RW noise in Zenith delay> <RW dH/dt noise>Gives the statistics for the atmospheric delays by site. A random walk process's moise increases as sqrt(time) where time is time in seconds. The valuesare the initial sigma in meters, RW changes in meters per sqrt(second) and (addedverion 1.2) a dH/dt variance term so that during rapid height changes more process noise can be added to zenith delay estimate. The process noise variance is(<RW dH/dt noise)*abs(dh/dt)>^2 per sqrt(sec) where dh/dt is m/s. The default valeis 0.00023 m^2/sqrt(s). (Term only needed for aircraft processing).e.g.,atm_statsall 0.10 0.0001 0.00000t39a 0.10 0.0001 0.00023base 0.00 0.0000 0.00000Sets the apriori sigma as 10cm and allows the delay to change 6 mm/sqrt(hr) and 2.3 mm/sqrt(s) when height is changing at 10 m/s (fast ascent or desent) for the site t39a (aircraft). In the case, here "base" is a reference station at whichthe atmospheric delay is not estimated (other stations absorbe the atmospheric delays at the base station). With kinematic positioning over few hundred kmsiteseparations, the atmospheric delay at one site is normally fixed.Ambiguity resolution/data control commands.-------------------------------------------trackRT uses a combination of the Melbourne-Wubbena widelane (MW-WL), the extrawidelane (EX-WL) and the floating point estimates of the ionospheric free ambiguity (LC)to resolve integer ambiguities. If we denote the number of integer cycle ambiguitiesat L1 and L2 by N1 and N2, The MW-WL is an estimate of N1-N2 based on phase and rangedata; EX-WL = N1 - f1/f2 N2 and is an integer for L1 cycles, but 1.283 N2 for L2 cycles;LC = 2.546 N1 - 1.984 N2. The EX-WL is unaffected by geometric ranges changes, but doesdepend on the ionospheric delay. For short baselines, the EX-WL should be near zero forcorrect choices of N1 and N2. The LC residual should also be near zero when N1 and N2are correct. The problem in ambiguity resolution is that difference choices ofN1 andN2 can make different linear combinations small. For example, errors in N1 and N2 of3 and4 cycles will change LC by 0.298 cycles (56.6 mm), MW-WL by 1 cycle and EX-WL by2.132 cycles (405 mm). On long baselines, at low elevation angles, ionospheric delaysof 400 mm are common. The most common error is a N1=N2=1 cycle error.For thiscombination, the MW-WL is unaffected and LC changes by 0.562 cyc (107 mm) and EX-WL by0.283 cyc (54 mm). Even on relative short baselines, 54 mm ionospheric delays are commonThe follow commands are use to control the ambiguity resolution.AMB_SET <RelRank> <FloatSigma (2)> <MWWL Fact> <EXWL Fact> <Min AmbSig> <MaxChi>Sets parameters for ambiguity resolution. The input parameters are:<RelRank><WL min> Minumum of values need to allow bias fixing<WL avN> Maximum number to be used in computing sigma ofmean MW-WL<FloatSigma (2)> Minimum sigma for LC and MW-WL for ambiguity fixinf<MWWL Fact> Weighting factor for MW-WL in chi**2<EXWL Fact> Weighting factor for extra-wide lange<Min AmbSig> Minimum sigma to assigned to float estimates of ambiguities.<MaxChi> Max chi**2 value allowed for ambiguity to be resolved.EXWL_SET <Jump> <Min Sigma> <Scale> <Elev Fact><Jump> -- magnitude of jump in EX-WL to have cycle slip added (default 0.10)<Min Sigma> -- minimum sigma for mean ex-wl (cycles, default 0.02 cycles)<Scale> -- Scaling factor for length. Scale 0.1 results in 0.1 cyclesover 100 km (default)<Elev Factor> -- Elevation angle factor that increases sigma as (1 +factor/sin(elev))MWWL_SET <Jump> <Min Sigma> <Max Averaging number> <Min number><Jump> -- magnitude of jump in MW-WL to have cycle slip added (default 5.0)<Min Sigma> -- minimum sigma for mean MW-WL (cycles, default 0.10 cycles) <Max Averaging number> -- Maximum number of values to use to compute mean sigma<Min number> -- Minimum number needed to resolve ambiguityDD_SET <Jump (cycle)> <Min Number>Sets parameters for double difference processing.<Jump (cycle)> -- Magnitude of jump in double differences on bias fixed data thatwill introduce a cycle slip<Min Number> -- Minumum number of double differences for epoch to be processed.(If too few than errors in the data can be detected than this cancause large position errors, default of 4 double differences allowsredunancy).RMS_EDIT_TOL <n-sigma Tolerance> <min sigma> <Reset number><n-sigma tolerance> is an n-sigma condition where sigma is based on data noise model.<min sigma> Minimum phase sigma to use to that no phase residual less than <<min sigma>*<n-sigma tolerance> are deleted<Reset number> number of sequential delete data, before ambigity and cycle slipsare reset (assumed missed cycle slip).Testing and evaluting commandsSTATUS <type> <# epoch>Writes status information to the current summary file at <# epoch> intervals The types of reports are given by "type"P -- Parameter estimatesA -- Ambiguity resolution report (shows resolved and unresolved)W -- Widelanes (Melbourne-Wubbena and Extra-widelanes)R -- Postfit residuals are current epochC -- Report current A and W entries onlytype PAWR will output all reports.NUM_EPOCHS <number><number> is the number of epochs of data to be processed before stopping. When shortduraton files are output (POS_ROOT command) the number of epochs may not be reachedbecuase the epoch counter is reset with each new file.START_TIME <yy mm dd hh min sec.>Used to set the start time. Useful with trackRTr to have the processing start at thesame time as the realtime stream.DEBUG <Up to 10 epoch numbers in pairs>Sets which epochs will report detailed debug and status information (see STATUS commandas well). The pairs are used as:1,2 -- Mostly model information. Useful if data are generating large residuals3,4 -- Parameter estimates, widelanes, residuals5,6 -- One-way OMC and single differences7,8 -- Not Used9,10 -- When 9 is non-zero, antenna model and SP3 information.RESET <ALL/list of sites>Reset command. Resets the filter state vector and resolved ambiquities for a list of sites or ALL sites (generally used in the update_file to fix problems)EXCLUDE_SVS <list of PRN numbers to be excluded>Excludes satellites from being processed. Useful when a satellite is notin an SP3 file or not available during the times there are data.Example: exclude_svs 26 22 13Output types------------The output position files have two header lines that explain the columns in the file. The output in the summary and output files contain the following basic blocks are written to the summary file. The time in these files is given byan epoch number which is a count of the number of epochs of data (at the referencesite) since the start of the current data files. The output position files givethis counter and the corresponding GPS times and deciminal days.CSLIP line are cycle slip detectionsCSLIP Ep 540 Site THMG G 29 DMW-WL/Tol 4.31 4.00DEX-WL/Tol 1.78 0.20 cyc, Elev 10.00 degIndicates the jump in the MW-WL and EW-WL (and the tolerance for the detection)that caused the slip to be detected.CSLIP Ep 961 Site THMG G 12 DD RESID 2.18 Tol 0.50 cyc, Elev 63.85 degIndicates a slip added due to repeated double difference residuals.When ambiguities are resolved, three lines are output giving the statisticsAMBFIX WMAP PRN29 EP 599 RelRank 101.76 FC ------ dL12 4 3 Dchi 6.54 691.02 AMB 36AMBWLS WMAP PRN29 EP 599 RG 133 599 FX 3 iL12 0 0 Means -0.16 -0.05 RMS 0.74 0.06 # 467 eN12 -0.54 -0.38 AzEl 274.32 12.35AMBCON WMAP PRN29 EP 599 NCont 3 MW Res -0.159 0.185 Chi20.74 EX Res -0.051 0.068 Chi2 0.55 LC Res -0.046 0.020 Chi2 5.25 AMBFIX gives the epoch number, relative rank, the Fix Code (always ------ when fixed), thechanges to the number of L1 and L2 cycles from the intially asigned values, the changes in chi**2 for the best and second best choices of ambiguities and the number of the ambiguity resolved. When the status feature is used, ambiguities not resolved yet are listed as AMBFREE entries.AMBWLS gives the mean values and RMS scatters of the MW-WL and EX-WL (cycles).The iL12 estimate here is the changes in L1 and L2 cycles implied by justthe widelanes with no contribution from the LC estimates. On short baselines, these values are expected to be zero.AMBCON gives the contributiond to chi**2 from each of the input type. Ncont is the number of contributions (3 for LC data type). Values given arethe residal (Res) with its sigmas and Chi**2 (Chi2) contribution,When STATUS output is requested the following output types are possible. STATUS REPORT Epoch 1800 Date 2010 05 06 08 36 47.000 Type PAWR PARAMETER ESTIMATES Epoch 1800 for NP 16 Type P STATUS REPORT ...PARAMETER Ep 1800 NP 1 MSCG dATM -0.0466 +- 0.0001 m Var 0.681E-08PARAMETER Ep 1800 NP 2 PSDM dATM -0.0001 +- 0.0001 m Var 0.436E-08...PARAMETER Ep 1800 NP 7 LJRN G 27 1 0.1281 +- 0.0253 cyc Var 0.639E-03PARAMETER Ep 1800 NP 8 LJRN G 05 1 6.1172 +- 0.0011 cyc Var 0.112E-05Parameter values for position adjustments (XYZ), atmospheric delay adjustementsand float ambiguities are given.Ambiguity report for both both fixed and free ambiguities. ResEpoch is the epoch whenthe ambiquity was fixed (if fixed in previous data file block, may be greater than current count).AMBIGUITY Report Ep 1800 Number of ambiquities 59AMBFIXD # 11 MSCG PRN10 RelRank 616.27 FCode ------ dL12-1 -1 Dchi 1.01 776.94 Elev 72.02 deg; ResEpoch 156 AMBFIXD # 12 MSCG PRN27 RelRank 135.91 FCode ------ dL12-1 -1 Dchi 5.39 766.54 Elev 29.86 deg; ResEpoch 156..AMBFREE # 41 LJRN PRN27 RelRank 11.37 FCode RW---- dL120 0 Dchi 25.69 294.89 Elev 30.60 deg; ResEpoch 0 AMBFREE # 42 LJRN PRN05 RelRank 11.76 FCode R----C dL12 7 6 Dchi 106.15 1251.17 Elev 36.16 deg; ResEpoch 0The FCode gives the reason for not being fixed the entries are: RWSSSCR - Relative rankW - Widelane sigma too large or not enough data yet (see MWWWL_set command)S - LC sigma too largeS - L1 sigma too large when L1+L2 data type used (not tested)S - L2 sigma too large when L1+L2 data type used (not tested)C - Chi**2 increment too large.Widelane report shows status of average values of widelanes and current LC estimate. X means fixed, R ambiguity is still free.WIDELANE Report Ep 1800 Number of ambiquities 59WIDELANE 11 MSCG PRN10 EP 1800 Range 2 1800 # 1799 NC 3 X MW Res 0.06 0.10 Chi2 0.30 EX Res -0.18 0.19 Chi2 0.82 LC Res -0.01 0.02 Chi2 0.10...WIDELANE 19 MSCG PRN29 EP 1800 Range 421 1800 # 1380 NC 3 R MW Res -0.09 0.15 Chi2 0.35 EX Res 0.07 0.26 Chi2 0.07 LC Res 0.17 0.02 Chi2 74.30Postfit residuals for data types being used.POSTFIT RESIDUAL Report Ep 1800 Number of DD 76POST 1 Ep 1800 MSCG PRN10 - CIT1 PRN02 Res -0.0945 +- 0.0381 LC AMB 11 Elev 37.31 degPOST 2 Ep 1800 MSCG PRN27 - CIT1 PRN02 Res -0.1040 +- 0.0441 LC AMB 12 Elev 22.17 deg...POST 41 Ep 1800 MSCG PRN04 - CIT1 PRN02 Res 2.9712 +-1.8988 PC AMB 13 Elev 30.60 degPOST 42 Ep 1800 MSCG PRN30 - CIT1 PRN02 Res -0.2485 +-1.8116 PC AMB 14 Elev 36.16 deg..++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++EXAMPLE COMMAND FILE:--------------------* FILE trackRT_pbo.cmd* Example TrackRT command file to process some PBO sites.* Run BNC with as the TCP/IP port (login information* can be obtained from .) and select sites P496, P497,* P498, P505 to be cast on specific port (port 3765 for example).* Run trackRT as* trackRT -p 3765 -r p497 -d p496 p498 p505 -f trackRT_pbo.cmd -n P497 & * (if BNC is run on another machine used -m <machine name> as well.* A relative directory ../sp3_files is assumed to contain up to date* sp3 file.sp3_dir ../sp3_files* These positions are post April 4, Baja earthquake (extracted from* globk output with sh_exglk -f <globk org file> -apr <apr file>* entries from apr_file* Nore: Additional sites can be given heresite_posP498_GGU -2313650.15153 -4835550.43945 3444474.55150 -0.00658 0.02008 0.01478 2010.259 0.0033 0.0059 0.0042 -1.0000 -1.0000 -1.0000P497_GGU -2315937.74208 -4838712.09671 3438545.10828 -0.00945 0.01927 0.01884 2010.259 0.0035 0.0063 0.0044 -1.0000 -1.0000 -1.0000P496_GGU -2319722.10265 -4842511.39340 3430709.79068 -0.01697 0.02370 0.02587 2010.259 0.0050 0.0089 0.0062 -1.0000 -1.0000 -1.0000P503_GGU -2325056.31018 -4826758.67413 3449209.09091 -0.01167 0.01951 0.01660 2010.259 0.0033 0.0059 0.0042 -1.0000 -1.0000 -1.0000P505_AGU -2309739.55681 -4802072.50955 3493258.58085 -0.00377 0.00608 0.00160 2010.404 0.0033 0.0057 0.0040* PBO sites antenna and receiver information (all are NetRS and thus C type) * Again more sites can given). (Note: If antenna information extracted from* gamit , one space must be removed before radome string).* Antenna/Radom combination must be in antmod_file for model to be used).ante_offp475 0.00 0.00 0.0083 TRM29659.00 SCIT Cp066 0.00 0.00 0.0083 TRM29659.00 SCIT C p472 0.00 0.00 0.0083 TRM29659.00 SCIT C p478 0.00 0.00 0.0083 TRM29659.00 SCIT C p494 0.00 0.00 0.0083 TRM41249.00 SCIT C p496 0.00 0.00 0.0083 TRM29659.00 SCIT C p497 0.00 0.00 0.0083 TRM29659.00 SCIT C p498 0.00 0.00 0.0083 TRM41249.00 SCIT C p500 0.00 0.00 0.0083 TRM29659.00 SCIT C p503 0.00 0.00 0.0083 TRM29659.00 SCIT C p505 0.00 0.00 0.0083 TRM29659.00 SCIT C p510 0.00 0.00 0.0083 TRM29659.00 SCIT C* Change path as needed.antmod_file /home/tah/gg/tables/antmod.datdcb_file /home/tah/gg/tables/dcb.datdata_type LCPCsite_statsall 0.1 0.1 0.1 0.025 0.025 0.025p497 0.0 0.0 0.0 0.0 0.0 0.0atm_statsall 0.20 0.00010 0.000 ! Unit m/sqrt(sec) -> 0.0001 = 0.03 m/sqrt(day) p497 0.00 0.00000 0.000pos_root ? 1h# Used at MIT for web output. Results can be viewed at# /kmeduna/# csv_root /net/chandler/var//kmeduna/trackrt_view/data/P497# Useful to have this option. Note file should be removed before trackRT# is run and created when needed.update_file upd_app.cmd++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++INSTALLATION------------trackRT requires the BKG NTRIP Client (BNC) and QT libraries and include files. Ntrip is the Networked Transport of RTCM via Internet Protocol and the client program BNC provides access to the realtime data streams and casts the date via an internet port GNSS raw data in a structure that is decoded by trackRT. Currently trackRT only uses the GPS data streams. TheNTRIP system is discussedat igs.bkg.bund.de/ntrip/. The BNC program can be downloaded fromigs.bkg.bund.de/ntrip/download. The BNC documentation constains the instructions for obtaining the QT libraries and includes needed for itsinstallation from source files. trackRT only needs the executable version of BNC. The QT software can be obtained directly from Noia at the siteqt.nokia./downloadsFor trackRT, the libraries and include files are needed and these are containedin Frameworks libraries 4.6.2 for your platform (linux or Mac). Earlier version should also work. For Linux and Max OSX, QT can be downloaded directly from qt.nokia..qt.nokia./downloads/linux-x11-cppqt.nokia./downloads/qt-for-open-source-cpp-development-on-mac-os-xQT takes several hours to install.On linux systems, other tools that can be needed are:sudo app-get install gfortransudo app-get install g++sudo app-get build-dep qt4-qmake(The last app-get may be need if the qt install can not find X11 libraries and includes.)For MacOSXThe basic qt install steps aresudo port install qt4-macIf you don't have MacPorts, it's a very easy way to get lots of *nix toolsfor the Mac. You can download MacPorts from .Once QT is installed, the makefile in the trackRT directory should be modifiedto include the directories for the libraries and includes. Example entries are given in the makefile. The actual entries are for the MIT QT installation (which is not quite standard).NOTE: There is a file makefile (not Makefile) for making trackRT series of programs.To run trackRT, bnc is run first and the desired stations added the mountpoints. Once bnc streaming is started, trackRT is started with the machine name (if different from the bnc machine) and port (given in the bnc FEED tab) given. From。