GAMIT学习资料整理
GAMBIT实例教程 5_引入和整理轿车IGES几何结构
5．引入和整理轿车几何结构在这个指南中，你将会引入一个包含轿车几何结构的IGES文件，整理改几何结构并用三角形和四面体对其划分网格。
在这个指南中你将学到如何：·导入IGES文件·指定给几何图形着色的方式·用手动和自动的方式接边·合并各面·生成一个三角形的表面网格·用四面体网格给体积划分网格·准备将这些网格读入FLUENT 55.1 前提这个指南假定你已经做完指南1的工作，因而熟悉GAMBIT GUI（图形用户界面）。
5.2 问题描述所要考虑的问题的示意图如5－1所示；它是一辆豪华轿车的外形，你要在轿车车体外面生成一套网格；因此，要生成一个围绕该几何结构的方体来代表流场。
5.3 解决方案在这个指南中，你会生成一套围绕作为IGES文件导入的汽车车体几何结构的完全非结构化四面体网格。
这个指南举例说明了要为一个导入的CAD几何结构划分网格所需遵循的典型步骤。
导入的几何结构是“不清晰的”－也就是说，有一些面之间有裂缝使得其不适用于生成一个CFD网格，你首先要用GAMBIT中提供的工具对几何结构进行整理。
大部分的缝隙可以在网格导入的过程中或者后续过程中通过“接边”指令自动修补。
原始的CAD几何结构在整理过程中没有改变；修补缝隙所需的更改要利用“虚拟的”的几何结构来进行，该结构在“真实的”的几何结构之上。
在原始的几何结构里有些边非常短且可以用“顶点连接”指令将其去除。
其它的边不能自动接合，因为它们离得比规定偏差远的多。
你要手动地连接这样的边。
导入的几何结构包括一定数量的小表面，这些面的边不必要地限制了网格地生成过程。
利用“合并表面”指令，GAMBIT可以让你毫不费力地在网格划分之前合并这些表面，然后你可以利用GAMBIT在汽车表面自动生成三角形网格。
因为导入的几何结构只包括车体，所以你必须围绕车体生成一个合适的区域以便于进行CFD分析（这大致等同于将汽车放置于一个风道内）。
GAMIT软件数据处理
概
述
长安大学 空间定位技术与信息研究所
GAMIT软件及其数据处理的主要内容
软件功能概述 GAMIT/GLOBK软件的安装 GAMIT软件的数据处理步骤
一、GAMIT/GLOBK软件功能概述
由MIT及SIO开发，编程语言为F77、C， 需要带X-Windows的UNIX/LINUX环境。 功能：处理GPS静态单时段观测数据，计 算测站坐标、卫星轨道等参数及协方差， 但不能计算钟差。 观测量：GPS单频、双频相位或伪距观测 数据及卫星星历。
二、GAMIT/GLOBK软件安装
设置运行环境参数
对 于 Unix System V 派 生 的 Unix 系 统 ， 通 常命令 环 境为 Cshell，设置文件为.cshrc，设置方法是增加以下路径： Set path=$PATH/home/temp/gamit1040/com /home/temp/gamit1040/gamit/bin /home/temp/gamit1040/kf/bin
leap.sec
luntab.2007.J2000 nutabl.2007 no rcvant.dat soltab.2007.J2000 svnav.dat no sestbl.regional sittbl. svs_exclude.dat
跳秒表
三、GAMIT软件的数据处理步骤
观测文件及星历文件准备
观测文件包含RINEX格式观测O文件及导航N文 件，星历文件指IGS精密星历。为了便于数据 组织，我们通常将观测文件及星历文件按天存 放于指定目录中。
数据处理目录
001 002 · · · · · · tables
GAMIT使用手册
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中就可以对任一参数强化约束。
GAMIT软件操作手册
GAMIT软件操作手册目录一GPS误差分析 (3)1.1 与GPS卫星有关的误差 (3)1.2 与信号传播有关的误差 (3)1.3 与接收设备有关的误差 (4)1.4 其他误差来源 (4)二GPS基线处理的几个关键问题 (5)2.1 星历 (5)2.2 对流层折射影响 (5)2.3 周跳是否修复是影响基线解算精度的因素之一 (6)2.4 基准点坐标的确定 (6)2.5 基线解算是否在地固系中进行 (6)2.6 整周未知数的确定 (7)三GPS应用软件介绍 (8)3.1 一般的商用软件 (8)3.2 高精度GPS软件 (8)四GAMIT软件简介 (9)4.1 概述 (9)4.2 主要模块介绍 (9)五GAMIT软件的安装 (11)六GAMIT软件的运行 (12)附录一、LINUX操作系统的安装: (15)附录二、GCC的安装: (18)附录三精密星历及相关表文件的获取 (20)附录四RINEX格式说明 (21)一GPS误差分析GPS是美国为了满足军事部门和民用对连续实时和三维导航的迫切要求于1973年开始研制的，至1994年整个系统全面建成。
这个系统的全称是“授时与测距导航系统/全球定位系统”（Navigation System Timing and Ranging/Global Positioning System—NAVSTAR/ GPS），通常称为“全球定位系统”（GPS）。
它能够在全球范围内提供全天候、高精度、连续实时的三维定位和测速，同时它还能够提供时间基准。
GPS是20世纪空间技术上的最大成就之一。
它的出现使大地测量产生了根本性的变革。
目前这一高新技术已广泛地应用于大地测量学、地球动力学、精密工程测量、地壳形变监测、石油勘探、资源调查、城市测量等领域。
影响GPS定位的误差按其主要来源可以分为如下几个部分：1.1 与GPS卫星有关的误差●星历误差与模型误差●卫星钟差与稳定性●卫星摄动●相位的不稳定性●卫星的相位中心1.2 与信号传播有关的误差●电离层折射●对流层折射●多路径效应1.3 与接收设备有关的误差●接收机钟差●天线的相位中心●观测误差（天线的整平与对中、量取天线高的误差）●接收机噪声1.4 其他误差来源●地球自转的影响（极移、UT1）●相对论效应的影响（信号传播与卫星钟）●地球潮汐（固体潮、海潮、大气负载潮）二 GPS 基线处理的几个关键问题在高精度GPS 测量中，影响定位精度的主要因素有：卫星的轨道精度、对流层折射的修正精度、多路径效应、相位中心的改正、接收机震荡器的稳定度、数据的后处理技术和起始点坐标的精度。
GAMIT简介
(2)文件准备(续)
kf . 10. 21. tar . Z kf 软件压缩包 libraries. 10. 21. tar . Z 库文件压缩包 templates. 10. 21. tar . Z 数据模板压缩包 map s. 10. 21. tar . Z 地图数据压缩包 将install_software设置为可执行文件 chmod +x install_software ./install_software //运行install_software
第一步： 数据编辑得出干净的观测数据文件——X-file. 数据编辑用快速解（Quick Solution）进行， 其目的并不是求解参数，而是生成干净的 Xfile. 流程见图1：
18
MAKEX MODEL AUTCLN SOLVE CFMRG CVIEW SCANDD
//将RINEX转换成 将 转换成GAMIT中的观测文件：X-file 中的观测文件： 转换成 中的观测文件
9
GAMIT安装 安装:(1)重新编译 重新编译GCC 安装 重新编译
从一下网站下载gcc-3.4.6.tar.bz2源代码： /gnu/gcc/gcc-3.4.6/gcc3.4.6.tar.bz2 在当前工作目录下解压： tar –jxvf gcc-3.4.6.tar.bz2 改设置：vi ~/gcc-3.4.6/libf2c/libI77/fio.h 将“ # define MAXUNIT100”的行更改为“ # define MAXUNIT 10000” 配置：./configure –prefix=/usr –enablethreads=posix –eanble—long-long –enablelanguages=c,c++,java,f77 重新编译和安装gcc/ g77 ： make make install
Gambit学习总结
容值：是指操作新老几何之间的最大缝合距离，可分为自动和手动，自动是指按照系统默认值进行操作，手动则可指定缝合距离（默认情况下是10-3）。
平滑示例：
修复示例：
3.6虚几何操作
反之可以利用Disconnect功能将连接好的实体断开。
3.4分割
分割操作可以利用一个实体将另外一个实体分割成一个或多个对象。以下图为例说明面的分割，蓝框表示图形A，绿框表示图形B。
以B来分割A，则出现以下两个面：
以A来分割B，则出现以下两个面：
如果选择互相分割，则会出现以下三个面：性是一个非常重要的概念，它可以使流体从一个面（或体）传递到另一个面（或体），但这两个面（或体）必须是连接的或者是非形定义的。
形网格：点在体的交界面处共享；
非形网格：点在体的交界面处不共享。
形网格非形网格
点、线、面都可以连接，该操作将消除所有重复的实体并重新连接拓扑结构，但只有实体才会被连接，在连接的同时，现有网格会被保存。
合并：将两个体合并成一个实体或虚体；
分割：将一个体分割成两个实体或虚体；
连接：将两个体连接成一个虚体。
3.7Unite、Merge和Connect的比较
Unite(Real):面必须有匹配的切线；
Unite（Virtual）：空隙和重叠是允许的，但不能连结边；
Merge：可以在实体或者虚体上操作，结果永远是虚体，面与面之间必须共享一条边，但不必相切：
Reset：删除所有的网格和几何；
Reset mesh：删除网格但保留几何。
2
Gambit的菜单栏主要包括三个部分，分别是文件菜单、编辑菜单和求解器菜单。
gamit笔记
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）（转帖）
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模块，精密定位还GAMIT、GLOBK两个模块都需要。
安装完成后的几个重要文件：gg/gamit（基线平差）和gg/kf（Kalman Filter）两个目录下到模块是用fortran编写的。
gg/com是cshell编写到脚本，重要用于gamit和kf目录下的模块的组织。
gg/tables是表文件。
年更新LUNTAB、SOLTAB 、NUTABL、LEAP.SEC周更新UT1、POLEsh_gamit批处理要求工程目录下至少有rinex brdc gfiles三个目录。
分别放O文件，N文件，卫星轨道文件g文件，这样做的目的是把文件分类，最后这些文件都会被link到单天的目录之下。
注意：需要将所有观测文件和表文件都link到单天目录下的，sh_gamit能自动完成link功能。
模型说明：1.otl 潮汐改正2.vmfl GMF 投影函数3.atml大气荷载模型，对高程影响较大，可消除周跳波动，可靠性需要进一步证实4.atl大气抄袭荷载模型和met气象模型星历文件：e/n, sp3, g,te/n为广播星历，主要用来你和卫星和接收机的种差g文件是根据sp3文件拟合的某天的圆形轨道参数t文件是根据观测文件和g文件求出的卫星位置，是gamit专用格式gcc编译器作用：将常见的编程语言转化为c语言。
安装gcc需要把原来到gcc覆盖。
在/usr目录下，具体怎么做，不是很清楚。
软件中的栅格文件：下载地址：ftp://1)海洋潮汐。
例如otl_FES2004.grid放在软件talbels目录下。
链接到otl.grid。
2)大气负荷。
例如atmldisp_cm.2006，每年更新一次。
连接到atmldisp_grid.20063)vmfl投影函数栅格，例如vmflgrd2006，连接到map.grid.2006。
每年更新一次。
以example为例作一个实例：1）在/media/Tool/TOOL/专业工具/GAMIT下新建文件夹10-05-18-EXAMPLE，在该目录下建立tables目录。
GAMIT-GLOBK入门介绍及应用(含sh_gamit和sh_glred批处理)
1、GAMIT/GLOBK介绍
主流高精度GPS处理软件(A、B级GPS网)：
GAMIT/GLOBK(美国麻省理工学院MIT) GIPSY/OASIS(美国喷气动力实验室JPL) Bernese(瑞士伯尔尼大学)
国内以GAMIT/GLOBK为主 GAMIT+ CosaGPS
I. 编辑Makefile.config
2012-4-6
国测二大队
8
3、 GAMIT/GLOBK安装
2012-4-6
国测二大队
9
3、 GAMIT/GLOBK安装
II. 修改安装脚本可执行属性 chmod +x install_software III. 执行./install_software(采用默认即可)
2012-4-6 国测二大队 26
4、 GAMIT/GLOBK处理流程
V. 基线解算 sh_gamit -s 2008 275 279 -expt c008 -noftp > sh_gamit.log
参数：
-d yr days 要处理的日期，年，年积日，例如： 1997 153 156 178 （处理1997年的153天、156 天和178天的数据） -s yr d1 d2 d1和d2是要处理的开始和结束的年 积日，例如：1997 153 178（处理1997年从 153天到178天的数据）
2012-4-6
国测二大队
22
4、 GAMIT/GLOBK处理流程
IV. 制作准备文件（ sittbl.）
测站的精度控制，高精度的已知坐标强约束， 待求点坐标松弛约束
2012-4-6
国测二大队
23
GAMIT简介pjh
GAMIT简介上海天文台程宗颐一，G AMIT的特点――用双差做基本观测量。
1，优点：消除了站钟和星钟的主要误差；轨道误差对测站相对位置不敏感，有利于精密定位；2，缺点：有效观测减少, 各站有效观测因其在网内的位置的变化而变化；不能用于单点定位；3, 在X-Windows 支持下的UNIX 操作系统或LINUX 操作系统下运行。
二，G AMIT的功能1, 估计卫星状态矢量（座标和速度）；2, 估计太阳光压参数（3－9个）；3, 估计y-bias；4, 估计erp 参数；5, 估计站坐标，可得站坐标的时间序列，确定测站位置随时间而变的规律，从而可研究板块运动和地壳形变以及其它地球物理现象；6, 估计对流层引起的天顶延迟，从而能开展GPS气象学方面的应用和研究。
三，G AMIT软件的使用简介。
1，建立以（积日）天命名的工作目录，tables目录(建立两个目录积日工作目录与tables 目录即可。
积日目录：比如，2003322是指2003年第322天的工作目录，里面存放o-files （相位与伪距观测文件），n-files（导航文件，卫星状态和星钟改正信息），以及该工作日的sp3-file （GPS卫星精密星历，从IGS网站下载）2 在tables目录中建立1〕测站坐标初始文件（L 文件lfile.），地心球坐标;2〕测站信息文件（ ），接收机和天线型号信息，天线高；3〕测段信息控制文件（sestbl. ）4〕测站信息控制文件（sittbl. ）5〕星号对照表（svnav.dat ）, 卫星质量，yaw rate,6〕接收机及其天线型号对照表（rcvant.dat ）7〕天线相位中心改正表（antmod.dat ）8〕地球形状参数表（gdetic.dat ）9〕太阳表( soltab. )10〕月亮表( luntab. )11〕章动表( nutabl. )12〕跳秒表( leap.sec )13〕周跳的自动探测和修改命令表( autcln.cmd )14〕海潮表( scherneck_grid, scherneck_stations.oct )15〕极移表（pole. ）16〕TAI－UT1 表（ut1. ）执行erp 命令即可产生这两个文件。
汇报—岳彩亚 - GAMIT
24
数据处理流程
（2）制作台站先验坐标 lfile. 执行：
% grep POSITION ../rinex/*o > lfile.rnx（提取o文件中的测站坐标） % rx2apr lfile.rnx $year $doy（转化成GAMIT识别的文件） % gapr_to_l lfile.rnx.apr lfile. $year $doy（将空间坐标转换成球面坐标） 这样，就得到了测站信息文件 lfile.（rinex/文件夹下），将其拷贝到tables/文 件夹即可。
二、IGS核心站数据的获取
广播星历： ftp:///pub/rinex/YYYY/DOY/ (以auto开头的文件，下载后改 为brdc开头) ftp:///pub/gps/data/daily/YYYY/brdc/ （以brdc开头） ftp://igs.ensg.ign.fr/pub/igs/data/YYYY/DOY/ （以brdc开头） sh_get_brdc -archive sopac -yr $year -doy $njr -ndays $ts ······
拷 贝 GAMIT/GLOBK 安 装 目 录 下 的 tables/ 文 件 夹 下 的 globk_comb.cmd + glorg_comb.cmd文件到解算目录下的gsoln/ 文件夹。
12
shell在GAMIT用途
GAMIT是由一些程序模块组成，并且这些程序模块可以实现调用，他们在 GAMIT数据处理过程中起到的作用互不相同。
使用shell脚本的编写可以将不同的程序模块组合起来实现不同的作用。 GAMIT 自身带了一些批处理命令，这些命令只是达到固定的结果，通过编写 shell脚本语言可以实现新的批处理过程，节约了时间，提高效率。
gambit 建模学习笔记
1、为什么更改后不显示边的颜色更改边的颜色命令允许用户更改与一条或者多条边相关的几何结构和/或者网格节点的显示颜色。
使用Modify Edge Color窗口要打开Modify Edge Color窗口（如下图），点击Geometry/Edge子工具框中的Modify Color 命令按钮即可。
Modify Edge Color窗口包含如下详细设定。
Edges 指定要更改颜色的一条或者多条边。
Color:--------------------------------------Geometry指定更改指定边的颜色。
Mesh指定更改与指定边相关的网格结点颜色。
体会：只能更改某点、边、面、体（没有组成更高级拓扑之前）的颜色，例如线在没有组成面之前是可以改变它的颜色的，而如果边组成了面它就不再可以改变颜色了。
2、不太懂2.3.2 连接/分离边Connect/Disconnect Edges命令按钮允许用户进行以下操作。
命令详细说明图标Connect Edges合并重合的实边或者生成代表一条或多条现有边的连接的虚边Disconnect About RealEdge分离具有公共实边的面和体以下部分将详细说明执行上面列举操作的步骤和要求的设定。
注意：Graphics/Windows Control工具框中的Specify Color Mode命令按钮允许用户显示基于实体连接而不是拓扑结构的模型颜色。
关于使用Specify Color Mode命令按钮的详细说明，请参阅GAMBIT User's Guide中的3.4.2部分。
合并边Connect Edges命令允许用户合并两条或者多条边。
要合并边，用户必须指定如下参数：∙要合并的两条或者多条边∙合并后的类型指定要合并的边要合并的边可以是实际的或者虚拟的，但是它们受到合并后类型的限制（见下面）。
指定合并后的类型GAMBIT允许以下类型的边合并操作：∙Real∙Virtual (Forced)∙Virtual (Tolerance)∙Real and Virtual (Tolerance)如果用户使用Virtual (Forced)，Virtual (Tolerance)或者Real and Virtual (Tolerance)操作合并一组边，GAMBIT允许用户通过Preserve first edge shape选项指定合并操作产生的边的位置和形状（见下面）。
gamit学习总结汇总
. 打开LINUX系统，登陆。
（选择安装语言包，关键更新)重启。
语言包的设置在system settings-language support- install languages 选择chinese（simplified）点击apply changes 等待系统自动安装，重启，重复步骤进入support 选择系统和菜单语言汉语ok（文件夹尽量保持英文状态，以便在终端中输入）；2.GAMIT安装包中没有包含海洋潮汐模型otl.grid 和map——vmf1函数模型需要自己下载连接。
otl.grid在ftp:///archive/garner/gamit/tables/ 中下载（约696M）vmf1在ftp:////pub/GRIDS/ 中下载（约220M）如没有下载otl.grid模型可以通过设置禁用该模型。
模型的禁用：打开/opt/gamit/tables/sestbl. 将use otl.grid的之改为N，表示禁用海洋潮汐模型，同时将Tides spplied 改为23(原始为31）。
独立下载模型的连接：进入/opt/gamit/tables 运行终端运行：ls -al 进行连接查看otl与map的连接情况rm otl.grid 移除已有的连接ln -s /opt/otlmap/otl_FES2004.grid otl.grid 进行连接（/opt/otlmap/otl_FES2004.grid为模型的路径）rm map.gridln -s /opt/otlmap/vmf1grd.2012_063 map.grid (map.grid 目前有vmf提供模型）ls -al 连接并查看是否已连接成功。
连接结束。
（连接完成后安装更新时会自动连接）其它otl模型：otl_FES99.grid 182.7 Mbotl_CSR4.grid 45.6 Mbotl_GOT00.grid 45.6 Mbotl_NAO99b.grid 87.1 MbAtmosphere模型也许年更新：atml files 185 Mb/year3.sites.defaults 准备介绍sites.defaults文件是站处理控制文件。
GAMIT使用手册
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中就可以对任一参数强化约束。
GAMIT使用手册 (1)
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中就可以对任一参数强化约束。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
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，查看程序是否已正确安装。
(注：在root用户下安装完后,回到自己用户不能使用，原因:需要修改软件的权限（chmod 777 软件名)数据准备阶段（个人理解：这一过程可以先在windows下准备好，到时候再拷到linux的工程文件夹）1.观测文件标准化：1）观测文件o文件的文件名要小写，观测文件里的点名一般要大写。
2）认真检查观测文件内接收机类型、天线类型并查看GAMIT的ant.dat/rcvant.dat中是否有相应的类型，否则需要更新tables中的文件。
3）天线高及天线测量方式。
2.需要从网上下载的数据文件1）/CORS/Gpscal.html查询观测0-files 文件所在的GPS周天，年积日与UTC 之间的关系。
2）导航电文文件brdcdoyn.11n、精密星历*.sp3文件以及用到的IGS跟踪站O-files数据。
(导航电文文件下载autodoyn.11n，或者brdcdoyn.12n；注意：文件名都要改成小写)3.一些常见表文件说明dcb.dat:码相关型接收机伪距改正参数统计表；(缺少这个文件基线处理最后一部会出错)antmod.dat：天线相位中心参数文件；rcvant.dat：接收机及天线名称对照表；guess_rcvant.dat: (如果rcvant.dat中没有rinex头文件中的接收机或天线型号，将从该文件中获取，用户可以自己修改)（该文件在处理过程中会自己产生，一般没用不需要准备）otl_FES2004.grid: GAMIT10.4中的潮汐改正文件，696M注：地球固体潮改正和海洋负荷潮改正：(详见《卫星导航定位新技术及高精度数据处理方法》p7-p8)地球固体潮改正：摄动天体（月亮、太阳）对弹性地球的引力使地球表面产生周期性的涨落，称为固体潮现象。
在小于100KM的GPS相对定位中，两个测站的固体潮影响几乎是相同的，在差分过程中可抵消，因此可不考虑此项改正。
海洋负荷潮改正：由海洋潮汐的周期性涨落所引起的。
与固体潮类似，但小一个数量级。
需要每周更新的：ut1.：地球自转参数；ln –s ***pole.：极移参数；ln –s ***需要每年更新的：leap.sec：跳秒表；soltab.：太阳星历；ln –s ***luntab.：月亮星力；ln –s ***nutabl.：章动参数表；ln –s ***gdetic.dat：各种大地坐标系参数；no: ln –s ***grid.otc: ln –s ***（这个好像有问题，GAMIT10.4中不是这个文件）？svs_exclude.dat：需要剔除卫星列表；(没有的话，处理过程中会提示waring，但也可以计算)（注：应该是每次计算都需要更新的，个人理解）vg_in：测站坐标速度表；autcln.cmd：（注：如果没有的话，数据处理的最后一步csh btest1.bat不能正常进行）需要重点准备的：sestbl：数据处理参数设置表；sittbl：测站约束表；：测站信息表；lfile.：站初始坐标（球坐标）；（自己制作）文件准备阶段（个人理解：在linux下为GAMIT基线解算准备文件）1.文件结构整理test：工程目录，在其下面包括：1）tables/存放各种表文件；2）igs/存放精密星历*.sp3文件；3）rinex/存放观测数据文件*doyn.11文件；4）brdc/存放广播星历文件brdcdoy0.11n文件；5）doy1/ doy2/ ……/存放各个时段解算的文件及链接文件。
2.需要重点准备的4个文件1) :记录各测站的站点、站点名称、开始观测时间、结束观测时间、天线高、天线量测方法、接收机类型、软硬件版本号、天线类型等信息。
严格按照已经给出的跟踪站格式编辑。
(注意：要现在当前目录下拷入文件，把内容删掉留下前面的几行，否则后面的步骤没法进行)具体操作方法有二：1、手动操作，按照相应文件格式对齐操作；可以从样本格式开始编辑。
2、自动提取观测文件的有关信息，生成 文件。
在各时段文件目录下操作：sh_upd_stnfo –files site1doy0.yyo site2doy0.yyo –orbt igsf –u –ref ../tables/运行上述命令后，就自动从观测文件0-files中提取重要信息，生成 文件。
我对文件的准备（武大郭博峰）：1、rinex文件夹下新建文件夹all_station(存放所有点)2、工程观测文件(o-files)复制到all_station，框架网（o-files）复制到all_station3、进入Terminalcd test/rinex/all_stationsh_upd_stnfo –files *.11o –expt test –orbt igsf –u –ref 4、查看文件是否正确gedit 5、将得到的复制到lfiles文件夹下2）l-file：记录各测站的近似坐标: Gamit要求是球面坐标，即l-file。
将所有o 文件整理在一起，然后分别运行如下命令:grep POSITION *.11o > lfile.rnxrx2apr lfile.rnx 2011 250 (注：根据所有的o文件中的x、y、z取平均值)gapr_to_l lfile.rnx lfile. “ “ 2011 250我对l-file文件的准备（武大郭博峰）：1、在rinex文件夹下新建文件夹all_data(存放所有点)2、进入该文件夹：cd test/rinex/all_data3、提取POSITIONgrep POSITION *.11o >lfile.xyz4、查看lfile.xyz文件是否正确gedit lfile.xyz5、rx2apr lfile.xyz 2011 250可得到lfile.xyz.apr文件（空间直角坐标）(注：根据所有的o文件中的x、y、z取平均值)6、lfile.xyz.apr转到lfile. (球面坐标)gapr_to_l lfile.xyz.apr lfile. “ ” 2011 2507、将得到的lfile.xyz.apr复制到lfile.文件夹下3）sestbl：（该表是gamit的关键，需要仔细学习）从tables文件夹下复制到lfiles中。
Gamit基线处理时的模型选择配置文件。
需要仔细参照技术文档和程序代码，一般情况下采用默认值，通常需要修改的地方有：Choice of Observable（观测值的选择类型）、Choice of Experiment（基线处理类型）两项。
这个模型参数的设置只要是对处理方式进行选取和误差改正模型的选取。
要是设定了误差模型改正，那在tables表里需要有这些模型改正文件，同时也需要将这些文件链接到各个观测时段文件里。
否则在处理的时候会出错的。
4）sittbl工程网和框架网需要分开生成，例如sittbl_kjw和sittbl_gcw。
也是从tables文件夹下复制到lfiles中。
各测站的精度控制指标。
一般情况下采用默认值。
对高精度的已知坐标强约束待求点坐标松弛约束。
是否需要igs跟踪站点的文件做约束吗？3.链接各个测站观测数据（包括brd c*.11n及igs*.sp3文件）在各时段文件夹下建立观测数据连接，在各时段文件夹下输入命令：ln –s ../rinex/*doy？.yyo ./ln –s ../igs/*igsweek.sp3 ./ln –s ../brdc/* doy？.yyn ./4.将tables表文件中需要的文件（包括luntab、soltab、nutabl、pole、ut1、pmu文件，还有一些误差改正模型文件等）链接到各时段文件夹。
输入以下命令建立链接：（注：ln –s 是软连接的意思）ln –s ../tables/luntab.2011.J2000 . /luntab.ln –s ../tables/soltab.2011.J2000 ./soltab.ln –s ../tables/nutabl.2011 . /nutabl.ln –s ../tables/no . /pole.ln –s ../tables/no . /ut1.ln –s ../tables/no . /pmu.另外，潮汐改正文件也是很容易出错的文件。
以10.4版本为例，默认的潮汐改正文件名为otl.grid（链接前是otl._FES2004.grid）。
一定要检查一下otl.grid 是否是真实链接，出现死链接的情况需要手动链接。