<div dir="ltr"><div><div><div><div>Hello,<br><br></div>I wanted to create 3d plots recently, bye some modification on the scripts here:<br><a href="https://www.ncl.ucar.edu/Applications/tdpack.shtml">https://www.ncl.ucar.edu/Applications/tdpack.shtml</a><br><br></div>I did manage to plot the surface plots by tdez2d function (plotted HGT from the WRF outputs), however my 3d variables cannot be plotted, getting always an error somehow although i check all the dimensions properly which they coincide.<br></div><div><br></div><div>I couldn't manage to upload an .nc outputs to the ftp server,<br></div><div><br></div><div><br></div>Here is my NCL script to create isosurface plot of graupel mixing ratio (called isosurf.ncl) different than 3d_hgt.ncl:<br><br><br>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"<br><br>begin<br><br>  a = addfile("wrfout_d01_2017-07-27_12:00:<a href="http://00.nc">00.nc</a>","r")<br><br>  qgraup = a->QGRAUP(49,:,:,:)                 ; Graupel mixing ratio at iteration #49<br>  lat = a->XLAT(0,:,0)                ; latitude<br>  lon = a->XLONG(0,0,:)               ; longitude <br><br>  xi = lat<br>  yi = lon<br>  zi = ispan(1,49,1)*1. ;vertical pressure level # 50 different pressure levels<br>  ui = 0.0001 ; graupel mixing ratio defined for isosurface plot<br><br><br>  ;dimlat = dimsizes(lat)<br>  ;dimlon = dimsizes(lon)<br>  dimz = dimsizes(zi)<br>  ;dimgraup = dimsizes(qgraup)<br> <br>  ;print(dimlat)<br>  ;print(dimlon)<br>  print(dimz)<br>  ;print(dimgraup)<br><br><br>  wks = gsn_open_wks("png","3d_isosurface_qgraup")<br><br>  tdez3d(wks,xi,yi,zi,qgraup,ui,1.8,-45.,58.,-4)<br><br>  frame(wks)<br><br>end<br><br><br></div>Armagan<br><div><div><div><br><br></div></div></div></div>