Re: Reading fortran unformatted data

[Steve Hankin <hankin@xxxxxxxxxxxxx>]

Hi Chris,

Check out the documentation at
http://ferret.wrc.noaa.gov/Ferret/Documentation/Users_Guide/v44/chapter2.html#_1_45

As long as your variable DOUT is single precision floating point you
should find that

   yes? FILE/FORMAT=UNF/VARIABLE=dout/COLUMNS=`nlg*nlt*nz`/GRID=gg
myfile.unf
   yes? CONTOUR dout[l=1,k=1]    ! first time step, first Z level 

does the trick, where "gg" is a 4-dimensional grid defined to match your
multiple time steps of DOUT. If the size of the grid is very large then
you may need to use the
command SET MEMORY/SIZE=xxx in advance of reading your data.

	Happy Ferreting - steve

========================================================================

Chris Picton wrote:
> 
> Hi
> 
> Im am trying to read 4-dimensional data from a file written by a fortran
> program.
> 
> The code which wrote the file looks something like this.
> 
>      WRITE (20)
>      +       (((DOUT(L,K,IZ,J,I,IP),L=1,NLG),K=NLT,1,-1),IZ=1,NZ)
> 
> IP is fixed, and there is a loop over J and I
> 
> ie, it writes J*I records (these are incremental in time), and each record
> is a block of 3-dimensional data.  This is written as an unformatted
> stream.
> 
> Is there any way to read this entire data file in ferret?  Any help would
> be appreciated
> 
> Thanks
> 
> Chris Picton
> chrisp@global.co.za

-- 
Steve Hankin
NOAA/PMEL, 7600 Sand Point Way NE, Seattle, WA 98115-0070
ph. (206) 526-6080 -- FAX (206) 526-6744