本文为小白笔记，如有错误请指正！

1，运行cpfxyz程序

将网站公布的cpf（ftp://cddis.gsfc.nasa.gov/pub/slr/）复制到文件中，第11行代表2014年01月02日0刻的坐标XYZ.

注：05年后是CPF，之前是IRV

Cpf与IRV不同，CPF只有XYZ没有速度，因此CPF需要运行CPFXYZ和matlab程序，而IRV数据只需要运行GST.f90程序即可，因为不需要插值速度。

下面为IRV格式（含有速度）。

运行CPF.F90程序，将生成的结果(即PX,PY,PZ)，放到matlab中，运行script1.m文件，生成相应的速度,见42行

作为下面GST程序中18行需要修改的速度值。速度文件

注：运行matlab时，要设置a值，根据卫星（怎么判断呢？）

2，运行gst2000.f90

从9行到18行都需要修改

根绝IERS网站公布的Bulletins，里面有相应的参数值。

https://www.iers.org/IERS/EN/Publications/Bulletins/bulletins.html

程序中的57388是MJD(儒略日)，需要相应改掉

跳秒，对应程序中13行的36，该文件中是37

UT1-UTC是用下面的还是上面的？（对应程序中10行斜杠前面的数值）

前面代表18年2月16日。。。

第四列代表儒略日，后面的X和Y为极移XP和YP（十四行和十五行斜杠前面的）

地心地固坐标系→准地固坐标系→瞬时真赤道坐标系（GEODYN需要的）→ J2000惯性坐标系

                 极移改正                 自转                                                    章动、岁差