How can I generate complex valued white noise in fortran?

[galaxy0625]

I have a complex function (or matrix in fortran) f(x,y,t) (where x and y are spacial variables) which is solution of dynamical system and I needed to add some noise which white noise in time and delta correlated in space. How could I generate this noise in fortran?