<div dir="ltr"><div class="gmail_default" style="font-size:small">Armagan,</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">In general, we recommend that people not use our 3D plotting package, as it is rather old and not anywhere near as developed as other free packages out there. </div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">If you are trying to plot WRF data in 3D, then I recommend VAPOR, which is tailored for plotting WRF data in 3D:</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default"><a href="https://www.vapor.ucar.edu">https://www.vapor.ucar.edu</a><br></div><div class="gmail_default"><br></div><div class="gmail_default">If you really need to use TDEZ3D for other reasons, then please include the full error message you're getting from NCL.</div><div class="gmail_default"><br></div><div class="gmail_default">Thanks,</div><div class="gmail_default"><br></div><div class="gmail_default">--Mary</div><div class="gmail_default"><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Sat, Dec 9, 2017 at 8:16 AM, Armağan Gözütok <span dir="ltr"><<a href="mailto:armagango@gmail.com" target="_blank">armagango@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><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" target="_blank">https://www.ncl.ucar.edu/<wbr>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/<wbr>nclscripts/csm/gsn_code.ncl"<br><br>begin<br><br> a = addfile("wrfout_d01_2017-07-<wbr>27_12:00:<a href="http://00.nc" target="_blank">00.nc</a>","r")<br><br> qgraup = a->QGRAUP(49,:,:,:) <wbr> ; 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_<wbr>isosurface_qgraup")<br><br> tdez3d(wks,xi,yi,zi,qgraup,ui,<wbr>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>
<br>______________________________<wbr>_________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/<wbr>mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>