This archive contains answers to questions sent to Unidata support through mid-2025. Note that the archive is no longer being updated. We provide the archive for reference; many of the answers presented here remain technically correct, even if somewhat outdated. For the most up-to-date information on the use of NSF Unidata software and data services, please consult the Software Documentation first.
Steve,
At the U.W. here we have some tiny modifications to dcgrib2 and
gribtonc code that I am no longer sure about. Included below are the
outputs of diff -c between the Unidata distributions (.clean
extension) and our code. Do you happen to know which of these is
correct, i.e. is that '- inrow' needed? Thanks for any help:
*** NAWIPS-5.8.1.clean/unidata/ldmbridge/dcgrib2/qlin.c Wed Feb 26 14:04:41 2003
--- NAWIPS-5.8.1/unidata/ldmbridge/dcgrib2/qlin.c Tue Apr 12 14:38:55 2005
***************
*** 116,122 ****
ni, /* number of output values to compute */
c); /* precomputed interpolation coefficients */
for (i=0; i < ni; i++) {
! double c1 = 1.0 - j*(nrows - 1.0)/(nj - 1.0);
double c2 = 1.0 - c1;
*outp = c1 * *outp + c2*row2[i];
outp++;
--- 116,122 ----
ni, /* number of output values to compute */
c); /* precomputed interpolation coefficients */
for (i=0; i < ni; i++) {
! double c1 = 1.0 - (j*(nrows - 1.0)/(nj - 1.0) - inrow);
double c2 = 1.0 - c1;
*outp = c1 * *outp + c2*row2[i];
outp++;
and for quasi.c:
*** NAWIPS-5.8.1.clean/unidata/ldmbridge/gribtonc_src/quasi.c Thu Jun 18
09:53:45 1998
--- NAWIPS-5.8.1/unidata/ldmbridge/gribtonc_src/quasi.c Tue Apr 12 14:38:56 2005
***************
*** 192,198 ****
ni, /* number of output values to compute */
c); /* precomputed interpolation coefficients */
for (i=0; i < ni; i++) {
! double c1 = 1.0 - j*(nrows - 1.0)/(nj - 1.0);
double c2 = 1.0 - c1;
*outp = c1 * *outp + c2*row2[i];
outp++;
--- 192,198 ----
ni, /* number of output values to compute */
c); /* precomputed interpolation coefficients */
for (i=0; i < ni; i++) {
! double c1 = 1.0 - (j*(nrows - 1.0)/(nj - 1.0) - inrow);
double c2 = 1.0 - c1;
*outp = c1 * *outp + c2*row2[i];
outp++;
--
David Ovens e-mail: address@hidden
Research Meteorologist phone: (206) 685-8108
Dept of Atm. Sciences plan: Real-time MM5 forecasting for the
Box 351640 Pacific Northwest
University of Washington http://www.atmos.washington.edu/mm5rt
Seattle, WA 98195 Weather Graphics and Loops
http://www.atmos.washington.edu/~ovens/loops