NCL大气经向流函数绘制 （三圈环流）

kongfeng0824

转自：http://blog.sina.com.cn/s/blog_4fe145e20101bubl.html

load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"

load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"

load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"

load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"

begin

ps_fi

=

addfile("ps.nc","r")

va_fi

=

addfile("va.nc","r")

ps

=

ps_fi->ps(0,:,:)

va

=

va_fi->va(0,::-1,:,:)

plev

=

va&plev

lat

=

va&lat

msf =

va(:,:,0)

msf =

zonal_mpsi(va,lat,plev,ps)

printVarSummary(msf)

plev_mb = plev/100.0

msf&plev = plev_mb

msf = where(ismissing(msf),msf,msf/1.0e+10)

wks =

gsn_open_wk

s("eps","msf")

printVarSummary(msf)

res =

True

res@vpHeightF

=

0.4

res@vpWidthF

=

0.75

res@trYReverse

=

True

res@gsnYAxisIrregular2Linear

= True

res@gsnXAxisIrregular2Linear

= True

res@gsnContourNegLineDashPat

tern = 1

res@cnLevelSelectionMode = "ExplicitLevels"

res@cnLevels

=

ispan(-20,20,1)

res@cnInfoLabelOn = False

res@tiMainFontHeightF

=

0.025

res@tiXAxisFontHeightF

=

0.025

res@tiYAxisFontHeightF

=

0.025

res@gsnLeftStringFontHeightF

=

0.025

res@gsnRightStringFontHeight

F

=

0.025

res@gsnCenterStringFontHeigh

tF

=

0.025

res@tmYLLabelFontHeightF

=

0.02

res@tmXBLabelFontHeightF

=

0.02

res@tiMainString

=

""

res@tiXAxisString

=

""

res@tiYAxisString

=

"pressure (mb)"

res@gsnLeftString

=

""

res@gsnRightString

=

""

res@gsnCenterString

=

""

plot = gsn_csm_contour(wks,msf,res)

end

zihan6688

我来抢个沙发

sun_shine_Xia

不错，顶一个~~

3870383

学习了学习了

Byby

请问版主NCL有没有求纬向质量流函数的功能

zongqing123

学习了学习了

zongqing123

学习了，学习了

杨小圆

{:eb502:}{:eb502:}{:eb502:}

猫的他

不错，画图的时候可以用gsn_csm_pressure_hegt画

nelflower

怎么控制垂直是否等距呢