1 /* -- translated by f2c (version 20100827).
2 You must link the resulting object file with libf2c:
3 on Microsoft Windows system, link with libf2c.lib;
4 on Linux or Unix systems, link with .../path/to/libf2c.a -lm
5 or, if you install libf2c.a in a standard place, with -lf2c -lm
6 -- in that order, at the end of the command line, as in
7 cc *.o -lf2c -lm
8 Source for libf2c is in /netlib/f2c/libf2c.zip, e.g.,
9
10 http://www.netlib.org/f2c/libf2c.zip
11 */
12
13 #include "libtinyf2c.h"
14
15 /* ----------------------------------------------------------------------- */
16 /* Copyright (C) 2000-2004 GFD Dennou Club. All rights reserved. */
17 /* ----------------------------------------------------------------------- */
vrinc1_(real * rx,real * ry,integer * n,integer * jx,integer * jy,real * rr)18 /* Subroutine */ int vrinc1_(real *rx, real *ry, integer *n, integer *jx,
19 integer *jy, real *rr)
20 {
21 /* System generated locals */
22 integer i__1;
23
24 /* Local variables */
25 static integer j, kx, ky;
26 static real rmiss;
27 extern /* Subroutine */ int glrget_(char *, real *, ftnlen);
28
29 /* Parameter adjustments */
30 --ry;
31 --rx;
32
33 /* Function Body */
34 glrget_("RMISS", &rmiss, (ftnlen)5);
35 kx = 1 - *jx;
36 ky = 1 - *jy;
37 i__1 = *n;
38 for (j = 1; j <= i__1; ++j) {
39 kx += *jx;
40 ky += *jy;
41 if (rx[kx] != rmiss) {
42 ry[ky] = rx[kx] + *rr;
43 } else {
44 ry[ky] = rmiss;
45 }
46 /* L10: */
47 }
48 return 0;
49 } /* vrinc1_ */
50
51