顺序存取数据的Fortran.pptVIP

顺序存取数据的Fortran生成几个变量的数据长度不一致 real VAR1(1000), VAR2(500) ... OPEN(UNIT=20,FORMAT=UNFORMATTED,ACCESS=SEQUENTIAL,FILE=MYFILE.DAT) ... WRITE (20) VAR1 WRITE (20) VAR2 .... DEFINE AXIS/X=1:500:1 xaxis DEFINE GRID/X=XAXIS mygrid FILE/FORMAT=stream/SKIP=1003/GRID=mygrid/VAR=var2 myfile.dat 1003表示纪录1中1000个数据与2个信息记录(2个为记录1的周边信息，另一个为记录2的开始信息)The argument 1003 is the sum of the 1000 data words in record 1, plus 2 words of record length information surrounding the data values in record 1 (variable var1), plus 1 word of record information preceding the data in record 2. 直接存取数据的生成 REAL VARI(10), VAR2(10), VAR3(10) ... OPEN(UNIT=20,FORMAT=UNFORMATTED,ACCESS=‘Direct,FILE=MYFILE.DAT) ... DO 10 I=1,10 WRITE (20,rec=i) VAR1(I), VAR2(I), VAR3(I) 10 CONTINUE ... file/format=stream/var=var1,var2,var3 myfile.dat 直接存取数据的生成 REAL*4 MYVAR(10,5) ... OPEN(UNIT=20, FILE=myfile.dat, ACCESS=DIRECT, RECL=10) ... DO 100 j = 1, 5 100 WRITE (20,REC=j) (MYVAR(i,j),i=1,10) .... DEFINE AXIS/X=1:10:1 x10 DEFINE AXIS/Y=1:5:1 y5 DEFINE GRID/X=x10/Y=y5 g10x5 FILE/VAR=MYVAR/GRID=g10x5/FORMAT=stream myfile.dat If the file consisted of a set of FORTRAN REAL*8 or C doubles, then the data would look like: dddddddd dddddddd dddddddd ... and the following Ferret commands would read the data into myvar: DEFINE AXIS/X=1:10:1 x10 DEFINE AXIS/Y=1:5:1 y5 DEFINE GRID/X=x10/Y=y5 g10x5 FILE/VAR=MYVAR/GRID=g10x5/FORMAT=stream/type=r8 myfile.dat 直接存取格点数据的加载 二进制格点数据,可以先定义x,y,z,t的轴,网格,再按照网格分布形式读入,大体形式如下 ... def axis/x=1:nx:1 xax def axis/y=1:ny:1 yax def axis/z=1:nz:1 zax def axis/t=1:nt:1 tax def grid/x=xax/y=yax/z=zax/t=tax mygrid file/format=stream/var=num,num1/grid=mygrid ($file) Netcdf数据的载入 use 文件目录\文件名 set data 文件目录\文件名 基本操作n部曲 确定数据 use 文件名(必须为Netcdf格式数据)set data 文件名 (或文件被赋予的标示号) 确定所需要的网格范围 set region/x/y/z/t/I/j/k/l 确定准备操作的变量或者变量表达式,绘图con (sst-273.16)*9/5 结果输出frame/file=aaa.gif 关闭数据cancel 文件名(或文件被赋予的标示号) * *