#include<stdio.h>
#include "H5IM.h" 

#define WIDTH      57                           /* dataset dimensions */
#define HEIGHT     57 
#define RANK   2

int main (void)
{
    hid_t       file;                           /* file handle */
    herr_t      status;                             
    unsigned char data[WIDTH][HEIGHT];          /* data to write */
    int         i, j, num, val;
    FILE *fp;

    /* Open ASCII file and read in values */
    fp = fopen ("storm110.txt", "r");

    for (i=0; i<WIDTH; i++)
        for (j=0; j<HEIGHT; j++) {
           num = fscanf (fp, "%d ", &val);
           data[i][j] = val;
        }

    /* Create new file */
    file = H5Fcreate ("storm.h5", H5F_ACC_TRUNC, H5P_DEFAULT, H5P_DEFAULT);

    /* Create and write image to file */
    status = H5IMmake_image_8bit (file, "Storm_Image", WIDTH, HEIGHT, 
                                   (const unsigned char *)data);

    /* Close file */
    status = H5Fclose (file);
}