On Tue, 22 Jul 2008 16:52:00 +0000 (UTC), Dr Ivan D. Reid
brunel.ac.uk>
wrote in loki.brunel.ac.uk>:> On Tue, 22 Jul 2008 17:28:50 +0200, fred
> wrote in <87tzeiuedp.fsf@free.fr>:
>
>> fread(&x, 1, sizeof(float), file_in);
>> fwrite(&x, 1, sizeof(float), file_out);
You seem to have got the arguments the wrong way around; you
are reading 4 1-byte items. From man fread:
size_t fread(void *restrict ptr, size_t size, size_t nitems,
FILE *restrict stream);
>> What am I doing wrong in the fortran version ?
>
> Have you tried something like:
> integer :: i, j, k
> integer :: nx, ny, nz
> real :: x(720)
> nx = 600; ny = 800; nz = 720
