OV-GSO-DC / dustem-wrapper_idl

Commit e1f224acc0d9151abe5d05b3d546b2e9eae7ed93

Authored by Ilyes Choubani
1 parent 9113e2b5
Exists in master

Not thoroughly tested: Added one line in write_isrf_plugin so that the ISRF form… 

…alism we discussed is implemented
Showing 1 changed file with 6 additions and 3 deletions   Show diff stats
Inline Side-by-side
src/idl/dustem_write_isrf_lv.pro
  Diff comments   View file @e1f224a
@@ -85,7 +85,7 @@ IF ctpop NE 0 or ctusrisrf NE 0 THEN BEGIN @@ -85,7 +85,7 @@ IF ctpop NE 0 or ctusrisrf NE 0 THEN BEGIN
85 ind_G0 = where(strupcase(strmid((*(*!dustem_fit).fixed_param_descs),5,/reverse)) EQ 'MS).G0', countg0) 85 ind_G0 = where(strupcase(strmid((*(*!dustem_fit).fixed_param_descs),5,/reverse)) EQ 'MS).G0', countg0)
86 IF countg0 NE 0 THEN G0_mathis = (*(*!dustem_fit).fixed_param_init_values)[ind_G0] 86 IF countg0 NE 0 THEN G0_mathis = (*(*!dustem_fit).fixed_param_init_values)[ind_G0]
87 ENDIF 87 ENDIF
88 - 88 + STOP
89 ;Locating gas.G0 in the free parameters vector 89 ;Locating gas.G0 in the free parameters vector
90 ind_GASG0 = where(strupcase(strmid((*(*!dustem_fit).param_descs),5,/reverse)) EQ 'GAS.G0', countgasg0) 90 ind_GASG0 = where(strupcase(strmid((*(*!dustem_fit).param_descs),5,/reverse)) EQ 'GAS.G0', countgasg0)
91 IF countgasg0 THEN BEGIN 91 IF countgasg0 THEN BEGIN
@@ -119,7 +119,10 @@ IF ctpop NE 0 or ctusrisrf NE 0 THEN BEGIN @@ -119,7 +119,10 @@ IF ctpop NE 0 or ctusrisrf NE 0 THEN BEGIN
119 c(1) = '# STELLAR POPULATION ISRF' 119 c(1) = '# STELLAR POPULATION ISRF'
120 c(2) = '# USER-DEFINED ISRF' 120 c(2) = '# USER-DEFINED ISRF'
121 ENDIF 121 ENDIF
122 - 122 +
  123 + IF G0_mathis NE 0 THEN G0_mathis=1
  124 +
  125 + ;the below condition will always be valid
123 IF G0_mathis NE 0 THEN BEGIN ;storing the mathis isrf in variable mathis_isrf if the user wants to use it 126 IF G0_mathis NE 0 THEN BEGIN ;storing the mathis isrf in variable mathis_isrf if the user wants to use it
124 ;This is reading the mathis isrf from the fortran release 127 ;This is reading the mathis isrf from the fortran release
125 ;Let's hope they keep releasing it otherwise we'll have to include it ourselves 128 ;Let's hope they keep releasing it otherwise we'll have to include it ourselves
@@ -141,7 +144,7 @@ IF ctpop NE 0 or ctusrisrf NE 0 THEN BEGIN @@ -141,7 +144,7 @@ IF ctpop NE 0 or ctusrisrf NE 0 THEN BEGIN
141 ENDIF ELSE c(Ncomments-3) = '# MATHIS ISRF' 144 ENDIF ELSE c(Ncomments-3) = '# MATHIS ISRF'
142 145
143 ENDIF 146 ENDIF
144 - ;Will this enable us to have the current 'updated' ISRF through the !dustem_params system variable? 147 +
145 st.isrf = final_isrf 148 st.isrf = final_isrf
146 149
147 ;Last .DAT header comments lines 150 ;Last .DAT header comments lines