NCL creates binary, read-in error in Fortran.

From: Yi Wang <wang29_at_nyahnyahspammersnyahnyah>
Date: Fri, 21 Apr 2006 12:34:47 -0500

I am sorry to post this again because there is no response in the
first round. Everyone is busy!?

Hi Everyone:

      I take the following example:

<<<Examples
The following will write ntimes records, each of length 64x128 float
words.

       ntimes = 100
      z = new ( (/ntimes,64,128/), float)
           .
      path = "/dummy/file"
      do nt=0,ntimes-1
         fbindirwrite(path,z(nt,:,:) )
      end do
The binary file will be the same as a one open via fortran's open
statement:
      open(..., form="unformatted", access="direct", recl= ) ; recl
is system dependent >>>

************************************************************************
*

    In my case, I use the NCL to create a binary data with one record
number. However, when I use ( my data is 14*myr in dimension,
float. I am using Mac OS f77, same thing for tempest, f77)

         open(32,file='inputdata/Single.OPREC.6498.dat',
      $ form='unformatted',access='direct',recl=4*14*myr)

         The data which is read in from above is wrong, except the
first one.

        Any clue about what is wrong here????

Yi Wang PhD
Center for Climatic Research
University of Wisconsin-Madison

_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Fri Apr 21 2006 - 11:34:47 MDT

This archive was generated by hypermail 2.2.0 : Mon Apr 24 2006 - 09:06:45 MDT