wabs*(vmekal+vmekal)+gaussmodel in which the temperatures are fixed and the NH and norms of the 2 vmekal components are allowed to vary.
Here is a plot of the recent changes in NH vs. phase. The previous cycle variations are plotted in red. Phases are calculated according to the following ephemeris:
HJD X-ray minimum = 1997.95 + 5.53609E
Here is a plot of the 2-10 keV observed flux (uncorrected for NH):