NCL绘制风场，位势高度场运行中aborted

肉粽子

脚本为：
begin
  a = addfile("./data/Trad_NCAROrb_OldLAI_cplpar_conv2.0_06_climo.nc","r")
  hyam=a->hyam
  hybm=a->hybm
  hyai=a->hyai
  hybi=a->hybi
  printVarSummary(hyam)
  printVarSummary(hybm)
  printVarSummary(hyai)
  printVarSummary(hybi)

  fils = systemfunc("ls ./data/Trad*.nc")
  f = addfiles(fils,"r")
  ListSetType (f,"cat")

  phis = f[:]->PHIS
  tv=f[:]->T*(1.+0.61*f[:]->Q)
  ps=f[:]->PS
  p0=100000
  tv!0="time"
  tv!1="mlev"
  tv!2="lat"
  tv!3="lon"
  time=f[:]->time
  mlev=f[0]->mlev
  lat=f[0]->lat         
  lon=f[0]->lon
  tv&time=time
  tv&mlev=mlev
  tv&lat=lat
  tv&lon=lon
  printVarSummary(tv)
  printVarSummary(ps)
  printVarSummary(phis)
  printVarSummary(p0)

  z = cz2ccm(ps,phis,tv,p0,hyam(::-1),hybm(::-1),hyai(::-1),hybi(::-1))
  u=f[:]->U
  v=f[:]->V
  plevo=(/500/)
  knew = dimsizes(plevo)
  P0mb=1000
  unew=vinth2p(u, hyam, hybm, plevo, ps, 1, P0mb, 1, True)
  vnew=vinth2p(v, hyam, hybm, plevo, ps, 1, P0mb, 1, True)
  znew=vinth2p(z, hyam, hybm, plevo, ps, 1, P0mb, 2, True)
  printVarSummary(unew)
  printVarSummary(vnew)
  printVarSummary(znew)

  utim  = dim_avg_n_Wrap(unew,0)
  vtim  = dim_avg_n_Wrap(vnew,0)
  ztim  = dim_avg_n_Wrap(znew,0)

  printVarSummary(utim)
  printVarSummary(vtim)
  printVarSummary(ztim)

  wks  = gsn_open_wks("png","uv500_wsgd")
  res  = True         
  res@gsnLeftString  = ""
  res@gsnRightString  = ""

  res@pmTickMarkDisplayMode = "Always"
  res@mpFillOn     = False
  res@mpCenterLonF =110
  res@mpMinLatF = -10         
  res@mpMaxLatF = 70
  res@mpMinLonF = 20
  res@mpMaxLonF = 200
  res@mpOutlineOn  = True

  res@vcGlyphStyle = "LineArrow"  
  res@vcMinMagnitudeF =1.0
  res@vcMinDistanceF =0.03
  res@vcLineArrowThicknessF = 2.5
  res@vcLineArrowHeadMaxSizeF =0.01

  res@vcRefAnnoOn = True  
  res@vcRefMagnitudeF = 1      
  res@vcRefAnnoString1On = True        
  res@vcRefAnnoString2On = False         
  res@vcRefAnnoString1 = "10m/s"
  res@vcRefAnnoSide = "Top"
  res@vcRefLengthF = 0.03
  res@vcRefAnnoOrthogonalPosF = -0.12
  res@vcRefAnnoParallelPosF = 0.999

  sres  = True         
  sres@cnLevelSelectionMode = "ManualLevels"
  sres@cnLevelSpacingF = 4
  sres@cnLinesOn = True
  sres@cnLineThicknessF = 3
  sres@gsnLeftString = ""
  sres@gsnRightString = ""         

  plot = gsn_csm_vector_map(wks, utim({500},:,:), vtim({500},:,:), res)
  plot_ov = gsn_csm_contour(wks,ztim({500},:,:),sres)

  overlay(plot,plot_ov)                     
  draw(plot)                              
  frame(wks)                             
end

肉粽子

能运行，但就是没办法出图！

戏子丶

好像是太大了，内存不够用，尽量弄小点，或者删除一些东西

肉粽子

戏子丶 发表于 2018-1-18 19:56
好像是太大了，内存不够用，尽量弄小点，或者删除一些东西

谢谢！其实我不知道怎么做下去了，感觉脚本没写错，可能还是z那一步有错误，因为我后面printVarSummary(z)没有显示，但是我又觉得没有错

戏子丶

肉粽子 发表于 2018-1-21 14:07
谢谢！其实我不知道怎么做下去了，感觉脚本没写错，可能还是z那一步有错误，因为我后面printVarSummary(z ...

脚本是没错。主要是数据太大导致的结果，我之前也是这样的，弄下一点后就好了

qtang1990

ddddddddddddddddddddddddddddddddd

肉粽子

qtang1990 发表于 2018-1-24 23:17
ddddddddddddddddddddddddddddddddd

您是想表达什么？？？？？？？？？？？