Hi there. I am trying to adjust to my particular requirements a Fast Fourier Transform package I found on the web. It works but I need to do some modifications. So, at one moment I place a few lines of code into a subroutine. I need to dump an array into a file to examine it subsequently with gnuplot. It is a very simple array with 2310 members. This is the code:
10 FORMAT (E15.5)
IOstatus = 0
counter = 0
do while (IOstatus.eq.0)
counter = counter + 1
if (MOD(counter,2).ne.0) then
write (2,10) r (counter)
What happens is the program breaks down. There are about 2600 lines of code and it is not immediately apparent at what point. If I comment the above lines out, everything works fine.
There are two ALLOCATE and two DEALLOCATE statements in that routine. I am sure if I rewrite the subroutine without them everything will work with and without this code block.
Why do I have this problem? Thanks, A.