[Thread Prev][Thread Next][Index]
question on reading multi variables
Dear ferret_users
There is an example in ferret's manual as listed below,
(Example 6—2 variables, 3 dimensions (time series))
yes? DEFINE AXIS/X=1:4:1 xplnt_tmyes
yes? DEFINE AXIS/Y=1:3:1 yplnt_tmyes
yes? DEFINE AXIS/T=1:12:1 tplnt_tmyes
yes? DEFINE GRID/X=xplnt_tm/Y=yplnt_tm/T=tplnt_tm gplant2
yes? FILE/VAR="height,leaf"/GRID=gplant2 green_time.dat
yes? PLOT/X=3/Y=2 height
i think the responding fortran program should be
do t=1:12
do y=1:3
do x=1:4
write(unit,formatted...)height(x,y,z),leaf(x,y,z)
end do
_______________________________________________________________
The question is if i meet the situation following, how should i reat the data in ferret?
certainly i could read all of them as a variable TEMP,then let each equal to the responing part of TEMP,but i don't think it is a good way, is there anybody has good idea?
On the other hand, if the format is UNFORMATTED and the access mode is SEQUENTIAL,I think the difference is just to design form='unformatted',is it right or not ?
do 10 t=1:12
do 20 y=1:3
do 20 x=1:4
write(unit,formatted/UNFORMATTED...)height(x,y,z)
20 continue
do 30 y=1:3
do 30 x=1:4
write(unit,formatted/UNFORMATTED)leaf(x,y,z)
30 continue
10 continue
by the way,is there any difference when considering the double quote""in the following routine,i.e.FILE/VAR="height,leaf"/GRID=gplant2 green_time.dat.
thanks a lot!
yours sincerely,gnzeng
gnzeng
led_rs@scsio.ac.cn
2003-09-16
[Thread Prev][Thread Next][Index]
Dept of Commerce /
NOAA /
OAR /
PMEL /
TMAP
Contact Us | Privacy Policy | Disclaimer | Accessibility Statement