<div dir="ltr">Hi Jared,<div><br></div><div>There is a trouble ticket for this issue and we'll be adding a fix in the upcoming 6.4.0 release.</div><div><br></div><div>Thanks,</div><div><br></div><div>Bill</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Fri, May 13, 2016 at 5:13 PM, Jared Lee <span dir="ltr"><<a href="mailto:jaredlee@ucar.edu" target="_blank">jaredlee@ucar.edu</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">Hi, I just started using the function wrf_user_vert_interp in one of my codes, and I came across an error message:<div><font face="monospace, monospace"><br></font></div><div><div><font face="monospace, monospace">fatal:wrf_vintrp: The pres and field arrays must have the same dimensionality</font></div><div><font face="monospace, monospace">fatal:["Execute.c":8575]:Execute: Error occurred at or near line 2886 in file $NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl</font></div></div><div><br clear="all"><div>Eventually I figured out that it was because the "field" variable that I was passing in didn't have the time dimension. I originally read in the original field (from a file with only one output time) and called the interp with these calls:</div><div><br></div><div><font face="monospace, monospace">TEMP_C = wrf_user_getvar(wrf_file, "tc", 0)</font><br></div><div><font face="monospace, monospace">TEMP_C_interp = wrf_user_vert_interp(wrf_file, TEMP_C, vert_coord, interp_levels, interp_opts)</font><br></div><div><br></div><div>But then when I got my temperature variable with the time dimension added (replacing the 0 with a -1), wrf_user_vert_interp worked successfully:</div><div><br><span style="font-family:monospace,monospace">TEMP_C = wrf_user_getvar(wrf_file, "tc", -1)</span><br></div><div><br></div><div>In the function wrf_user_vert_interp, when the pressure variable is read in, it has the time dimension. Therefore, I think the NCL documentation page for wrf_user_vert_interp, which currently reads:</div><div><br></div><div>"<span style="color:rgb(51,51,51);font-family:verdana,sans-serif;font-size:13.3333px;line-height:16px">The field variable to be interpolated. This must be a 3D array dimensioned bottom_top x south_north x west_east, or a 4D array dimensioned time x bottom_top x south_north x west_east.</span>"</div><div><br></div><div>should be changed to reflect that the function won't work if the field variable is only 3D. The field variable *must* be a 4D array dimensioned time x bottom_top x south_north x west_east. Either that, or the function wrf_vintrp should be updated to allow for "field" variables that don't have a time dimension.</div><div><br></div><div>Thanks,</div><div>Jared</div><span class="HOEnZb"><font color="#888888"><div><br></div>-- <br><div><div dir="ltr"><div><div><div><div><div><font face="courier new, monospace">===============================<br>Jared A. Lee, Ph.D.<br></font></div><font face="courier new, monospace">Project Scientist I<br></font></div><font face="courier new, monospace">Research Applications Laboratory<br></font></div><font face="courier new, monospace">National Center for Atmospheric Research<br>Boulder, Colorado, USA<br><br></font></div><font face="courier new, monospace">Email: <a href="mailto:jaredlee@ucar.edu" target="_blank">jaredlee@ucar.edu</a> (w)<br></font></div><font face="courier new, monospace">Phone: <a href="tel:303.497.8485" value="+13034978485" target="_blank">303.497.8485</a> (w)</font><div><font face="courier new, monospace">Web: <a href="https://staff.ucar.edu/users/jaredlee" target="_blank">https://staff.ucar.edu/users/jaredlee</a><br>===============================</font><br></div></div></div>
</font></span></div></div>
<br>_______________________________________________<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/mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>