Pro get_amie_data,data_1,data_2,data_3,iday1,iday2

;e.g., get_amie_data,'superd_dmsp_mags','binnedSuperD.binnedDMSP.mags.weimer','amie.wmsd',13,16
;e.g., get_amie_data,'superd_mags','binnedSuperD.mags.weimer','amie.wms',13,16
;e.g., get_amie_data,'mags','mags.weimer','amie.wm',13,16
;e.g., get_amie_data,'weimer','weimer','amie.w',13,16
;data_1='superd_dmsp_mags'
;data_2='binnedSuperD.binnedDMSP.mags.weimer'
;data_3='amie.wmsd'


;outpath = '/home/jsshim/AMIE_geoff/idl/20061213-15/'
;openw,11,outpath+'20061213-15.'+data_3
openw,11,'20061213-15.'+data_3+'.emery'
printf,11,'year  mon day hh  mm  ss     mlts   lats     epot(kV)  aurora_eflux(mW/m^2) Qjh  aurora_kev'
for iday=iday1,iday2 do begin
print,'iday = ',iday
ciday=string(iday,format='(i2.2)')
date='0612'+ciday

;afile='../'+data_1+'/12.13-15/amieout.'+date+'.0000.'+data_2+'.hires.244.afile_sond'
afile='amieout.'+date+'.0000.'+data_2+'.hires.244.afile_sond'
epot=extract_amie_data(afile, 'efpot')
au_eflux=extract_amie_data(afile, 'wpm2')
jh=extract_amie_data(afile, 'smqj')
au_kev=extract_amie_data(afile, 'emkev')

;afile='../'+data_1+'/12.13-15/amieout.'+date+'.0000.SH.'+data_2+'.hires.244.afile_sond'
afile='amieout.'+date+'.0000.SH.'+data_2+'.hires.244.afile_sond'
epot_S=extract_amie_data(afile, 'efpot')
au_eflux_S=extract_amie_data(afile, 'wpm2')
jh_S=extract_amie_data(afile, 'smqj')
au_kev_S=extract_amie_data(afile, 'emkev')

  for i=0,288-1 do begin
  time=epot.juls[i]
;   print,i,epot.juls[i]
  caldat,time,month,day,year,hour,minute,second
;  print,time,year,month,day,hour,minute,second,format="(f12.3,x,I4.4,x,5(I2.2,x))"

   for j=0,37-1 do begin
;  northern part
    for k=31-1,0,-1 do begin
     printf,11, year, month, day, hour, minute, second,epot.mlts[j],epot.mlat[k],$
     epot.amiedata[i,k,j]/1e3,au_eflux.amiedata[i,k,j]*1e3, jh.amiedata[i,k,j]*1e3,$
     au_kev.amiedata[i,k,j], $
     format="(i4.4,2x,5(i2.2,2x),2(f6.1,2x),2x,f9.4,2x,3(f10.5,2x))"
     endfor ; k-loop
; southern part
    for k=31-1,0,-1 do begin
     printf,11, year, month, day, hour, minute, second,epot_S.mlts[j],epot_S.mlat[k],$
     epot_S.amiedata[i,k,j]/1e3,au_eflux_S.amiedata[i,k,j]*1e3, jh_S.amiedata[i,k,j]*1e3, $
     au_kev_S.amiedata[i,k,j], $
     format="(i4.4,2x,5(i2.2,2x),2(f6.1,2x),2x,f9.4,2x,3(f10.5,2x))"
    endfor ; k-loop
    endfor ; j-loop
  endfor  ; i-loop
endfor  ; day
  
  close,11


end