- 积分
- 9
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2016-9-23
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
本帖最后由 dongfang16 于 2017-6-17 10:06 编辑
求助:
画1971-2000年北极海冰在9,10,11个月的平均趋势。
用ncl画图,在海岸线附近图不连续,可以看到很明显的拐角,是不是画错了?
lat=192,lon=256
;==========================================================================
;==========================================================================
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/wrf/WRF_contributed.ncl"
begin
MODEL = 18
start = 1850
model = "MIROC-ESM-CHEM"
AD = 8
season = "SON"
f1 = addfile("sic_OImon_MIROC-ESM-CHEM_historical_r1i1p1_185001-200512.nc","r")
xa = f1->sic(0,:,:)
xb = f1->sic(0,:,:)
xc = f1->sic(0,:,:)
;==========================================================================
aa = 1971
dd = 1985
x1 = f1->sic((aa-start)*12+AD:(dd-start)*12+AD:12,:,:)
x2 = f1->sic((aa-start)*12+1+AD:(dd-start)*12+1+AD:12,:,:)
x3 = f1->sic((aa-start)*12+2+AD:(dd-start)*12+2+AD:12,:,:)
aa = 1986
dd = 2000
x4 = f1->sic((aa-start)*12+AD:(dd-start)*12+AD:12,:,:)
x5 = f1->sic((aa-start)*12+1+AD:(dd-start)*12+1+AD:12,:,:)
x6 = f1->sic((aa-start)*12+2+AD:(dd-start)*12+2+AD:12,:,:)
;==========================================================================
xa = (dim_avg_n_Wrap(x1,0)+dim_avg_n_Wrap(x2,0)+dim_avg_n_Wrap(x3,0))/45
xb = (dim_avg_n_Wrap(x4,0)+dim_avg_n_Wrap(x5,0)+dim_avg_n_Wrap(x6,0))/45
xc = xb - xa
;==========================================================================
res = True
res@gsnPolar = "NH"
res@mpMinLatF = 60.
res@mpCenterLonF = 90.
res@mpGridLonSpacingF = 30.
res@mpGridLatSpacingF = 10.
res@mpOutlineBoundarySets ="National"
res@mpDataBaseVersion = "MediumRes"
res@cnFillOn = True
res@cnLinesOn = False
res@gsnLeftString = season
res@gsnCenterString = ""
res@gsnRightString = "Trend"
res@tiMainString = model
res@pmLabelBarHeightF = 0.1
res@pmLabelBarWidthF = 0.5
res@cnLevelSelectionMode= "ManualLevels"
res@cnMinLevelValF = -1.5
res@cnMaxLevelValF = 1.5
res@cnLevelSpacingF = 0.15
res@lbLabelFontHeightF = 0.015
res@gsnMaximize = True
;==========================================================================
wks = gsn_open_wks("png",MODEL+"-S1T1-trd-sic-SON-test")
gsn_define_colormap(wks,"BlueWhiteOrangeRed")
gsn_reverse_colormap(wks)
plot = new(1,graphic)
plot = gsn_csm_contour_map_polar(wks,xc,res)
end
|
-
|