地球物理大地测量
大型科学
计算平台

CASM.jpg

地球自转参数ERP潮汐激发效应预报算法Fortran代码

ERP潮汐效应预报算法Fortran代码

 二维码 141
发表时间:2024-12-26 10:54作者:章传银网址:http://www.zcyphygeodesy.com
文章附图

[计算目标]

    给定时间跨度和采样间隔,计算地球自转参数(自转极移和日长变化)的长周期和短周期(周日/半日)潮汐效应预报时间序列。

    此功能调用IERS协议2010中收录的模型和Fortran代码模块实现,这些代码可从IERS网站获得。

幻灯片30.JPG幻灯片31.JPG

[测试入口程序]

    ERPeffectforcastexcitedbytide.f90

    输入参数knd-计算的EPR潮汐效应类型。knd=0长周期潮汐效应,knd=1短周期(周日/半日)潮汐效应。

    输入参数bgntm, endtm, tmdlt-预报时间序列开始日期、结束日期(ETideLoad长整型格式)和时间间隔(天小数)。

    输出预报时间序列文件reslt.txt。

    当knd=0时记录格式:预报日期(约定长整型),相对于第一个历元的天数,自转极移m_1(μas),m_2(μas),日长变化(μs/day),UT1变化(ms),有效角动量变化χ_1(mas),有效角动量变化χ_2(mas),自转速率变化(10-14rad/s)。

    当knd=1时记录格式:预报日期(约定长整型),相对于第一个历元的天数,自转极移m_1(μas),m_2(μas),日长变化(μs/day),UT1变化(μs)。

[调用模块]

(1)自转极移长周期海潮激发效应计算模块

    polarlongtide(jd,pole)

    输入jd-自JD2000.0历元时刻起算的计算历元时刻的世纪数(实数)。

    返回pole(7)-计算时刻pole(1:2)自转极移参数的长周期海潮效应(as,角秒),pole(5:6)自转极移有效角动量的长周期海潮效应(as,角秒)。

(2)日长及自转速率带谐潮汐效应计算模块

    RG_ZONT2(jd, cor_ut1, cor_lod, cor_w)

    输入jd-自JD2000.0历元时刻起算的计算历元时刻的世纪数(实数)。

    返回cor_ut1-计算时刻UT1带谐潮汐效应(s,秒)。

    返回cor_lod-计算时刻日长LOD带谐潮汐效应(s/day,秒/天)。

    返回cor_w-计算时刻地球自转速率带谐潮汐效应(e-14rad/s)。

(3)地球自转参数周日半日海潮效应计算模块

    PMUT1_OCEANS (mjd,cor_x,cor_y,cor_ut1,cor_lod)

    输入mjd-自JD2000.0历元时刻起算的计算历元时刻的天数(实数)。

    返回cor_x, cor_y-计算时刻自转极移参数短周期海潮效应(as,角秒)。

    返回cor_ut1-计算时刻UT1短周期海潮效应(s,秒)。

    返回cor_lod-计算时刻日长LOD短周期海潮效应(s/day,秒/天)。

(4)地球自转参数的日月潮汐效应计算模块

    PM_GRAVI(mjd,cor_x,cor_y)

    输入mjd-自JD2000.0历元时刻起算的计算历元时刻的天数(实数)。

    返回cor_x, cor_y-计算时刻自转极移参数日月效应(as,角秒)。

(5)时间系统转换包

    CAL2JD (IY0, IM0, ID0, DJM, J);JD2CAL(DJ1, DJ2, IY, IM, ID, FD, J)

    tmcnt(tm, iyr, imo, idy, ihr, imn, sec)

    mjdtotm(mjd, ltm); tmtostr(tm, tmstr)

(6)其他辅助模块

    EanesERP(mjd,X,Y,UT1,N,rjd_int,x_int,y_int,ut1_int)

    FUNDARG(jd,L,LP,F,D,OM); LAGINT (X,Y,n,xint,yout)

[编译连接]

    Fortran固定格式代码,任何fortran编译器,无需任何外部连接库。

[算法公式]ETideLoad4.5说明书

    8.6.3地球自转参数潮汐效应及其计算

附件zip压缩包:visual studio_intel fortran 集成环境测试项目、DOS可执行测试程序、地球物理模型和全部测试输入输出数据。