- 积分
- 751
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2016-12-5
- 最后登录
- 1970-1-1
|
楼主 |
发表于 2016-12-19 10:52:10
|
显示全部楼层
程序如下
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
coords = asciiread("./lonlat1404_c.txt",(/1502,2/),"float")
lon = coords(:,0)
lat = coords(:,1)
epsi1 = asciiread("./TB1404_c.txt",-1,"float")
wks = gsn_open_wks ("pdf","1404-TB-m") ; open workstation
res = True
res@gsnMaximize = True
res@gsnDraw = False
res@gsnFrame = False
res@cnLevelSelectionMode = "ManualLevels"
res@cnMinLevelValF = 150
res@cnMaxLevelValF = 270
res@cnLevelSpacingF = 5
res@mpMinLatF = 39
res@mpMaxLatF = 42
res@mpMinLonF = 99
res@mpMaxLonF = 104
res@sfDataArray = epsi1
res@sfXArray = lon
res@sfYArray = lat
gsn_define_colormap(wks,"rainbow")
res@cnFillOn = True
res@cnFillMode = "RasterFill" ; Turn on raster color
res@cnLinesOn = False
res@cnLineLabelsOn = False
plot = gsn_csm_contour_map(wks,epsi1,res)
draw(plot)
frame(wks)
end
|
|