Inconsistency in GW instantaneous frequency computation for strain

Hi,

I noticed the following inconsistent result of the GW strain analysis: Given an output file Rpsi4mode_r7.l1, parameters M=2.7 solar masses and M omega_22 = 0.038, watpy can generate an output file using the following script:

#!/usr/bin/env python3

import math, os
import numpy as np
from watpy.wave.wave import wave

thisdir = os.path.abspath(os.path.dirname(__file__))
rpsi4_filename = "Rpsi4mode22_r7.l1"

# initial data, usually given in tables of NR papers
Momega22 = 0.038 # M_solar
M = 2.7 # M_solar
omega22 = Momega22 / M

# analyze h, output data and generate minimal plot
gw = wave(thisdir, code="bam", filename=rpsi4_filename,
          mass=M, f0=omega22/(math.pi*2))
gw.write_to_txt("h", thisdir)
gw.show_strain(os.path.join(thisdir, "strain.pdf"))

The output is written to Rh_l2_m2_00550.txt. The data is organized as

u/M:0 Reh/M:1 Imh/M:2 Momega:3 A/M:4 phi:5 t:6

Around line 1579-1580 you find the following data

-1.310089570425269079e-01 -3.323043395883315299e+01 -1.813291466060879387e+01 9.900920017518063942e-02 3.785583621030307455e+01 8.925263724628601736e+00 5.745600000000000591e+02
4.102154068564063498e-03 -3.392523472187944122e+01 -1.793441314586844371e+01 9.662303427120182142e-02 3.837401107287734447e+01 8.938477652926124506e+00 5.749248000000000047e+02

The fourth column contains the value M omega_22 which is around 0.096 around u=0. However, I would expect to find a value of 0.038 here. Interestingly, the difference between the numbers is roughly a factor of M, e.g. 0.096 / 2.7 is approx 0.036.

I am suspecting a missing factor of M somewhere in watpys internals, but I could not find any. Also: redid the computations myself independently of watpy and I come to the same (wrong) result.

@agonzalez @mbreschi @sbernuzzi Anybody got a clue about what could be the problem?

Here is the data file (taken from Tullio) for the above script: Rpsi4mode22_r7.l1