Read unformatted data which generated from fortran

From: 程军 <chengjun_at_nyahnyahspammersnyahnyah>
Date: Sat Apr 03 2010 - 12:24:03 MDT

Dear NCL users,

I'm trying to read a "unformatted" FORTRAN data file with NCL, but the my NCL script seems like not work, could someone help me to resolve this problem? Thanks at first.

The FORTRAN program involving the data is as following:

   parameter(nlon=100,nlat=116)

    open(17,file='mask_map_pb.dat',status='unknown',
   & form='unformatted',access='direct',recl=nlon*nlat*4)

    write(17,rec=1) ((mask_map(ii,jj),ii=1,nlon),jj=1,nlat)
    write(32,*) mask_map


The NCL script to read this data file:

fili = "mask_map_pb.dat" ; input file
nrec = 0
dims = (/100,116/)
t = fbindirread(fili, nrec, dims, "float")

Regards,

Jun

夗z{S蕳瓄嗁
夗z{S靰昴澥妜j呛<span class="Apple-style-span" style="font-family: Helvetica; ">Dear NCL users,<br><br>I'm trying to read a "unformatted" FORTRAN data file with NCL, but the my NCL script seems like not work, could someone help me to resolve this problem? Thanks at first.<br><br>The FORTRAN program involving the data is as following:<br><br>&nbsp;&nbsp;&nbsp;parameter(nlon=100,nlat=116)<br><br>&nbsp;&nbsp;&nbsp;&nbsp;open(17,file='mask_map_pb.dat',status='unknown',<br>&nbsp;&nbsp;&nbsp;&amp; &nbsp;form='unformatted',access='direct',recl=nlon*nlat*4)<br><br>&nbsp;&nbsp;&nbsp;&nbsp;write(17,rec=1) ((mask_map(ii,jj),ii=1,nlon),jj=1,nlat)<br>&nbsp;&nbsp;&nbsp;&nbsp;write(32,*) mask_map<br><br><br>The NCL script to read this data file:<br><br>fili = "mask_map_pb.dat" &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;; input file<br>nrec = 0<br>dims = (/100,116/)<br>t &nbsp;&nbsp;&nbsp;= fbindirread(fili, nrec, dims, "float")<br><br>Regards,<br><br>Jun</span>
夗z{S蕳瓄萙0U觎妷
夗z{S靰昴澥妜;
夗z{C娛h+b)逎蒻jY&j)b b.+-妠--妷觳骒r笡z)*'查旃&迒踚鶜奩殲nfX)撸栔
Received on Sat Apr 3 13:02:14 2010

This archive was generated by hypermail 2.1.8 : Wed Apr 07 2010 - 07:12:50 MDT