Foil dynamics formula

[chuanhui-hu]

Hi, thank you for you nice work! I am trying to read the code and learn to use the simulation.
Here I find something confusing to me in the foil dynamics plugin: usv_dynamics/foil_dynamics_plugin/src/foil_dynamics_plugin.cpp

In the file, the lift and drag coefficients are defined as:
cl = this->mult_lift * sin (2*this->alpha);
cd = this->mult_drag * (1 - cos (2 * this->alpha));

But in your paper (Paravisi, Marcelo, et al. "Unmanned surface vehicle simulator with realistic environmental disturbances." Sensors 19.5 (2019): 1068.), the coefficients are:
CL = 1.5*(1-cos(2alpha))
CD = sin(2alpha)

It looks like interchanged.
Should I make any adjustment to the plugin before I simulate with it?