- 积分
- 718
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2012-7-10
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
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
; read in fnl file
files = systemfunc("ls ./fnl2013/fnl*")+".grb"
do i=0,1;,dimsizes(files)-1
fi = addfile(files(i),"r")
; read in zonal winds
hgt = fi->HGT_P0_L100_GLL0
tim1=str_split(files(i),"_")
;tim = hgt@initial_time
tim=tim1(1)+tim1(2)
file_name="500hPa"+tim
; create plot
wks_type = "png"
wks_type@wkWidth = 2500
wks_type@wkHeight = 2500
wks = gsn_open_wks(wks_type,file_name)
gsn_define_colormap(wks,"gui_default")
res = True
res@cnFillOn = True
res@cnLinesOn = True
res@cnLineLabelsOn = True
res@cnInfoLabelOn = False
res@gsnSpreadColors = True ; 6.1.0
res@mpMinLatF = 10
res@mpMaxLatF = 65
res@mpMinLonF = 65
res@mpMaxLonF = 165
res@cnLevelSelectionMode = "ManualLevels"
res@cnMinLevelValF = 5000
res@cnMaxLevelValF = 6000
res@cnLevelSpacingF = 20
res@lbBoxLinesOn = False ; Label bar
res@lbLabelAutoStride = True ; let NCL determine label spacing
res@tiMainString = tim+":"+"500hPa"
plot = gsn_csm_contour_map(wks, hgt(13,:,:), res)
clear(wks)
end do
end
上面是我程序,这样画出来图,为什么文件名不是我想要的。
500hPa2013010100.000001.png
500hPa2013010100.000002.png
500hPa2013010106.000001.png
500hPa2013010106.000002.png
自动加了000001,怎么处理一下呢
|
|