Path: csiph.com!newsfeed.hal-mli.net!feeder3.hal-mli.net!newsfeed.hal-mli.net!feeder1.hal-mli.net!news.tele.dk!feed118.news.tele.dk!news.tele.dk!small.news.tele.dk!newsgate.cistron.nl!newsgate.news.xs4all.nl!post.news.xs4all.nl!not-for-mail
Content-Type: text/plain; charset="utf-8"
MIME-Version: 1.0
Content-Transfer-Encoding: quoted-printable
From: David Froger <david.froger@inria.fr>
User-Agent: alot/0.3.5
To: python-list@python.org
References: <ac6f6659-0053-482c-b587-454066a4a515@googlegroups.com>
In-Reply-To: <ac6f6659-0053-482c-b587-454066a4a515@googlegroups.com>
Subject: Re: writing fortran equivalent binary file using python
Date: Thu, 14 Nov 2013 07:57:39 +0100
Precedence: list
Newsgroups: comp.lang.python
Message-ID: <mailman.2585.1384415276.18130.python-list@python.org>
Lines: 77
NNTP-Posting-Host: 2001:888:2000:d::a6
Xref: csiph.com comp.lang.python:59416

Hi Sudheer,

Fortan binary format is not portable, and it's hard to work with it
with different langages, compilers, architectures...

In you're Python code, you also use NetCDF, which solve all these problems.

I would suggest to use only NetCDF files in both Python and Fortran code. (=
So
never use the buildit open and write of Python and Fortran, but always the
NetCDF library).

Hope it helps,
David

Quoting Sudheer Joseph (2013-11-14 01:53:42)
> Hi,
>          I need to write a binary file exactly as written by fortran code=
 below to be read by another code which is part of a model which is not adv=
isable to edit.I would like to use python for this purpose as python has mo=
de flexibility and easy coding methods.
>   =

>   character(40) :: TITLE=3D"122322242"
>   integer :: IWI,JWI
>   real :: XFIN,YFIN,DXIN=3D0.5,DYIN=3D0.5,WDAY(6000)
>   XFIN=3D0.0,YFIN=3D-90.0,NREC=3D1461,DXIN=3D0.5;DYIN=3D0.5;IWI=3D720;JWI=
=3D361 =

>   real,allocatable,dimension(:,:,:) :: VAR1_VAL
>   real,allocatable,dimension(:,:,:) :: VAR2_VAL
> =

>   open(11,file=3Doutf,form=3D'UNFORMATTED')
>   WRITE(11) TITLE
>   WRITE(11) NX,NY,XFIN,YFIN,DXIN,DYIN,NREC,WDAY
>   write(*,'(A10,2f10.3)') "START=3D",VAR1_VAL(1,1,1),VAR2_VAL(1,1,1)
>   write(*,'(A10,2f10.3)') "END=3D",VAR1_VAL(nx,ny,nrec),VAR2_VAL(nx,ny,nr=
ec)
>   do i=3D1,NREC
>   WRITE(11) VAR1_VAL(:,:,i),VAR2_VAL(:,:,i)
>   WRITE(*,'(2I10,f10.3)') NX,NY,WDAY(i)
>   enddo
> =

> My trial code with Python (data is read from file here)
> =

> from netCDF4 import Dataset as nc
> import numpy as np
> XFIN=3D0.0,YFIN=3D-90.0,NREC=3D1461,DXIN=3D0.5;DYIN=3D0.5
> TITLE=3D"NCMRWF 6HOURLY FORCING MKS"
> nf=3Dnc('ncmrwf_uv.nc')
> ncv=3Dnf.variables.keys()
> IWI=3Dlen(nf.variables[ncv[0]])
> JWI=3Dlen(nf.variables[ncv[1]])
> WDAY=3Dnf.varlables[ncv[2]][0:NREC]
> U=3Dnf.variables[ncv[3]][0:NREC,:,:]
> V=3Dnf.variables[ncv[4]][0:NREC,:,:]
> bf=3Dopen('ncmrwf_uv.bin',"wb")
> f.write(TITLE)
> f.write(IWI,JWI,XFIN,YFIN,DXIN,DYIN,NREC,WDAY)
> for i in np.arange(0,NREC):
>     f.write(U[i,:,:],V[i,:,:])
> f.close()
> =

> But the issue is that f.write do not allow multiple values( it allows one=
 by one so throws an error with above code ) on same write statement like i=
n the fortran code. experts may please advice if there a solution for this?
> =

> with best regards,
> Sudheer
> -- =

> https://mail.python.org/mailman/listinfo/python-list