implement color correction in d12

Close #128

pysm3/data/presets.cfg

@@ -180,6 +180,7 @@ freq_ref = "353 GHz"
 map_mbb_index = "mkd_dust/2048/thermaldust_specind{layer}.fits"
 map_mbb_temperature = "mkd_dust/2048/thermaldust_temp{layer}.fits"
 num_layers = 6
+color_correction = 0.911
 unit_layers = "MJy/sr"
 unit_mbb_temperature = "K"
 max_nside = 2048

pysm3/models/dust_layers.py

@@ -18,6 +18,7 @@ def __init__(
         map_mbb_index,
         map_mbb_temperature,
         nside,
+        color_correction=1,
         max_nside=None,
         num_layers=1,
         unit_layers=None,
@@ -51,6 +52,9 @@ def __init__(
             Alternatively an array of shape (num_layers, num_pix)
         nside: int
             Resolution parameter at which this model is to be calculated (with `ud_grade`)
+        color_correction: float
+            Scalar correction factor multiplied to the maps, implemented to add
+            a color correction factor to Planck HFI 353 GHz maps
         """
         super().__init__(nside=nside, max_nside=max_nside, map_dist=map_dist)
         num_pix = hp.nside2npix(nside)
@@ -66,6 +70,9 @@ def __init__(
         else:
             self.layers = u.Quantity(map_layers, unit_layers)
 
+        self.color_correction = color_correction
+        self.layers *= self.color_correction
+
         self.freq_ref = u.Quantity(freq_ref).to(u.GHz)
 
         with u.set_enabled_equivalencies(u.cmb_equivalencies(self.freq_ref)):