code: def _get_phi_vector(self, phi_mean, phi_std, imt_list):
          
          p_bar = {}
          p_std = {}
  
          for imt in imt_list:
              p_bar[imt] = []
              p_std[imt] = []
              for mag in self.magnitude_limits:
                  phi_ss_mean = get_phi_ss(imt, mag, phi_mean)
                  phi_ss_std = get_phi_ss(imt, mag, phi_std)
                  if self.ergodic:
                      # Add on the phi_s2ss term according to Eqs. 5.15 and 5.16
                      # of Al Atik (2015)
                      phi_ss_mean = np.sqrt(
                          phi_ss_mean ** 2. +
                          PHI_S2SS_MODEL[self.phi_s2ss_model][imt]["mean"] ** 2.
                          )
                      phi_ss_std = np.sqrt(
                          phi_ss_std ** 2. +
                          PHI_S2SS_MODEL[self.phi_s2ss_model][imt]["var"] ** 2.
                          )
                  p_bar[imt].append(phi_ss_mean)
                  p_std[imt].append(phi_ss_std)
          return p_bar, p_std
docs: Gets the vector of mean and variance of phi values corresponding to
          the specific model and returns them as dictionaries