三维读取

云中月影

'reinit'
'sdfopen d:\nah\hgt.mon.mean.nc'
'set gxout fwrite'
'set fwrite d:\nah\hgt6.grd'
'set lev 100'
'set lat 10 50'
'set lon 10 140'
i=6
while(i<=756)
'define h=hgt(t='i')'
'd h'
i=i+12
endwhile
'disable fwrite'
这个是我准备导出每年6月的纬度10~50，经度10~140的100hPa高度场的gs，导出后应该是个三维吧。
然后再用fortran读的时侯老师出错，不知道三维的怎么个读法

topmad

你仔细安排一下数据，按照什么样的方式  三维数组其实很好理解的，你把每一层的高度场想成一张纸，长宽分别对应经纬度，有几层就是几张纸 然后就组织一下好了

lqouc

如果你读取的时候用form='unformatted'的话，空间隐式循环，时间显示循环就好了。

mofangbao

'set lat 10 50'
'set lon 10 140'
改成set x set y的形式，要不然你无法准确确定提取出来的x是多少点 y是多少点 当然读取会出问题