GrADS站点资料作图详细解决方案

勇敢的悲伤

我有一个问题，想问下。对于多变量，多层次的情况里面有一段：
！接下来写高空的10层
do k=1,10
write(2) lev(k),z(i,j,k)
enddo
这个里面的lev(k)是指925,850,700,500,400,300,250,200,150,100吗?这个可以直接用吗？z(i,j,k)又是代表什么呢？前面程序好像没有提及到呀

mofangbao

勇敢的悲伤 发表于 2014-7-25 20:07
我有一个问题，想问下。对于多变量，多层次的情况里面有一段：
！接下来写高空的10层
do k=1,10

lev(k)你可以自己按照自己的理解试试看哈~
z当然代表的是变量 某曾某个位置

勇敢的悲伤

mofangbao 发表于 2014-7-25 21:39
lev(k)你可以自己按照自己的理解试试看哈~
z当然代表的是变量 某曾某个位置

如果temp既有地面资料又有高空资料，那下面这个程序temp（i,j）是不是应该换成地面temp的资料，k应该是do k=2,11，temp(i,j)要改成temp（i,j,k）吗？
do j=1,sta
Write(2)stid(j),lat(j),lon(j),tim,nlev,flag,rain(i,j),temp(i,j)
!接下来写高空的10层
do k=1,10
   write(2)lev(k),z(i,j,k)   !这里的lev和 z也都要定义为 real型
enddo
enddo

mofangbao

勇敢的悲伤 发表于 2014-7-25 22:54
如果temp既有地面资料又有高空资料，那下面这个程序temp（i,j）是不是应该换成地面temp的资料，k应该是do ...

多试试

勇敢的悲伤

mofangbao 发表于 2014-7-26 09:14
多试试

好的。
还有如果没有地面资料的话，有11层高空（/1000,925,850,700,500,400,300,250,200,150,100/)的话，是把前面一句给删除掉，直接开始循环写吗？像这样：

time=0.0
nlev=11 （这个nlev=11对吗？还用加地面的那一层吗？）
flag=1
do j=1,sta
do k=1,11
write(2)  sta(j),lat(j),lon(j),tim,nlev,nflag,lev(k),temp(j,k) (你给的程序说明里面是和地面资料一起写的，如果是高空的话，前面颜色部分还要写进去吗？)
enddo
enddo
nlev=0
write(2) sta,lat,lon,tim,nlev,nflag

mofangbao

勇敢的悲伤 发表于 2014-7-26 14:03
好的。
还有如果没有地面资料的话，有11层高空（/1000,925,850,700,500,400,300,250,200,150,100/)的话 ...

你没有明白我的意思哈  我是说当你遇到不确定的时候 你可以都试试 就知道哪个结果正确了  这样比问别人要记得牢 也能明白原理

sunground

这个是个好东西！

勇敢的悲伤

mofangbao 发表于 2014-7-26 15:26
你没有明白我的意思哈  我是说当你遇到不确定的时候 你可以都试试 就知道哪个结果正确了  这样比问别人要 ...

恩恩。我试验了几次终于可以了，而且也更懂这个方法了。
谢谢啦~~

四叶草

超级大坏蛋 发表于 2014-7-14 08:52
急求帮助，我做到这一步怎么没有数据显示，同样是用的你给的160站点数据，数据格式就像图上这样，请高手指 ...

刚开始确实不好理解，但是你自己揣摩之后，会记忆尤深。

小二

终于下到了