changed default behaviour when CC method is unknown

Annie Hughes
1 parent 008ab3a3
Showing 1 changed file with 7 additions and 3 deletions Show diff stats
src/idl/dustem_cc.pro
@@ -214,9 +214,13 @@ FOR i=0L,Nbands-1 DO BEGIN
          cc[i]=num/(spec0*den*(st.central_wavelengths[ii])^beta)
       END
       ELSE:BEGIN
-        message,strupcase(fluxconvs[i])+' not recognized',/info
-        cc[i]=1.
-        ;stop
+        message,strupcase(fluxconvs[i])+' not recognized, assuming NUINU=CSTE',/info
+        ;cc[i]=1. ; changed default behaviour in version 4.3
+         beta=-1.
+         spec_int=interpol(spec2,wavein2,*(st.use_wavelengths(ii)))
+         num=integral(*(st.use_wavelengths(ii)),spec_int*(*(st.use_transmissions(ii)))/(*(st.use_wavelengths(ii)))^2.,st.use_wmin(ii),st.use_wmax(ii),/double)
+         den=integral(*(st.use_wavelengths(ii)),((*(st.use_transmissions(ii)))/((*(st.use_wavelengths(ii)))^(2+beta))),st.use_wmin(ii),st.use_wmax(ii),/double)               
+         cc[i]=num/(spec0*den*(st.central_wavelengths[ii])^beta)
       END
     ENDCASE
     !dustem_previous_cc=ptr_new(cc)