Commit e1600af70a547a56acbcf96ee9b2e7ff25ac1f7f

Authored by Annie Hughes
1 parent 84cabbe5
Exists in master

first commit

Showing 1 changed file with 76 additions and 0 deletions   Show diff stats
src/idl/dustem_fill_ext_dependent_columns.pro 0 → 100644
... ... @@ -0,0 +1,76 @@
  1 +FUNCTION dustem_fill_ext_dependent_columns,ext,help=help
  2 +
  3 +;+
  4 +; NAME:
  5 +; dustem_fill_ext_dependent_columns
  6 +;
  7 +; PURPOSE:
  8 +; fill up columns of a DustEMWrap EXT that are dependent on I,Q,U and associated (co-)variances
  9 +;
  10 +; CATEGORY:
  11 +; DustEMWrap, Distributed, High-Level, User-convenience
  12 +;
  13 +; CALLING SEQUENCE:
  14 +; output=dustem_fill_ext_dependent_columns(ext[,/help])
  15 +;
  16 +; INPUTS:
  17 +; ext : a valid DustEMWrap ext
  18 +;
  19 +; OPTIONAL INPUT PARAMETERS:
  20 +; None
  21 +;
  22 +; OUTPUTS:
  23 +; None
  24 +;
  25 +; OPTIONAL OUTPUT PARAMETERS:
  26 +; None
  27 +;
  28 +; ACCEPTED KEY-WORDS:
  29 +; help = If set print this help
  30 +;
  31 +; COMMON BLOCKS:
  32 +; None
  33 +;
  34 +; SIDE EFFECTS:
  35 +; None
  36 +;
  37 +; RESTRICTIONS:
  38 +; The DustEMWrap IDL code must be installed
  39 +;
  40 +; PROCEDURE:
  41 +; None
  42 +;
  43 +; EXAMPLES
  44 +; filled_ext=dustem_fill_ext_dependent_columns(ext)
  45 +;
  46 +; MODIFICATION HISTORY:
  47 +; Written by AH Nov-2022
  48 +; Evolution details on the DustEMWrap gitlab.
  49 +; See http://dustemwrap.irap.omp.eu/ for FAQ and help.
  50 +;-
  51 +
  52 +
  53 +IF keyword_set(help) THEN BEGIN
  54 + doc_library,'dustem_fill_ext_dependent_columns'
  55 + output=0
  56 + goto,the_end
  57 +ENDIF
  58 +
  59 +output=ext
  60 +
  61 +polar_iqu2ippsi,ext.ext_I,ext.ext_Q,ext.ext_U,p,psi,p2=p2,/lac
  62 +output.ext_smallp=p
  63 +output.psi=psi
  64 +output.ext_p=la_mul(output.ext_smallp,output.ext_I)
  65 +
  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 + variance_smallp,variance_psi,variance_largep,/lac
  68 +
  69 +output.SIGEXTp=variance_largep
  70 +output.SIGEXTsmallp=variance_smallp
  71 +output.SIGEXT_psi=variance_psi
  72 +
  73 +the_end:
  74 +
  75 +RETURN,output
  76 +END
... ...