IRI output arrays
Output file OUTF(1:20,1:1000):*********************************************ALL ANGLES ARE IN DEGREEALL DENSITIES ARE IN M-3
ALL ALTITUDES ARE IN KM
ALL TEMPERATURES ARE IN KELVINALL TIMES ARE IN DECIMAL HOURS**********************************************************************

OUTF(1,*) electron density in m-3
OUTF(2,*) neutral temperature in K
OUTF(3,*) ion temperature in K
OUTF(4,*) electron temperature in K
OUTF(5,*) O+ ion density in % or in m-3 if JF(22)=false
OUTF(6,*) H+ ion density in % or in m-3 if JF(22)=false
OUTF(7,*) He+ ion density in % or in m-3 if JF(22)=false
OUTF(8,*) O2+ ion density in % or in m-3 if JF(22)=false
OUTF(9,*) NO+ ion density in % or in m-3 if JF(22)=false
and, if JF(6)=false:
OUTF(10,*) Cluster ion density in % or in m-3 if JF(22)=false
OUTF(11,*) N+ ion density in % or in m-3 if JF(22)=false
OUTF(12,*) free
OUTF(13,*) free
If JF(24)=true
OUTF(14,1:11) standard IRI electron density at 60, 65, 70, ... , 110 km
If JF(24)=false
OUTF(14,12:22) Friedrich (FIRI) model at these heights
OUTF(14,23:33) standard Danilov model (SW=0, WA=0)
OUTF(14,34:44) Danilov model for minor Stratospheric Warming (SW=0.5)
OUTF(14,45:55) Danilov model for major Stratospheric Warming (SW=1)
OUTF(14,56:66) Danilov model for weak Winter Anomaly (WA=0.5) conditions
OUTF(14,67:77) Danilov model for strong Winter Anomaly (WA=1) conditions
OUTF(15-20,*) free

Output file OARR(1:100):
Several parameters (marked #) are used for input (user input) as well as output (IRI computed).
OARR(1)=NMF2S#OARR(2)=HMF2#OARR(3)=NMF1#OARR(4)=XHMF1#OARR(5)=NMES#OARR(6)=HME#F2-peak density in m-3
F2-peak height in km
F1-peak density in m-3
F1-peak height in km
E-peak density in m-3
E-peak height in km

OARR(7)=NMD
OARR(8)=HMD
OARR(9)=HHALF#OARR(10)=B0
OARR(11)=VNER
OARR(12)=HEF
OARR(13)=ATE(2)OARR(14)=AHH(2)#OARR(15)=ATE(3)#OARR(16)=ATE(4)OARR(17)=ATE(5)OARR(18)=ATE(6)OARR(19)=ATE(7)OARR(20)=ATE(1)OARR(21)=TI1
OARR(22)=XTETIOARR(23)=XHI3
OARR(24)=SUNDECOARR(25)=DIP
OARR(26)=MAGBROARR(27)=MODIPOARR(28)=LATIOARR(29)=SAX200
OARR(30)=SUX200OARR(31)=SEASONOARR(32)=LONGI#OARR(33)=rssn
OARR(34)=COV#OARR(35)=B1
OARR(36)=xm3000
OARR(37)=TEC
OARR(38)=TEC-top%#OARR(39)=gind
OARR(40)=f1pb
OARR(41)=f107d
OARR(42)=c1
OARR(43)=daynr
OARR(44)=drift
OARR(45)=stormcorr#OARR(46)=f10781
OARR(47)=estormcor
OARR(48)=spreadf
OARR(49)=MLAT
OARR(50)=MLONG

density in m-3 of D-region inflection point
height in km of D-region inflection point
height in km used by Gulyaeva B0 model
bottomside thickness parameter in km
density in m-3 at E-valley bottom
height in km of E-valley top (Ne(HEF)=NmE)
electron temperature Te in K at AHH(2)intermediate height between 120km and 300/350km
Te at 300km/350km for BIL-1995/TBT2012+SA model
Te at 400km/550km for BIL-1995/TBT2012+SA model
Te at 600km/850km for BIL-1995/TBT2012+SA model
Te at 1400km/1400km for BIL-1995/TBT2012+SA model
Te at 3000km/2000km for BIL-1995/TBT2012+SA model
Te at 120km = neutral temperature from CIRA
ion temperature in K at 430km
altitude in km where Te=Ti
solar zenith angle at 200 km in degree
sun declination angle in degree
IGRF magnetic inclination (dip) in degree
IGRF dip latitude in degree
modified dip latitude in degree
geographic latitude in degree
time of sunrise at 200 km in decimal hours
time of sunset at 200 km in decimal hours=1 spring, =2 summer, =3 fall =4 winterassumes equal length seasons (92 days) with spring
starting at day-of-year=45
geographic longitude in degree
12-month running mean of sunspot number
12-month running mean of F10.7
Bottomside shape parameter
Propagation factor M(3000)F2
Total Electron Content in m-2 (placeholder for for IRIWeb)
TEC-top/TEC * 100(placeholder for for IRIWeb)
12-month running mean of IG indexprobability for an F1 layer
daily solar radio flux at 10.7cm (F10.7)
shape parameter for F1 layer
day of year
vertical ion drift at equator in m/s
ratio foF2_storm/foF2_quiet
81-day average of F10.7
ratio foE_storm/foE_quiet
probability of spread-F occurrence
IGRF magnetic latitude in degree
IGRF magnetic longitude in degree

OARR(51)=index_3h_ap
OARR(52)=IAP_daily
OARR(53)=invdip
OARR(54)=XMLTOARR(55)=cgm_lat
OARR(56)=cgm_lon
OARR(57)=cgm_mlt
OARR(58)=cgmlat
OARR(59)=cl(MLT=0)
OARR(60)=cl(MLT=1)……
…..
OARR(81)=cl(MLT=22)
OARR(82)=clt(MLT=23)
OARR(83)=xkp
OARR(84)=dec
OARR(85)=fl
OARR(86)=dimo
OARR(87)=SAX300
OARR(88)=SUX300#OARR(89)=HNEA#OARR(90)=HNEE#ap index for current UTdaily ap index
invariant dip latitude in degrees
Magnetic Local Time (MLT) in decimal hours
Corrected Geomagnetic (CGM) latitude&
Corrected Geomagnetic (CGM) longitude &Magnetic Local Time (MLT) for CGM coordinates&
CGM latitude& of equatorward boundary
CGM latitude& at MLT=0
CGM latitude& at MLT=1….
CGM latitude& at MLT=0
CGM latitude& at MLT=23
Kp at the time specified by the usermagnetic declination in degrees
L-value
Earth's dipole moment
sunrise at 300km in decimal hours
sunset at 300km in decimal hours
lower boundary in km of IRI profile
upper boundary in km of IRI profile

Parameter is used for input (user input) as well as output (IRI computed).

&Please check subroutine GEOCGM01 in file IGRF.FOR for more information on the Corrected
Geomagnetic (CGM) coordinates. CGM coordinates are only calculated if you select JF(33)=.true.
(auroral boundaries are calculated).