ncl做批量画图时改进问题

上帝也厕所

do while(date_start .lt. "20101231 23" )
year=systemfunc("date +%Y  -d '"+date_start+"'")
month=systemfunc("date +%m -d '"+date_start+"'")
day=systemfunc("date +%d   -d '"+date_start+"'")
hour=systemfunc("date +%H  -d '"+date_start+"'")
  precip=cbinread("FY2E_FDI_ALL_NOM_"+year+""+month+""+day+"_"+hour+"00.rad",(/601,1001/),"double")
res@tiMainString = "GuangXi Radiation 1km Data"
res@gsnLeftString=date_start
res@gsnRightString="W/m~S~2~N~"
plot=gsn_csm_contour_map(wks,precip,res)   
draw(plot)
system("rm -f 1km_gx*.png")
frame(wks)
system("mv 1km_gx*.png "+year+month+day+hour+".png")
date_start=systemfunc("date +%Y%m%d' '%H -d ' "+date_start+" 1 hour'")
print(date_start+"")
end do
这个已经可以画图，但是因为有的数据没有，中间就会停止画图，还得人工改时间才能进行画图。有没有一个更好的方法来进行批量画图，希望大神能知道一下，谢谢

likailing2008

FILES = systemfunc (" ls -1 "+"*.nc ")        ;从这里开始到下面都是读取文件的命令
  numFILES = dimsizes(FILES)
  print("numFILES = " + numFILES)
  print(FILES)
  print (" ")
do i = 0,numFILES-1            
    a = addfile(FILES(i),"r")                
.......................后面的都一样
        print(FILES(i))                               
    tos=a->tos(:,:,:)   

destiny小鱼

我也遇到了同样的问题，能否留下联系方式qq125165813

看天吃饭

我批量画图的脚本，不知道能不能给你帮助。我也是才学，菜鸟级，大家一起讨论。

白兰

正想做批量画图这件事，学习一下！

freekiller

推荐likailing2008的方法，我一直在用，感觉还行，顺求lz的完整脚本绘制风云卫星

andrewsoong

likailing2008 发表于 2013-9-26 11:20

                               
登录/注册后可看大图

FILES = systemfunc (" ls -1 "+"*.nc ")        ;从这里开始到下面都是读取文件的命令
  numFILES = dims ...

求完整的脚本啊~~~

likailing2008

andrewsoong 发表于 2014-2-28 19:20

                               
登录/注册后可看大图

求完整的脚本啊~~~

begin

  FILES = systemfunc (" ls -1 "+"*.nc ")        ;从这里开始到下面都是读取文件的命令
  numFILES = dimsizes(FILES)
  print("numFILES = " + numFILES)
  print(FILES)
  print (" ")


  type = "pdf"     ;开始指定文件输出的类型

  wks = gsn_open_wks(type,"sst_100year_avg")  ;指定文件输出的类型和名称
  gsn_define_colormap(wks,"gui_default")
  plot = new(24,graphic)                                    ;括号中的24表示要画24张图片

        res=True
       
     这里就是res设置
       
;=====================================================================================
    do i = 0,numFILES-1            ; TIME LOOP这个是最重要的，开始循环读入文件，ncl是从0开始，到文件夹所在文件数数量减一结束   
    a = addfile(FILES(i),"r")         ;把所有读入的文件写入一个叫做a的文件中       
        print(FILES(i))                                ;把读入的文件打印出来       
       
        tos=a->tos(:,:,:)
        tos_ave=dim_avg_Wrap( tos(lat|:, lon|:, time|:) )
         plot(i)= gsn_csm_contour_map_ce(wks,tos1(:,:),res)  
  end do      
这个是求海温平均的程序，循环结构就是这样，然后中间的自己设置

andrewsoong

likailing2008 发表于 2014-3-1 22:42

                               
登录/注册后可看大图

begin

  FILES = systemfunc (" ls -1 "+"*.nc ")        ;从这里开始到下面都是读取文件的命令

非常感谢！

zhuhongxia

你好，请问批量读取数据再批量画图（daily资料txt文件），一张图放几张图片，如何出图