/*
 * Update derived fields.
 */

#include "assign/assign_protein_type.h"
#include "error/check_h5_error.h"
#include <stdio.h>
#include <signal.h>

#define FILE "influenza.h5"

int
main ()
{
  /*
   * Open the HDF5 file.
   */
  hid_t file_id = H5Fopen (FILE, H5F_ACC_RDWR, H5P_DEFAULT);
  if (file_id < 0)
    check_h5_error (file_id, __FILE__, __LINE__);

  /*
   * Assign protein type values to the sequence records.
   */
  assign_protein_type (file_id);

  /*
   * Close the HDF5 file.
   */
  herr_t status = H5Fclose (file_id);
  if (status < 0)
    check_h5_error (status, __FILE__, __LINE__);

  return 0;
}