<div class="gmail_quote">On Sat, Jun 28, 2008 at 1:44 AM, lantao sun <<a href="mailto:lantaosun@gmail.com">lantaosun@gmail.com</a>> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
Dear all,<br><br>I write a simple file in FORTRAN by:<br><span></span><span><span> open (51,file='uwind.dat',form='unformatted, access='direct',recl=2001*2001*4)</span></span></blockquote><div><br>
</div><div>Which Fortran compiler are you using? The units for specifying the RECL was system dependent before Fortran 90 (could be in bytes or words); with f90 the units are words (4 bytes). I suspect you are using f90, in which case the RECL above is wrong (try recl=2001*2001), having a "gap" (read:junk values) right after the rad variable. You can verify this by checking the size of your resulting file; the correct size is 2001*2001*4 bytes.</div>
<div><br></div><div> Arlindo</div><div><br></div><div><br></div></div>-- <br>Arlindo da Silva<br><a href="mailto:dasilva@alum.mit.edu">dasilva@alum.mit.edu</a>