- 积分
- 65
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2015-10-14
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
求助,为什么画涡度场的时候写了高度场就报错,不写上500hpa就好了(标红部分),我想画500hPa涡度场的图。始终没找到问题所在呀,小白求助!!!
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"
begin
;************************************************
; variable and file handling
;************************************************
fu = addfile("E:/hyd/uwnd.2020.nc","r") ; open netcdf file
fv = addfile("E:/hyd/vwnd.2020.nc", "r")
u = fu->uwnd(203:203,{500},:,:) ; pull u off file
v = fv->vwnd(203:203,{500},:,:) ; pull v off file
;************************************************
; calculate vorticity on a Gaussian Grid
; scale for visual convenience
;************************************************
scale = 1.e05
vrt = u ; retain coordinates
vrt = uv2vrG_Wrap(u,v) * scale
vrt@long_name = "vorticity"
vrt@units = "scaled"
;************************************************
; create plot
;************************************************
wks = gsn_open_wks("eps","E:/hyd/vort") ; send graphics to PNG file
res = True ; plot mods desired
res@cnFillOn = True ; turn on color
res@cnFillPalette = "BlWhRe" ; set color map
res@cnLinesOn = False ; no contour lines
res@mpMinLatF = 0 ; range to zoom in on
res@mpMaxLatF = 60
res@mpMinLonF = 60
res@mpMaxLonF = 150
res@cnLevelSelectionMode = "ManualLevels" ; set manual contour levels
res@cnMinLevelValF = -4 ; set min contour level
res@cnMaxLevelValF = 4 ; set max contour level
res@cnLevelSpacingF = 0.5 ; set contour spacing
plot = gsn_csm_contour_map(wks,vrt(0,2,:,:),res) ; create plot
;************************************************
end
|
|