diff --git a/CHANGELOG.md b/CHANGELOG.md
index e04e02f30e93b56fe1d5176db7f6db39f6792e25..61f3354d44114fbb47dec17c849c0ea13188bdd7 100644
--- a/CHANGELOG.md
+++ b/CHANGELOG.md
@@ -19,6 +19,7 @@
 - In the absence of a nebular component `restframe\_parameters` would crash when attempting to compute the equivalent widths of the lines listed in `EW_lines`. Now they are simply ignored. (Médéric Boquien)
 - The luminosity spectrum of the best fit was saved assuming the distance corresponding to the redshift rounded to two decimals. This was an issue in particular at very low redshift as a difference of 0.005 in redshift can translate to a large difference on the luminosity distance. Now the exact luminosity distance of the object is used to compute the spectrum luminosity. (Médéric Boquien, reported by Jorge Melnick)
 - When using the `parameters\_file` option, the indices of the models now correspond to the line number of the input file. (Médéric Boquien)
+- When using the `parameters\_file` option, the list of modules is read from `sed\_modules` rather than being inferred from the input file. (Médéric Boquien)
 
 ### Optimised
 - The cache architecture has been simplified, making it somewhat faster. It speeds up the model generation by ~1%. (Médéric Boquien)
diff --git a/pcigale/managers/parameters.py b/pcigale/managers/parameters.py
index 44df308c4695d2bf21443c2a4236d299fc7a7405..14a282ea4251e6816d2d0ee0825fed10d24f2109 100644
--- a/pcigale/managers/parameters.py
+++ b/pcigale/managers/parameters.py
@@ -178,12 +178,7 @@ class ParametersManagerFile(object):
         table = read_table(conf['parameters_file'])
 
         self.size = len(table)
-
-        self.modules = []
-        for colname in table.colnames:
-            module = colname.split('.', 1)[0]
-            if module not in self.modules:
-                self.modules.append(module)
+        self.modules = conf['sed_modules']
         self.blocks = self._split(range(self.size),
                                   conf['analysis_params']['blocks'])