code: def _get_stddevs(self, C, stddev_types, num_sites):
          
          # standard deviation is converted from log10 to ln
          std_total = np.log(10 ** C['sigma'])
          stddevs = []
          for _ in stddev_types:
              stddevs.append(np.zeros(num_sites) + std_total)
  
          return stddevs
docs: Return total standard deviation.