爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 6765|回复: 5

[经验总结] 由画Ekman螺线想到的一点移动默认坐标轴的小技巧

[复制链接]

新浪微博达人勋

发表于 2014-12-7 21:32:13 | 显示全部楼层 |阅读模式

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

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

x
在画Ekman螺线时,ncl默认的图X轴总是处于y=-1的地方,无法穿过最后的那个小弯弯。最后通过分为两步在同一帧叠加两个不同的图,可实现调整X轴往上到y=0出穿过小弯弯。脚本及最后结果图(粗略)如下:load"$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load"$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
begin
u=new(1256,float)
v=new(1256,float)
i=fspan(0,12.55,1256)
  u=1-exp(-i)*cos(i)
  v=exp(-i)*sin(i)
wks=gsn_open_wks("pdf","ekman")
res = True
res@tiMainString = "EkmanSpiral"
;----first step----
res@trYMaxF = 1.0
res@trYMinF = 0.0
;-------Turn off left,right,top axes---
res@tmXTBorderOn = False
res@tmXTOn       = False
res@tmYRBorderOn = False
res@tmYROn       = False
res@tmYLBorderOn = False
res@tmYLOn       = False
res@gsnFrame     = False
;------plot x-axes and set to right position----
u1=fspan(0,0,100)
v1=fspan(0,1.2,100)
res@vpYF = 0.855     ;critical
res@tmXBMinorPerMajor = 3
plot1=gsn_xy(wks,v1,u1,res)
;--------second step---------
u@long_name = "u"
v@long_name = "v"
res1 = True
res1@trYMaxF = 1.0
;-------Turn off right,top,bottom axes----
res1@tmXTBorderOn = False
res1@tmXTOn       = False
res1@tmXBBorderOn = False
res1@tmXBOn       = False
res1@tmYRBorderOn = False
res1@tmYROn       = False
plot2=gsn_xy(wks,u,v,res1)
frame(wks)
end


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

新浪微博达人勋

发表于 2014-12-8 06:46:25 | 显示全部楼层
不错不错
密码修改失败请联系微信:mofangbao
回复

使用道具 举报

新浪微博达人勋

发表于 2016-4-11 18:32:34 | 显示全部楼层
请问一下,脚本中i的取值范围(0,12.55)就是风速大小的范围吧?如比说我有一天每个时刻每个高度的垂直风速大小,如何去画Ekman图呢?是直接在数据中观测最大最小值范围就能定义i吗?
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2016-4-11 18:48:16 | 显示全部楼层
dylan-du 发表于 2016-4-11 18:32
请问一下,脚本中i的取值范围(0,12.55)就是风速大小的范围吧?如比说我有一天每个时刻每个高度的垂直风 ...

那是 4 pi,,,,
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2016-4-11 22:51:11 | 显示全部楼层

那我计算的时候是否可以拿i当成观测的水平速度U,然后通过i(U)计算水平风速的u,v分量?
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2016-4-12 08:27:58 | 显示全部楼层
不错不错
密码修改失败请联系微信:mofangbao
回复

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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