diff --git a/src/idl/dustem_fill_ext_dependent_columns.pro b/src/idl/dustem_fill_ext_dependent_columns.pro
index ad87f48..9101b75 100644
--- a/src/idl/dustem_fill_ext_dependent_columns.pro
+++ b/src/idl/dustem_fill_ext_dependent_columns.pro
@@ -66,6 +66,13 @@ output.ext_p=la_mul(output.ext_smallp,output.ext_I)
 polar_variance_iqu2ippsi,ext.ext_I,ext.ext_Q,ext.ext_U,ext.SIGEXTII,ext.SIGEXTQQ,ext.SIGEXTUU,ext.SIGEXTIQ,ext.SIGEXTIU,ext.SIGEXTQU, $
 						variance_smallp,variance_psi,variance_largep,/lac
 
+xx=where(finite(variance_largep) eq 0, nanct)
+if nanct gt 0 then variance_largep[xx]=la_undef()
+xx=where(finite(variance_smallp) eq 0, nanct)
+if nanct gt 0 then variance_smallp[xx]=la_undef()
+xx=where(finite(variance_psi) eq 0, nanct)
+if nanct gt 0 then variance_psi[xx]=la_undef()
+
 output.SIGEXTp=variance_largep
 output.SIGEXTsmallp=variance_smallp
 output.SIGEXT_psi=variance_psi
diff --git a/src/idl/dustem_fill_sed_dependent_columns.pro b/src/idl/dustem_fill_sed_dependent_columns.pro
index 6f92149..9ad34ed 100644
--- a/src/idl/dustem_fill_sed_dependent_columns.pro
+++ b/src/idl/dustem_fill_sed_dependent_columns.pro
@@ -66,10 +66,18 @@ output.largep=la_mul(output.smallp,output.stokesI)
 polar_variance_iqu2ippsi,sed.stokesI,sed.stokesQ,sed.stokesU,sed.sigmaII,sed.sigmaQQ,sed.sigmaUU,sed.sigmaIQ,sed.sigmaIU,sed.sigmaQU, $
 						variance_smallp,variance_psi,variance_largep,/lac
 
+xx=where(finite(variance_largep) eq 0, nanct)
+if nanct gt 0 then variance_largep[xx]=la_undef()
+xx=where(finite(variance_smallp) eq 0, nanct)
+if nanct gt 0 then variance_smallp[xx]=la_undef()
+xx=where(finite(variance_psi) eq 0, nanct)
+if nanct gt 0 then variance_psi[xx]=la_undef()
+
 output.sigma_largep=variance_largep
 output.sigma_smallp=variance_smallp
 output.sigma_psi=variance_psi
 
+
 the_end:
 
 RETURN,output
--
libgit2 0.21.2