u1 = u1 + PORT_RAND_NORM()*flt_noise
v1 = v1 + PORT_RAND_NORM()*flt_noise
w1 = w1 + PORT_RAND_NORM()*flt_noise
Real*8 PORT_RAND_NORM
u2 = u2 + PORT_RAND_NORM()*flt_noise
v2 = v2 + PORT_RAND_NORM()*flt_noise
w2 = w2 + PORT_RAND_NORM()*flt_noise
u3 = u3 + PORT_RAND_NORM()*flt_noise
v3 = v3 + PORT_RAND_NORM()*flt_noise
w3 = w3 + PORT_RAND_NORM()*flt_noise
u4 = u4 + PORT_RAND_NORM()*flt_noise
v4 = v4 + PORT_RAND_NORM()*flt_noise
w4 = w4 + PORT_RAND_NORM()*flt_noise
EXTERNAL , PORT_RAND_NORM
real*8 FUNCTION PORT_RAND_NORM()