三维波作用通量ncl脚本分享

fengyy789

谢谢楼主分享~

fengyy789

戴眼镜的猫 发表于 2015-12-3 16:33
在使用过程中发现有一个变量有问题bsmon未定义，目前我看完程序初步推测这个变量应该是smon。如有更正，另 ...

求问楼主，.nc文件里提到的两个hgt文件, 一个是monthly-climate hgt；另一个是和air,hgt,uwind,vind一波的用来做计算的hgt， 它俩是一样的嘛？ 相关代码如下：
;  direcotry name
diri = "/e4b/ncep/plev/monthly/nc/"
dirci = "/e4b/ncep/plev/monthly-climate/nc/"
;
; monthly mean data
; modify direcotry name
filename = systemfunc("ls "+diri+"hgt.mon.mean.nc")
zfile = addfile(filename,"r")

; monthly climatology data
filename = systemfunc("ls "+dirci+"air.mon.ltm.nc")
btfile = addfile(filename,"r")
filename = systemfunc("ls "+dirci+"hgt.mon.ltm.nc")
bzfile = addfile(filename,"r")
filename = systemfunc("ls "+dirci+"uwnd.mon.ltm.1968-1996.nc")
bufile = addfile(filename,"r")
filename = systemfunc("ls "+dirci+"vwnd.mon.ltm.1968-1996.nc")
bvfile = addfile(filename,"r")

fengyy789

戴眼镜的猫 发表于 2015-12-3 16:33
在使用过程中发现有一个变量有问题bsmon未定义，目前我看完程序初步推测这个变量应该是smon。如有更正，另 ...

求问楼主，.nc文件里提到的两个hgt文件, 一个是monthly-climate hgt；另一个是和air,hgt,uwind,vind一波的用来做计算的hgt， 它俩是一样的嘛？ 相关代码如下：
;  direcotry name
diri = "/e4b/ncep/plev/monthly/nc/"
dirci = "/e4b/ncep/plev/monthly-climate/nc/"
;
; monthly mean data
; modify direcotry name
filename = systemfunc("ls "+diri+"hgt.mon.mean.nc")
zfile = addfile(filename,"r")

; monthly climatology data
filename = systemfunc("ls "+dirci+"air.mon.ltm.nc")
btfile = addfile(filename,"r")
filename = systemfunc("ls "+dirci+"hgt.mon.ltm.nc")
bzfile = addfile(filename,"r")
filename = systemfunc("ls "+dirci+"uwnd.mon.ltm.1968-1996.nc")
bufile = addfile(filename,"r")
filename = systemfunc("ls "+dirci+"vwnd.mon.ltm.1968-1996.nc")
bvfile = addfile(filename,"r")

面南者

楼主，我想计算平面的波活动通量，可以吗

不会游泳的fish

想问下楼主，如果我的量都是异常场，可以直接用异常场来算吗还是换成原场呢。我看程序里好像高度场是用的距平，其他量是气候场？

20141170086

不会游泳的fish 发表于 2019-12-18 15:06
想问下楼主，如果我的量都是异常场，可以直接用异常场来算吗还是换成原场呢。我看程序里好像高度场是用的距 ...

你好，我最近也在计算三维波作用量，我理解的是应该不能用异常场放进去哦，里面有些计算用到了uvtde气候平均的量哇，想请问一下您弄懂没有哦，我不懂为什么最后要把u小于5的地方设为缺测呢？
Mask out where westerlies is small or negative (less than 5 m/s).
by using mask

;x-component of (38)
Fx = mask( levtmp/1000./(2.*cumag*re*re)*( cuvar/coslattmp * xuterm + cvvar * xvterm), cuvar.lt.5,False)

;y-component
Fy = mask(levtmp/1000./(2.*cumag*re*re)*( cuvar*xvterm + coslattmp*cvvar*yvterm),  cuvar.lt.5,False)

;z-component
Fz = mask(levtmp/1000.*ftmp*ftmp/(2.*cumag*NN*re)*(cuvar*zuterm + coslattmp*cvvar*zvterm), cuvar.lt.5, False)

a123456

Dannyqiu 发表于 2019-5-7 13:42
我也是....要怎么合并画图呢?请教楼主，谢谢~

你好，请问你解决了吗？我要画500hpa高度场上的，是不是就是用x分量，y分量的nc数据，确定高度时间之后画矢量场就可以呢

Dannyqiu

a123456 发表于 2020-6-6 18:16
你好，请问你解决了吗？我要画500hpa高度场上的，是不是就是用x分量，y分量的nc数据，确定高度时间之后画 ...

是的，不用z分量就行

a123456

Dannyqiu 发表于 2020-6-6 21:25
是的，不用z分量就行

谢谢了

wsx

fengyy789 发表于 2019-11-12 20:17
求问楼主，.nc文件里提到的两个hgt文件, 一个是monthly-climate hgt；另一个是和air,hgt,uwind,vind一波 ...

这两种nc资料的差别和第二种ltm的资料应该如何下载呢？