GAMIT使用笔记

3.2.1 mesh edgessoft link :建立边与边的一个联接●Form－建立soft link，一条边只能在一个soft link 链中●Maintain－●BreakReverse反转边的方向。

Pick with links选择：对存在于同一个soft link中的边同时进行各种（如reverse，mesh）操作。

不选：可以对soft link中的某一个边单独进行操作。

此时如果要保持softlink 请选择Maintain。

Use first edge setting例如你已经对边a1，a2，a3进行了不同的mesh操作然后你又依次选择了a1，a2，a3如果选择了Use first edge setting那么gambit自动将a1的grading参数应用到其他边，在grading子面板中显示的的是a1的参数如果不选择那么grading子面板呈非激活状态，显示最近选择的边的参数Specifying the Grading Scheme非对称梯度结构（前六种）此时gamit保证相邻的线段的比R是常数］各种方法的不同之处只是在于确定R方法不同如下：Double－Sided Grading如果n是奇数，两边共用一个节点，且该与节点相邻的左右两线段长度相等。

如果n是偶数，两边共用一条线段。

对称梯度结构（后两种）Bi-exponent（双对数）Bell shaped正态分布参数输入与节点分布情况同上表specifying Node SpacingInerval Countm=n+1m是节点数n是间隔数Interval sized是Interval size，如果n是不是整数，取最近的整数。

Shortest edge（％）取Interval size为d＝x％L minLmin为最短边的长度,3.2.2SET Edge Element Type3节点是在每个网格元素的中间添加一个节点3.2.3Link／Unlink Edge Meshs（连接／断开边网格）link Edge Endpoint Vertices当建立一系列edge的硬连接时,gambit是把这些edge的起点都连接在一起,把终点也连接在一起，如果两条edge的端点已经被某个hark－link连接在一起，那么这两条edge之间就不能再建立hark－linkReversing the Grading Orientation在建立一系列edge的hardlink的同时gambit允许其中某些edge的grading orientateon（梯度方向）与其他edge相反，具体做法是shift－middle－click这些edge。

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中就可以对任一参数强化约束。

GAMIT10.4安装（基于ubuntu10.04）1、安装虚拟机vmware、ubuntu10.04；2、进入终端输入：sudo passwd root为root用户创建密码，并以root用户登录，或sudo su 回车；3、系统更新、汉化；4、安装gcc：# apt-get install gcc ；5、安装csh：# apt-get install csh 修改bash为csh，重启；6、安装gfortran ：# apt-get install gfortran ；7、安装libx11-dev库支持# apt-get install libx11-dev ；8、修改shall为bash，重启，并设置路径：回到用户根目录，打开.bashrc ，在最后加上如下代码即可export PATH=$PATH:/opt/gamit/gamit/bin:/opt/gamit/com:/opt/gamit/kf/binexport setenv HELP_DIR=/opt/gamit/help/9、将gamit安装包放在目录opt/gamit/下进入目录：# cd /opt/gamit10、修改install_software文件内容：# gedit install_software ，打开install_software文件，在文件的中下部修改“usr -name libX11.a”为"usr -name libX11.so" 。

（动态共享库）11、运行install_software ：# ./install_software，按提示输入两次Y后，修改makefile.config ，在/opt/gamit/libraries里，修改Makefile.config中的一组参数（1）MAXSIT 55 、MAXSAT 32 、MAXATM 25 、MAXEPC 5760（2）# Specific to FC5（F6,F7,F8 ）然后，在输入两次Y完成安装；12、安装完后，打开终端输入：doy，查看程序是否已正确安装。

Glorg_cmd:No specific station coordinates are tightly constrainedGlobk_cmdCom_file包含运行信息的common fileSrt_fileHfile时间排列list的二进制文件Sol_file结果和协方差矩阵的二进制文件To use glorg as stand-alone module, the samecom_file command must appear in both thecommand file for globk and the command file for the(subsequent) glorg run为使用glorg作为独立模块，globk和glorg控制文件运行都必须包含同样的com_file命令APR_XXX This command will free (estimate) a parameter.如果这个命令缺失，相应的参数将被固定到它的先验值如果先验sigma值定为0，将不估计参数（effectively left unconstrained)）强制一个参数到它的先验值，使用F作为sigmaApr_neu测站位置估计Apr_neu site sigN、sigEsigUsigVNsigVEsigVUSite是站名，可以应用到所有站，sigNsigEsigU位置的先验sigma（in meters）sigVNsigVEsigVU先验速度sigma（in m/yr）apr_neu SITE 10 10 10 F 0.1 0设置每个测站NEU方向先前sigma为10米，强制N速度为先验值，设置E速度先前sigma为0.1m/yr，不估计垂直向速度Mar_neu Estimating site positionstochasticallymar_neuMarkov随机游走过程（random walk process）Mar_neu site rwNrwErwUrwVNrwVErwVUSite 站名，可以all，rwNrwErwU位置的随机游走（m^2/yr）rwVNrwVErwVU速度的随机游走（in (m/yr)^2/yr）通常使用0 0 0apr_svsEstimating satellite orbitparameters设置航天飞行器（space vehicle）参数的先验sigmaApr_svs PRN X Y Z VxVyYzrad_parmsPRN 卫星号PRN_nn或者ALLXYZ 位置sigma（3个值）（m）VxVyVz速度sigma（3个值）（mm/s）rad_parms辐射参数（radiation parameter）sigma 高达11个值，0-1的小数全球测站估计轨道apr_svs all 100 100 100 10 10 10 1R只有当地坐标，约束轨道到先验值apr_svs all 0.1 0.1 0.1 0.01 0.01 0.01 0.01R最后的R表示剩下的radiation 参数EOP参数先验不确定度Apr_wob X Y XdotYdotX Y 极位置（pole position ）（mas）毫角秒XdotYdot变化率（mas/day）Apr_ut1 <ut1><ut1 dot>Ut1 in mas，ut1 dot in mas/day噪声处理一样形式Mar_wob RWX RWY RWXdRWYd单位都是（mas2）/yr（mas/day）2/yr Mar_ut1 for ut1(mas2)and (mas/day)2/yr通常使用以下两个形式全球网– apr_wob 10 10 1 1– apr_ut1 10 1区域网(constrained)– apr_wob 0.2 0.2 0.02 0.02– apr_ut1 0.2 0.02UT1问题需要小心使用UT1，因为它不能从卫星轨道节点分离如果UT1约束的话，UT1表和轨道必须使用同一框架，如果处理当地，是可以的If local processing, then should be OK如果使用不同的EOP表来确定轨道，UT1应该松弛Frame estmation commands 框架估计命令Apr_tran< X Y Z><VxVyVz>估计坐标系统的平移translationapr_scale<ppb><ppb/yr>估计缩放和缩放率scale scale ratemar_tran和mar_scale指定处理噪声通常全球网下使用File 命令一些命令来控制GLOBK读的文件Apr_file<name >使用一些先验坐标文件格式是：• Site_name X Y Z VxVyVz epoch• Site names are 8 characters, for GPS ABCD_GPS• X Y Z are geocenter Cartesian (m)• Velocities are m/yr• Epoch is decimal （十进位的、小数）yearsApriori position files先验位置文件Apr_files命令can be issued multiple times withthe latest values taken不过没指定，使用GAMIT先验坐标（没有速度）EXTENDED选项允许更复杂的behaviors （见globk.hlp）Orbit files轨道文件处理单天解时，默认使用GAMIT结果里的轨道多天解，需要用GLOBK,make_svs命令产生轨道文件格式：make_svs<file name>• <file name> will be overwritten• Command must be near top of command file.EOP files EOP文件默认使用Gamit的EOP，但可以通过用in_pmu命令改变格式：– in_pmu<file name>文件格式：– yy mm ddhr min Xpole +- Ypole +- UT1-AT +-– Must be uniformly spaced– Pole in arc-seconds, UT1 in time-secondsEarthquake files 地震文件使用地震文件可以自动说明地震引起的位移主要功能，用2-char code 说明地震，重新命名测站重命名命令：• Rename <old><new> [HFile code] [Epoch ranges][Position shift]Earthquakes in eq_file• For earthquake:• eq_def<code>Lat Long Radius Depth Epoch（精度纬度半径深度历元）• eq_cosei<code><Static Variance><distance dependentvariances>• eq_post<code><dur><Static RW><distance dependentRW>• eq_rename<code>• The rename option cause site to renamed fromxxxx_GPS to xxxx_G[Code] Glorg运行通过应用明确的旋转和平移，而不是一些站的紧约束来定义坐标系统可以单独程序调用或者在globk时运行Glorg imbedded commands嵌入的glorg命令• The globk控制文件嵌入的GLORG命令有：• org_cmd<glorg command file name>• org_opt<Options for output>• org_out<output file name>如果org_out没给，输出文件名将用.org替代Miscellaneous commands其他五花八门的命令max_chi<max chi**2 Increment><max prefitdifference><max rotation> • 允许自动删除坏的H文件和坏的坐标• app_ptid允许应用极潮改正，如果gamit运行里没有包括的话因为如果应用改正，这里没有指定，小心使用SINEX文件• crt_opt, prt_opt, org_opt指定屏幕输出、打印和org文件选项• glorg help 给出所有选项，主要有：– ERAS –写之前檫掉文件(通常追加文件)– NOPR –不写输出文件– BLEN –基线长度– BRAT -- baseline rates when velocities estimated估计速度时，基线率– VSUM –经纬度速度的summary (需要画速度)– PSUM –经纬度位置的summary– GDLF –包括H文件的列表和和运行时的卡方值chi**2 increments– CMDS –在输出文件里重复globk控制文件这些选项的主要作用是控制输出文件小点Command file layout控制文件布局–通常，globk中的命令可以无序安排–也可以重复，后边的将代替前面的•如apr_neu all 10 10 10 0 0 0• apr_neuiisc_gps 0.01 0.01 0.01 0 0 0• 将用10m约束所有站，除过IISC站用10mm约束–下面这些命令必须在前面出现: com_file, srt_file, make_svs, eq_fileGlobk velocity command file• 如果多年的数据，尝试解算速度• globk控制文件必须包括– apr_neu all 10 10 10 1 1 1–其中1 1 1 是速度不确定度in m/yr–如果速度和指定站相关，添加：– apr_neuiisc_gps 0.01 0.01 0.01 F FF–将固定IISC的速度–通常不这么做，而是用glorg来广义Glred为了检查H文件的独立位置的质量，我们可以只是用H文件列表生成gdl文件（we could create gdl files with just single hfiles listed），每一个用包含下面内容的控制文件运行globk•apr_neu all 10 10 10 0 0 0• apr_neuiisc_gps 0.01 0.01 0.01 0 0 0• iisc再次被约束and normal process• 用GLOBK运行一系列的独立的H文件是冗繁的（Running globk on a set of individual hfiles is tedious），所以我们用可以自动做这些的程序：glredglred和globk同样的运行方式，例如：– glred<crt><prt><log><gdl file><command file>–而不是联合gdl文件里的所有数据,每个文件用globk单独处理。

Gambit学习笔记（2）（转帖）⼀.边界层边界层是指定与边或者⾯相邻区域的⽹格接点的距离，⽬的是控制⽹格密度，从⽽控制感兴趣区域计算模型的有⽤的信息量。

例如：在⼀个液体流管中，我们知道靠近壁⾯处的速度剃度⼤，⽽中⼼处的速度剃度⼩，为了使得壁⾯处的⽹格密⽽中⼼处的⽹格稀疏，我们就在壁⾯处加⼀边界层。

这样我们就能控制⽹个密度。

要定义⼀个边界层，你要定义以下参数：1）边界曾依附的边或者⾯2）指定边界层⽅向的⾯或者体3）第⼀列⽹格的⾼度4）相邻列之间的⽐例因⼦5）总列数，指定了边界层的深度同时，你也可以指定⼀个过度边界层。

要指定⼀个过度边界层，你需要定义以下参数（过度模式，过度的列数）1)⽣成边界层需要定义以下参数：i)size:包括指定第⼀列的⾼度和相邻列的⽐例因⼦ii)internal continuity ：当在体的某个⾯上施加边界层时，gambit会把边界层印在与这个⾯相邻的所有⾯上，如果在体的两个或者更多的⾯上施加边界层，那么边界层就有可能重叠，internal continuity 这个参数就决定了边界层如何重叠当选择internal continuity 时，gambit不会在相邻的⾯上互相施加边界层否则就会在相邻的⾯上互相施加边界层具体可以看guide的图⽰：同时这个参数还影响施加了边界层的体可以采⽤何种⽅式划分⽹格corner shape ：gambit 允许你控制conner（即边界层依附的的两条边的连接点处）附近区域⽹格的形状iii)Transition Characteristics需要定义以下两个参数：Transition pattern指的是边界层远离依附边或者⾯那⼀侧的节点的排布情况Number of transition rows这个列数肯定要⼩于前⾯指定的那个列数。

iv)Attachment Entity and Direction指定⽅向⾮常重要，可以通过⿏标和list 对话框来完成。

GAMIT使用
GAMIT的使用需要依赖于多个输入数据源。

首先，需要GPS观测数据
来确定地球表面的位置。

这些观测数据通常通过全球分布的GPS测量站获得。

其次，需要卫星轨道测量数据来确定卫星位置和运动的精度。

最后，
还需要一个大地坐标系以及与卫星测量数据相关的各种参数。

这些数据最
终用于计算地壳形变的模型估计。

GAMIT的工作流程主要分为数据处理和模型估计两个部分。

在数据处
理方面，数据首先通过一个预处理阶段进行数据检查和修正，以去除可能
的误差和异常。

接下来，数据经过一个精确的观测数据处理流程，包括时
钟和轨道参数的建模和估计。

最后，通过使用所有这些数据和参数，计算
出地壳的运动模型。

在模型估计方面，GAMIT使用单差估计方法，将所有的观测数据与参
考站数据进行比较。

通过不同站点之间的差异，可以进一步估计出地壳运
动的速率和方向。

此外，GAMIT还可以提供地壳变形模型的参数，如弹性
滞后和地震震级。

这些估计结果对于地球物理学家来说非常重要，可以帮
助他们更好地理解地壳运动的动力学过程。

GAMIT具有一些突出的特点和优势。

首先，它是一个开源的软件工具，可以免费使用。

其次，它具有较高的计算精度和稳定性，可以处理大规模
的观测数据。

此外，GAMIT还具有较好的可扩展性，可以与其他地球物理
研究工具进行集成和扩展。

最后，GAMIT提供了用户友好的图形界面和详
细的文档，以帮助用户更好地理解和使用该工具。

GLOBK: Global Kalman filter analysis programGLOBK Ver 5.18: Global Kalman filter analysis program.5.18版本：Runstring:运行字符串：% GLOBK <std out> <print file> <log file> <exper. list> <command file> <OPTION>where <std out> is a numerical value (if 6 is typed then output will besent to current window, any other numerical value willsend output to a file fort.nn)<std out>是一个数值（如果是6，则结果将会发送至当前窗口，如果是其它数字，结果将会输出到fort.nn文件）<print file> is the name for the output print file with thesolution in it. If the print file already exists, thenthe new solution will be appended to it.<print file>是输出文件的名称，解决方案保存在输出文件中。

如果该输出文件已经存在，那么新的解决方案将会附加其中。

<log file> is a log file which contains the running time for theprogram and the pre-fit chi**2 value for each inputcovariance matrix file. If the log file already exists,then the new solution will be appended to it.<log file>是一个日志文件，记录了程序运行时间和每个输入的协方差矩阵文件的pre-fit chi**2值。

gamit操作步骤Gamit是一种在全球定位系统（GPS）应用中用于数据处理和分析的软件，可以用于精确测量和分析地球上的运动、形变和地震活动。

下面是使用Gamit进行数据处理和分析的详细操作步骤。

1.数据准备：首先，您需要准备GPS观测数据和相关的参考桩，以便进行后续的数据处理和分析。

确保数据的准确性和完整性，并将其存储在计算机上的合适位置。

2.数据导入：打开Gamit软件，并从菜单中选择“数据导入”选项。

在弹出的对话框中，选择您准备好的GPS观测数据文件，并导入到Gamit中。

4.固定点选择：根据您的需要选择一些固定点，这些点的坐标已知且稳定。

这些固定点将用于基线解算和数据校正。

5.数据质量控制：在进行后续的数据处理和分析之前，您需要进行数据质量控制。

这可能包括检查数据的完整性、纠正观测误差、排除异常值等。

6.基线解算：在数据准备和质量控制完成后，进行基线解算以获取各个GPS观测站的坐标。

使用Gamit中的基线解算工具，输入固定点和待解算的观测点，然后选择合适的解算方法和参数。

7.高斯坐标转换：如果您需要将GPS观测数据的坐标转换为其他坐标系统（如高斯坐标系），则可以使用Gamit中的坐标转换工具。

根据项目需求输入相关参数并执行坐标转换。

8.形变分析：如果您希望通过GPS观测数据进行形变分析，可以使用Gamit中的形变分析工具。

输入相关数据和参数，并运行形变分析以获取形变测量和分析结果。

9.地震活动分析：若想利用GPS观测数据进行地震活动分析，可以使用Gamit中的地震活动分析工具。

输入相关数据和参数，并运行地震活动分析以获取地震活动监测结果。

10.结果输出：在完成数据处理和分析后，您可以选择将结果输出到相关的文件或报告中。

Gamit提供了多种输出格式和选项，以满足不同需求。

以上是使用Gamit进行数据处理和分析的一般操作步骤。

然而，具体的操作步骤和流程可能根据项目需求和数据类型而有所不同。

在使用Gamit之前，建议先熟悉软件的基本功能和操作方法，并参考软件的用户手册和指南。