- 积分
- 751
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2016-12-5
- 最后登录
- 1970-1-1
|

楼主 |
发表于 2016-12-7 22:11: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
coords = asciiread("./testlonlat.txt",(/437896,2/),"float")
lon = coords(:,0)
lat = coords(:,1)
epsi1 = asciiread("./tb.txt",-1,"float")
;epsi1 = epsi(::10)
;printVarSummary(epsi1)
;print(epsi(10))
;print(epsi1(1))
wks = gsn_open_wks ("pdf","1546-tb-m") ; open workstation
res = True
res@gsnMaximize = True
res@gsnDraw = False
res@gsnFrame = False
res@cnLevelSelectionMode = "ManualLevels"
res@cnMinLevelValF = 160
res@cnMaxLevelValF = 290
res@cnLevelSpacingF = 10
;res@mpProjection = "Satellite" ; choose map projection
;res@mpCenterLonF = 95 ; choose center lon
;res@mpCenterLatF = 40 ; choose center lat
;res@mpSatelliteDistF = 3.0 ; choose satellite view
;res@mpGridLineDashPattern = 2 ; make lat/lon lines dash
;res@mpGridAndLimbOn = True
;res@mpPerimOn = False
;res@mpMinLatF = 35
;res@mpMaxLatF = 40
;res@mpMinLonF = 75
;res@mpMaxLonF = 90
res@sfDataArray = epsi1
res@sfXArray = lon
res@sfYArray = lat
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
求教各位大神,应该怎么进行修改 |
|