; read in Vikram's data file:
; The columns are:
;   Radius (cms), Density (g/cc), Pressure, Velocity (cms/sec), Temperature
;   0.86983E+18  0.74111E-23  0.61268E-11  0.38385E+08  0.10000E+05
;   0.87374E+18  0.74111E-23  0.61268E-11  0.38558E+08  0.10000E+05

readcol, "/nfs/spectra/d1/dd/Science/Vikram/wrbub.1067", vik_rad, vik_dens, $
  vik_pres, vik_vel, vik_temp

amu = 1.660e-24
parsec = 3.1e18  ; cm

!p.multi = [0,1,3]



; Density
; amu/cm^3
plot_io, vik_rad/parsec, vik_dens/amu, PSYM=3 

;Pressure 
;;plot_io, vik_rad/parsec, vik_pres, PSYM=3 
; this is ~ the same as pres:
; p ~ number_dens *  kT 
;;oplot, vik_rad/parsec, vik_temp * vik_dens * 0.8e8, PSYM=3

; Velocity
; in km/s
plot_io, vik_rad/parsec, vik_vel*1.e-5, PSYM=3 

; Temperature
plot_io, vik_rad/parsec, vik_temp, PSYM=3 

!p.multi=0