workaround now that fortran uses negative numbers for gas G0

Annie Hughes
1 parent 52f17b34
Showing 1 changed file with 9 additions and 6 deletions Show diff stats
src/idl/dustem_write_isrf_release.pro
@@ -78,8 +78,11 @@ IF ctpop NE 0 or ctusrisrf NE 0 THEN BEGIN
  
     G0_mathis = 1. ;Default value
  
-    IF ((*!dustem_params).g0) NE 1. then G0_mathis = ((*!dustem_params).g0) 
-    IF ((*!dustem_params).gas.g0) NE 1. then G0_mathis = float((*!dustem_params).gas.g0)
+    IF ((*!dustem_params).g0) NE 1. then G0_mathis = ((*!dustem_params).g0)
+; Fortran now uses negative numbers to specify Gas.dat G0
+; If positive G0 present in Gas.dat, it takes precedence over grain.dat G0   
+;   IF ((*!dustem_params).gas.g0) NE 1. then G0_mathis = float((*!dustem_params).gas.g0)
+    IF ((*!dustem_params).gas.g0) gt 0. then G0_mathis = float((*!dustem_params).gas.g0)
  
     ;Getting the ISRF default wavelengths
     ;using this structure
@@ -107,9 +110,9 @@ IF ctpop NE 0 or ctusrisrf NE 0 THEN BEGIN
     IF G0_mathis GT 1E-10 THEN BEGIN  ;discussed with JP
         ;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
-        ma_isrf_dir=!dustem_soft_dir+'data/ISRF_MATHIS.DAT'
-        ma_isrf=dustem_read_isrf(ma_isrf_dir)
-        final_isrf=final_isrf/G0_mathis+ma_isrf.isrf ;[0]
+        mathis_isrf_datfile=!dustem_soft_dir+'data/ISRF_MATHIS.DAT'
+        mathis_isrf=dustem_read_isrf(mathis_isrf_datfile)
+        final_isrf=final_isrf/G0_mathis+mathis_isrf.isrf ;[0]
         Ncomments=4
         c=strarr(Ncomments)
         c(0)='# DUSTEM_WRAP: exciting radiation field featuring'
@@ -147,7 +150,7 @@ IF ctpop NE 0 or ctusrisrf NE 0 THEN BEGIN
  
  
 ;WRITING THE NEW ISRF
-    
+
     file=!dustem_dat+'data/ISRF.DAT'  
     openw,unit,file,/get_lun