How to convert this code to scan random datas instead of binned datas?

Question

I am new to fortran and I am trying to write a code using random datas instead of binned datas in x, y, z as shown in my sample code.

         implicit real*8