From e1f224acc0d9151abe5d05b3d546b2e9eae7ed93 Mon Sep 17 00:00:00 2001 From: ichoubani Date: Tue, 22 Nov 2022 14:48:10 +0100 Subject: [PATCH] Not thoroughly tested: Added one line in write_isrf_plugin so that the ISRF formalism we discussed is implemented --- src/idl/dustem_write_isrf_lv.pro | 9 ++++++--- 1 file changed, 6 insertions(+), 3 deletions(-) diff --git a/src/idl/dustem_write_isrf_lv.pro b/src/idl/dustem_write_isrf_lv.pro index d8474d9..c1aaeee 100755 --- a/src/idl/dustem_write_isrf_lv.pro +++ b/src/idl/dustem_write_isrf_lv.pro @@ -85,7 +85,7 @@ IF ctpop NE 0 or ctusrisrf NE 0 THEN BEGIN ind_G0 = where(strupcase(strmid((*(*!dustem_fit).fixed_param_descs),5,/reverse)) EQ 'MS).G0', countg0) IF countg0 NE 0 THEN G0_mathis = (*(*!dustem_fit).fixed_param_init_values)[ind_G0] ENDIF - + STOP ;Locating gas.G0 in the free parameters vector ind_GASG0 = where(strupcase(strmid((*(*!dustem_fit).param_descs),5,/reverse)) EQ 'GAS.G0', countgasg0) IF countgasg0 THEN BEGIN @@ -119,7 +119,10 @@ IF ctpop NE 0 or ctusrisrf NE 0 THEN BEGIN c(1) = '# STELLAR POPULATION ISRF' c(2) = '# USER-DEFINED ISRF' ENDIF - + + IF G0_mathis NE 0 THEN G0_mathis=1 + + ;the below condition will always be valid IF G0_mathis NE 0 THEN BEGIN ;storing the mathis isrf in variable mathis_isrf if the user wants to use it ;This is reading the mathis isrf from the fortran release ;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 ENDIF ELSE c(Ncomments-3) = '# MATHIS ISRF' ENDIF - ;Will this enable us to have the current 'updated' ISRF through the !dustem_params system variable? + st.isrf = final_isrf ;Last .DAT header comments lines -- libgit2 0.21.2