FUNCTION dustem_fill_sed_dependent_columns,sed,help=help

;+
; NAME:
;    dustem_fill_sed_dependent_columns

;

; PURPOSE:

;    fill up columns of a DustEMWrap SED that are depednent on I,Q,U and associated (co-)variances

;

; CATEGORY:

;    DustEMWrap, Distributed, High-Level, User-convenience

;

; CALLING SEQUENCE:
;    output=dustem_fill_sed_dependent_columns(sed[,/help])

;

; INPUTS:
;    sed    : a valid dustemwrap sed

;

; OPTIONAL INPUT PARAMETERS:
;    None

;

; OUTPUTS:
;    None

;

; OPTIONAL OUTPUT PARAMETERS:
;    None

;

; ACCEPTED KEY-WORDS:
;    help      = If set print this help

;

; COMMON BLOCKS:
;    None

;

; SIDE EFFECTS:
;    None

;

; RESTRICTIONS:

;    The DustEMWrap IDL code must be installed

;

; PROCEDURE:
;    None

;

; EXAMPLES
;    filled_sed=dustem_fill_sed_dependent_columns(sed)

;

; MODIFICATION HISTORY:

;    Written by JPB June-2022
;    Evolution details on the DustEMWrap gitlab.
;    See http://dustemwrap.irap.omp.eu/ for FAQ and help.  

;-

IF keyword_set(help) THEN BEGIN
  doc_library,'dustem_fill_sed_dependent_columns'
  output=0
  goto,the_end
ENDIF

output=sed

polar_iqu2ippsi,sed.stokesI,sed.stokesQ,sed.stokesU,p,psi,p2=p2,/lac
output.smallp=p
output.psi=psi
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

END