- 积分
- 7130
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2016-7-4
- 最后登录
- 1970-1-1
|
发表于 2016-12-21 14:23:55
|
显示全部楼层
begin
data = asciiread("/home/fc/datafile/baoyu.txt",156,"float")
sta = asciiread("/home/fc/datafile/station156.txt",(/156,3/),"float")
baoyu = data(:)
lon = sta(:,1)
lat = sta(:,2)
olon = fspan(97,109,48)
olat = fspan(25,35,40)
olon!0 = "lon"
olon@long_name = "lon"
olon@units = "degrees-east"
olon&lon = olon
olat!0 = "lat"
olat@long_name = "lat"
olat@units = "degrees_north"
olat&lat = olat
baoyu@_FillValue = -999.000000
rscan = (/1,0.5,0.25/)
final = obj_anal_ic_deprecated(lon,lat,baoyu,olon,olat,rscan,False)
wks = gsn_open_wks ("pdf","baoyu")
gsn_define_colormap(wks,"rainbow")
res = True
res@gsnMaximize = True
res@gsnDraw = False
res@gsnFrame = False
res@tiMainString =""
;>--------------------------------------------<
; set for the map
;>--------------------------------------------<
res@mpMinLatF = 25.
res@mpMaxLatF = 35.
res@mpMinLonF = 97.
res@mpMaxLonF = 109.
res@gsnMajorLonSpacing=3
res@gsnMajorLatSpacing=2
res@mpFillOn = True
res@mpOutlineOn = True
res@cnFillDrawOrder = "PreDraw"
res@mpDataBaseVersion = "MediumRes"
res@mpDataSetName = "Earth..4"
res@mpAreaMaskingOn = True
res@mpMaskAreaSpecifiers =(/"sichuan"/)
; res@mpMaskAreaSpecifiers = (/"China","Taiwan","Disputed area between India and China","India:Arunachal Pradesh"/)
res@mpLandFillColor = "black"
res@mpInlandWaterFillColor = "white"
res@mpOceanFillColor = "white"
res@mpOutlineBoundarySets = "NoBoundaries"
;<================================>;
;>--------------------------------------------<
; set for the plot
;>--------------------------------------------<
res@cnFillOn = True
res@cnLinesOn = False
res@cnLevelSpacingF = 0.5
res@gsnSpreadColors = True
res@lbLabelAutoStride = True
res@gsnAddCyclic = False
res@cnSmoothingOn =True ;--平滑线条
map = gsn_csm_contour_map(wks,final,res)
station = asciiread("/home/fc/datafile/station156.txt",(/156,3/),"float")
res2 = True
res2@gsMarkerIndex = 16
res2@gsMarkerSizeF = 6.
res2@gsMarkerColor = "black"
res2@tfPolyDrawOrder = "PostDraw"
res2@cnFillDrawOrder = "PostDraw"
plots=gsn_add_polymarker(wks,map,station(:,1),station(:,2),res2)
delete(res2)
draw(map)
frame(wks)
end
|
-
|