This file of magnetic field model input parameters (link at bottom of page) has been derived from OMNIWeb data using the method described in Qin et al. [Solar Wind Parameters for Magnetospheric Magnetic Field Modelling, Space Weather Journal, 2007]. Please note that this database is not intended to replace the OMNIWeb database and recent changes to the OMNIWeb database may or may not be encorporated into this file. Beginning in January, 2008, this data will be available through the Virtual Radiation Belt Observatory (VIRBO, http://virbo.org). Our hope is that more recent data (after May 2007) will also be regularly made available there.
Appropriate acknowledgement would depend on the use. For minor use, you can add a sentence to the acknowledgements section such as "Magnetic field input parameters (derived from data on OMNIWeb) were provided by Richard Denton and Zhengui Qin." Please also acknowledge OMNIWeb as indicated on its website (see http://omniweb.gsfc.nasa.gov/html/ow_data.html#contact). For a more substantial use (when the input parameters are crucial for the study) please also acknowledge a grant number that Richard Denton will give you when you contact him (see below). Then possible acknowledgement could look like "Magnetic field input parameters (derived from data on OMNIWeb) were provided by Richard Denton and Zhengui Qin, supported by so-and-so grant." If this data set is a major part of the paper (like a statistical study of the magnetic field for a large portion of the time period covered), we would appreciate an offer for us to collaborate as coauthors. Please feel free to make the decision yourself as to which form of acknowledgement is appropriate, but if possible let Richard Denton know when you publish a paper acknowledging the use of this database. You can contact Richard Denton at
As described in Qin et al. , the status variables indicate the quality of the data. The "8 status" field gives 8 status digits corresponding to the 8 preceding quantities, ByIMF through G3. The "12 status" field gives 12 status digits corresponding to the 12 preceding quantities, Bz1 through W6. If the status digit is 2, the corresponding quantity has been measured, or calculated from measured quantities (G parameters), or the average value of the status of the quantities used to get it was 2 (W parameters). If the status digit is 1, the corresponding quantity was interpolated within about 2 correlation times (so that the neighboring measured quantities had an effect on the values), or interpolated values were used to get the value (G parameters), or the average value of the status of the quanties used to get it was 1 (W parameters). Similarly, if the status digit is 0, the quantity came from the 20 day average, or the values used to get the value came from a 20 day average (G parameters), or the average value of the status variables used to get it was 0 (W parameters).
Note carefully: If the status variable is 2, the quantity you are using is fairly well determined. If it is 1, the value has some connection to measured values, but is not directly measured. These values are still better than just using an average value, but not as good as those with the status variable equal to 2. If the status variable is 0, the quantity is based on average quantities, and the values listed are no better than an average value. The lower the status variable, the less confident you should be in the value.
Items in the file:
Year = 4 digit year
Day = day of year
Hr = hour
ByIMF = IMF By (GSM) in nT
BzIMF = IMF Bz (GSM) in nT
V_SW = solar wind speed in km/s
Den_P = solar wind proton density in particles / cm^3
Pdyn = solar wind dynamic pressure in nPa
G1,G2,G3 = Tsyganenko's G parameters for the T01 and T01s models
8 status = field of 8 status variables corresponding to the 8 preceding quantities
kp = Kp
akp3 = average value of Kp at time t using Kp values at preceding time t_p averaged wit the weighting function exp( - (t - t_p)/(3 days) )
dst = Dst - for recent dates, may be provisional or quicklook values from Kyoto
Bz1-Bz6 = Bzi values averaged for use in calculating Wi as described by Qin et al.
W1-W6 = W parameters used in the T04s model
12 status = status variables for the preceding 12 quantities
In Fortran, a line of data can be read with the following code:
dimension istatus8(8),istatus12(12),Bzi(6),W(6) read(iunitnumber,100) iYear,iDay,iHr,ByIMF,BzIMF,V_SW,Den_P,Pdyn,G1,G2,G3 & ,(istatus8(i),i=1,8),akp,akp3,dst,(Bzi(i),i=1,6),(W(i),i=1,6) & ,(istatus12(i),i=1,12) 100 format(i5,i4,i3,2f7.2,f7.1,2f7.2,3f7.2,i3,7i1,f7.2,f7.2,f7.1 & ,6f7.2,6f9.3,i3,11i1)
The solar wind parameters file may be freely distributed as long as the header is kept. Here's the link to the file. (On Windows, it's best to use "Save Target As".)
Database for input parameters for all Tsyganenko magnetic field models from 1963 to May 31, 2007 (80 MB Ascii with header)
The same database gzipped (14.4 MB)
The development of this database was supported by NSF grants ATM-0632740 and ATM-0120950 (Center for Integrated Space Weather Modeling, CISM,
funded by the NSF Science and Technology Centers Program) and NASA grants NNG05GJ02G and NNX06AD81G.