The file format is as follows: The bulk of the file is made up of 2-D arrays of model output, dimensioned by (altitude,time). (There are 2 sets of these arrays, one for Arequipa and one for Jicamarca.) The time array has 25 times, running from 0 UT to 0 UT by hour The altitude array has 26 altitudes, from 100 to 400 km with finer spacing in the lower thermosphere. The actual values are found in the file. The 2-D arrays are Tn, neutral temperature N2 nitrogen number densities in cm-3 O2 molecular oxygen number densities in cm-3 O atomic oxygen number densities in cm-3 ne electron densities in cm-3 un neutral zonal wind, +E vn neutral meridional wind, +N ui zonal ion drift, +E wi vertical ion drift +up -- all winds/drifts are in m/s MSIS and HWM output for the same gridpoints and conditions as the TIEGCM are also in the file. In detail, the file format is: - first 3 lines: altitude array - next 3 lines: UT array - then blocks of 2-D arrays. - For the temperature and density file: - first come all temperature and density output for Jicamarca in this order: Tn O O2 N2 Ne MSIS Tn MSIS O MSIS O2 MSIS N2 - then the same set of output for Arequipa - For the winds file: - first all the winds and drifts at Jicamarca in this order un vn ui wi HWM un HWM vn - then the same set of output for Arequipa The two ascii files can be saved and then read by the idl code "read.model" The temperature and density file is called (in the idl code) "tn_dens"; the wind file is called "winds". If you don't use these names, you'll have to change the names of the files to be opened and read by the idl code. (Look for the lines "openr,1,....")