These are the files i used to make the plots in the paper.  
the actual HIAMCM data is higher resolution (1/2 degree) and is finer in z.
You can request that data from Dr. Erich Becker if you'd like to have it.

The hiamcm data params file is ascii.
fileparams='params.txt'
you can read it in idl with this code fragment:
openr,1,fileparams
readf,1, nlon,nlat,nz,nt
nlon=fix(nlon+0.01) & nlat=fix(nlat+0.01) & nz=fix(nz+0.01) & nt=fix(nt+0.01)
lon=fltarr(nlon) & lat= fltarr(nlat) & zzkm=fltarr(nz) & timehrs=fltarr(nt)
readf,1,lon
readf,1,lat
readf,1,zzkm
readf,1,timehrs
close,1

Here, nlon-# longitude pts, nlat=# lat pts, 
nz=# z levels, nt=# time pts in file.
lon=arry with longitudes
lat=arry with latitudes
zzkm=arry with heights (in km)
timehrs=arry with times on 15 jan 2022 (in hrs)


the HIAMCM data files are binary files
filename='JAN12_lonmap.dat'
you can read them in idl with this code fragment:
    openr,1,filename

  ; perturbations quantities (waves with horizontal wavelengths>2000 km)
      duvel= array(0)  ; zonal velocity (m/s) 
      dvvel= array(1)  ; meridional velocity  (m/s)
      dwvel= array(2)  ; vertical velocity  (m/s)
      dtemp= array(3)  ; temperature (in T)
      drho= array(4)  ; density (in g/m^3)
  ; background quantities (waves with horizontal wavelengths>2000 km)
      uvelb= array(5)  ; zonal veocity (m/s)
      vvelb= array(6)  ; meridional veocity (m/s)
      wvelb= array(7)  ; vertical veocity (m/s)
      tempb= array(8)  ; temperature perturbation (in K)
      rhob= array(9)  ; density perturbation (in g/m^3)
    close,1