登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
本帖最后由 wcyycw 于 2021-4-20 14:44 编辑
代码: 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 num_st=37 loc =asciiread("./chun.csv",(/num_st,3/),"float") lon = loc(:,0) lat = loc(:,1) R = loc(:,2)
print(lat) print(lon) print(R) ………………等等 R@_FillValue = 999999.000000 rscan = (/10,5,3/)
data1 =obj_anal_ic_deprecated(lon,lat,R,olon,olat,rscan,False)
print(data1)
wks_type = "pdf" wks =gsn_open_wks(wks_type,"chun") ...............等等 res = True res@gsnAddCyclic = False res@mpDataSetName ="/home/ncl/Earth..4" res@mpDataBaseVersion = "MediumRes" res@mpOutlineOn = True res@mpOutlineSpecifiers =(/"China","China:Provinces"/) ……………等
map=gsn_csm_contour_map_ce(wks,data1(:,:),res)
dgres = True dgres@gsLineThicknessF = 2.0 dgres@gsLineColor ="Black"
plot11=gsn_add_shapefile_polylines(wks,map,"/home/boufiles/bou2_4l_cn50.shp",dgres)
draw(map) frame(wks) end
最后图片这样:
种种疑惑。。。
|