u1 = u1 + port_rand_norm()*flt_noise
v1 = v1 + port_rand_norm()*flt_noise
w1 = w1 + 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