OV-GSO-DC / dustem-wrapper_idl

Commit e3e586f0f8f10df9908a7abadc01bf1fed383704

Authored by Annie Hughes
1 parent a22cb625
Exists in master

handling nans

Showing 2 changed files with 15 additions and 0 deletions   Show diff stats
Inline Side-by-side
src/idl/dustem_fill_ext_dependent_columns.pro
  Diff comments   View file @e3e586f
... ... @@ -66,6 +66,13 @@ output.ext_p=la_mul(output.ext_smallp,output.ext_I)
66 66 polar_variance_iqu2ippsi,ext.ext_I,ext.ext_Q,ext.ext_U,ext.SIGEXTII,ext.SIGEXTQQ,ext.SIGEXTUU,ext.SIGEXTIQ,ext.SIGEXTIU,ext.SIGEXTQU, $
67 67 variance_smallp,variance_psi,variance_largep,/lac
68 68  
  69 +xx=where(finite(variance_largep) eq 0, nanct)
  70 +if nanct gt 0 then variance_largep[xx]=la_undef()
  71 +xx=where(finite(variance_smallp) eq 0, nanct)
  72 +if nanct gt 0 then variance_smallp[xx]=la_undef()
  73 +xx=where(finite(variance_psi) eq 0, nanct)
  74 +if nanct gt 0 then variance_psi[xx]=la_undef()
  75 +
69 76 output.SIGEXTp=variance_largep
70 77 output.SIGEXTsmallp=variance_smallp
71 78 output.SIGEXT_psi=variance_psi
... ...
src/idl/dustem_fill_sed_dependent_columns.pro
  Diff comments   View file @e3e586f
... ... @@ -66,10 +66,18 @@ output.largep=la_mul(output.smallp,output.stokesI)
66 66 polar_variance_iqu2ippsi,sed.stokesI,sed.stokesQ,sed.stokesU,sed.sigmaII,sed.sigmaQQ,sed.sigmaUU,sed.sigmaIQ,sed.sigmaIU,sed.sigmaQU, $
67 67 variance_smallp,variance_psi,variance_largep,/lac
68 68  
  69 +xx=where(finite(variance_largep) eq 0, nanct)
  70 +if nanct gt 0 then variance_largep[xx]=la_undef()
  71 +xx=where(finite(variance_smallp) eq 0, nanct)
  72 +if nanct gt 0 then variance_smallp[xx]=la_undef()
  73 +xx=where(finite(variance_psi) eq 0, nanct)
  74 +if nanct gt 0 then variance_psi[xx]=la_undef()
  75 +
69 76 output.sigma_largep=variance_largep
70 77 output.sigma_smallp=variance_smallp
71 78 output.sigma_psi=variance_psi
72 79  
  80 +
73 81 the_end:
74 82  
75 83 RETURN,output
... ...