求助ncep数据问题。

CDH

'reinit'
'sdfopen k:\ncep\hgt.2012.nc'
'set lon 100 130'
'set lat 30 55'
'set lev 500'
t1=1
'set gxout fwrite'
'set fwrite K:\ncep\fenxing\2012.grd'
WHILE(t1<1464)
'set t 't1''
'd hgt'
t1=t1+4
ENDWHILE
'disable fwrite'
我用上述脚本解出的格点数据，在用他计算的时候我在Fortran里定义hgt(m,n)来读入刚才这个格点数据，那hgt(m,n)这个数组里的m是不是要对应的是经度，n对应的是纬度？不能颠倒过来对吧？
不知道我描述清楚了没有，我现在用hgt（m,n）这个数组去计算，我怀疑我把这个经纬度给搞反了，希望理解我意思的高手给指点一下，谢谢。

river

这个高手低手关系不大吧，你的数据明显除了经纬度变化之外时间也在变化，怎么可能简单的用一个二维数组就能读出来呢

维克多·蓝晶

你要看你fortran是怎么写读取循环语句，先循环m还是先循环n。最好把你的fortran也贴出来。

CDH

river 发表于 2015-1-9 17:03
这个高手低手关系不大吧，你的数据明显除了经纬度变化之外时间也在变化，怎么可能简单的用一个二维数组就能 ...

我是用的循环读入的，读入一个时次计算一次。假如我就读入一个hgt(m,n)那我上边的问题能回答吗？

CDH

维克多·蓝晶 发表于 2015-1-9 17:16
你要看你fortran是怎么写读取循环语句，先循环m还是先循环n。最好把你的fortran也贴出来。

parameter m=11,n=13  ！m代表纬向格点数(30-55N),n代表经向格点数(100-130)
real*4 hgt(m,n)
open(1,file='K:\ncep\fenxing\'//year(yr)//'.grd',form='unformatted',access='direct',recl=m*n)
read(1,rec=dy)hgt
上面是我读入数据的部分语句，下边计算的时候是先循环的m然后循环的n。一般的数据排列不都是行是纬度，然后列是经度吗。参考了计算结果以后我怀疑我计算的时候把经纬度搞反了。

这次计算结束之后我把这个格点数据给输出成txt格式的了，这样可以看见数据，我后续处理起来方便很多，下边是我输出的语句，m和n还是之前定义的大小。
write(3,*)dy  ！日期
do p=1,m
write(3,'(13f8.0)')(hgt(p,q),q=1,n)
enddo
我的终极目的就是想知道我上述的处理方式是不是把NCEP数据的经纬度给弄反了。希望我表达清楚我的问题了。

river

CDH 发表于 2015-1-11 11:18
我是用的循环读入的，读入一个时次计算一次。假如我就读入一个hgt(m,n)那我上边的问题能回答吗？

m是经度，n是纬度

CDH

river 发表于 2015-1-11 17:19
m是经度，n是纬度

好的，谢谢指导。

CDH

river 发表于 2015-1-11 17:19
m是经度，n是纬度

能讲解一下是如何判断的吗？跟着学习一下，谢谢。

river

CDH 发表于 2015-1-12 09:22
能讲解一下是如何判断的吗？跟着学习一下，谢谢。

这个不需要怎么判断啊，你了解grads的资料排列方式的话就应该知道

CDH

river 发表于 2015-1-12 12:23
这个不需要怎么判断啊，你了解grads的资料排列方式的话就应该知道

我使用的其他数据习惯的排列方式都是行是纬度，列是经度，您能给我介绍一下grads的资料排列方式吗？没找到介绍相关问题的资料。谢谢指导。