[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