HLSL >> Gerstner Wave

float4 vec = float4(0,0,0,0);

float pW = 2 / waveLength;
float wSpd = (time * phase * pW);
float wDot = (dot(wPos.xy, normalize(windVec).xy) * pW) + wSpd;
float calcSteepness = amplitude * (steepness / (amplitude * pW * 36));

vec.x = cos(wDot) * windVec.x * calcSteepness;
vec.y = cos(wDot) * windVec.y * calcSteepness;
vec.z = sin(wDot) * amplitude;
vec.w = amplitude * pW * pW * phase;    //Amplitude Mask

return vec;