Hello Costas
What set the default value is the dictionnary (tomawac.dico) and Tomawac_cata_auto.py is automatically built from the dictionnary so he is the one that own the truth.
I am in the documentation process now so it is the moment to correct the manual
when we say that we use a phillip constant of 0.0081 it is only for the option of of limiting waves with strong current (see the routine limite.f) We use a phillip constant of 0.0081 only locally to set a special Spectrum.
Sorry for the confusion
thierry