<html>
<body>
Folks,<br><br>
In many Fortran implementations, unformatted read/write operated on space
or sometimes comma delimited files in ASCII or EBCDIC and were often
human readable. If formats were specified, the input&nbsp; (output) could
be completely run together (3I3 would read 001002003 as three integers 1,
2, and 3). <br><br>
Brandon<br><br>
At 02:41 PM 3/8/2007, you wrote:<br>
<blockquote type=cite class=cite cite>It's been a very long time since I
used Fortran, but as I remember, the &quot;unformatted&quot; I/O was
&quot;system dependent.&quot;&nbsp; A quick Google search turned up some
software for &quot;Fortran Unformatted File Utilities.&quot;<br><br>
<a href="http://www.engineers.auckland.ac.nz/~snor007/software.html#uread">http://www.engineers.auckland.ac.nz/~snor007/software.html#uread</a><br><br>
If you have a user's manual for your Fortran, it might define the format
for &quot;unformatted.&quot;<br>
john<br><br>
<br>
WolfgangZ wrote: <br>
<blockquote type=cite class=cite cite><br>
<pre>Lorenzo Isella schrieb:
&nbsp;
</pre><font face="Courier New, Courier"></font><blockquote type=cite class=cite cite><br>
<pre>Dear All,
I hope this is not too off-topic. I have been given an old but
reliable fortran code for fluid dynamic simulations. It saves a a lot
of data using 3D arrays (q1,q2,q3,pr) as a binary file.
I cut and paste the part of the fortran code saving the data into a
binary file:
(in the following iav=1 and iwrq2=1).


if(iav.eq.1) then
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; namfil='field'//ipfi//'.dat'
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; pnamfil='field'//ipfi//'.dat'
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; open(13,file=namfil,form='unformatted')
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
else

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; pnamfil=filcnw
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; print*, &quot;filcnw is&quot;, filcnw
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; open(13,file=filcnw,form='unformatted')
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
endif
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; print*, &quot;iav is&quot;, iav

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; write(6,*) pnamfil,'written at t=',
&nbsp;&nbsp;&nbsp;&nbsp;
1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; time, ' prma
mi=',prma,prmi
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; nfil=13
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; rewind(nfil)
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; write(nfil) n1,n2,n3
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; write(nfil) ros,alx3d,ren,time
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; if(iwrq2.eq.1) then
c
c&nbsp;&nbsp;&nbsp; large memory occupancy
c
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; print*,&quot;I am saving the
extended results&quot;
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; write(nfil)
(((q1(i,j,k),i=1,n1),j=1,n2),k=1,n3),
&nbsp;&nbsp;&nbsp;&nbsp;
1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
(((q2(i,j,k),i=1,n1),j=1,n2),k=1,n3),
&nbsp;&nbsp;&nbsp;&nbsp;
1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
(((q3(i,j,k),i=1,n1),j=1,n2),k=1,n3),
&nbsp;&nbsp;&nbsp;&nbsp;
1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
(((pr(i,j,k),i=1,n1),j=1,n2),k=1,n3)
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
else
c
c&nbsp;&nbsp;&nbsp; reduced memory occupancy
c&nbsp;&nbsp;&nbsp; the pressure is not necessary for restarting files
c&nbsp;&nbsp;&nbsp; even for post processing can be saved but then
c&nbsp;&nbsp;&nbsp; the advancement of a time step should be performed
c
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; write(nfil)
(((q1(i,j,k),i=1,n1),j=1,n2),k=1,n3),
&nbsp;&nbsp;&nbsp;&nbsp;
1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
(((q3(i,j,k),i=1,n1),j=1,n2),k=1,n3),
&nbsp;&nbsp;&nbsp;&nbsp;
1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
(((pr(i,j,k),i=1,n1),j=1,n2),k=1,n3)
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
endif
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; write(nfil) ntime,ntt,nav
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; close(nfil)



The results is for instance file field0010.dat, which I try reading in
Python by using pylab and the statement:

s = file( './field0010.dat','rb' ).read( )

newarr = fromstring(s ,Float)

but the content of newarr seems absolutely wrong (number of the order
of 1e+309 which are not produced or saved in the simulations...).
Am I doing something wrong?

Kind Regards

Lorenzo
&nbsp;&nbsp;&nbsp;
</pre><font face="Courier New, Courier"></font></blockquote><br>
<pre>

are you sure that fortran saves binary files? Have you&nbsp; checked the 
files with a text editor?

_______________________________________________
SciPy-user mailing list
<a href="mailto:SciPy-user@scipy.org">SciPy-user@scipy.org</a>
<a href="http://projects.scipy.org/mailman/listinfo/scipy-user" eudora="autourl">http://projects.scipy.org/mailman/listinfo/scipy-user</a>


&nbsp;
</pre><font face="Courier New, Courier"></font></blockquote>_______________________________________________<br>
SciPy-user mailing list<br>
SciPy-user@scipy.org<br>
<a href="http://projects.scipy.org/mailman/listinfo/scipy-user" eudora="autourl">http://projects.scipy.org/mailman/listinfo/scipy-user</a></blockquote>
<x-sigsep><p></x-sigsep>
<br>
Brandon C. Nuttall<br><br>
BNUTTALL@UKY.EDU&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
Kentucky Geological Survey<br>
(859)
257-5500&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
University of Kentucky<br>
(859) 257-1147
(fax)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
228 Mining &amp; Mineral Resources Bldg<br>
<a href="http://www.uky.edu/KGS/home.htm" eudora="autourl">http://www.uky.edu/KGS/home.htm</a>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
Lexington, Kentucky 40506-0107 </body>
</html>