idl_scripts/petroglyph.pro at master · cnshingledecker/idl_scripts · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
PRO petroglyph
file = 'abundance.csv'


limit = ','

big_data = READ_CSV(file, N_TABLE_HEADER=2, DELIMITER=limit)

x = 1.0e-5 ; cm
y = 1.0e-5 ; cm
z = 1.0e-5 ; cm

volume = x*y*z

time = big_data.FIELD1
fluence = big_data.FIELD2
sp1 = big_data.FIELD3
sp2 = big_data.FIELD4

n1 = (1E20*sp1)/(volume) ;number density of O atoms
n2 = (1E20*sp2)/(volume) ;number density of O3 atoms
 ; fit = curvefit(fluence,n2)
;print, fit


PRINT, time[0],fluence[0],sp1[0],sp2[0]

aspect_ratio=1.5
xsize=9
ysize=xsize/aspect_ratio
set_plot, 'x'
;device,filename='O3_new.eps',encapsulated=1
;device, xsize=xsize,ysize=ysize

PLOT, fluence,sp2,/XLOG, XTITLE='Fluence (H$^{+}$/cm$^{2}$)',YTITLE='[O$_{3}$](10$^{20}$ cm$^{-3}$)'
;plot2 = PLOT(SMOOTH(n2,7),/OVERPLOT)
;plot2 = PLOT(time,sp1, /XLOG, /OVERPLOT)

;leg = LEGEND(TARGET=[plot1] , /AUTO_TEXT_COLOR,label='O$_{3}$', /RELATIVE, POSITION=[0.3,0.9],FONT_NAME="Hershey 3",FONT_SIZE=30)
;device, /close
set_plot, 'x'

;PRINT, sp1[0]


END