/******************************************************************
 * FUNCTION: filmat - fill in one row of the matrix of output     *
 *                    variables                                   *
 * FILENAME filmat.c                                              *
 * PROGRAMMER I.Urieli    (FORTRAN)                               *
 *            E. Malroy   (C Translation)                         *
 * DATE:      12/01/94  (FORTRAN), 03/13/96 (C Translation)       *
 * LAST MODIFIED: 08/17/96 by I.Urieli                            *
 * INCLUDE:   "adiabatic.h"                                       *
 * GLOBAL VARIABLES: none                                         *
 * SYMBOLIC CONSTANTS: ROWV, ROWD, COL                            *
 * PROTOTYPES:                                                    *
 *    void filmat(int, double y[], double dy[],                   *
 *                double var[][COL], double dvar[][COL]);         *
 * PRE:   vectors y & dy have values                              *
 * POST:  the j_th column of the var & dvar matrices have values  *
 ******************************************************************/
#include "adiabatic.h"

void filmat(int j,     /* matrix column index */
            double y[ROWV],  /* dependant variable vector */
            double dy[ROWD], /* derivative of dependant variable vector */
            double var[ROWV][COL],      /* variable matrix */
            double dvar[ROWD][COL])     /* derivative matrix */
{
   int i;

   for(i = 0; i<ROWV; i++)
      var[i][j] = y[i];

   for(i = 0; i<ROWD; i++)
      dvar[i][j] = dy[i];
}