diff options
Diffstat (limited to 'pkg/xtools/icfit/iclistr.x')
| -rw-r--r-- | pkg/xtools/icfit/iclistr.x | 45 |
1 files changed, 45 insertions, 0 deletions
diff --git a/pkg/xtools/icfit/iclistr.x b/pkg/xtools/icfit/iclistr.x new file mode 100644 index 00000000..4e2b2c14 --- /dev/null +++ b/pkg/xtools/icfit/iclistr.x @@ -0,0 +1,45 @@ +# Copyright(c) 1986 Association of Universities for Research in Astronomy Inc. + +include "icfit.h" +include "names.h" + +# IC_LIST -- List X, Y, FIT, W. + +procedure ic_listr (ic, cv, x, y, wts, npts, file) + +pointer ic # ICFIT pointer +pointer cv # Curfit pointer +real x[ARB] # Ordinates +real y[ARB] # Abscissas +real wts[ARB] # Weights +int npts # Number of data points +char file[ARB] # Output file + +int i, fd, open() +real rcveval() +errchk open() + +begin + # Open the output file. + fd = open (file, APPEND, TEXT_FILE) + + if (npts == IC_NFIT(ic)) { + do i = 1, npts { + call fprintf (fd, "%8g %8g %8g %8g\n") + call pargr (x[i]) + call pargr (y[i]) + call pargr (rcveval (cv, x[i])) + call pargr (wts[i]) + } + } else { + do i = 1, IC_NFIT(ic) { + call fprintf (fd, "%8g %8g %8g %8g\n") + call pargr (Memr[IC_XFIT(ic)+i-1]) + call pargr (Memr[IC_YFIT(ic)+i-1]) + call pargr (rcveval (cv, Memr[IC_XFIT(ic)+i-1])) + call pargr (Memr[IC_WTSFIT(ic)+i-1]) + } + } + + call close (fd) +end |