<div dir="ltr">When you use wrf_user_getvar with the diagnostic "eth", that triggers a calculation to compute equivalent potential temperature, which should have values near 300-350 K at the lower levels. None of the numbers you quoted above make sense for equivalent potential temperature, so do you have an already created 'eth' variable in your files that isn't equivalent potential temperature? If so, then don't use wrf_user_getvar to get the variable out, since this is going to trigger the computation (and it appears there may be a unit issue causing those weird numbers). You can use the regular netcdf variable extraction:<div><br></div><div>For example:</div><div><br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">eth = aa[:]->eth</blockquote><div><br></div><div>Hope this helps,</div><div><br></div><div>Bill </div><div><br></div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr">On Tue, Nov 27, 2018 at 6:15 AM Komal Shukla <<a href="mailto:komalshukla1992@gmail.com">komalshukla1992@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div dir="ltr"><div dir="ltr">Hi,</div><div dir="ltr"><br></div><div>I have being using ncl scripts to plot wrf_chem model outputs. It has been working well so far, Recently I needed to plot for a variable named "eth" . I have already plotted many variables of same category, and basically all I had to do was to change the variable name to eth and get the plot.</div><div>When I ran the below mentioned script, my plot is getting values of a range 332000 to 440000, wherein the rreal range of data is in decimals, ranging somewhere around 0.06-0.12.</div><div>I opened the netcdf data files and checked the vales, they were around 0.06-0.12 range only.</div><div><br></div><div>I am not understanding why my plot output is showing such axis values.</div><div><br></div><div>Can you look into this? The script is pasted below.</div><div><br></div><div dir="ltr">begin<div> a = addfile("../wrfout_d02_2016-05-01_00:00:<a href="http://00.nc" target="_blank">00.nc</a>","r")</div><div> it = 0 ; first time step</div><div> hgt = wrf_user_getvar(a,"HGT_M",it) ; Terrain elevation</div><div> hgt@lat2d = wrf_user_getvar(a,"XLAT",it) ; latitude/longitude</div><div> hgt@lon2d = wrf_user_getvar(a,"XLONG",it) ; required for plotting</div><div> wks = gsn_open_wks("png","eth")</div><div> FILES = systemfunc ("ls ../wrfout_d02_2016-05-01*") ; file paths</div><div>numFILES = dimsizes(FILES)<br></div><div> aa = addfiles (FILES+".nc", "r") ;add multiple files</div><div> eth = wrf_user_getvar(aa,"eth",-1) ; 3D tc<br></div><div> eth_avg=dim_avg_n(eth,0)</div><div> eth_avg2D=eth_avg(0,:,:)</div><div> eth_avg2D@lat2d = wrf_user_getvar(a,"XLAT",it) ; latitude/longitude</div><div> eth_avg2D@lon2d = wrf_user_getvar(a,"XLONG",it) ; required for plotting</div><div>;---Set some basic plot options</div><div> res = True</div><div> res@gsnMaximize = True ; maximize plot in frame</div><div> res@cnFillOn = True</div><div>res@cnLinesOn = False<br></div><div> res@mpProjection = "CylindricalEquidistant" ; The default<br></div><div> res@gsnAddCyclic = False</div><div>;---Additional resources desired</div><div> res@pmTickMarkDisplayMode = "Always" ; nicer tickmarks</div><div>res@mpDataBaseVersion = "MediumRes" ; better and more map outlines<br></div><div> res@mpDataSetName = "Earth..4"</div><div> res@mpOutlineBoundarySets = "AllBoundaries"</div><div><div> res@mpOutlineOn = True</div><div>res@lbOrientation = "Vertical"<br></div><div> res@tiMainOffsetYF = -0.03 ; Move the title down</div><div>contour = gsn_csm_contour_map(wks,eth_avg2D,res)<br></div><div>draw(contour)</div><div> frame(wks)</div></div><div>---</div><div><br></div><div>Thanks</div><div><br clear="all"><div><div dir="ltr" class="m_-906682365817601325gmail_signature"><div dir="ltr"><br><div>Komal Shukla</div><div>PhD Research Scholar</div><div>Indian Institute of Technology Delhi (IIT Delhi)</div><div>Hauz Khas, Delhi-110016</div></div></div></div></div></div></div></div>
_______________________________________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu" target="_blank">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/mailman/listinfo/ncl-talk</a><br>
</blockquote></div>