ncl批量画轨迹!

对方的凤飞飞

第一次批量读取行数39,列13的txt
begin
   diri="/media/wjj/data/FLEXPART/"
   files = systemfunc ("ls " + diri + "*.txt")
   do i=0, dimsizes(files)-1
   data= asciiread(files(i),(/39,13/),"float")
   end do
   wks = gsn_open_wks("png","traj")               ; send graphics to PNG file
  
   do i=0, dimsizes(files)-1
   想要提取每个txt的所有行,其中的2和3列,这里我有点迷,估计是最近python写多了,想用上[],对ncl的数组和矩阵不是很熟悉,就卡这一步了
      lat(i)=data(i).(:,3)
      lon(i)=data(i).(:,2)
   
      res             = True                              ; plot mods desired
      res@gsnFrame    = False                             ; don't advance yet
      res@gsnMaximize = True
      res@gsnDraw          = False   ; so we can add poly stuff
      res@mpDataBaseVersion     =  "MediumRes" ;
      res@mpPerimOn             = False          ;-- don't draw the box around the plot
      res@gsnPolar="SH"
      res@mpCenterLonF          = 0              ;-- center at lon=10
      res@mpMaxLatF          = -30   
      res@mpGridLineColor       = "grey30"         ;-- set grid line color
      res@tiMainString    = ""  ;
      res@gsnCenterString = ""   ; center string
      res@mpFillOn               = False    ;
      res@mpOutlineOn            = True    ; turn on continental outline


      map = gsn_csm_map(wks,res)                        ; gsn_csm_map_ce
      pres                  = True               ; polyline resources
      pres@gsLineThicknessF = 2.0                ; line thickness
  
      pres@gsLineColor = "red"                        
      gsn_polyline(wks,map,lon(i),lat(i),pres)      ; draw the traj
   end do
   draw(map)                                          
   frame(wks)                                         
  
end
想要批量读取txt数据的经纬度数据,画几条轨迹,
begin
   diri="/media/wjj/data/FLEXPART/"
   files = systemfunc ("ls " + diri + "*.txt")
   do i=0, dimsizes(files)-1
   data= asciiread(files(i),(/39,13/),"float")
   end do
   wks = gsn_open_wks("png","traj")               ; send graphics to PNG file
  
   do i=0, dimsizes(files)-1
   想要提取每个txt的所有行,其中的2和3列,这里我有点迷,估计是最近python写多了,想用上[],对ncl的数组和矩阵不是很熟悉,就卡这一步了
      lat(i)=data(i).(:,3)
      lon(i)=data(i).(:,2)
还请各位看官指点一二!

对方的凤飞飞

自己解决了

wsx

可以分享下解决后的程序吗？谢谢