;+
; NAME:
;     plasma_temp_plot
; PURPOSE:
;     plot height vs temperature
; CATEGORY:
; CALLING SEQUENCE: 
; INPUTS:    
; OPTIONAL (KEYWORD) INPUT PARAMETERS:
; OUTPUTS:    
; COMMON BLOCKS:
; SIDE EFFECTS:
; RESTRICTIONS:
;     
; MODIFICATION HISTORY:
;     (SC) Fri Mar 5 2010
;     (SC) Fri Apr 30 2010: Added the output table
;-
 

thick = 1 & xthick = 1 & ythick = 1 & charthick = 1
if !d.name eq 'PS' then !p.font = 0 else !p.font = -1
if !D.name eq 'PS' then begin
  thick=5 & xthick=5 & ythick=5 & charthick = 3
  device,/encapsulate,bits=8,file='plasma_temperature.eps'
  device,/color
endif

model1 = 'models/Model1001.tab.txt'
model2 = 'models/Model1005.tab.txt'
model3 = 'models/Model1006.tab.txt'

aa = read_font(model1)
bb = read_font(model2)
cc = read_font(model3)

plot_io, aa.height, aa.temp,$
  xtit = 'Height, cm', ytit='Temperature, K',$
  title = 'Plasma Temperature',$
  thick=thick, xthick=xthick, ythick=ythick, charthick=charthick,$
  yr=[1e3,1e6],yst=1,xr=[-2e7, 2e8],xst=1, psym=-2, linestyle=1
hsi_linecolors

oplot, aa.height, aa.temp, color=4, thick=thick, psym=-2, linestyle=1
oplot, bb.height, bb.temp, color=6, thick=thick, psym=-2, linestyle=1
oplot, cc.height, cc.temp, color=3, thick=thick, psym=-2, linestyle=1


get_lun, lun
openw, lun, 'plasma_temp.dat'
printf, lun, 'Densities in Fontenla Model 1001', format='(A-32)'
printf, lun, 'Height, cm', 'Temperature, K', 'Electron density',$ 
    'Proton density', 'Hydrogen density',$ 
    format='(A-17,x,A-17,x,A-17,x,A-17,x,A-17)'
format='(e-18.2, e-18.2, e-18.2, e-18.2, e-18.2)'
FOR k = 0,90 DO BEGIN
  heighta = aa.height[k]
  tempa = aa.temp[k]
  n_ea = aa.n_e[k]
  n_pa = aa.n_p[k]
  n_ha = aa.n_h[k]
  printf, lun, heighta, tempa, n_ea, n_pa, n_ha, format=format
ENDFOR

printf, lun, 'Densities in Fontenla Model 1005', format='(A-32)'
printf, lun, 'Height, cm', 'Temperature, K', 'Electron density',$
    'Proton density', 'Hydrogen density',$
    format='(A-17,x,A-17,x,A-17,x,A-17,x,A-17)'
format='(e-18.2, e-18.2, e-18.2, e-18.2, e-18.2)'
FOR k = 0,100 DO BEGIN
  heightb = bb.height[k]
  tempb = bb.temp[k]
  n_eb = bb.n_e[k]
  n_pb = bb.n_p[k]
  n_hb = bb.n_h[k]
  printf, lun, heightb, tempb, n_eb, n_pb, n_hb, format=format
ENDFOR

printf, lun, 'Densities in Fontenla Model 1006', format='(A-32)'
printf, lun, 'Height, cm', 'Temperature, K', 'Electron density',$
    'Proton density', 'Hydrogen density',$
    format='(A-17,x,A-17,x,A-17,x,A-17,x,A-17)'
format='(e-18.2, e-18.2, e-18.2, e-18.2, e-18.2)'
FOR k = 0,79 DO BEGIN
  heightc = cc.height[k]
  tempc = cc.temp[k]
  n_ec = cc.n_e[k]
  n_pc = cc.n_p[k]
  n_hc = cc.n_h[k]
  printf, lun, heightc, tempc, n_ec, n_pc, n_hc, format=format
ENDFOR

close,lun

;loadct,0

if !D.name eq 'PS' then device,/close 
end