add functions to use polarization parameters for fitting rather than shape parameters (not quite a success)