[Thread Prev][Thread Next][Index]

[ferret_users] Reading fortran stream file



Hi ferret users,

I wrote a stream file of two variable with these fortran commands:

Real*4 myvar1(8,5); Real*4 myvar2(8,5)

open (unit=20, file="test.dat", access="direct", recl=16)

... ... ... ....

myvar1(i,j)= Real(i) + Real(j) * 10.0
myvar2(i,j)=-(Real(i) + Real(j) * 10.0)

... ... ... ....
do 120 j=1, 5

write (20,rec=j), (myvar1(i,j),myvar2(i,j), i=1,8)

120 continue
end


myvar1(i,j) and myvar2(i,j) have same values but with opposite signs.I read the
output file with the ferret commands:

define axis/x=1:8:1 x8
define axis/y=1:5:1 y5

define grid/x=x8/y=y5 g8x5
file/var=myvar1,myvar2/grid=g8x5/format=stream test.dat
show data
list myvar1
list myvar2

I expect to see all positive values in myvar1 and all negative in myvar2. But what I get is


yes? list myvar1
VARIABLE : MYVAR1
FILENAME : test.dat
SUBSET : 8 by 5 points (X-Y)
1 2 3 4 5 6 7 8
1 2 3 4 5 6 7 8

1 / 1: 11.00 -11.00 12.00 -12.00 13.00 -13.00 14.00 -14.00
2 / 2: 15.00 -15.00 16.00 -16.00 17.00 -17.00 18.00 -18.00
3 / 3: 21.00 -21.00 22.00 -22.00 23.00 -23.00 24.00 -24.00
4 / 4:
25.00 -25.00 26.00 -26.00 27.00 -27.00 28.00 -28.00
5 / 5: 31.00 -31.00 32.00 -32.00 33.00 -33.00 34.00 -34.00
yes? list myvar2
VARIABLE : MYVAR2
FILENAME : test.dat
SUBSET : 8 by 5 points (X-Y)

1 2 3 4 5 6 7 8
1 2 3 4 5 6 7 8
1 / 1: 35.00 -35.00 36.00 -36.00 37.00 -37.00 38.00 -38.00
2 / 2: 41.00
-41.00 42.00 -42.00 43.00 -43.00 44.00 -44.00
3 / 3: 45.00 -45.00 46.00 -46.00 47.00 -47.00 48.00 -48.00
4 / 4: 51.00 -51.00 52.00 -52.00 53.00 -53.00 54.00 -54.00
5 / 5: 55.00 -55.00 56.00
-56.00 57.00 -57.00 58.00 -58.00


What is wrong? (test.dat is attached)

Nadeem
*********************************************************************************

--
Imran Nadeem
PhD Student
Institute of Meteorology
Department of Water, Atmosphere and Environment
Uni. of Natural Resources and Applied Life Sciences (BOKU)

Peter-Jordan Strasse 82
1190 Vienna, Austria

Mobile: +43 699 1194 3044
Tel.: +43 1 47654 5614
Fax: +43 1 47654 5610

Attachment: test.dat
Description: MPEG movie


[Thread Prev][Thread Next][Index]

Dept of Commerce / NOAA / OAR / PMEL / TMAP

Contact Us | Privacy Policy | Disclaimer | Accessibility Statement