如何利用循环批量生成Overlay图？

河南大学环规院

因为绘图数量较多，因此采用一个循环批量出图。
在overlay位势高度场和温度场时，虽然利用循环进行overlay，且在"x11"模式下可以看到每次overlay的结果，但是生成png文件时仅有最后一张overlay显示出来，请问是什么原因，该如何修改？请大神指教！
代码如下：

begin
a = addfile("Basicmean_field.nc","r")
gh = a->PH
t = a->T

wks = gsn_open_wks("png","test_L41_")

=============位势高度场设置====================
res = True
res@cnFillOn = True
res@cnLinesOn = False
res@lbLabelBarOn = True
res@cnFillPalette         = "BlAqGrYeOrReVi200"

res@gsnLeftString = "Geopotential height & Potential temperature"
res@gsnRightString = ""
res@gsnMaximize = True
res@gsnDraw = False
res@gsnFrame = False

nplots = 2
level = 41
plots = new(nplots,graphic)

do i=0,nplots-1
  plots(i) = gsn_csm_contour(wks,gh(i,level,:,:),res)
end do

=========温度场设置==================
tres = True

tres@cnFillOn = False
tres@lbLabelBarOn = False
tres@cnLevelSpacingF = 4.
tres@cnLineLabelsOn = True
tres@cnLineThicknessF = 4.
tres@cnLineColors = "red"
tres@gsnLeftString = " "
tres@gsnRightString = " "
tres@cnMonoLineDashPattern = True
tres@cnLineDashPattern = 1
tres@gsnDraw = False

tplots = new(nplots,graphic)

do i=0,nplots-1
  tplots(i) = gsn_csm_contour(wks,t(i,level,:,:),tres)
end do

============循环overlay===================
do i=0,nplots-1
  overlay(plots(i),tplots(i))
end do

draw(plots)

frame(wks)

end

小其其格

FILES = systemfunc (" ls -1 " +"./.nc") ;目标文件夹下多个.nc资料
numFILES = dimsizes(FILES)
print("numFILES = " + numFILES)
print(FILES)

do it = 0,numFILES-1               ; TIME LOOP
    a = addfile(FILES(it),"r")
    gh = a->PH
    t = a->T
    wks = gsn_open_wks("png","test_L41_")

   

小其其格

或者你把wks = gsn_open_wks("png","test_L41_")放到overlay循环里面去

河南大学环规院

小其其格 发表于 2019-7-25 10:50
FILES = systemfunc (" ls -1 " +"./.nc") ;目标文件夹下多个.nc资料
numFILES = dimsizes(FILES)
print( ...

这个结果显示fatal:file (a) isn't defined

河南大学环规院

小其其格 发表于 2019-7-25 10:53
或者你把wks = gsn_open_wks("png","test_L41_")放到overlay循环里面去

这个放到循环后显示wks变量为定义，不知道是哪一点出问题了，还望大神进一步解答，谢谢！

小其其格

河南大学环规院 发表于 2019-7-26 16:20
这个结果显示fatal:file (a) isn't defined

你把完整代码发出来

看云淡风清

小其其格 发表于 2019-7-25 10:50
FILES = systemfunc (" ls -1 " +"./.nc") ;目标文件夹下多个.nc资料
numFILES = dimsizes(FILES)
print( ...

您好，我想问下批量出图怎么实现，我用您的代码尝试后发现后面的图会覆盖掉之前的图，问题和楼主的很相似，您知道怎么解决吗？