- 积分
- 541
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2015-10-19
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
想画东亚夏季风的降水气候态,但是一直报错。见图片。month_to_season这个函数用错了吗
prlr(time,lat,lon)
;************************************
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
;************************************
in = addfile("/lustre/home/moshi/prlr_CanCM3.nc","r")
time=in->time
print(time)
lat=in->lat
lon=in->lon
prlr=in->prlr ; get rid of cyclic points (req to plt)
;print(prlr)
prlr_mm=prlr*(365*24*60*60*1000)
;print(prlr_mm)
;************************************
wks = gsn_open_wks("x11","prlr") ; open a x11 file
gsn_define_colormap(wks,"BlueDarkRed18")
res = True
res@cnFillOn = True ; turn on color fill
res@mpOutlineOn = True ; turn on continental outline(地图轮廓)
res@gsnDraw = False ; do not draw picture
res@gsnFrame = False ; do not advance frame
res@cnLinesOn =False ; turn on contour lines
res@mpMinLatF = -30
res@mpMaxLatF = 65
res@mpMinLonF = -70
res@mpMaxLonF = 110
res@mpCenterLonF = 180
prlr_JJA = month_to_season(prlr_mm,"JJA")
prlr_mmm=dim_avg_n_Wrap(prlr_JJA,0)
print(prlr_JJA)
print(time)
plot = gsn_csm_contour_map_ce(wks,prlr_JJA,res) ; create plot
draw(plot)
frame(wks)
;************************************
end
|
-
|