Hello,
In TOMAWAC, the peak enhancement factor is a parameter you can modify in the steering file, via the keyword BOUNDARY PEAK FACTOR and INITIAL PEAK FACTOR. It is used to set the JONSWAP spectrum at the boundaries and for the inital conditions.
The value of 3.3 for the peak enhancement factor is a mean value, which is generally used ad default value. In the literature, generally, the peak enhancement factor varies between 1 and 7.
As far as I know this "gamma" factor is used to take into account the fetch limited wind sea : the higher the peak factor is, the "younger" the sea-state is. With gamma=1 you get the Pierson-Moscowitz spectrum, which represents a fully developed sea-state.
If you are looking for a mathematical formulation of "gamma", you can refer to :
Donelan, M. A., J. Hamilton and W. H. Hui (1985): Directional spectra of wind-generated waves. Phil. Trans. R. Soc. Lond., A315, 509–562.
In this model, to calculate "gamma" you need to provide the wind speed U10 and the fetch length.
Hope it helps,
Giovanni