fixed bug due to change of geometry in the GAS.DAT file

Jean-Philippe Bernard
1 parent 19c45d4a
Showing 1 changed file with 4 additions and 0 deletions Show diff stats
src/idl/dustem_read_gas.pro
@@ -91,7 +91,11 @@ CASE !dustem_which OF
     st.nH=double(strv[ii]) & ii=ii+1
     st.nH2=double(strv[ii]) & ii=ii+1
     st.CR_rate=double(strv[ii]) & ii=ii+1
+    ;modified because the file geometry did change !!
+    ;st.G0=double(strv[ii]) & ii=ii+1
     st.G0=double(strv[ii]) & ii=ii+1
+    readf,unit,str
+    strv=str_sep(strcompress(strtrim(str,2)),' ') & ii=0L
     st.n_charges=long(strv[ii])
     charges_st=replicate(one_charge_st,st.n_charges)
     FOR i=0L,st.n_charges-1 DO BEGIN