<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 class=""><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><div><br></div>-- <br><div class="gmail_signature"><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: 303.497.8485 (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>
</div></div>