Hi, > do i=1,lx > do j=1,ly > do k=1,lz Maybe you meant lx1, ly1, lz1 here. Most probably (unless you defined them somewhere else), lx, ly, lz are 0, which would explain this: > but the example.dat file is empty. -- Ismael