new file mode 100644
--- /dev/null
+++ b/progs/mincmakevector/Makefile
@@ -0,0 +1,20 @@
+# --------------------------------------------------------------------
+#
+# MINC Makefile
+#
+
+ROOT = ../..
+include $(ROOT)/Make_machine_specific
+include $(ROOT)/Make_configuration
+
+# Executable names
+PROGS    = mincmakevector
+EXTRA_OBJS = 
+HEADERS  = 
+CDEFINES = -DDEBUG                 # cpp defines
+LDOPT    = $(PROG_LDOPT)
+MANSECT  = 1
+MANPAGES = $(PROGS).$(MANSECT)
+
+include $(ROOT)/progs/Make_progs
+

new file mode 100644
--- /dev/null
+++ b/progs/mincmakevector/mincmakevector.c
@@ -0,0 +1,247 @@
+/* ----------------------------- MNI Header -----------------------------------
+@NAME       : mincmakevector
+@INPUT      : argc, argv - command line arguments
+@OUTPUT     : (none)
+@RETURNS    : status
+@DESCRIPTION: Program to make a vector minc file from a series of scalar 
+              minc files.
+@METHOD     : 
+@GLOBALS    : 
+@CALLS      : 
+@CREATED    : August 11, 1997 (Peter Neelin)
+@MODIFIED   : $Log: mincmakevector.c,v $
+@MODIFIED   : Revision 1.1  1997-08-12 19:03:46  neelin
+@MODIFIED   : Initial revision
+@MODIFIED   :
+---------------------------------------------------------------------------- */
+
+#ifndef lint
+static char rcsid[]="$Header: /private-cvsroot/minc/progs/mincmakevector/mincmakevector.c,v 1.1 1997-08-12 19:03:46 neelin Exp $";
+#endif
+
+#include <stdlib.h>
+#include <stdio.h>
+#include <string.h>
+#include <ctype.h>
+#include <math.h>
+#include <float.h>
+#include <minc.h>
+#include <ParseArgv.h>
+#include <time_stamp.h>
+#include <minc_def.h>
+#include <voxel_loop.h>
+
+#ifndef public
+#  define public
+#endif
+
+#ifndef TRUE
+#  define TRUE 1
+#  define FALSE 0
+#endif
+
+#define DEFAULT_RANGE DBL_MAX
+#define NCOPTS_DEFAULT NC_VERBOSE | NC_FATAL
+#define INVALID_DATA (-DBL_MAX)
+
+/* Structure for program info */
+typedef struct {
+   int dummy;
+} Program_Data;
+
+/* Function prototypes */
+public void do_makevector(void *caller_data, long num_voxels,
+                          int input_num_buffers, int input_vector_length,
+                          double *input_data[],
+                          int output_num_buffers, int output_vector_length,
+                          double *output_data[], Loop_Info *loop_info);
+public long get_vector_length(int mincid);
+
+/* Argument variables */
+int clobber = FALSE;
+int verbose = TRUE;
+nc_type datatype = NC_UNSPECIFIED;
+int is_signed = FALSE;
+double valid_range[2] = {0.0, 0.0};
+int buffer_size = 10 * 1024;
+
+/* Argument table */
+ArgvInfo argTable[] = {
+   {"-clobber", ARGV_CONSTANT, (char *) TRUE, (char *) &clobber,
+       "Overwrite existing file."},
+   {"-noclobber", ARGV_CONSTANT, (char *) FALSE, (char *) &clobber,
+       "Don't overwrite existing file (default)."},
+   {"-verbose", ARGV_CONSTANT, (char *) TRUE, (char *) &verbose,
+       "Print out log messages (default)."},
+   {"-quiet", ARGV_CONSTANT, (char *) FALSE, (char *) &verbose,
+       "Do not print out log messages."},
+   {"-buffer_size", ARGV_INT, (char *) 1, (char *) &buffer_size,
+       "Set the internal buffer size (in kb)."},
+   {"-filetype", ARGV_CONSTANT, (char *) NC_UNSPECIFIED, (char *) &datatype,
+       "Use data type of first file (default)."},
+   {"-byte", ARGV_CONSTANT, (char *) NC_BYTE, (char *) &datatype,
+       "Write out byte data."},
+   {"-short", ARGV_CONSTANT, (char *) NC_SHORT, (char *) &datatype,
+       "Write out short integer data."},
+   {"-long", ARGV_CONSTANT, (char *) NC_LONG, (char *) &datatype,
+       "Write out long integer data."},
+   {"-float", ARGV_CONSTANT, (char *) NC_FLOAT, (char *) &datatype,
+       "Write out single-precision floating-point data."},
+   {"-double", ARGV_CONSTANT, (char *) NC_DOUBLE, (char *) &datatype,
+       "Write out double-precision floating-point data."},
+   {"-signed", ARGV_CONSTANT, (char *) TRUE, (char *) &is_signed,
+       "Write signed integer data."},
+   {"-unsigned", ARGV_CONSTANT, (char *) FALSE, (char *) &is_signed,
+       "Write unsigned integer data (default if type specified)."},
+   {"-valid_range", ARGV_FLOAT, (char *) 2, (char *) valid_range,
+       "Valid range for output data."},
+   {NULL, ARGV_END, NULL, NULL, NULL}
+};
+
+/* Main program */
+
+int main(int argc, char *argv[])
+{
+   char **input_files;
+   char *output_file;
+   char *arg_string;
+   int num_input_files;
+   int inmincid;
+   Loop_Options *loop_options;
+   Program_Data program_data;
+
+   /* Save time stamp and args */
+   arg_string = time_stamp(argc, argv);
+
+   /* Get arguments */
+   if (ParseArgv(&argc, argv, argTable, 0) || (argc < 3)) {
+      (void) fprintf(stderr, 
+                     "\nUsage: %s [options] <in1.mnc> [...] <out.mnc>\n",
+                     argv[0]);
+      (void) fprintf(stderr,
+                     "       %s -help\n\n",
+                     argv[0]);
+      exit(EXIT_FAILURE);
+   }
+   input_files = &argv[1];
+   num_input_files = argc - 2;
+   output_file = argv[argc-1];
+
+   /* Open the first input file and get the vector length */
+   inmincid = miopen(input_files[0], NC_NOWRITE);
+   if (get_vector_length(inmincid) > 1) {
+      (void) fprintf(stderr, "Input file %s is not a scalar file\n",
+                     input_files[0]);
+      exit(EXIT_FAILURE);
+   }
+
+   /* Set up looping options */
+   loop_options = create_loop_options();
+   set_loop_clobber(loop_options, clobber);
+   set_loop_verbose(loop_options, verbose);
+   set_loop_datatype(loop_options, datatype, is_signed, 
+                     valid_range[0], valid_range[1]);
+   set_loop_output_vector_size(loop_options, num_input_files);
+   set_loop_buffer_size(loop_options, (long) buffer_size * 1024);
+   set_loop_first_input_mincid(loop_options, inmincid);
+   set_loop_accumulate(loop_options, TRUE, 0, NULL, NULL);
+
+   /* Do loop */
+   voxel_loop(num_input_files, input_files, 1, &output_file, 
+              arg_string, loop_options,
+              do_makevector, (void *) &program_data);
+
+   exit(EXIT_SUCCESS);
+}
+
+/* ----------------------------- MNI Header -----------------------------------
+@NAME       : do_makevector
+@INPUT      : caller_data - pointer to structure containing lookup info
+              num_voxels - number of voxels to work on
+              input_num_buffers - number of input buffers
+              input_vector_length - length of input vector dimension
+              input_data - vector of pointers to input buffer data
+              output_num_buffers - number of output buffers
+              output_vector_length - length of output vector dimension
+              start - vector specifying start of hyperslab (not used)
+              count - vector specifying count of hyperslab (not used)
+@OUTPUT     : output_data - vector of pointers to output buffer data
+@RETURNS    : (nothing)
+@DESCRIPTION: Routine to copy values from scalar files into a vector.
+@METHOD     : 
+@GLOBALS    : 
+@CALLS      : 
+@CREATED    : August 11, 1997 (Peter Neelin)
+@MODIFIED   : 
+---------------------------------------------------------------------------- */
+public void do_makevector(void *caller_data, long num_voxels,
+                          int input_num_buffers, int input_vector_length,
+                          double *input_data[],
+                          int output_num_buffers, int output_vector_length,
+                          double *output_data[], Loop_Info *loop_info)
+     /* ARGSUSED */
+{
+   Program_Data *program_data;
+   long ivoxel, ovoxel;
+   int current_input_file;
+
+   /* Get pointer to lookup info */
+   program_data = (Program_Data *) caller_data;
+
+   /* Check that values correspond */
+   if ((input_num_buffers != 1) || (output_num_buffers != 1) ||
+       (input_vector_length != 1)) {
+      (void) fprintf(stderr, "Bad internal values.\n");
+      exit(EXIT_FAILURE);
+   }
+
+   /* Figure out which file we are looking at */
+   current_input_file = get_info_current_file(loop_info);
+
+   /* Loop through the voxels */
+   for (ivoxel=0, ovoxel = current_input_file; 
+        ivoxel < num_voxels; 
+        ivoxel++, ovoxel += output_vector_length) {
+
+      output_data[0][ovoxel] = input_data[0][ivoxel];
+
+   }
+
+   return;
+}
+
+/* ----------------------------- MNI Header -----------------------------------
+@NAME       : get_vector_length
+@INPUT      : mincid - minc file id
+@OUTPUT     : (none)
+@RETURNS    : Length of vector dimension or zero if no such dimension.
+@DESCRIPTION: Routine to get the length of the vector dimension in a minc file.
+@METHOD     : 
+@GLOBALS    : 
+@CALLS      : 
+@CREATED    : November 30, 1994 (Peter Neelin)
+@MODIFIED   : 
+---------------------------------------------------------------------------- */
+public long get_vector_length(int mincid)
+{
+   int imgid;
+   int ndims;
+   int dim[MAX_VAR_DIMS];
+   char dimname[MAX_NC_NAME];
+   long vector_length;
+
+   /* Get image variable id */
+   imgid = ncvarid(mincid, MIimage);
+
+   /* Get the image dimension info */
+   (void) ncvarinq(mincid, imgid, NULL, NULL, &ndims, dim, NULL);
+
+   /* Check for vector dimension */
+   (void) ncdiminq(mincid, dim[ndims-1], dimname, &vector_length);
+   if ((strcmp(dimname, MIvector_dimension) != 0) || (ndims <= 2)) {
+      vector_length = 0;
+   }
+
+   return vector_length;
+}
+

new file mode 100644
--- /dev/null
+++ b/progs/mincmakevector/mincmakevector.man1
@@ -0,0 +1,111 @@
+.\" Copyright 1997 Peter Neelin, McConnell Brain Imaging Centre,
+.\" Montreal Neurological Institute, McGill University.
+.\" Permission to use, copy, modify, and distribute this
+.\" software and its documentation for any purpose and without
+.\" fee is hereby granted, provided that the above copyright
+.\" notice appear in all copies.  The author and McGill University
+.\" make no representations about the suitability of this
+.\" software for any purpose.  It is provided "as is" without
+.\" express or implied warranty.
+.\"
+.\" $Header: /private-cvsroot/minc/progs/mincmakevector/mincmakevector.man1,v 1.1 1997-08-12 19:03:46 neelin Exp $
+.\"
+.TH MINCMAKEVECTOR 1
+
+.SH NAME
+mincmakevector - convert a list of scalar minc files into one vector file
+
+.SH SYNOPSIS
+.B mincmakevector
+[<options>] <in>.mnc [...] <out>.mnc
+
+.SH DESCRIPTION 
+.I Mincmakevector
+converts a list of scalar minc files into one vector minc file. A vector 
+minc file is one that contains the dimension 
+.I vector_dimension
+as the fastest varying dimension of the image data and represents vector
+data at each voxel such as RGB images or gradient volumes. A scalar
+minc file does not contain the vector dimension and represents
+grayscale or intensity data.
+
+.SH OPTIONS
+Note that options can be specified in abbreviated form (as long as
+they are unique) and can be given anywhere on the command line.
+
+.SH General options
+.P
+.I -clobber:
+Overwrite an existing file.
+.P
+.I -noclobber:
+Don't overwrite an existing file (default).
+.P
+.I -no_clobber:
+Synonym for -noclobber.
+.P
+.I -verbose:
+Print out progress information for each chunk of data copied (default).
+.P
+.I -quiet:
+Do not print out progress information.
+.P
+.P
+.I -buffer_size
+<size>:
+Specify the maximum size of the internal buffers (in kbytes). Default
+is 10 MB.
+.P
+.I -filetype:
+Create an output file with the same type as the first input file
+(default).
+.P
+.I -byte:
+Create an output file of type byte.
+.P
+.I -short:
+Create an output file of type short.
+.P
+.I -long:
+Create an output file of type long.
+.P
+.I -float:
+Create an output file of type float.
+.P
+.I -double:
+Create an output file of type double.
+.P
+.I -signed:
+Create an output file with data stored in a signed type. This only
+has an effect if the one of the -byte, -short or -long options is
+specified.
+.P
+.I -unsigned:
+Create an output file with data stored in an unsigned type. This only
+has an effect if the one of the -byte, -short or -long options is
+specified.
+.P
+.I -valid_range
+<min> <max>:
+Create an output file with integer data stored in the specified
+valid range. This only has an effect if the one of the -byte, -short
+or -long options is specified.
+
+.SH Generic options for all commands:
+.P
+.I -help:
+Print summary of command-line options and abort.
+
+.SH EXAMPLES
+
+To convert files containing red, green and blue colour components into
+an RGB file:
+
+   mincmakevector red.mnc green.mnc blue.mnc rgb.mnc
+
+.SH AUTHOR
+Peter Neelin
+
+.SH COPYRIGHTS
+.ps 18
+Copyright 1997 by Peter Neelin