请选择 进入手机版 | 继续访问电脑版
爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 4302|回复: 9

[分享资料] 对做任意两点的时间斜剖面的一些想法

[复制链接]

新浪微博达人勋

发表于 2017-4-22 19:28:31 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
首先是基于先找点,再写入新的文件,最后画图的思想。
首先一个时次,计算两点之间多个点的经纬度,然后再写入二进制文件,在画图。下面是我的gs和ctl文件以及一张实验图,不知道对不对。
  1. 'reinit'
  2. 'sdfopen E:\uwnd.nc'
  3. 'set gxout fwrite '
  4. 'set fwrite E:\uwnd.dat'
  5. 'set grads off'
  6. 'set zlog on'
  7. 'set csmooth on'

  8. 'set t 1 20'
  9. lon1 = 70.0
  10. lon2 = 120.0
  11. lat1 = 50.0
  12. lat2 = 25.0
  13. lon = lon1

  14. while (lon <= lon2)
  15.   lat = lat1 + (lat2-lat1)*(lon-lon1) / (lon2-lon1)
  16. 'set lat 'lat''
  17. 'set lon 'lon''
  18. 'a=u'
  19. 'set y 1'
  20. 'set lon 'lon''
  21. 'd a'
  22.   lon = lon + 0.5

  23. endwhile
  24. 'disable  fwrite '
  25. 'c'
  26. 'reinit '
  27. ;
复制代码
  1. dset F:\Obliqueprofile\uwnd.dat
  2. undef 9999
  3. title grid data
  4. xdef 101 linear 60 0.5
  5. ydef 1 linear 20 0.5
  6. zdef 1 linear 1000 10
  7. tdef 20 linear jan2003 1yr
  8. vars 1
  9. xie 0 99 xie data
  10. endvars
复制代码

test.png
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2017-4-23 09:12:28 | 显示全部楼层
本帖最后由 river 于 2017-4-23 09:14 编辑

'set y 1'
'set lon 'lon''
纬向为什么是固定到1,'set y 1'为什么就代表了北纬20度?
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2017-4-23 15:27:07 | 显示全部楼层
river 发表于 2017-4-23 09:12
'set y 1'
'set lon 'lon''
纬向为什么是固定到1,'set y 1'为什么就代表了北纬20度?

读取的时候是读取的任意两点的一条斜线,写入时我就想写成一条水平线,这样ctl我才想的出怎么写,所以你明白我的那种意思吗,这样想能实现画图,但是我不知道对不对
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2017-4-26 08:25:05 | 显示全部楼层
chnlhttom 发表于 2017-4-23 15:27
读取的时候是读取的任意两点的一条斜线,写入时我就想写成一条水平线,这样ctl我才想的出怎么写,所以你 ...

我觉得你的思路是对的,但是set y 1不能代表set lat 20 吧,你直接用set lat 20看看结果一样不?
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2017-4-26 15:16:03 | 显示全部楼层
先占个座,看看大神怎么说!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2017-5-1 21:10:42 | 显示全部楼层
期待大神回复
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2017-6-8 14:40:57 | 显示全部楼层
楼主在吗?是沿一条斜线?那横坐标应该同时标注经度和纬度啊?不然图看上去就是沿经或纬线的传播。不知道脚本是否正确 想用用
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2017-6-8 16:27:21 | 显示全部楼层
while (lon <= lon2)

  lat = lat1 + (lat2-lat1)*(lon-lon1) / (lon2-lon1)

'set lat 'lat''

'set lon 'lon''(斜线上等分点的经纬度)

'a=u'

'set y 1'

'set lon 'lon''(可是又set  y 1为什么呢?这个y是否是上一步的lat值呢?如果不是的话那相当于提取沿着y 1的纬线提取数值了啊)

'd a'

  lon = lon + 0.5
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2017-6-14 20:53:22 | 显示全部楼层

不好意思,这几天忙着没上。这个坐标轴我还没有解决,你如果需要私信我你qq号我们可以交流一下
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2017-6-15 12:42:35 | 显示全部楼层
chnlhttom 发表于 2017-6-14 20:53
不好意思,这几天忙着没上。这个坐标轴我还没有解决,你如果需要私信我你qq号我们可以交流一下

我后来用插值法  得到这条直线上的点 然后画出来了 用你的方法也画过 感觉是错的
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表