Introduction to GAMIT/GLOBK Applications of GLOBK. Lecture 11 OVERVIEW

Transcription

1 Itroductio to GAMIT/GLOBK Applicatios of GLOBK Lecture 11 GAMIT/GLOBK Lec11 1 OVERVIEW o I this lecture we cover: o Basic types of aalyses with globk l Velocity ad repeatability rus o GLOBK acillary programs l GLRED l GLIST l GLSAVE l HFUPD l GLBTOSNX GAMIT/GLOBK Lec11 2

2 Basic types of globk aalyses o Geerally globk is used i two ways: Combiig may years of data to geerate velocities solutios Aalyzig groups of idividual experimets to assess positio quality o For example: After may years of processig, coordiate systems are likely to chage (eg., ITRF94 to ITRF2000) ad globk ca brig all early results ito a cosistet frame Reaso we write "loose" hfiles from gamit GAMIT/GLOBK Lec11 3 Velocity estimatio o A gdl file is created by makig a list of biary-hfiles. If ruig from directory gsol, the ls../glbf/h*.glx > idia.gdl o (.glr files are bias free solutios,.glx are bias fixed. Siex files are.gls) o If multiple years of data, the try velocity solutio GAMIT/GLOBK Lec11 4

3 Globk velocity commad file o Mai feature of i globk commad file would be l apr_eu all l where the is velocity ucertaities i m/yr l If velocity relative to specific site the add l apr_eu iisc_gps F F F l would fix velocity of iisc l Usually do ot do this, but rather use glorg for geeralized defiitio GAMIT/GLOBK Lec11 5 GLRED o To check the quality of the idividual positios i the hfiles we could create gdl files with just sigle hfiles listed ad ru globk o each oe usig a commad file with apr_eu all apr_eu iisc_gps Agai iisc is costraied ad ormal process o Ruig globk o a set of idividual hfiles is tedious so we have a program that automatically does this: glred GAMIT/GLOBK Lec11 6

4 GLRED o glred is ru the same way as globk ie. l glred <crt> <prt> <log> <gdl file> <commad file> l but rather tha combiig all the data i the gdl files, each file is processed separately by globk o glred is a coveiet way of makig multiple globk rus GAMIT/GLOBK Lec11 7 Geeral format of gdl files o Geeral format of a gdl file is: [biary hfile ame] <scale> <diagoal scale> <+> o all except the file ame are optioal o scale is a multiplier applied to the covariace matrix o diagoal scale scales the diagoal (oly used for poorly coditioed siex files) GAMIT/GLOBK Lec11 8

5 GLRED gdl file o There are some optios i the gdl file that add flexibility to glred o a + ca be added at the eds of hfile ame ames so that groups of hfiles will be combied eg.../glbf/h991001_gps1.glx +../glbf/h991002_gps1.glx../glbf/h000601_gps1.glx +../glbf/h000602_gps1.glx o Would combie the first ad secod ad third ad fourth files. GAMIT/GLOBK Lec11 9 GLIST o To determie what is cotaied i a set of biary files use glist glist <gdl file> <output_file> <sort_directio> [eq_file] [out gdl] [apr file] Oly first argumet required addig eq_file will apply earthquakes ad reames out_gdl is a time sorted gdl file which ca be used i globk apr_file is a optioal apr_file which the apr positios i the biary hfiles will be compared to. Differeces greater tha 0.3 m are reported. GAMIT/GLOBK Lec11 10

6 GLIST o Durig tutorial we will use this program to check cotets of local hfiles o A grep 'P$' will extract a list of positios of the sites ad the duratio of their use GAMIT/GLOBK Lec11 11 GLSAVE o Allows solutios geerated by globk to be saved as biary hfiles o Normally ivoked i globk commad file with the out_glb <ame/structure> where ame/structure is either: Expliciit ame of file wild card costructor: Forms H@.GLX will use the gdl file ame replacig the first letter ad extet H------_SIO.GLX will use the 2-7 characters from the gdl file ame GAMIT/GLOBK Lec11 12

7 GLSAVE l if the com_file commad was used i globk, the glsave ca be ru as separate program l glsave <com file> [out global ame] [descriptio] <sol_file ame> l where com file is the ame of the file give i com_file commad i globk l out global ame is ame for output biary hfile l descriptio is strig descriptio eclosed i ' ' GAMIT/GLOBK Lec11 13 GLSAVE sol_file o sol file ame is a optioal ew solutio file (ot usually used) o Geerated i glorg where the glorg solutio file ca be output o For example: with coordiate system defiitio applied ad costraits o earby sites o Advaced feature which should be carefully used. GAMIT/GLOBK Lec11 14

8 HFUPD o hfupd allow updates to be made to the biary hfile files o hfupd <optios> <list of hfile ames> o Optios iclude: o -s <siex file/statio.ifo> ame of ew statio.ifo or siex header file that cotais ew iformatio about sites (eg. correct atea height) o -e <edit file> file cotaiig edited ad reamed sites GAMIT/GLOBK Lec11 15 HFUPD optios o -r report cotets of statio.ifo/siex file o -u update the hfiles (must give to actually update files, without this optio program reports chages to be made but does chage files) o -p <hf/pmu_file> applied pole tide correctio if eeded. uses either polar motio values i hfile or i pmu_file o other optios are advaced features GAMIT/GLOBK Lec11 16

9 HFUPD edit file o Edit file cotais lies of form EDIT <NAME> <HF code> <Epoch rage> o deletes site NAME i hfile cotaiig <HF code> i ame for period give by epoch rage o RESTORE <NAME> <HF code> <Epoch rage> o will restore a site provided ot a combied biary file GAMIT/GLOBK Lec11 17 HFUPD edit file o RENAME <OLD NAME> <NEW NAME> <HF code> <Epoch rage> o will reame a site from OLD to NEW o Epoch rage i form start (yr, mth, day, hr, mi) ed (yr, mth, day, hr, mi) o start must before start of data i hfile ad ed must after ed of data i hfile GAMIT/GLOBK Lec11 18

10 GLBTOSNX SINEX (Solutio idepedet exchage format) is iteratioal stadard for exchagig solutios of the types cotaied i biary hfiles glbtosx <dir> [-h] <commets file> <iput biary hfile> <output file ame> dir is directory for ame -h will write oly headers ad ot full covariace matrix commets files is file cotaiig commets (see $HELP_DIR/head.sx output file ame cotrols structure of geerated file ame GAMIT/GLOBK Lec11 19 Geeral items o Whe globk program are ru usig ' ' for a iput will ofte use a default value o Called a place holder whe specific iput required o eg. glist file.gdl file.glist '' '' out.gdl o will use default sort directio (+1) ad o earthquake file GAMIT/GLOBK Lec11 20

11 Creatig apr_file o apr_files are supplied i gg/templates o They ca be created from globk output with l grep '^Uc. ' <prt/org file> awk '{prit substr($0,5)}' >! ew.apr o Satellite ephemerides ca be obtaied with grep '^Eph. ' o Program plate ad corcom ca also be used to geerate apr files GAMIT/GLOBK Lec11 21 Summary o Examied basic types of globk rus o The exteral programs used with globk to l List the cotets of h-files l Save combied h-files l Ru separate globk aalyses l Update hfiles with ew iformatio l Covert h-files to SINEX files GAMIT/GLOBK Lec11 22