[GAMIT/GLOBK學習筆記]globk_comb.cmd/glorg_comb.cmd文件詳解

在看了陳超老師的視頻教學後總結出來的,在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


發表評論
所有評論
還沒有人評論,想成為第一個評論的人麼? 請在上方評論欄輸入並且點擊發布.
相關文章