set_plot,'ps'
      device,/inches,yoffset=0.5,ysize=10.
      !TYPE=28
      !FANCY=1
      !GRID=0

      !P.THICK=2.
      !X.THICK=3.
      !Y.THICK=3

      alt=[100.,105.,110.,115.,120.,125.,130.,135.,140.,$
             145.,150.,160.,170.,180.,190.,200.,220.,240.,260.,$
             280.,300.,320.,340.,360.,380.,400.]

      ut=findgen(25)

      nalt=26
      ntime=25
      nlat=2


; msis

      msis_tn_areq=fltarr(ntime,nalt)
      msis_n2_areq=fltarr(ntime,nalt)
      msis_o2_areq=fltarr(ntime,nalt)
      msis_ox_areq=fltarr(ntime,nalt)
      msis_tn_jica=fltarr(ntime,nalt)
      msis_n2_jica=fltarr(ntime,nalt)
      msis_o2_jica=fltarr(ntime,nalt)
      msis_ox_jica=fltarr(ntime,nalt)

; hwm 

      hwm_un_areq=fltarr(ntime,nalt)
      hwm_vn_areq=fltarr(ntime,nalt)
      hwm_un_jica=fltarr(ntime,nalt)
      hwm_vn_jica=fltarr(ntime,nalt)


; tiegcm

tn_areq=fltarr(ntime,nalt)
un_areq=fltarr(ntime,nalt)
vn_areq=fltarr(ntime,nalt)
ui_areq=fltarr(ntime,nalt)
vi_areq=fltarr(ntime,nalt)
wi_areq=fltarr(ntime,nalt)
n2_areq=fltarr(ntime,nalt)
o2_areq=fltarr(ntime,nalt)
ox_areq=fltarr(ntime,nalt)
xne_areq=fltarr(ntime,nalt)

tn_jica=fltarr(ntime,nalt)
un_jica=fltarr(ntime,nalt)
vn_jica=fltarr(ntime,nalt)
ui_jica=fltarr(ntime,nalt)
vi_jica=fltarr(ntime,nalt)
wi_jica=fltarr(ntime,nalt)
n2_jica=fltarr(ntime,nalt)
o2_jica=fltarr(ntime,nalt)
ox_jica=fltarr(ntime,nalt)
xne_jica=fltarr(ntime,nalt)

close,1
file=' '

; temp and densities

openr,1,'tn_dens'                  ; this is the file containing the
                                   ; ascii model output for the 
                                   ; densities and temperatures

a=' '

dum=fltarr(ntime,nalt)

   readf,1,alt
   readf,1,ut

   work=fltarr(ntime,nalt)

; jicamarca read first

   readf,1,work
   tn_jica=work
   readf,1,work
   ox_jica=work
   readf,1,work
   o2_jica=work
   readf,1,work
   n2_jica=work
   readf,1,work
   xne_jica=work
   readf,1,work
   msis_tn_jica=work
   readf,1,work
   msis_ox_jica=work
   readf,1,work
   msis_o2_jica=work
   readf,1,work
   msis_n2_jica=work

; then arequipa

   readf,1,work
   tn_areq=work
   readf,1,work
   ox_areq=work
   readf,1,work
   o2_areq=work
   readf,1,work
   n2_areq=work
   readf,1,work
   xne_areq=work
   readf,1,work
   msis_tn_areq=work
   readf,1,work
   msis_ox_areq=work
   readf,1,work
   msis_o2_areq=work
   readf,1,work
   msis_n2_areq=work


; winds

close,1
openr,1,'winds'                  ; this is the file containing
                                 ; the ascii model output for the winds
a=' '

dum=fltarr(ntime,nalt)

   readf,1,alt
   readf,1,ut

; jicamarca read first

   readf,1,work
   un_jica=work
   readf,1,work
   vn_jica=work
   readf,1,work
   ui_jica=work
   readf,1,work
   wi_jica=work
   readf,1,work
   hwm_un_jica=work
   readf,1,work
   hwm_vn_jica=work

; then arequipa 

   readf,1,work
   un_areq=work
   readf,1,work
   vn_areq=work
   readf,1,work
   ui_areq=work
   readf,1,work
   wi_areq=work
   readf,1,work
   hwm_un_areq=work
   readf,1,work
   hwm_vn_areq=work



; plots

      !p.multi=[0,2,4,0,1]
;     window,1
      !x.range=[0,24]
      !xticks=6


      work=fltarr(ntime)
      time=fltarr(ntime)
      tickn=['06','10','14','18','22','02','06']
      !ytitle=' '
      !xtitle=' '
      !y.range=[650,900]
      !yticks=5
      !mtitle='Tn - Arequipa - 300 km'
      !linetype=0
      for it=11,23 do work(it-11)=msis_tn_areq(it,20)
      for it=0,11 do work(it+13)=msis_tn_areq(it,20)
      plot,work,font=0,charsize=1.3,xticklen=0.06,xtickv=time,$
         xticks=6,xtickname=tickn,yticklen=0.04
      !linetype=2
      for it=11,23 do work(it-11)=tn_areq(it,20)
      for it=0,11 do work(it+13)=tn_areq(it,20)
      oplot,work
      !linetype=0
      oplot,[15,17],[852,852]
      !linetype=2
      oplot,[15,17],[818,818]
      xyouts,17.5,850,'MSIS',font=0,charsize=1.1
      xyouts,17.5,815,'TIEGCM',font=0,charsize=1.1

      !y.range=[1.e8,1.e9]
      !yticks=1
      !mtitle='N2 density (cm^3) - Arequipa - 220 km'
      !linetype=0
      for it=11,23 do work(it-11)=msis_n2_areq(it,17)
      for it=0,11 do work(it+13)=msis_n2_areq(it,17)
      plot,work,font=0,charsize=1.3,xticklen=0.06,xtickv=time,$
         xticks=6,xtickname=tickn,/ylog,yticklen=0.04
      !linetype=2
      for it=11,23 do work(it-11)=n2_areq(it,17)
      for it=0,11 do work(it+13)=n2_areq(it,17)
      oplot,work


      !y.range=[-160,160]
      !yticks=4
      !ytitle=' '
      !xtitle=' '
      !mtitle='zonal wind (+E) - Arequipa - 300 km'
      !xtitle=' '
      !linetype=0
      for it=11,23 do work(it-11)=hwm_un_areq(it,20)
      for it=0,11 do work(it+13)=hwm_un_areq(it,20)
      plot,work,font=0,charsize=1.3,xticklen=0.06,xtickv=time,$
         xticks=6,xtickname=tickn,yminor=2,yticklen=0.04
      for it=11,23 do work(it-11)=un_areq(it,20)
      for it=0,11 do work(it+13)=un_areq(it,20)
      !linetype=2
      oplot,work
      !linetype=1
      oplot,[0,24],[0,0]
      !linetype=0
      oplot,[1,3],[110,110]
      !linetype=2
      oplot,[1,3],[70,70]
      xyouts,3.5,108,'HWM',font=0,charsize=1
      xyouts,3.5,68,'TIEGCM',font=0,charsize=1

end