任意两点间斜剖面做法总结【利用GRADS自带函数】

river

Ivy 发表于 2013-6-18 15:04

                               
登录/注册后可看大图

没有正负号，我觉得还是得改一下吧？

为什么没有呢？如果没有应该也没有问题的啊，因为那是一个比例的关系，分子分母都是同号的。你再仔细想一下，你投影的方向虽然和我的是反过来的，那选取坐标的原点的时候也是反过来的，所以是一样的。

Ivy

river 发表于 2013-6-18 16:14

                               
登录/注册后可看大图

为什么没有呢？如果没有应该也没有问题的啊，因为那是一个比例的关系，分子分母都是同号的。你再仔细想一 ...

好的，我再仔细想想，再次感谢！

river

Ivy 发表于 2013-6-18 20:04

                               
登录/注册后可看大图

好的，我再仔细想想，再次感谢！

不客气，提出问题说明你在思考，这样也会促使我自己思考。这是好事儿，欢迎多提出问题

Ivy

river 发表于 2013-6-18 21:20

                               
登录/注册后可看大图

不客气，提出问题说明你在思考，这样也会促使我自己思考。这是好事儿，欢迎多提出问题

嗯，好的。

小船

好东西，顶一个，谢楼主了，正需要

风儿飘飘

非常好，学习中

小小reed

楼主，我用你的代码出不来结果呢？觉得没用错哇~
'reinit'
'open F:\Products\5min_3\mref\ref\rad.ctl'
'set grads off'
'set  time 17:00Z19MAR2013'
'set x 1'
'set y 1'
'set lev 0.3 20'
lon1 = 109
lon2 = 110
lat1 = 29
lat2 = 25.5
lon = lon1
'collect 1 free'
while (lon <= lon2)
  lat = lat1 +(lat2-lat1)*(lon-lon1) / (lon2-lon1)
  'collect 1 gr2stn(mdbz,'lon','lat')'
     lon=lon+0.01                                                                                                                                                  
endwhile
'set map 15 1 9'
*'set parea 1 10.2 1.5 8'
*'set font 0'
'set lon 109 110'
*'set xaxis 'lon1' 'lon2
'set clab on'
'set gxout shaded'
'set ccols 14 9 4 11 5 13 3'
'd coll2gr(1)'
'cbar'

我的数据是雷达数据，Z坐标，0.01°的，ctl定没错，因为我能用它画出某经圈或纬圈的剖面
xdef 1499 linear 103.91 0.01
ydef 1081 linear 21.397 0.01
zdef 29 levels  
   0.30000
   0.50000
   0.80000
   1.00000
   1.20000
   1.50000
   2.00000
   2.50000
   3.00000
   3.50000
   4.00000
   4.50000
   5.00000
   5.50000
   6.00000
   7.00000
   8.00000
   9.00000
  10.00000
  11.00000
  12.00000
  13.00000
  14.00000
  15.00000
  16.00000
  17.00000
  18.00000
  19.00000
  20.00000
tdef 12 linear 13:00Z19MAR2013 30mn
vars 1
mdbz 29 99 qiya
endvars

river

小小reed 发表于 2013-7-10 17:51

                               
登录/注册后可看大图

楼主，我用你的代码出不来结果呢？觉得没用错哇~
'reinit'
'open F:\Products\5min_3\mref\ref\rad.ctl'
...

提问得说清楚问题啊，你说的出不来是怎么个意思，运行提示有错？那错误提示呢。还是说运行没报错，但也没出图，总得说清楚么。还有就是你没搞清楚coll2gr的用法吧，后面还有个参数呢，那个应该是不可以省略的吧。楼主先自己解决一下吧，接下来一个周我都出差，没时间来回答你的问题了。

Lorenz混沌

真是好东东

Lorenz混沌

但是好像这个gs文件中hgt文件没有用到啊？？我看 动力论坛上有一句'define dz=orog.5(t=1)-hgtprs.4(t=1)'？？？ http://bbs.lasg.ac.cn/bbs/viewth ... 2%E2%2B%C6%CA%C3%E6，river，此orog是不是不一样？