- 积分
- 8280
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2014-3-11
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
程序如下:
;************************************************
; pop2lat_1.ncl
;************************************************
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"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/popRemap.ncl"
;************************************************
begin
in = addfile("case01.pop.h.nmonth.1908-01.nc","r")
ice = in->SST
printVarSummary(ice)
print(in)
;************************************************
; convert pop to a 1x1 degree grid
;************************************************
remap = PopLatLon(ice,"gx1v6","1x1d","bilin","da","696451")
printVarSummary(remap)
;************************************************
; create plots
;************************************************
wks = gsn_open_wks("x11","pop2lat") ; open ps file
gsn_define_colormap(wks,"WhBlGrYeRe") ; choose colormap
;
; This will not be necessary in V6.1.0 and later. Named colors can
; be used without having to first add them to the color map.
;
i = NhlNewColor(wks,0.7,0.7,0.7) ; add gray to colormap
res = True
res@gsnPolar = "NH" ; choose hemisphere
res@cnFillOn = True ; turn on color
res@cnLinesOn = False ; no contour lines
res@mpLandFillColor = "grey" ; color of land
res@mpMinLatF = 50 ; adjust map limit
res@cnLevelSelectionMode = "ManualLevels" ; set manual contour levels
res@cnMinLevelValF = 0. ; set min contour level
res@cnMaxLevelValF = 5. ; set max contour level
res@cnLevelSpacingF = 0.5 ; set contour spacing
res@gsnSpreadColors = True ; use full colormap
res@gsnSpreadColorStart = 2
res@gsnSpreadColorEnd = 97
res@tiMainString = "gx1v3 to 1x1d"
; because this is ice data, which has a gap in the tropics, we need to
; explicitly pass the range of the data to plot. Since we are coloring
; the vectors, this range should also match the MinLatF above, since the
; range for the colors is chosen over the full data passed, and not the
; map limits.
plot = gsn_csm_contour_map_polar(wks,remap(0,{50:90},:),res)
end
报错如下:
(0)
ls: cannot access map_*: No such file or directory
(0) create_rmpPopFileName: file not recognized/found: map_gx1v6_to_1x1d_bilin_da_696451.nc
(0) create_rmpPopFileName: dir searched: ./
求大神指教啊!
|
|