- 积分
- 392
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2016-2-15
- 最后登录
- 1970-1-1
|
楼主 |
发表于 2018-5-4 16:22:14
|
显示全部楼层
这是我的脚本
begin
nt=1860
nlat=72
nlon=128
preci=fbindirread("/public1/home/hehuanran/ECMWF/data/ld0_20_precip.grd", 0, (/nt,nlat,nlon/),"float")
lat=fspan(53.5,18,72)
lon=fspan(72,135.5,128)
lat@units="degrees_north"
lon@units="degrees_east"
preci@units="mm"
preci!0="t"
preci!1="lat"
preci!2="lon"
preci&lat=lat
preci&lon=lon
preci=where(preci.lt.0, -32767, preci)
preci@_FillValue= -32767
wks=gsn_open_wks("eps", "fffff")
; gsn_define_colormap(wks,"matlab_hsv")
res=True
res@gsnDraw = False
res@gsnFrame = False
res@gsnAddCyclic = False
res@mpMinLatF = 15
res@mpMaxLatF = 55
res@mpMinLonF = 73
res@mpMaxLonF = 135
res@mpGridLatSpacingF=5
res@mpGridLonSpacingF=10
res@pmTickMarkDisplayMode="Always"
res@mpDataBaseVersion="MediumRes" ;中等分辨率
res@mpDataSetName="Earth..4" ;第4版地图,有中国边界数据
res@mpGeophysicalLineColor="black"
res@mpGeophysicalLineDashPattern=16
res@mpGeophysicalLineDashSegLenF=0.2
res@mpGeophysicalLineThicknessF=0.5
res@mpOutlineOn = True
res@mpOutlineSpecifiers=(/"China:states","Taiwan"/) ;中国边界有问题,缺藏南、台湾
res@mpOutlineBoundarySets ="NoBoundaries"
res@mpFillOn = True
res@cnFillDrawOrder = "PreDraw"
;res@mpDataBaseVersion = "MediumRes"
;res@mpDataSetName = "Earth..4"
res@mpAreaMaskingOn = True
res@mpMaskAreaSpecifiers = (/"China","Taiwan","Disputed area between India and China","India:Arunachal Pradesh"/)
res@mpLandFillColor = "white"
res@mpInlandWaterFillColor = "white"
res@mpOceanFillColor = "white"
res@cnLinesOn=False
res@cnFillOn=True
plot = gsn_csm_contour_map(wks,preci(500,:,:),res)
draw(plot)
frame(wks)
end |
|