1subroutine noisegen(d4,nmax)
2
3  real*4 d4(4,nmax)
4
5  do i=1,nmax
6     d4(1,i)=gran()
7     d4(2,i)=gran()
8     d4(3,i)=gran()
9     d4(4,i)=gran()
10  enddo
11
12  return
13end subroutine noisegen
14