<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">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="gmail_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>