ncl做批量画图时改进问题

597265865

  time = fa->time
  a=2017012200
  b=2017021400
  t=cd_calendar(time, -3)
  c=ind(a.eq.t)
  d=ind(b.eq.t)
  do c=c,d
  h1=h(c,:,:,:)
  air1=air(c,:,:,:)
  air1=air1-273.15
  copy_VarCoords(air(0,:,:,:),air1)
  wks  = gsn_open_wks("png",t(c)+"air")     
  res=True
  res@gsnDraw = False
  res@gsnFrame = False
  res@cnFillOn = True
  res@cnLinesOn = True
  res@cnLineLabelsOn =False
  res@cnFillPalette = "BlAqGrYeOrRe"
  res@gsnLeftString=""
  res@gsnRightString=""
  res@cnLevelSelectionMode="ManualLevels"
  res@cnMinLevelValF=-50
  res@cnMaxLevelValF=0
  res@cnLevelSpacingF=2
  res@tiMainString = t(c)+"  500hPa"
  plot=gsn_csm_contour_map(wks, air1(5,:,:), res)
  draw(plot)
  frame(wks)
  end do

597265865

597265865 发表于 2019-3-25 11:17
time = fa->time
  a=2017012200
  b=2017021400

http://bbs.06climate.com/forum.p ... &extra=page%3D2