国际参考电离层

一. 源码获取

源码地址见 irimodel, 请按图1所示下载所需文件。

图1. IRI-2020 Fortran源码下载

模型在线运行: CCMC Instant Run System

您也可以通过构造POST请求来获取数据

REM computer system: Windows 10
@echo off

set URL=https://kauai.ccmc.gsfc.nasa.gov/instantrun/api/iri/2020/
set OUTPUT_FILE=temp.txt
set POST_DATA={...}(json格式)

wget --header="Content-Type: application/json" --limit-rate=600k --post-data="%POST_DATA%" %URL% -O %OUTPUT_FILE%

二. 程序参数

1. 时间信息

# iy(year),imd(mmdd or -ddd),iut(=0/1,LT/UT),hour(hour with LT/UT)
iy = 2023
imd = -001  # mmdd = imd (mmdd in subroutine iri_web, 下同)
iut = 1
hour = 12  # phour = hour

2. 位置信息

# jm(=0/1,geog/geom),xlat(latitude/deg),xlon(longitude/deg)
jm = 0  # jmag = jm
xlat = 0.0
xlon = 0.0

# hx(height/km)
hx = 300 # hxx = hx

# lower and upper height [km] for TEC integration (valid range: 65-30,000 km)(0,0 for no TEC) 
htec_min = 65
htec_max = 2000

3. 循环选项

# ivar=1: height
# ivar=2: latitude
# ivar=3: longitude
# ivar=4: year
# ivar=5: month
# ivar=6: day of month
# ivar=7: dat of year
# ivar=8: hour
ivar = 3
vbeg = 0
vend = 360
vstp = 5

4. 输入选项

# Enter 0 to use standard or 1 to enter your own
jchoice = 0

5. 输出选项

# Enter 0 for standard table of IRI parameters
# Enter 1 for list of peak heights and densities
# Enter 2 for plasma frequencies, B0, M3000, valley, width and depth
# Enter 3 for 6 parameters of your choice
# Enter 4 for D-region models at 60,65,..,110 km
# Enter 5 special test output
# if ivar=1 then choose 0, 3, 4, or 5
piktab = 0

提示

1. Fortran源码采用F77固定格式编写, 为便于修改, 可将iritest.for转换为F90格式
2. 如果只计算TEC数据，可以设置jf(2)=.false. & jf(3)=.false.以提高运算速度！

三. 数据演示

图2. IRI-2020 数据演示