sh_gamit 批处理总结

首先建立必要的目录结构，在工程名下建立3个文件名，分别是：rinex,igs,brdc 然后拷贝gamit目录下的tables 文件夹，如下：
linux-zhao% mkdir rinex igs
linux-zhao% mkdir brdc
linux-zhao% cp -R /home/zhao/gg/tables/ ./
cp: 无法打开 “/home/zhao/gg/tables/igs_01_antex11.pcv” 读取数据: 权限不够
linux-zhao%
（说明：我是在工程目录下操作的。）(提醒：要检查tables文件夹中是否缺文件，文件中是否有数据，例如ut1.文件。我在拷贝后，竟然这个文件是空的！原因不明。所以检查必要表文件很重要。还有就是海洋潮汐负载文件，要自己下，然后自己连接到tables中，具体方法见其他文章，这个连接最好是绝对连接，保证随时可用！)
2、下载数据，为了熟悉流程，只处理中国的IGS站其中3个，有上海shao,昆明kunm，乌鲁木齐南山站guao。
观测文件、广播星历，精密星历都到cddis.gsfc.nasa.gov下载，其中精密星历到gps/products/年下找，所需GPS周，（利用doy 程序：doy 年 年纪日 ），广播星历在/pub/gps/data/brdc下。观测o文件很好找了，自己摸索一下吧。
3、修改sites.defaults:添加所需要的测站，把不学要的＃掉。
修该process.defaults:我的没有该，等熟练后在修改吧。
修改sestbl. 文件：时段处理的设置，例如所选模型等。我的没有修俄该。
修改 sittbl. 文件：把新添加的站加进去，约束站点。我这里采用0.05,0.05,0.05。

sites.defaults如下：
# individual station commands will override all_sites if they follow it      
# mate_gps tubi ftprnx xstinfo glrepu glreps  
# ankr_gps tubi ftprnx xstinfo glrepu glreps glts
# kit3_gps tubi ftprnx xstinfo glrepu glreps
# nssp_gps tubi ftprnx xstinfo glrepu glreps
# graz_gps tubi ftprnx xstinfo glrepu glreps  
# noto_gps tubi ftprnx xstinfo glrepu glreps  
# nico_gps tubi ftprnx xstinfo glrepu
# sofi_gps tubi ftprnx xstinfo glrepu
# zeck_gps tubi ftprnx xstinfo glrepu
# zwen_gps tubi ftprnx xstinfo glrepu glreps
# bahr_gps tubi ftprnx xstinfo glrepu
# tela_gps tubi ftprnx xstinfo glrepu       
# tubi_gps tubi ftpraw xstinfo glrepu        glts
# mert_gps tubi ftpraw         glrepu        glts
# kant_gps tubi ftpraw         glrepu        glts
# ttth_gps tubi xstinfo xsite:1999_256-1999_278 glreps xsite:1999_300-1999_365
# thht_gps tubi xstinfo xsite glreps
shao_gps   2008
kunm_gps 2008
guao_gps 2008

sittbl.如下：
SITE              FIX    --COORD.CONSTR.--
      << default for regional stations >>
ALL               NNN    100. 100.   100.    
       << IGS china stations >>
shao shao_gps     NNN    0.050 0.050 0.05
kunm kunm_GPS     NNN    0.050 0.050 0.05
guao guao_GPS     NNN    0.050 0.050 0.05
        << IGS core stations >>
VILL VILL_GPS     NNN    0.050 0.050 0.05
MAS1 MAS1_GPS     NNN    0.050 0.050 0.05
FORT FORT_GPS     NNN    0.050 0.050 0.05
BRAZ BRAZ_GPS     NNN    0.050 0.050 0.05
KOUR KOUR_GPS     NNN    0.050 0.050 0.05
OHIG OHIG_GPS     NNN    0.050 0.050 0.05
BRMU BRMU_GPS     NNN    0.050 0.050 0.05
THU1 THU1_GPS     NNN    0.050 0.050 0.05
SANT SANT_GPS     NNN    0.050 0.050 0.05
AREQ AREQ_GPS     NNN    0.050 0.050 0.05
WES2 WES2_GPS     NNN    0.050 0.050 0.05
GODE GODE_GPS     NNN    0.050 0.050 0.05
ALGO ALGO_GPS     NNN    0.050 0.050 0.05
NLIB NLIB_GPS     NNN    0.050 0.050 0.05
MDO1 MDO1_GPS     NNN    0.050 0.050 0.05
PIE1 PIE1_GPS     NNN    0.050 0.050 0.05
YELL YELL_GPS     NNN    0.050 0.050 0.05

4、把igs目录下的精密星历解压，好像不手动的话，会报报错。

然后运行如下命令：
linux-zhao% sh_gamit -expt 2008 -d 2008 201 202 203 -orbit IGSF -copt x k p -dopt c ao -noftp -remakes Y > & ! sh_gamit.log

5、运行：
linux-zhao% sh_glred -s 2008 201 2008 203 -expt 2008 -opt H G E
Input options -s 2008 201 2008 203 -expt 2008 -opt H G E
glrepuse:
glrepstab:
glveluse:
glvelstab:
gltimeser:
Checking and making required directories
~/project/2008/glbf ~/project/2008
Htoglb converting GAMIT ascii h-files to GLOBK binary h-files
htoglb /home/zhao/project/2008/glbf /home/zhao/project/2008/tables/svt /home/zhao/project/2008/glbf/h????a.08201
htoglb /home/zhao/project/2008/glbf /home/zhao/project/2008/tables/svt /home/zhao/project/2008/201/h????a.08201
~/project/2008
Running the GLRED combination
~/project/2008/gsoln ~/project/2008
rm: 无法删除 “binhlist.glr”: 没有那个文件或目录
rm: 无法删除 “binhlist.glx”: 没有那个文件或目录
rm: 无法删除 “binhlist.GLX”: 没有那个文件或目录
ls: No match.
ls: No match.
H-files included in this combination: /home/zhao/project/2008/glbf/h0807191200_2008.glx
Input options -eqf none -aprf /home/zhao/project/2008/tables/itrf00.apr -pmuf /home/zhao/project/2008/tables/pmu.usno -globkf /home/zhao/project/2008/gsoln/globk_comb.cmd -glorgf /home/zhao/project/2008/gsoln/glorg_comb.cmd -hfname 2008 -netlist ALL -repuse -repstab
glred 6 globk_2008_08201.prt globk_2008_08201.log globk_2008_08201.gdl globk_comb.cmd
~/project/2008
~/project/2008/glbf ~/project/2008
Htoglb converting GAMIT ascii h-files to GLOBK binary h-files
htoglb /home/zhao/project/2008/glbf /home/zhao/project/2008/tables/svt /home/zhao/project/2008/glbf/h????a.08202
htoglb /home/zhao/project/2008/glbf /home/zhao/project/2008/tables/svt /home/zhao/project/2008/202/h????a.08202
~/project/2008
Running the GLRED combination
~/project/2008/gsoln ~/project/2008
ls: No match.
ls: No match.
H-files included in this combination: /home/zhao/project/2008/glbf/h0807201200_2008.glx
Input options -eqf none -aprf /home/zhao/project/2008/tables/itrf00.apr -pmuf /home/zhao/project/2008/tables/pmu.usno -globkf /home/zhao/project/2008/gsoln/globk_comb.cmd -glorgf /home/zhao/project/2008/gsoln/glorg_comb.cmd -hfname 2008 -netlist ALL -repuse -repstab
/home/zhao/project/2008/gsoln/globk_comb.cmd already exists. NOT overwritten
/home/zhao/project/2008/gsoln/glorg_comb.cmd already exists. NOT overwritten
glred 6 globk_2008_08202.prt globk_2008_08202.log globk_2008_08202.gdl globk_comb.cmd
~/project/2008
~/project/2008/glbf ~/project/2008
Htoglb converting GAMIT ascii h-files to GLOBK binary h-files
htoglb /home/zhao/project/2008/glbf /home/zhao/project/2008/tables/svt /home/zhao/project/2008/glbf/h????a.08203
htoglb /home/zhao/project/2008/glbf /home/zhao/project/2008/tables/svt /home/zhao/project/2008/203/h????a.08203
~/project/2008
Running the GLRED combination
~/project/2008/gsoln ~/project/2008
ls: No match.
ls: No match.
H-files included in this combination: /home/zhao/project/2008/glbf/h0807211200_2008.glx
Input options -eqf none -aprf /home/zhao/project/2008/tables/itrf00.apr -pmuf /home/zhao/project/2008/tables/pmu.usno -globkf /home/zhao/project/2008/gsoln/globk_comb.cmd -glorgf /home/zhao/project/2008/gsoln/glorg_comb.cmd -hfname 2008 -netlist ALL -repuse -repstab
/home/zhao/project/2008/gsoln/globk_comb.cmd already exists. NOT overwritten
/home/zhao/project/2008/gsoln/glorg_comb.cmd already exists. NOT overwritten
glred 6 globk_2008_08203.prt globk_2008_08203.log globk_2008_08203.gdl globk_comb.cmd
~/project/2008
Creating time series between: 2008 201 --> 2008 203
~/project/2008/gsoln ~/project/2008
sh_plotcrd -f tmp.ensum_input -s short -cols 1 -e 2008

ENSUM: Program to extract ENU components from GLOBK solution
------------------------------------------------------------

Reading stations from tmp.ensum_input                                                                                                                
Solution: Combination of ALL networks
There are       27 baseline entries between    3 sites
Generating solution summary
Processing data from site SHAO_GPS
Processing data from site KUNM_GPS
Processing data from site GUAO_GPS
Checking for bad records in input file SUM.2008
Bad records found, removing the following records from histogram:
SHAO_GPS-       N 1    3    3461996.3266 0.1134 38.1 0.19   ******** 50809.69 31.5 0.16 0.01 2008.55       1-      6
SHAO_GPS-       E 1    3   11552739.0090 1.8431 606.0 0.19   ******** ******** 603.6 0.19 0.01 2008.55      11-     16
SHAO_GPS-       U 1    3         21.7975 1.3055 421.5 0.19   ******** ******** 348.8 0.15 0.01 2008.55      21-     26
KUNM_GPS-       N 1    3    2786275.4130 0.5547 94.0 0.10   ******** ******** 120.8 0.13 0.01 2008.55      31-     36
KUNM_GPS-       E 1    3   10368693.2635 1.7886 466.7 0.15   ******** ******** 596.7 0.19 0.01 2008.55      41-     46
KUNM_GPS-       U 1    3       1985.8331 1.3581 505.1 0.21   ******** ******** 17.5 0.01 0.01 2008.55      51-     56
GUAO_GPS-       N 1    3    4839181.7230 1.0904 291.3 0.15   ******** ******** 244.8 0.13 0.01 2008.55      61-     66
GUAO_GPS-       E 1    3    7042905.2398 1.6815 468.9 0.16   ******** ******** 652.2 0.22 0.01 2008.55      71-     76
GUAO_GPS-       U 1    3       2028.5929 1.3685 317.6 0.13   ******** ******** 112.9 0.05 0.01 2008.55      81-     86

Soft kill LENGTH only :
Soft kill EAST   only :
Soft kill NORTH only :
Soft kill UP     only :
GMT Warning: Old GMT 3 .gmtdefaults file found. May not be fully compatible with GMT 4.
GMT Warning: It is recommended that you migrate your GMT 3 settings to GMT 4 settings.

Created postscript file : pshist_wrms.SUM.2008

To view      : gs pshist_wrms.SUM.2008

To print     : lpr -s pshist_wrms.SUM.2008

Checking for bad records in input file SUM.2008
Bad records found, removing the following records from histogram:
SHAO_GPS-       N 1    3    3461996.3266 0.1134 38.1 0.19   ******** 50809.69 31.5 0.16 0.01 2008.55       1-      6
SHAO_GPS-       E 1    3   11552739.0090 1.8431 606.0 0.19   ******** ******** 603.6 0.19 0.01 2008.55      11-     16
SHAO_GPS-       U 1    3         21.7975 1.3055 421.5 0.19   ******** ******** 348.8 0.15 0.01 2008.55      21-     26
KUNM_GPS-       N 1    3    2786275.4130 0.5547 94.0 0.10   ******** ******** 120.8 0.13 0.01 2008.55      31-     36
KUNM_GPS-       E 1    3   10368693.2635 1.7886 466.7 0.15   ******** ******** 596.7 0.19 0.01 2008.55      41-     46
KUNM_GPS-       U 1    3       1985.8331 1.3581 505.1 0.21   ******** ******** 17.5 0.01 0.01 2008.55      51-     56
GUAO_GPS-       N 1    3    4839181.7230 1.0904 291.3 0.15   ******** ******** 244.8 0.13 0.01 2008.55      61-     66
GUAO_GPS-       E 1    3    7042905.2398 1.6815 468.9 0.16   ******** ******** 652.2 0.22 0.01 2008.55      71-     76
GUAO_GPS-       U 1    3       2028.5929 1.3685 317.6 0.13   ******** ******** 112.9 0.05 0.01 2008.55      81-     86

Soft kill LENGTH only :
Soft kill EAST   only :
Soft kill NORTH only :
Soft kill UP     only :
GMT Warning: Old GMT 3 .gmtdefaults file found. May not be fully compatible with GMT 4.
GMT Warning: It is recommended that you migrate your GMT 3 settings to GMT 4 settings.

Created postscript file : pshist_nrms.SUM.2008

To view      : gs pshist_nrms.SUM.2008

To print     : lpr -s pshist_nrms.SUM.2008

select ALL
Limiting output sigmas to   0.0000 m Horizontal and   0.0000 m height
Created : mb_SHAO_GPS.dat1
Created : mb_SHAO_GPS.dat2
Created : mb_SHAO_GPS.dat3
Created : mb_KUNM_GPS.dat1
Created : mb_KUNM_GPS.dat2
Created : mb_KUNM_GPS.dat3
Created : mb_GUAO_GPS.dat1
Created : mb_GUAO_GPS.dat2
Created : mb_GUAO_GPS.dat3
sh_baseline -n 3 -cols 1 -o 0 -ps 2008 -u 1    -f mb_*

Running sh_baseline to create GMT plots of time series

>>>>>> scaling uncertainties by : 1
COMMENT :
creating psbase_2008.*

GMT Warning: Old GMT 3 .gmtdefaults file found. May not be fully compatible with GMT 4.
GMT Warning: It is recommended that you migrate your GMT 3 settings to GMT 4 settings.
=+=+=+=+=      9 frames      3 pages =+=+=+=+=

Plotting: GUAO
Created: psbase_2008.GUAO
Plotting: KUNM
Created: psbase_2008.KUNM
Plotting: SHAO
Created: psbase_2008.SHAO

Created postript files    : psbase_2008.*

Created baseline statistics : base_stats_psbase_2008

To view postript files    : gs psbase_2008.*

Hardcopy   postript files : lpr -s psbase_2008.*

~/project/2008
linux-zhao%