[Thread Prev][Thread Next][Index]

Several variables in one unformatted file


Using a Fortran program, I wrote an unformatted file containing four
two-dimensional fields.

The file was written like this:

   real A(64,32),B(64,32),C(64,32),D(64,32)

   write(10) A,B,C,D

I try to read it with ferret using

! Definition of gaussian grid 64x32
file /var=lats/columns=1 lats.gauss.32
define axis/x=0:354.375:5.625/unit=degrees x64
define axis/y/from_data/name=y32/unit=degrees lats
define grid/x=x64/y=y32 g64x32

file /var="T650,T1000,Tdiff,r"/grid=g64x32/format=UNFORMATTED/COLUMNS=`64*32*4` tempmoy.dat

Unfortunately, when I try to display the fields with like this:

fill /line/key/title="Average temperature at 650 mb (1980) (Kelvins)"/i=1:64/j=1:32 T650
fill /line/key/title="Average temperature at 1000 mb (1980) (Kelvins)"/i=1:64/j=1:32 T1000
fill /lin/key/title="T1000 - T650 (1980) (Kelvins)"/i=1:64/j=1:32 Tdiff
fill /line/key/title="r(T650,T1000) (1980)"/i=1:64/j=1:32 r

I get strange results: the bottom part of each plot looks like
temperature data (values range from 240 to 300 approximately) and the
top part looks like it comes from either the third or fourth field.

What am I doing wrong ?  I suspect that I have to change something
either in the way I write the fields or the way I use "file"...

Michel Beland                 beland@CERCA.UMontreal.CA
professionnel de recherche    tel: (514)369-5223  fax: (514)369-3880
CERCA (CEntre de Recherche en Calcul Applique)
5160, boul. Decarie, bureau 400(423), Montreal (Quebec), Canada, H3X 2H9

[Thread Prev][Thread Next][Index]

Dept of Commerce / NOAA / OAR / ERL / PMEL / TMAP

Contact Us | Privacy Policy | Disclaimer | Accessibility Statement