在看了陳超老師的視頻教學後總結出來的,在globk_comb.cmd和glorg_comb.cmd文件主體上附加上解釋。
* GLOBK command file to generate daily time series and to combine
%命令控制文件 用於產生時間序列和多日解合併
* h-files over 2 to 30 days.
* For combination, set COMB as a globk command-line option to
* invoke the saving of the output h-file
* Last edited by rwk 150618
* << column 1 must be blank if not comment >>
* This group of commands must appear before any others:
srt_file @.srt
srt_dir +1
%對時間進行排序,-1表示降序
eq_file ~/gg/tables/itrf08_comb.eq %地震文件
# Optionally add a second eq_file for analysis-specific renames
* End commands that must appear first
* ITRF2008 augmented by now-defunct sites and recent IGS solutions;
# matched to itrf08_comb.eq
apr_file ~/gg/tables/itrf08_comb.apr %近似座標文件採用itrf08 如有igs站則起作用
# Optionally add additional apr files for other sites
x apr_file ../../tables/regional.apr
* Set maximum chi2 and prefit coordinate difference (m),for an h-file to be used;
%卡方增量參數:組合新的文件中卡方增量的最大值默認值100;座標的擬合殘差最大限差默認值100000m(先驗座標比較大時此值可以很小)單位是m;最後一個默認值爲10000,慣用30 50 2000;
max_chii 13 3
# increase tolerances to include all files for diagnostics
x max_chi 100 5.0
# Not necessary unless combining h-files with different a priori EOP
%地球自轉參數先驗值文件 需要手動連接(非必須)
in_pmu ../tables/pmu.usno
* Invoke glorg %調用glorg命令文件 合併多日解時不需要 可以註釋此行
org_cmd glorg_comb.cmd
* Print file options %打印文件
crt_opt NOPR(不輸出該文件,一般不用 可註釋此行)
prt_opt NOPR GDLF(輸出卡方增量列表) CMDS(結果文件中加上命令文件內容) MIDP
org_opt PSUM(顯示座標值得到點位座標,VSUM表示得到速度信息一般時間跨度在2年及以上 BLEN輸出基線長度BRAT輸出基線變化率) CMDS GDLF MIDP FIXA RNRP
# sh_glred will name the glorg print files %glorg文件名
x org_out globk_comb.org
* Coordinate parameters to be estimated and a priori constraints%對測站座標進行約束
apr_neu all 10 10 10 0 0 0 分別爲n e u vn ve vu 0爲不輸出(如果上面要求輸出速度 則此處必須設置非0)10爲鬆約束 一般igs站 先驗座標精確 可以進行緊約束
可以對空間直角座標進行約束:apr_site urum 10 10 10 0 0 0 後面參數與上述相同
* Rotation parameters to be estimated and a priori constraints
apr_wob 10 10 0 0 %全球站點解算時採用10 10 1 1 給鬆弛約束 如果是區域網則強約束例如:0.2 0.2 0.02 0.02 見下
apr_ut1 10 0
%全球站點解算時採用10 1
# If combining with global h-files, allow EOPS to change
# between days
x mar_wob 3650 3650 365 365
x mar_ut1 365 365
# EOP tight if translation-only stabilization in glorg %區域網則爲緊約束
x apr_wob .25 .25 .1 .1
x apr ut1 .25 .1
* Write out a combined H-file
# Can substitute your analysis name for 'COMB' in the file name below %合併多日解文件
COMB out_glb H------_COMB.GLX
* Optionally put a uselist and/or sig_neu and mar_neu reweight in a source file
x source ../tables/uselist
x source ../tables/daily_reweights
* Turn off quake log estimates if in the eq_file
%是否在地震文件中關閉估計值
free_log -1
* Remove scratch files for repeatability runs
%是否刪除臨時文件
del_scra yes
* Correct the pole tide when not compatible with GAMIT
app_ptid all
* If orbits free in GAMIT (RELAX) and you want them fixed, use:
x apr_svs all F F F F F F FR
* but if you are combining with globk h-files, better to leave them
* on but, if the models are incompatible, turn off radiation-pressure parameters,
x apr_svs all 100 100 100 10 10 10 0R
* When using MIT GLX files which have satellite phase center positions
* estimated use:
apr_svan all F F F ! Fix antenna offset to IGS apriori values.
* Glorg command file for daily repeatabilities or combinations
* Last edited by rwk 130701
* Parameters to be estimated %選用的參考框架轉換的方式是平移加旋轉
pos_org xtran ytran ztran xrot yrot zrot
# or if translation-only %或者是隻有平移
x pos_org xtran ytran ztran
* Downweight of height relative to horizontal (default is 10)%控制高度的相對權重
# Heavy downweight if reference frame robust and heights suspect
x
cnd_hgt 1000
* Controls for removing sites from the stabilization
# Vary these to make the stabilization more robust or more precise%迭代次數
stab_it 4 0.8 3.0
x stab_it 4 0.5 4.0
* A priori coordinates
# ITRF2008 may be replaced by an apr file from a priori velocity solution
apr_file ~/gg/tables/itrf08_comb.apr
x apr_file ../../tables/regional.apr
* List of stabilization sites
# This should match the well-determined sites in the apr_file
%穩定站列表
stab_site clear先清除所有穩定站
stab_site xxxx xxxx可以自己設置穩定站
source ~/gg/tables/igb08_hierarchy.stab_site%默認將改文件裏面的站點設爲穩定站
x source ../../tables/regional_stab_site