怎么样用fortran固定层次

201883301147

dset D:\shuju\hgt.mon.mean.grd
UNDEF  -9.96921e+36f
title monthly mean hgt from the NCEP Reanalysis
xdef 144 linear 0.000 2.500
ydef 73 linear -90.000 2.500
zdef 17 levels 1000 925 850 700 600 500  400  300  250  200  150  100  70  50  30  20  10
tdef 900 linear jan1948 1mo
vars 1
h 17 -999 monthly mean  
endvars

下载了一个NCEP的高度场逐月平均资料描述文件为上，想用fortran固定层次为850hpa，但在fortran里设置数组读取的时候读不出来，想知道一般这种数据是怎么存放的

do k=1,17
         do it=1,72
        do imo=1,12
        do j=1,73
        do i=1,144
        read(3)(hgt1(i,j,imo,it,k))
        enddo;enddo;enddo;enddo;enddo
这是我设置的读取程序

清宵逝水

grd文件由外到内依次是T, 变量V，Z, Y, X，这也对应循环由外到内，fortran变量的索引正好反过来（i, j, k, v, t）

201883301147

清宵逝水 发表于 2022-4-4 19:49
grd文件由外到内依次是T, 变量V，Z, Y, X，这也对应循环由外到内，fortran变量的索引正好反过来（i, j, k,  ...

想问下necp的月平均资料是怎么样排放的呢 我在Fortran一直读取不出来好像陷入死循环了

清宵逝水

201883301147 发表于 2022-4-4 23:47
想问下necp的月平均资料是怎么样排放的呢 我在Fortran一直读取不出来好像陷入死循环了

如果你下载的资料就是grd，而且没有改过ctl，那么资料就是按我说的那个规则排放的呀，如果GrADS读取画图都没有问题，那么就是你的fortran读取程序有问题。请检查你的fortran循环顺序，循环变量范围，变量维度定义是否正确，是不是和你要读取的数据范围一致