COSMIC-PAH / photoelectric-heating

Commit 0917f1d1c8ce735f6c509324057651d59cdbbbd8

Authored by Jalabert
1 parent 8f494b35
Exists in master

feat add commentaries

Showing 2 changed files with 7 additions and 2 deletions   Show diff stats
Inline Side-by-side
@@ -77,4 +77,8 @@ Here you can find : @@ -77,4 +77,8 @@ Here you can find :
77 2. Then you have the files containing the data of the stars or the interstellar medium for different approximations. Each file has a column with wavelengths (in nm) associated with an intensity (in erg cm-2 s-1 nm-1 sr-1). 77 2. Then you have the files containing the data of the stars or the interstellar medium for different approximations. Each file has a column with wavelengths (in nm) associated with an intensity (in erg cm-2 s-1 nm-1 sr-1).
78 Here you can find : 78 Here you can find :
79 1. ISRF file 79 1. ISRF file
80 - 2. stars file  
81 \ No newline at end of file 80 \ No newline at end of file
  81 + 2. stars file
  82 +
  83 +## Authors
  84 +
  85 +The authors of the programs are O. Berné, S. Foschino and F. Jalabert
82 \ No newline at end of file 86 \ No newline at end of file
radiation_fields.py
  Diff comments   View file @0917f1d
@@ -84,7 +84,7 @@ def radiation_field(filename, star_radius , parsec, ISRF=False, RF_list=False): @@ -84,7 +84,7 @@ def radiation_field(filename, star_radius , parsec, ISRF=False, RF_list=False):
84 wavelength_intensity = wave_intensity['col2'] #in erg cm-2 s-1 nm-1 sr-1 84 wavelength_intensity = wave_intensity['col2'] #in erg cm-2 s-1 nm-1 sr-1
85 #=> intensity per nm per sr 85 #=> intensity per nm per sr
86 86
87 - ''' Usefull parameters ''' 87 + ''' useful parameters '''
88 radiation_field = [] 88 radiation_field = []
89 distance = [] 89 distance = []
90 energy = ( ( (h/ev) * (c*1e2) ) / (wavelength*1e-7) )[::-1] #in ev 90 energy = ( ( (h/ev) * (c*1e2) ) / (wavelength*1e-7) )[::-1] #in ev
@@ -157,6 +157,7 @@ def radiation_field(filename, star_radius , parsec, ISRF=False, RF_list=False): @@ -157,6 +157,7 @@ def radiation_field(filename, star_radius , parsec, ISRF=False, RF_list=False):
157 distance.append(d) 157 distance.append(d)
158 radiation_field.append(g_0) 158 radiation_field.append(g_0)
159 159
  160 + ''' creation of data files '''
160 document_1 = Table([wavelength , wavelength_intensity] , names =('col1','col2')) 161 document_1 = Table([wavelength , wavelength_intensity] , names =('col1','col2'))
161 document_1.meta['comments'] = ['wavelength (nm)' , 'intensity (erg s-1 cm-2 nm-1 sr-1)'] 162 document_1.meta['comments'] = ['wavelength (nm)' , 'intensity (erg s-1 cm-2 nm-1 sr-1)']
162 document_1.write('intensity_per_wavelength.txt', format = 'ascii', overwrite=True) 163 document_1.write('intensity_per_wavelength.txt', format = 'ascii', overwrite=True)