Re: fbinread can't read binary file created by Fortran correctly

Xitain,

It seems there is a compiler option to do the trick:

 gfortran -frecord-marker=4 your-fortran-code

It seems on newer compilers, the default is 4,
but on old compiler, it is 8.

Compile your code with 4, it will work.

Wei

huangwei@ucar.edu
VETS/CISL
National Center for Atmospheric Research
P.O. Box 3000 (1850 Table Mesa Dr.)
Boulder, CO 80307-3000 USA
(303) 497-8924

On Jun 13, 2012, at 11:07 AM, Xitian Cai wrote:

> Hi Wei,
>
> Thank you very much for your quick reply. Unfortunately, all our machines only have gfortran 4.1.2. It is better for me switch to ifort.
>
> Best,
>
> Xitian
>
> On Wed, Jun 13, 2012 at 12:23 PM, Wei Huang <huangwei@ucar.edu> wrote:
> Xitian,
>
> Check your gfortran version. You may have a very old gfortran compiler.
>
> Mine is:
>
> gfortran -v
> Using built-in specs.
> Target: x86_64-apple-darwin10
> Configured with: ../gcc-4.5-20090910/configure --enable-languages=fortran,c++ --host=x86_64-apple-darwin10 --build=x86_64-apple-darwin10
> Thread model: posix
> gcc version 4.5.0 20090910 (experimental) (GCC)
>
>
> I checked on a machine with gfortran 4.1.2, which has the similar problem.
>
> Wei
>
> huangwei@ucar.edu
> VETS/CISL
> National Center for Atmospheric Research
> P.O. Box 3000 (1850 Table Mesa Dr.)
> Boulder, CO 80307-3000 USA
> (303) 497-8924
>
>
>
>
>
>
> On Jun 13, 2012, at 9:48 AM, Xitian Cai wrote:
>
>> Hi All,
>>
>> I'm writing a NCL script to read regular binary files created by FORTRAN. For example, when I compiled the FORTRAN code using ifort compiler and write the binary file, the NCL script is working fine to read it, however when I use gfortran, fbinread got the first value as 0 and neglect the last value. The attached are the example FORTRAN code, NCL script and the binary file created by the FORTRAN code using FORTRAN compiler.
>>
>> (1) Compile code: gfortran write_binary_test.f
>> (2) Run a.out, which creates the binary file Test_WrittenByF.bin
>> (3) Run "ncl ReadWrite_Binary_NCL.ncl" which will print the data.
>>
>> The result should be: 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15
>> But I got:
>> (0,0) 0
>> (0,1) 1
>> (0,2) 2
>> (0,3) 3
>> (0,4) 4
>> (1,0) 5
>> (1,1) 6
>> (1,2) 7
>> (1,3) 8
>> (1,4) 9
>> (2,0) 10
>> (2,1) 11
>> (2,2) 12
>> (2,3) 13
>> (2,4) 14
>>
>> Does anyone know the problem?
>>
>> Thanks,
>>
>> Xitian
>> <ReadWrite_Binary_NCL.ncl><write_binary_test.f><Test_WrittenByF.bin>_______________________________________________
>> ncl-talk mailing list
>> List instructions, subscriber options, unsubscribe:
>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed Jun 13 15:14:26 2012