-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathREADME
50 lines (29 loc) · 1.62 KB
/
README
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
49
Infrared absorption bands of interstellar dust
Full theory of work with links contains in report.
For changing input parameters you should edit constants.f95
Data about scattering index of ice I got from
Hudgins et. al., Mid and Far infrared spectrscopy of ices Astrophys.J.Suppl. 86:713-870, 1993 June
The head of "esp_ice.dat" is
lambda(cm^-1) --- n ---- k
Data about astronomical silicate I get from Draine (1985) Astrophys.J.Suppl. 57, 587.
The head of "eps_sil.dat" is
lamda(mkm) ------ esp1 ------ eps2 ------ n-1---- k eps = dielectic coefficients
After that I interpolated this data to same GRID with interpolation module writed on octave and put it on data2.txt
The head of "data2.txt"
lambda(mkm) n_ice, k _ice , n_sil, k_sil
You can change data about ice and silicate, but after that you should change interpolation rules.
It`s not difficult with octave
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!! After change num you should rebuild data2.txt with octave !!!!
!! you should put in command line this string "octave interpolation.m"!!!!
!! !!!!
!! then you should replace NA to nearest value in data2.txt !!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
First of all you should run 'make' in you command line
For run it you should type 'make run'
For make graphics from result data you should run 'make plot'
You'll get result for prolate particle, you can change just several functions in 'subfun.f95'
with agree to theory for get rusult for oblate case.
Get fun ^^
Aleksey Artamonov