地球自转参数ERP潮汐激发效应预报算法Fortran代码ERP潮汐效应预报算法Fortran代码 二维码
141
发表时间:2024-12-26 10:54网址:http://www.zcyphygeodesy.com ![]() [计算目标] 给定时间跨度和采样间隔,计算地球自转参数(自转极移和日长变化)的长周期和短周期(周日/半日)潮汐效应预报时间序列。 此功能调用IERS协议2010中收录的模型和Fortran代码模块实现,这些代码可从IERS网站获得。 [测试入口程序] 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可执行测试程序、地球物理模型和全部测试输入输出数据。 |