uu = uu + uu*(PORT_RAND(seed)-0.5)*flt_noise
vv = vv + vv*(PORT_RAND(seed)-0.5)*flt_noise
ww = ww + ww*(PORT_RAND(seed)-0.5)*flt_noise
u1 = u1 + u1*(PORT_RAND(seed)-0.5)*flt_noise
v1 = v1 + v1*(PORT_RAND(seed)-0.5)*flt_noise
w1 = w1 + w1*(PORT_RAND(seed)-0.5)*flt_noise
Real*8 PORT_RAND
u2 = u2 + u2*(PORT_RAND(seed)-0.5)*flt_noise
v2 = v2 + v2*(PORT_RAND(seed)-0.5)*flt_noise
w2 = w2 + w2*(PORT_RAND(seed)-0.5)*flt_noise
u3 = u3 + u3*(PORT_RAND(seed)-0.5)*flt_noise
v3 = v3 + v3*(PORT_RAND(seed)-0.5)*flt_noise
w3 = w3 + w3*(PORT_RAND(seed)-0.5)*flt_noise
u4 = u4 + u4*(PORT_RAND(seed)-0.5)*flt_noise
v4 = v4 + v4*(PORT_RAND(seed)-0.5)*flt_noise
w4 = w4 + w4*(PORT_RAND(seed)-0.5)*flt_noise
real*8 FUNCTION PORT_RAND(seed)