<div dir="ltr"><div>The <b>*wrf_user* </b>functions and documentation were created by WRF developers not NCL developers.</div><div>NCL provides an interface to to provided derivation code [usually fortran].</div><div>Hence, we can not support the underlying codes. Please contact WRF.</div><div>========</div><div>That said:<br></div><div><b></b></div><div><br></div><div> NCDF = addfile("<a href="http://wrfout_d02_2017-10-21_00_00_00.nc">wrfout_d02_2017-10-21_00_00_00.nc</a>","r")</div> ht = wrf_user_getvar(NCDF,"z",0)<br><div> printVarSummary(ht) ; ?? [bottom_top | ..] x [south_north | ..] x [west_east | ..]?? height</div><div> tk = wrf_user_getvar(NCDF,"tk",0) <br> printVarSummary(tk) ; ; ?? [bottom_top | ..] x [south_north | ..] x [west_east | ..]?? temperature</div><div><br></div><div> J = 45 ; arbitrary latitude index<br></div><div> I = 50 ; longitude<br></div><div><br></div><div> print( ht(:,J,I)+" "+tk(:,J,I)) ; print all vertical values at a specified location<br></div><div><br></div><div>===</div><div><a href="https://www.ncl.ucar.edu/Document/Functions/WRF_arw/wrf_user_interp_level.shtml"><b>wrf_user_interp_level</b></a> was introduced in the NCL 6.6 release. I do not have NCL 6.6 on my MAC</div><div><br></div><div> temp_1km=<b>wrf_user_interp_level</b>(tk,ht,1000,0)</div><div> printVarSummary( temp_1km )<br></div><div><br></div><div>However, you could print the derived 1000m temperature at the J,K<br></div><div><br></div><div>print( temp_1km(J,K) )</div><div><br></div><div>This should indicate to you if the derived value is correct.</div><div><br></div><div>===<br></div><div><br></div><div><br></div><div><br></div><div><div><br></div></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, Nov 25, 2021 at 10:27 AM Israt Jahan via ncl-talk <<a href="mailto:ncl-talk@mailman.ucar.edu">ncl-talk@mailman.ucar.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div><div><br></div><div dir="auto">Hello,</div><div dir="auto">It would be much appreciated if someone could shed some light on this. Does wrf_user_getvar(NCDF_File,"z",0) measure z from ground level or sea level?</div></div><div dir="auto">I need to find temperature at 1000 m above ground level. I am doing the following:</div><div dir="auto"><div dir="auto">1. ht=wrf_user_getvar(NCDF_File,"z",0) ;</div><div dir="auto">where NCDF_File is a list of WRF netCDF files</div><div dir="auto">2. tk=wrf_user_getvar(NCDF_File,"tk",0) ;get temperature in K</div><div dir="auto">2. temp_1km=wrf_user_interp_level(tk,ht,1000,0) ;</div><div dir="auto">get temperature at 1 km above ground level</div><div dir="auto">Is my approach correct?</div><div dir="auto">Best,</div><div dir="auto">Israt</div></div><div dir="auto"><br></div><div><div dir="auto"><br><div class="gmail_quote" dir="auto"><div dir="ltr" class="gmail_attr">On Tue, Nov 16, 2021 at 4:26 PM Israt Jahan <<a href="mailto:israt.jahan@uconn.edu" target="_blank">israt.jahan@uconn.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr">Dear NCL users,<div>I have a question about the NCL function "wrf_user_getvar". I see this function can extract data from WRF files and calculate the requested diagnostic. One of the available diagnostics is "z/height" which is labeled as the "full model height in m" in this site <a href="https://www.ncl.ucar.edu/Document/Functions/WRF_arw/wrf_user_getvar.shtml" target="_blank">wrf_user_getvar (ucar.edu)</a>. From WRF output files, I need to find temperature at 1000 m above ground level. I am doing the following:</div><div>1. ht=wrf_user_getvar(NCDF_File,"z",0) where NCDF_File is a list of WRF netCDF files</div><div>2. tk=wrf_user_getvar(NCDF_File,"tk",0) to get temperature in K</div><div>2. temp_1km=wrf_user_interp_level(tk,ht,1000,0) to get temperature at 1 km </div><div><br></div><div>My question is: is the diagnostic "z/height" that we get from "wrf_user_getvar" measured from ground level or sea level?</div><div>If it is measured from sea level, what should I do to get z from ground level?</div><div><br></div><div>Thank you so much.</div><div>Best,</div><div>Israt</div></div></div></div></div></div></div></div></div></div></div></div><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div><br><div><br></div>-- <br><div dir="ltr"><div dir="ltr"><div><font size="2" face="arial, sans-serif" color="#000000">Regards,</font></div><div><font size="2" face="arial, sans-serif" color="#000000">Israt Jahan</font></div><div><font size="2" face="arial, sans-serif" color="#000000">Graduate Assistant</font></div><div><font size="2" face="arial, sans-serif" color="#000000">Department of Civil and Environmental Engineering</font></div><div><font size="2" face="arial, sans-serif" color="#000000">School of Engineering | University of Connecticut<br></font></div><div><font size="2" face="arial, sans-serif" color="#000000">Email: <a href="mailto:israt.jahan@uconn.edu" target="_blank">israt.jahan@uconn.edu</a></font><br></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div>
</blockquote></div></div>
</div>-- <br><div dir="ltr"><div dir="ltr"><div><font size="2" face="arial, sans-serif" color="#000000">Regards,</font></div><div><font size="2" face="arial, sans-serif" color="#000000">Israt Jahan</font></div><div><font size="2" face="arial, sans-serif" color="#000000">Graduate Assistant</font></div><div><font size="2" face="arial, sans-serif" color="#000000">Department of Civil and Environmental Engineering</font></div><div><font size="2" face="arial, sans-serif" color="#000000">School of Engineering | University of Connecticut<br></font></div><div><font size="2" face="arial, sans-serif" color="#000000">Email: <a href="mailto:israt.jahan@uconn.edu" target="_blank">israt.jahan@uconn.edu</a></font><br></div></div></div>
_______________________________________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@mailman.ucar.edu" target="_blank">ncl-talk@mailman.ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="https://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">https://mailman.ucar.edu/mailman/listinfo/ncl-talk</a></blockquote></div>