u1 = u1 + port_rand_norm()*flt_noise
v1 = v1 + port_rand_norm()*flt_noise
w1 = w1 + port_rand_norm()*flt_noise
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
smoothTmpFld(i,j,bi,bj)=port_rand_norm()
_RL port_rand,port_rand_norm
smoothTmpFld(i,j,k,bi,bj)=port_rand_norm()
_RL port_rand, port_rand_norm
port_rand_norm = t*x1