The SAMI3 data for January 15, 2022 is in this directory.
The data is recorded every 5 minutes.
Note: the electron density, dene0, is included in the file, although we did not use
this quantity in this paper.


to read the data in IDL, use the following code fragment:


dir=''
dir='sami3_5min_data/'  ; the directory where you have the sami3 5-min data
fileinarry=dir + 'sami3_15JAN2022_5min_'+['tonga.nc','notonga.nc']

filein=fileinarry(0)  ; tonga run on 15 January 2022
filein=fileinarry(1)  ; no-tonga run on 15 January 2022 (MESORAC not included)

; first get lon, lat and dimensions of arrays.
;NCDF_CLOSE, Id
filename=filein
Id      = NCDF_OPEN(filename,/NOWRITE)
print,'         #################          '
print,'opening file=',filename

;dimensions:
; grid in native sami3 grid, use lon,lat,alt
;        nz = 304 ; (grid along B field)
;        nf = 124 ; ('number' of field lines)
;        nl = 720 ;  (grid in longitude)
;        nt = UNLIMITED ; // (289 currently)  (number of time steps)
; grid interpolated 
;        nx = 100 ; (latitude in interpolated grid)
;        ny = 100 ; (altitude in interpolated grid)
; tec/nmf2/hmf2 on interpolated grid, use lon0,lat0 with ny=0 
ncdf_varget, Id, 'time', timehr  ; time(nt) in hrs on 15 january 2022. 
ncdf_varget, Id, 'lon0', lon0  ; (nx,ny,nl) in deg. plot tec/nmf2/hmf2 with lat0/lon0 with ny=0 
ncdf_varget, Id, 'lat0', lat0  ; (nx,ny,nl) in deg.
ncdf_varget, Id, 'dene0', dene0  ; (nx,ny,nl,nt) electron density / grid interpolated
ncdf_varget, Id, 'tec', tec  ; (nx,nl,nt) in m^{-2}
ncdf_varget, Id, 'nmf2', NmF2  ; (nx,nl,nt) in cm^{-3}
ncdf_varget, Id, 'hmf2', hmF2  ; (nx,nl,nt) in km
NCDF_CLOSE, Id


nt= n_elements(timehr) & sz = size(lon0)

; i know that the data is every 5 min, so ill put that in instead.
;timehr=fltarr(nt)
;dt=5./60.
;for it=0,nt-1 do timehr(it)= float(it)*dt

print,'timehr',timehr(0),timehr(1),timehr(nt-1)
print,'max(tec)', max(tec)

nlon=nl & nlat=nx & k=0
print,nlon,nlat,nx,k
lonarry=fltarr(nlon,nlat) & latarry=fltarr(nlon,nlat)
tecarry=fltarr(nlon,nlat,nt) & NmF2arry=fltarr(nlon,nlat,nt) & hmF2arry=fltarr(nlon,nlat,nt)
for i=0,nlon-1 do begin
  for j=0,nlat-1 do begin
    lonarry(i,j)=lon0(j,k,i)
    latarry(i,j)=lat0(j,k,i)
    for it=0,nt-1 do begin
      tecarry(i,j,it)= tec(j,i,it)
      NmF2arry(i,j,it)= NmF2(j,i,it)
      hmF2arry(i,j,it)=hmF2(j,i,it)
    endfor
  endfor
endfor