/************************************************************

  This example shows how to create intermediate groups with
  a single call to H5Gcreate.

  This file is intended for use with HDF5 Library version 1.8

 ************************************************************/

#include "hdf5.h"
//#include "/usr/local/hdf5/include/H5Dpublic.h"
#include <stdio.h>
#include <stdlib.h>

//#define FILE            "myh5ex_g_pwa_interface.h5"
#define DATASET         "title"
#define DATASET2	"primeID"
#define	DATASET3	"/elementData/elementNames"
#define	DATASET4	"/elementData/elementPrimeID"
#define DATASET5        "/phaseData/eqState"
#define DATASET6        "/phaseData/phase"
#define DATASET7        "/reactionData/reactionEquations" 
#define DATASET8        "/reactionData/reactionPrimeID" 
#define DATASET9        "/reactionData/reactionRateType"
#define DATASET10        "/reactionData/arrheniusData/arrheniusPrimeID"
#define DATASET11        "/reactionData/falloffData/falloffType"
#define DATASET12        "/reactionData/falloffData/falloffNames"
#define DATASET13        "/reactionData/rateCoefLinks/r00009454"
#define DATASET14        "/reactionData/rateCoefLinks/r00009484"  
#define DATASET15        "/reactionData/rateCoefLinks/r00009509"             
#define DATASET16        "/reactionData/rateCoefLinks/r00011821"
#define DATASET17        "/reactionData/rateCoefLinks/r00011842"            
#define DATASET18        "/reactionData/rateCoefLinks/r00011866"                                                                                                                
#define DATASET19        "/reactionData/rateCoefLinks/r00012023"                                                                                                                
#define DATASET20        "/reactionData/rateCoefLinks/r00012138"                                                                                                                
#define DATASET21        "/reactionData/rateCoefLinks/r00012186"                                                                                                                
#define DATASET22        "/reactionData/rateCoefLinks/r00012255"                                                                                                                
#define DATASET23        "/reactionData/rateCoefLinks/r00013690"
#define DATASET24        "/reactionData/rateCoefLinks/r00013692"                                                                                                                 
#define DATASET25        "/reactionData/rateCoefLinks/r00013701"                                                                                                                 
#define DATASET26        "/reactionData/rateCoefLinks/r00013764"                                                                                                                 
#define DATASET27        "/reactionData/rateCoefLinks/r00013767"                                                                                                                 
#define DATASET28        "/reactionData/rateCoefLinks/r00013781"                                                                                                                 
#define DATASET29        "/reactionData/rateCoefLinks/r00013804"                                                                                                                 
#define DATASET30        "/reactionData/rateCoefLinks/r00013805"                                                                                                                 
#define DATASET31        "/reactionData/rateCoefLinks/r00013869"                                                                                                                 
#define DATASET32        "/reactionData/rateCoefLinks/r00015825"
#define DATASET33        "/speciesData/speciesNames"
#define DATASET34        "/speciesData/speciesPrimeID"                                                                                                                 
#define DATASET35        "/speciesData/thermoPolynomials/thermoPrimeID"
#define DATASET36        "/speciesData/thermoPolynomials/thermoType"

#define DATASET37	"/reactionData/arrheniusData/arrheniusParam"


#define NX2		43
#define NY2		3
#define RANK2		3

#define DIM0            1
#define DIM1		1
#define DIM2		4
#define DIM3            4
#define DIM4            1
#define DIM5            1
#define DIM6            20
#define DIM7            20
#define DIM8            20
#define DIM9            43
#define DIM10           1
#define DIM11           4
#define DIM12           1
#define DIM13           1
#define DIM14           1
#define DIM15           1
#define DIM16           1
#define DIM17           1
#define DIM18           1
#define DIM19           1
#define DIM20           1
#define DIM21           1
#define DIM22           1                                                                                                                                                                                 
#define DIM23           1
#define DIM24           1
#define DIM25           1
#define DIM26           1
#define DIM27           1
#define DIM28           1
#define DIM29           1
#define DIM30           1
#define DIM31           1
#define DIM32           10
#define DIM33           10
#define DIM34           10
#define DIM35           1
#define DIM36		43

/*
 * Operator function to be called by H5Ovisit.
 */
herr_t op_func (hid_t loc_id, const char *name, const H5O_info_t *info,
            void *operator_data);
//int myDataSets(void);

int
main(void)
{
    hid_t       file, group1, group2, group3, group4, group5,group6, group7, group8, gcpl1, gcpl2, gcpl3, gcpl4, gcpl5, gcpl6, gcpl7, gcpl8;      /* Handles */
//    hid_t	file2, dset1, dset2, dspace;
    hid_t	filetype, filetype2, filetype3, filetype4, memtype, memtype2, memtype3, memtype4, space, space2,  space3, space4, dset, dset2, dset3, dset4;
    hid_t       filetype5, filetype6, memtype5, memtype6, space5, space6, dset5, dset6, filetype7, filetype8, filetype9, filetype10, memtype7, memtype8; 
    hid_t	memtype9, memtype10, space7, space8, space9, space10, dset7, dset8, dset9, dset10;
    hid_t       filetype11, filetype12, filetype13, filetype14, memtype11, memtype12, memtype13, memtype14, space11, space12,  space13, space14, dset11, dset12, dset13, dset14;
    hid_t       filetype15, filetype16, filetype17, filetype18, memtype15, memtype16, memtype17, memtype18, space15, space16, space17, space18, dset15, dset16, dset17, dset18;
    hid_t       filetype19, filetype20, filetype21, filetype22, memtype19, memtype20, memtype21, memtype22, space19, space20, space21, space22, dset19, dset20, dset21, dset22;
    hid_t       filetype23, filetype24, filetype25, filetype26, memtype23, memtype24, memtype25, memtype26, space23, space24, space25, space26, dset23, dset24, dset25, dset26;
    hid_t       filetype27, filetype28, filetype29, filetype30, memtype27, memtype28, memtype29, memtype30, space27, space28, space29, space30, dset27, dset28, dset29, dset30; 
    hid_t       filetype31, filetype32,  filetype33, filetype34, memtype31, memtype32, memtype33, memtype34, space31, space32, space33, space34, dset31, dset32, dset33, dset34; 
    hid_t       filetype35, filetype36, memtype35, memtype36, space35, space36, dset35, dset36;

    hid_t	filetype37, memtype37, space37, dset37;  
    hsize_t	dimsf2[2];
//    float	wdata2[NX2][NY2];
//    float	rdata2[NX2][NY2];
//    float	myval2;
//    int 	myi, myy, mynum;
    
//    FILE 	*fp2;
    

    //hsize_t	dims[2];
    hsize_t     dims[1] = {DIM0};
    hsize_t	dims2[1] = {DIM1};
    hsize_t	dims3[4] = {DIM2};
    hsize_t     dims4[4] = {DIM3};
    hsize_t     dims5[1] = {DIM4};
    hsize_t     dims6[1] = {DIM5};
    hsize_t     dims7[20] = {DIM6};
    hsize_t     dims8[20] = {DIM7};
    hsize_t     dims9[20] = {DIM8};
    hsize_t     dims10[43] = {DIM9};
    hsize_t     dims11[1] = {DIM10};
    hsize_t     dims12[4] = {DIM11};
    hsize_t     dims13[1] = {DIM12};
    hsize_t     dims14[1] = {DIM13};
    hsize_t     dims15[1] = {DIM14};
    hsize_t     dims16[1] = {DIM15};
    hsize_t     dims17[1] = {DIM16};
    hsize_t     dims18[1] = {DIM17};
    hsize_t     dims19[1] = {DIM18};
    hsize_t     dims20[1] = {DIM19};
    hsize_t     dims21[1] = {DIM20};
    hsize_t     dims22[1] = {DIM21};
    hsize_t     dims23[1] = {DIM22};
    hsize_t     dims24[1] = {DIM23};
    hsize_t     dims25[1] = {DIM24};                                                                                                                                                                      
    hsize_t     dims26[1] = {DIM25};                                                                                                                                                                      
    hsize_t     dims27[1] = {DIM26};                                                                                                                                                                      
    hsize_t     dims28[1] = {DIM27};                                                                                                                                                                      
    hsize_t     dims29[1] = {DIM28};
    hsize_t     dims30[1] = {DIM29};
    hsize_t     dims31[1] = {DIM30};
    hsize_t     dims32[1] = {DIM31};
    hsize_t     dims33[10] = {DIM32};
    hsize_t     dims34[10] = {DIM33};
    hsize_t     dims35[10] = {DIM34};
    hsize_t     dims36[1] = {DIM35};
   // hsize_t	dims37[43][3] = {DIM36,NY2};

    float       wdataf2[NX2][NY2];
    float       rdataf2[NX2][NY2];
    float       myval2;
    int         myi, myj, mynum;

//    FILE *fp2;





    char        *wdata[DIM0] = {"PrIMe H2-Mech 1.0"},
                                            /* Write buffer */
                **rdata;                    /* Read buffer */
    int 	ndims, i;

    char        *wdata2[DIM1] = {"m00000004"},                                                                                                                                                   
                                            /* Write buffer */
                **rdata2;                    /* Read buffer */
    int         ndims2, i2;

    char        *wdata3[DIM2] = {"Ar", "H", "N", "O"},
                                            /* Write buffer */
                **rdata3;                    /* Read buffer */
    int         ndims3, i3;

    char        *wdata4[DIM3] = {"e00000018", "e00000001", "e00000007", "e00000008"},
                                            /* Write buffer */
                **rdata4;                    /* Read buffer */
    int         ndims4, i4;

    char        *wdata5[DIM4] = {"ideal gas"},
                                            /* Write buffer */
                **rdata5;                    /* Read buffer */
    int         ndims5, i5;

    char        *wdata6[DIM5] = {"gas"},
                                            /* Write buffer */
                **rdata6;                    /* Read buffer */
    int         ndims6, i6;

   char        *wdata7[DIM6] = {"O + O <=> O2", "H + O <=> OH", "H2 + O <=> OH + H", "HO2 + O <=> OH + O2", "O + H2O2 <=> OH + HO2", "O2 + H <=> HO2", "O2 + H <=> OH + O", "2 H <=> H2", "OH + H <=> H2O", "HO2 + H <=> H2O + O", "HO2 + H <=> H2 + O2", "HO2 + H <=> 2 OH", "H2O2 + H <=> H2 + HO2", "H2O2 + H <=> OH + H2O", "H2 + OH <=> H2O + H", "2 OH <=> H2O2", "2 OH <=> H2O + O", "HO2 + OH <=> H2O + O2", "OH + H2O2 <=> HO2 + H2O", "2 HO2 <=> H2O2 + O2"},
                                            /* Write buffer */
                **rdata7;                    /* Read buffer */
    int         ndims7, i7;

    char        *wdata8[DIM7] = {"r00013869", "r00012023", "r00012186", "r00012138", "r00015825", "r00011821", "r00012255", "r00009454", "r00013764", "r00013804", "r00009509", "r00011866", "r00009484", "r00013701", "r00013781", "r00013690", "r00013805", "r00013767", "r00011842", "r00013692"},
                                            /* Write buffer */
                **rdata8;                    /* Read buffer */
    int         ndims8, i8;

    char        *wdata9[DIM8] = {"third body", "third body", "mass action", "mass action", "mass action", "third body", "mass action", "third body", "third body", "mass action", "mass action", "mass action", "mass action", "mass action", "mass action", "unimolecular", "mass action", "sum", "sum", "sum"},
                                            /* Write buffer */
                **rdata9;                    /* Read buffer */
    int         ndims9, i9; 

    char        *wdata10[DIM9] = {"rk00000016", "rk00000017", "rk00000022", "rk00000015", "rk00000003", "rk00000004", "rk00000009", "rk00000002", "rk00000005", "rk00000002", "rk00000003", "rk00000006", "rk00000007", "rk00000009", "rk00000011", "rk00000015", "rk00000042", "rk00000019", "rk00000020", "rk00000023", "rk00000024", "rk00000010", "rk00000011", "rk00000014", "rk00000009", "rk00000002", "rk00000007", "rk00000020", "rk00000008", "rk00000020", "rk00000015", "rk00000030", "rk00000023", "rk00000024", "rk00000029", "rk00000022", "rk00000012", "rk00000002", "rk00000003", "rk00000003", "rk00000004", "rk00000002", "rk00000003"},
                                            /* Write buffer */
                **rdata10;                    /* Read buffer */
    int         ndims10, i10;

    char        *wdata11[DIM10] = {"troe"},
                                            /* Write buffer */
                **rdata11;                    /* Read buffer */
    int         ndims11, i11;

    char        *wdata12[DIM11] = {"a", "T***", "T*", "T**"},
                                            /* Write buffer */
                **rdata12;                    /* Read buffer */
    int         ndims12, i12;




    char        *wdata13[DIM12] = {"rk00000028"},
                                            /* Write buffer */
                **rdata13;                    /* Read buffer */
    int         ndims13, i13;



    char        *wdata14[DIM13] = {"rk00000008"},
                                            /* Write buffer */
                **rdata14;                    /* Read buffer */
    int         ndims14, i14;

    char        *wdata15[DIM14] = {"rk00000007"},
                                            /* Write buffer */
                **rdata15;                    /* Read buffer */
    int         ndims15, i15;

    char        *wdata16[DIM15] = {"rk00000025"},
                                            /* Write buffer */
                **rdata16;                    /* Read buffer */
    int         ndims16, i16;

    char        *wdata17[DIM16] = {"rk00000043"},
                                            /* Write buffer */
                **rdata17;                    /* Read buffer */
    int         ndims17, i17;



    char        *wdata18[DIM17] = {"rk00000020"},                                                                                                                                                         
                                            /* Write buffer */
                **rdata18;                    /* Read buffer */                                                                                                                                           
    int         ndims18, i18;                          

    char        *wdata19[DIM18] = {"rk00000011"},                                                                                                                                                         
                                            /* Write buffer */
                **rdata19;                    /* Read buffer */                                                                                                                                           
    int         ndims19, i19;                          

    char        *wdata20[DIM19] = {"rk00000002"},                                                                                                                                                         
                                            /* Write buffer */
                **rdata20;                    /* Read buffer */                                                                                                                                           
    int         ndims20, i20;                          

    char        *wdata21[DIM20] = {"rk00000005"},                                                                                                                                                         
                                            /* Write buffer */
                **rdata21;                    /* Read buffer */                                                                                                                                           
    int         ndims21, i21;                          

    char        *wdata22[DIM21] = {"rk00000042"},                                                                                                                                                         
                                            /* Write buffer */
                **rdata22;                    /* Read buffer */                                                                                                                                           
    int         ndims22, i22;                          

    char        *wdata23[DIM22] = {"rk00000042"},                                                                                                                                                         
                                            /* Write buffer */
                **rdata23;                    /* Read buffer */                                                                                                                                           
    int         ndims23, i23;                          

    char        *wdata24[DIM23] = {"rk00000048"},                                                                                                                                                         
                                            /* Write buffer */
                **rdata24;                    /* Read buffer */                                                                                                                                           
    int         ndims24, i24;                          


    char        *wdata25[DIM24] = {"rk00000020"},
                                            /* Write buffer */
                **rdata25;                    /* Read buffer */
    int         ndims25, i25;

    char        *wdata26[DIM25] = {"rk00000017"},
                                            /* Write buffer */
                **rdata26;                    /* Read buffer */
    int         ndims26, i26;

    char        *wdata27[DIM26] = {"rk00000050"},
                                            /* Write buffer */
                **rdata27;                    /* Read buffer */
    int         ndims27, i27;

    char        *wdata28[DIM27] = {"rk00000015"},
                                            /* Write buffer */
                **rdata28;                    /* Read buffer */
    int         ndims28, i28;

    char        *wdata29[DIM28] = {"rk00000002"},
                                            /* Write buffer */
                **rdata29;                    /* Read buffer */
    int         ndims29, i29;

    char        *wdata30[DIM29] = {"rk00000012"},
                                            /* Write buffer */
                **rdata30;                    /* Read buffer */
    int         ndims30, i30;

    char        *wdata31[DIM30] = {"rk00000026"},
                                            /* Write buffer */
                **rdata31;                    /* Read buffer */
    int         ndims31, i31;


    char        *wdata32[DIM31] = {"rk00000003"},
                                            /* Write buffer */
                **rdata32;                  /* Read buffer */                   
    int         ndims32, i32;

    char        *wdata33[DIM32] = {"Ar", "N2", "H", "H2", "O", "O2", "OH", "H2O", "HO2", "H2O2"},
                                            /* Write buffer */
                **rdata33;                   /* Read buffer */
    int         ndims33, i33;

    char        *wdata34[DIM33] = {"s00000049", "s00010231", "s00009800", "s00009809", "s00010285", "s00010295", "s00010102", "s00009881", "s00010103", "s00009882"},
                                            /* Write buffer */
                **rdata34;                   /* Read buffer */
    int         ndims34, i34;


    char        *wdata35[DIM34] = {"thp00000001", "thp00000001", "thp00000001", "thp00000001", "thp00000001", "thp00000001", "thp00000001", "thp00000001", "thp00000001", "thp00000001"},
                                            /* Write buffer */
                **rdata35;                   /* Read buffer */
    int         ndims35, i35;


    char        *wdata36[DIM35] = {"NASA7"},
                                            /* Write buffer */
                **rdata36;                   /* Read buffer */
    int         ndims36, i36;

    float        *wdata37[NX2][NY2],
                                            /* Write buffer */
                **rdata37;                   /* Read buffer */
    int         ndims37, i37;



    //FILE *fp2;


    herr_t      status;

    FILE *fp;


    fp = fopen ("arrheniusParam.csv", "r+");

    for (myi=0; myi<NX2; myi++)
        for (myj=0; myj<NY2; myj++)
        {
           mynum = fscanf (fp, "%f", &myval2);
           wdataf2[myi][myj] = myval2;
        }

/*
    printf ("\n");
    for (myi=0; myi<NX2; myi++)
    {
        for (myj=0; myj<NY2; myj++)
        {
         printf (" %f", rdataf2[myi][myj]);
        }
        printf ("\n");
   }

*/

//     myDataSets();    

    /*
     * Create a new file using the default properties.
     */
    //pid  = H5Pcreate (H5P_FILE_ACCESS);
    file = H5Fcreate ("myh5ex_g_pwa_interface.h5", H5F_ACC_TRUNC, H5P_DEFAULT, H5P_DEFAULT);
     




    /*
     * Describe the size of the array and create the data space for fixed
     * size dataset. 
     */
    dimsf2[0] = NX2;
    dimsf2[1] = NY2;
    space37 = H5Screate_simple(2, dimsf2, NULL);

//    filetype37 = H5Tcopy (H5T_FORTRAN_S1);
//    status = H5Tset_size (filetype37, H5T_VARIABLE);

    //filetype37 = H5Tcopy (H5T_NATIVE_FLOAT);                                                                                                                                                                
    //status = H5Tset_size (filetype37, H5T_VARIABLE); 

    dset37 = H5Dcreate(file, DATASET37, H5T_IEEE_F32BE, space37, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT);

    //dset37 = H5Dcreate(file, DATASET37, filetype37, space37, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT);
//     dset5 = H5Dcreate (file, DATASET5, filetype5, space5, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT);

    /*
     * Write the data to the dataset using default transfer properties.
     */
   

   status = H5Dwrite(dset37, H5T_NATIVE_FLOAT, H5S_ALL, H5S_ALL, H5P_DEFAULT, wdataf2);

    status = H5Dread(dset37, H5T_NATIVE_FLOAT, H5S_ALL, H5S_ALL, H5P_DEFAULT, rdataf2); 



    printf ("\n");
    for (myi=0; myi<NX2; myi++)
    {
        for (myj=0; myj<NY2; myj++)
        {
         printf (" %s[%d][%d]: %f", DATASET37, myi, myj, rdataf2[myi][myj]);
        }
        printf ("\n");
   }

      /*
     * Close/release resources.
     */
    H5Dclose(dset37);
    H5Sclose(space37);
    //fclose(fp);
    
 


 //    myDataSets();
    
     /*
     * Create file and memory datatypes.  For this example we will save
     * the strings as FORTRAN strings.
     */
    filetype = H5Tcopy (H5T_FORTRAN_S1);
    status = H5Tset_size (filetype, H5T_VARIABLE);
    filetype2 = H5Tcopy (H5T_FORTRAN_S1);
    status = H5Tset_size (filetype2, H5T_VARIABLE);
    filetype3 = H5Tcopy (H5T_FORTRAN_S1);
    status = H5Tset_size (filetype3, H5T_VARIABLE);
    filetype4 = H5Tcopy (H5T_FORTRAN_S1);
    status = H5Tset_size (filetype4, H5T_VARIABLE);
    filetype5 = H5Tcopy (H5T_FORTRAN_S1);
    status = H5Tset_size (filetype5, H5T_VARIABLE);
    filetype6 = H5Tcopy (H5T_FORTRAN_S1);
    status = H5Tset_size (filetype6, H5T_VARIABLE);
    filetype7 = H5Tcopy (H5T_FORTRAN_S1);
    status = H5Tset_size (filetype7, H5T_VARIABLE);
    filetype8 = H5Tcopy (H5T_FORTRAN_S1);
    status = H5Tset_size (filetype8, H5T_VARIABLE);
    filetype9 = H5Tcopy (H5T_FORTRAN_S1);
    status = H5Tset_size (filetype9, H5T_VARIABLE);
    filetype10 = H5Tcopy (H5T_FORTRAN_S1);
    status = H5Tset_size (filetype10, H5T_VARIABLE);
    filetype11 = H5Tcopy (H5T_FORTRAN_S1);
    status = H5Tset_size (filetype11, H5T_VARIABLE);
    filetype12 = H5Tcopy (H5T_FORTRAN_S1);
    status = H5Tset_size (filetype12, H5T_VARIABLE);
    filetype13 = H5Tcopy (H5T_FORTRAN_S1);
    status = H5Tset_size (filetype13, H5T_VARIABLE);
    filetype14 = H5Tcopy (H5T_FORTRAN_S1);
    status = H5Tset_size (filetype14, H5T_VARIABLE);
    filetype15 = H5Tcopy (H5T_FORTRAN_S1);
    status = H5Tset_size (filetype15, H5T_VARIABLE);
    filetype16 = H5Tcopy (H5T_FORTRAN_S1);
    status = H5Tset_size (filetype16, H5T_VARIABLE);    
    filetype17 = H5Tcopy (H5T_FORTRAN_S1);
    status = H5Tset_size (filetype17, H5T_VARIABLE);
    filetype18 = H5Tcopy (H5T_FORTRAN_S1);
    status = H5Tset_size (filetype18, H5T_VARIABLE);
    filetype19 = H5Tcopy (H5T_FORTRAN_S1);
    status = H5Tset_size (filetype19, H5T_VARIABLE);
    filetype20 = H5Tcopy (H5T_FORTRAN_S1);
    status = H5Tset_size (filetype20, H5T_VARIABLE);
    filetype21 = H5Tcopy (H5T_FORTRAN_S1);
    status = H5Tset_size (filetype21, H5T_VARIABLE);
    filetype22 = H5Tcopy (H5T_FORTRAN_S1);
    status = H5Tset_size (filetype22, H5T_VARIABLE);
    filetype23 = H5Tcopy (H5T_FORTRAN_S1);
    status = H5Tset_size (filetype23, H5T_VARIABLE);
    filetype24 = H5Tcopy (H5T_FORTRAN_S1);
    status = H5Tset_size (filetype24, H5T_VARIABLE); 
    filetype25 = H5Tcopy (H5T_FORTRAN_S1);
    status = H5Tset_size (filetype25, H5T_VARIABLE);
    filetype26 = H5Tcopy (H5T_FORTRAN_S1);
    status = H5Tset_size (filetype26, H5T_VARIABLE);
    filetype27 = H5Tcopy (H5T_FORTRAN_S1);
    status = H5Tset_size (filetype27, H5T_VARIABLE);
    filetype28 = H5Tcopy (H5T_FORTRAN_S1);
    status = H5Tset_size (filetype28, H5T_VARIABLE);
    filetype29 = H5Tcopy (H5T_FORTRAN_S1);
    status = H5Tset_size (filetype29, H5T_VARIABLE);
    filetype30 = H5Tcopy (H5T_FORTRAN_S1);
    status = H5Tset_size (filetype30, H5T_VARIABLE);
    filetype31 = H5Tcopy (H5T_FORTRAN_S1);
    status = H5Tset_size (filetype31, H5T_VARIABLE);
    filetype32 = H5Tcopy (H5T_FORTRAN_S1);
    status = H5Tset_size (filetype32, H5T_VARIABLE);
    filetype33 = H5Tcopy (H5T_FORTRAN_S1);
    status = H5Tset_size (filetype33, H5T_VARIABLE);
    filetype34 = H5Tcopy (H5T_FORTRAN_S1);                                                                                                                                                                
    status = H5Tset_size (filetype34, H5T_VARIABLE);                                                                                                                                                      
    filetype35 = H5Tcopy (H5T_FORTRAN_S1);                                                                                                                                                                
    status = H5Tset_size (filetype35, H5T_VARIABLE);                                                                                                                                                      
    filetype36 = H5Tcopy (H5T_FORTRAN_S1);                                                                                                                                                                
    status = H5Tset_size (filetype36, H5T_VARIABLE);  

    memtype = H5Tcopy (H5T_C_S1);
    status = H5Tset_size (memtype, H5T_VARIABLE);
    memtype2 = H5Tcopy (H5T_C_S1);
    status = H5Tset_size (memtype2, H5T_VARIABLE);
    memtype3 = H5Tcopy (H5T_C_S1);
    status = H5Tset_size (memtype3, H5T_VARIABLE);
    memtype4 = H5Tcopy (H5T_C_S1);
    status = H5Tset_size (memtype4, H5T_VARIABLE);
    memtype5 = H5Tcopy (H5T_C_S1);
    status = H5Tset_size (memtype5, H5T_VARIABLE);
    memtype6 = H5Tcopy (H5T_C_S1);
    status = H5Tset_size (memtype6, H5T_VARIABLE);
    memtype7 = H5Tcopy (H5T_C_S1);
    status = H5Tset_size (memtype7, H5T_VARIABLE);
    memtype8 = H5Tcopy (H5T_C_S1);
    status = H5Tset_size (memtype8, H5T_VARIABLE);
    memtype9 = H5Tcopy (H5T_C_S1);
    status = H5Tset_size (memtype9, H5T_VARIABLE);
    memtype10 = H5Tcopy (H5T_C_S1);
    status = H5Tset_size (memtype10, H5T_VARIABLE);
    memtype11 = H5Tcopy (H5T_C_S1);
    status = H5Tset_size (memtype11, H5T_VARIABLE);
    memtype12 = H5Tcopy (H5T_C_S1);
    status = H5Tset_size (memtype12, H5T_VARIABLE);
    memtype13 = H5Tcopy (H5T_C_S1);
    status = H5Tset_size (memtype13, H5T_VARIABLE);
    memtype14 = H5Tcopy (H5T_C_S1);
    status = H5Tset_size (memtype14, H5T_VARIABLE);
    memtype15 = H5Tcopy (H5T_C_S1);
    status = H5Tset_size (memtype15, H5T_VARIABLE);
    memtype16 = H5Tcopy (H5T_C_S1);
    status = H5Tset_size (memtype16, H5T_VARIABLE);
    memtype17 = H5Tcopy (H5T_C_S1);
    status = H5Tset_size (memtype17, H5T_VARIABLE);
    memtype18 = H5Tcopy (H5T_C_S1);
    status = H5Tset_size (memtype18, H5T_VARIABLE);
    memtype19 = H5Tcopy (H5T_C_S1);
    status = H5Tset_size (memtype19, H5T_VARIABLE);
    memtype20 = H5Tcopy (H5T_C_S1);
    status = H5Tset_size (memtype20, H5T_VARIABLE);
    memtype21 = H5Tcopy (H5T_C_S1);
    status = H5Tset_size (memtype21, H5T_VARIABLE);
    memtype22 = H5Tcopy (H5T_C_S1);
    status = H5Tset_size (memtype22, H5T_VARIABLE);
    memtype23 = H5Tcopy (H5T_C_S1);
    status = H5Tset_size (memtype23, H5T_VARIABLE);
    memtype24 = H5Tcopy (H5T_C_S1);
    status = H5Tset_size (memtype24, H5T_VARIABLE);
    memtype25 = H5Tcopy (H5T_C_S1);
    status = H5Tset_size (memtype25, H5T_VARIABLE);
    memtype26 = H5Tcopy (H5T_C_S1);
    status = H5Tset_size (memtype26, H5T_VARIABLE);
    memtype27 = H5Tcopy (H5T_C_S1);
    status = H5Tset_size (memtype27, H5T_VARIABLE);
    memtype28 = H5Tcopy (H5T_C_S1);
    status = H5Tset_size (memtype28, H5T_VARIABLE);
    memtype29 = H5Tcopy (H5T_C_S1);
    status = H5Tset_size (memtype29, H5T_VARIABLE);
    memtype30 = H5Tcopy (H5T_C_S1);
    status = H5Tset_size (memtype30, H5T_VARIABLE);
    memtype31 = H5Tcopy (H5T_C_S1);
    status = H5Tset_size (memtype31, H5T_VARIABLE);
    memtype32 = H5Tcopy (H5T_C_S1);
    status = H5Tset_size (memtype32, H5T_VARIABLE);
    memtype33 = H5Tcopy (H5T_C_S1);
    status = H5Tset_size (memtype33, H5T_VARIABLE);
    memtype34 = H5Tcopy (H5T_C_S1);
    status = H5Tset_size (memtype34, H5T_VARIABLE);
    memtype35 = H5Tcopy (H5T_C_S1);
    status = H5Tset_size (memtype35, H5T_VARIABLE);
    memtype36 = H5Tcopy (H5T_C_S1);
    status = H5Tset_size (memtype36, H5T_VARIABLE);

    //memtype37 = H5Tcopy (H5T_C_S1);
    //status = H5Tset_size (memtype37, H5T_VARIABLE);


    /*
     * Create dataspace.  Setting maximum size to NULL sets the maximum
     * size to be the current size.
     */
    space = H5Screate_simple (1, dims, NULL);
    space2 = H5Screate_simple (1, dims2, NULL);
    space3 = H5Screate_simple (1, dims3, NULL);
    space4 = H5Screate_simple (1, dims4, NULL);
    space5 = H5Screate_simple (1, dims5, NULL);
    space6 = H5Screate_simple (1, dims6, NULL);
    space7 = H5Screate_simple (1, dims7, NULL);
    space8 = H5Screate_simple (1, dims8, NULL);
    space9 = H5Screate_simple (1, dims9, NULL);
    space10 = H5Screate_simple (1, dims10, NULL);
    space11 = H5Screate_simple (1, dims11, NULL);
    space12 = H5Screate_simple (1, dims12, NULL);
    space13 = H5Screate_simple (1, dims13, NULL);
    space14 = H5Screate_simple (1, dims14, NULL);
    space15 = H5Screate_simple (1, dims15, NULL);
    space16 = H5Screate_simple (1, dims16, NULL);
    space17 = H5Screate_simple (1, dims17, NULL);
    space18 = H5Screate_simple (1, dims18, NULL);
    space19 = H5Screate_simple (1, dims19, NULL);
    space20 = H5Screate_simple (1, dims20, NULL);
    space21 = H5Screate_simple (1, dims21, NULL);
    space22 = H5Screate_simple (1, dims22, NULL);
    space23 = H5Screate_simple (1, dims23, NULL);
    space24 = H5Screate_simple (1, dims24, NULL);
    space25 = H5Screate_simple (1, dims25, NULL);
    space26 = H5Screate_simple (1, dims26, NULL);
    space27 = H5Screate_simple (1, dims27, NULL);
    space28 = H5Screate_simple (1, dims28, NULL);
    space29 = H5Screate_simple (1, dims29, NULL);
    space30 = H5Screate_simple (1, dims30, NULL);
    space31 = H5Screate_simple (1, dims31, NULL);
    space32 = H5Screate_simple (1, dims32, NULL);
    space33 = H5Screate_simple (1, dims33, NULL);
    space34 = H5Screate_simple (1, dims34, NULL);
    space35 = H5Screate_simple (1, dims35, NULL);
    space36 = H5Screate_simple (1, dims36, NULL);

    //space37 = H5Screate_simple (3, dims37, NULL);

     space37 = H5Screate_simple(RANK2, dimsf2, NULL);

    /*
     * Create the dataset and write the variable-length string data to
     * it.
     */
    //dset = H5Dcreate (file, DATASET, filetype, space, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT);
    //status = H5Dwrite (dset, memtype, H5S_ALL, H5S_ALL, H5P_DEFAULT, wdata);

    /*
     * Create group creation property list and set it to allow creation
     * of intermediate groups.
     */
    gcpl1 = H5Pcreate (H5P_LINK_CREATE);
    status = H5Pset_create_intermediate_group (gcpl1, 1);
    gcpl2 = H5Pcreate (H5P_LINK_CREATE);
    status = H5Pset_create_intermediate_group (gcpl2, 1);
    gcpl3 = H5Pcreate (H5P_LINK_CREATE);
    status = H5Pset_create_intermediate_group (gcpl3, 1);
    gcpl4 = H5Pcreate (H5P_LINK_CREATE);
    status = H5Pset_create_intermediate_group (gcpl4, 1);
    gcpl5 = H5Pcreate (H5P_LINK_CREATE);
    status = H5Pset_create_intermediate_group (gcpl5, 1);
    gcpl6 = H5Pcreate (H5P_LINK_CREATE);
    status = H5Pset_create_intermediate_group (gcpl6, 1);
    gcpl7 = H5Pcreate (H5P_LINK_CREATE);
    status = H5Pset_create_intermediate_group (gcpl7, 1);
    gcpl8 = H5Pcreate (H5P_LINK_CREATE);
    status = H5Pset_create_intermediate_group (gcpl8, 1);
    //gcpl9 = H5Pcreate (H5P_LINK_CREATE);
    //status = H5Pset_create_intermediate_group (gcpl9, 1);


    //dset = H5Dcreate (file, DATASET, filetype, space, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT);
    //status = H5Dwrite (dset, memtype, H5S_ALL, H5S_ALL, H5P_DEFAULT, wdata);


    /*
     * Create the group /G1/G2/G3.  Note that /G1 and /G1/G2 do not
     * exist yet.  This call would cause an error if we did not use the
     * previously created property list.
     */
    group1 = H5Gcreate (file, "/elementData", gcpl1, H5P_DEFAULT, H5P_DEFAULT);
    group2 = H5Gcreate (file, "/phaseData", gcpl2, H5P_DEFAULT, H5P_DEFAULT);
    group3 = H5Gcreate (file, "/reactionData", gcpl3, H5P_DEFAULT, H5P_DEFAULT);
    group4 = H5Gcreate (file, "/speciesData", gcpl4, H5P_DEFAULT, H5P_DEFAULT);
    group5 = H5Gcreate (file, "/reactionData/arrheniusData", gcpl5, H5P_DEFAULT, H5P_DEFAULT);
    group6 = H5Gcreate (file, "/reactionData/falloffData", gcpl6, H5P_DEFAULT, H5P_DEFAULT);
    group7 = H5Gcreate (file, "/reactionData/rateCoefLinks", gcpl7, H5P_DEFAULT, H5P_DEFAULT);
    group8 = H5Gcreate (file, "/speciesData/thermoPolynomials", gcpl8, H5P_DEFAULT, H5P_DEFAULT); 



     dset = H5Dcreate (file, DATASET, filetype, space, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT);
     status = H5Dwrite (dset, memtype, H5S_ALL, H5S_ALL, H5P_DEFAULT, wdata);
     dset2 = H5Dcreate (file, DATASET2, filetype2, space2, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT);
     status = H5Dwrite (dset2, memtype2, H5S_ALL, H5S_ALL, H5P_DEFAULT, wdata2);
     dset3 = H5Dcreate (file, DATASET3, filetype3, space3, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT);
     status = H5Dwrite (dset3, memtype3, H5S_ALL, H5S_ALL, H5P_DEFAULT, wdata3);                                                                                                                      
     dset4 = H5Dcreate (file, DATASET4, filetype4, space4, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT);
     status = H5Dwrite (dset4, memtype4, H5S_ALL, H5S_ALL, H5P_DEFAULT, wdata4);
     dset5 = H5Dcreate (file, DATASET5, filetype5, space5, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT);
     status = H5Dwrite (dset5, memtype5, H5S_ALL, H5S_ALL, H5P_DEFAULT, wdata5);
     dset6 = H5Dcreate (file, DATASET6, filetype6, space6, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT);
     status = H5Dwrite (dset6, memtype6, H5S_ALL, H5S_ALL, H5P_DEFAULT, wdata6);
     dset7 = H5Dcreate (file, DATASET7, filetype7, space7, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT);
     status = H5Dwrite (dset7, memtype7, H5S_ALL, H5S_ALL, H5P_DEFAULT, wdata7);
     dset8 = H5Dcreate (file, DATASET8, filetype8, space8, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT);
     status = H5Dwrite (dset8, memtype8, H5S_ALL, H5S_ALL, H5P_DEFAULT, wdata8);
     dset9 = H5Dcreate (file, DATASET9, filetype9, space9, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT);
     status = H5Dwrite (dset9, memtype9, H5S_ALL, H5S_ALL, H5P_DEFAULT, wdata9);
     dset10 = H5Dcreate (file, DATASET10, filetype10, space10, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT);
     status = H5Dwrite (dset10, memtype10, H5S_ALL, H5S_ALL, H5P_DEFAULT, wdata10);
     dset11 = H5Dcreate (file, DATASET11, filetype11, space11, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT);
     status = H5Dwrite (dset11, memtype11, H5S_ALL, H5S_ALL, H5P_DEFAULT, wdata11);
     dset12 = H5Dcreate (file, DATASET12, filetype12, space12, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT);
     status = H5Dwrite (dset12, memtype12, H5S_ALL, H5S_ALL, H5P_DEFAULT, wdata12);
     dset13 = H5Dcreate (file, DATASET13, filetype13, space13, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT);
     status = H5Dwrite (dset13, memtype13, H5S_ALL, H5S_ALL, H5P_DEFAULT, wdata13);
     dset14 = H5Dcreate (file, DATASET14, filetype14, space14, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT);
     status = H5Dwrite (dset14, memtype14, H5S_ALL, H5S_ALL, H5P_DEFAULT, wdata14);
     dset15 = H5Dcreate (file, DATASET15, filetype15, space15, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT);
     status = H5Dwrite (dset15, memtype15, H5S_ALL, H5S_ALL, H5P_DEFAULT, wdata15);
     dset16 = H5Dcreate (file, DATASET16, filetype16, space16, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT);
     status = H5Dwrite (dset16, memtype16, H5S_ALL, H5S_ALL, H5P_DEFAULT, wdata16);
     dset17 = H5Dcreate (file, DATASET17, filetype17, space17, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT);
     status = H5Dwrite (dset17, memtype17, H5S_ALL, H5S_ALL, H5P_DEFAULT, wdata17);
     dset18 = H5Dcreate (file, DATASET18, filetype18, space18, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT);
     status = H5Dwrite (dset18, memtype18, H5S_ALL, H5S_ALL, H5P_DEFAULT, wdata18);
     dset19 = H5Dcreate (file, DATASET19, filetype19, space19, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT);
     status = H5Dwrite (dset19, memtype19, H5S_ALL, H5S_ALL, H5P_DEFAULT, wdata19);
     dset20 = H5Dcreate (file, DATASET20, filetype20, space20, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT);
     status = H5Dwrite (dset20, memtype20, H5S_ALL, H5S_ALL, H5P_DEFAULT, wdata20);
     dset21 = H5Dcreate (file, DATASET21, filetype20, space21, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT);
     status = H5Dwrite (dset21, memtype21, H5S_ALL, H5S_ALL, H5P_DEFAULT, wdata21);
     dset22 = H5Dcreate (file, DATASET22, filetype22, space22, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT);
     status = H5Dwrite (dset22, memtype22, H5S_ALL, H5S_ALL, H5P_DEFAULT, wdata22);
     dset23 = H5Dcreate (file, DATASET23, filetype23, space23, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT);
     status = H5Dwrite (dset23, memtype23, H5S_ALL, H5S_ALL, H5P_DEFAULT, wdata23);
     dset24 = H5Dcreate (file, DATASET24, filetype24, space24, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT);
     status = H5Dwrite (dset24, memtype24, H5S_ALL, H5S_ALL, H5P_DEFAULT, wdata24);
     dset25 = H5Dcreate (file, DATASET25, filetype25, space25, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT);
     status = H5Dwrite (dset25, memtype25, H5S_ALL, H5S_ALL, H5P_DEFAULT, wdata25);
     dset26 = H5Dcreate (file, DATASET26, filetype26, space26, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT);
     status = H5Dwrite (dset26, memtype26, H5S_ALL, H5S_ALL, H5P_DEFAULT, wdata26);                                                                                                                   
     dset27 = H5Dcreate (file, DATASET27, filetype27, space27, H5P_DEFAULT, H5P_DEFAULT,H5P_DEFAULT); 
     status = H5Dwrite (dset27, memtype27, H5S_ALL, H5S_ALL, H5P_DEFAULT, wdata27);
     dset28 = H5Dcreate (file, DATASET28, filetype28, space28, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT);
     status = H5Dwrite (dset28, memtype28, H5S_ALL, H5S_ALL, H5P_DEFAULT, wdata28);                                                                                                                   
     dset29 = H5Dcreate (file, DATASET29, filetype29, space29, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT);                                                                                                
     status = H5Dwrite (dset29, memtype29, H5S_ALL, H5S_ALL, H5P_DEFAULT, wdata29);
     dset30 = H5Dcreate (file, DATASET30, filetype30, space30, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT);
     status = H5Dwrite (dset30, memtype30, H5S_ALL, H5S_ALL, H5P_DEFAULT, wdata30);  
     dset31 = H5Dcreate (file, DATASET31, filetype31, space31, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT);
     status = H5Dwrite (dset31, memtype31, H5S_ALL, H5S_ALL, H5P_DEFAULT, wdata31);
     dset32 = H5Dcreate (file, DATASET32, filetype32, space32, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT);
     status = H5Dwrite (dset32, memtype32, H5S_ALL, H5S_ALL, H5P_DEFAULT, wdata32);
     dset33 = H5Dcreate (file, DATASET33, filetype33, space33, H5P_DEFAULT, H5P_DEFAULT,H5P_DEFAULT);
     status = H5Dwrite (dset33, memtype33, H5S_ALL, H5S_ALL, H5P_DEFAULT, wdata33);
     dset34 = H5Dcreate (file, DATASET34, filetype34, space34, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT);                                                                                                
     status = H5Dwrite (dset34, memtype34, H5S_ALL, H5S_ALL, H5P_DEFAULT, wdata34);
     dset35 = H5Dcreate (file, DATASET35, filetype35, space35, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT);
     status = H5Dwrite (dset35, memtype35, H5S_ALL, H5S_ALL, H5P_DEFAULT, wdata35);
     dset36 = H5Dcreate (file, DATASET36, filetype36, space36, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT);    
     status = H5Dwrite (dset36, memtype36, H5S_ALL, H5S_ALL, H5P_DEFAULT, wdata36);                                                                                  



    // dset37 = H5Dcreate (file, DATASET37, filetype37, space37, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT);
     //status = H5Dwrite (dset37, memtype37, H5S_ALL, H5S_ALL, H5P_DEFAULT, wdata37);


                                                  
    /*
     * Print all the objects in the files to show that intermediate
     * groups have been created.  See h5ex_g_visit for more information
     * on how to use H5Ovisit.
     */
    printf ("Objects in the file:\n");
    status = H5Ovisit (file, H5_INDEX_NAME, H5_ITER_NATIVE, op_func, NULL);

   /*
     * Now we begin the read section of this example.  Here we assume
     * the dataset has the same name and rank, but can have any size.
     * Therefore we must allocate a new array to read in data using
     * malloc().
     */

    /*
     * Open file and dataset.
     */
    //file = H5Fopen (FILE, H5F_ACC_RDONLY, H5P_DEFAULT);
    //dset = H5Dopen (file, DATASET);

    /*
     * Get the datatype.
     */
    filetype = H5Dget_type (dset);
    filetype2 = H5Dget_type (dset2);
    filetype3 = H5Dget_type (dset3);
    filetype4 = H5Dget_type (dset4);
    filetype5 = H5Dget_type (dset5);
    filetype6 = H5Dget_type (dset6);
    filetype7 = H5Dget_type (dset7);
    filetype8 = H5Dget_type (dset8);
    filetype9 = H5Dget_type (dset9);
    filetype10 = H5Dget_type (dset10);
    filetype11 = H5Dget_type (dset11);
    filetype12 = H5Dget_type (dset12);
    filetype13 = H5Dget_type (dset13);
    filetype14 = H5Dget_type (dset14);
    filetype15 = H5Dget_type (dset15);
    filetype16 = H5Dget_type (dset16);
    filetype17 = H5Dget_type (dset17);
    filetype18 = H5Dget_type (dset18);
    filetype19 = H5Dget_type (dset19);
    filetype20 = H5Dget_type (dset20);
    filetype21 = H5Dget_type (dset21);
    filetype22 = H5Dget_type (dset22);
    filetype23 = H5Dget_type (dset23);
    filetype24 = H5Dget_type (dset24);
    filetype25 = H5Dget_type (dset25);
    filetype26 = H5Dget_type (dset26);
    filetype27 = H5Dget_type (dset27);
    filetype28 = H5Dget_type (dset28);
    filetype29 = H5Dget_type (dset29);
    filetype30 = H5Dget_type (dset30);
    filetype31 = H5Dget_type (dset31);
    filetype32 = H5Dget_type (dset32);
    filetype33 = H5Dget_type (dset33);
    filetype34 = H5Dget_type (dset34);
    filetype35 = H5Dget_type (dset35);
    filetype36 = H5Dget_type (dset36);
    //filetype37 = H5Dget_type (dset37);

    /*
     * Get dataspace and allocate memory for read buffer.
     */
    space = H5Dget_space (dset);
    ndims = H5Sget_simple_extent_dims (space, dims, NULL);
    rdata = (char **) malloc (dims[0] * sizeof (char *));
    space2 = H5Dget_space (dset2);
    ndims2 = H5Sget_simple_extent_dims (space2, dims2, NULL);
    rdata2 = (char **) malloc (dims2[0] * sizeof (char *));
    space3 = H5Dget_space (dset3);
    ndims3 = H5Sget_simple_extent_dims (space3, dims3, NULL);
    rdata3 = (char **) malloc (dims3[0] * sizeof (char *));
    space4 = H5Dget_space (dset4);
    ndims4 = H5Sget_simple_extent_dims (space4, dims4, NULL);
    rdata4 = (char **) malloc (dims4[0] * sizeof (char *));
    space5 = H5Dget_space (dset5);
    ndims5 = H5Sget_simple_extent_dims (space5, dims5, NULL);
    rdata5 = (char **) malloc (dims5[0] * sizeof (char *));
    space6 = H5Dget_space (dset6);
    ndims6 = H5Sget_simple_extent_dims (space6, dims6, NULL);
    rdata6 = (char **) malloc (dims6[0] * sizeof (char *));
    space7 = H5Dget_space (dset7);
    ndims7 = H5Sget_simple_extent_dims (space7, dims7, NULL);
    rdata7 = (char **) malloc (dims7[0] * sizeof (char *));
    space8 = H5Dget_space (dset8);
    ndims8 = H5Sget_simple_extent_dims (space8, dims8, NULL);
    rdata8 = (char **) malloc (dims8[0] * sizeof (char *));
    space9 = H5Dget_space (dset9);
    ndims9 = H5Sget_simple_extent_dims (space9, dims9, NULL);
    rdata9 = (char **) malloc (dims9[0] * sizeof (char *));
    space10 = H5Dget_space (dset10);
    ndims10 = H5Sget_simple_extent_dims (space10, dims10, NULL);
    rdata10 = (char **) malloc (dims10[0] * sizeof (char *));
    space11 = H5Dget_space (dset11);
    ndims11 = H5Sget_simple_extent_dims (space11, dims11, NULL);
    rdata11 = (char **) malloc (dims11[0] * sizeof (char *));
    space12 = H5Dget_space (dset12);
    ndims12 = H5Sget_simple_extent_dims (space12, dims12, NULL);
    rdata12 = (char **) malloc (dims12[0] * sizeof (char *));
    space13 = H5Dget_space (dset13);
    ndims13 = H5Sget_simple_extent_dims (space13, dims13, NULL);
    rdata13 = (char **) malloc (dims13[0] * sizeof (char *));
    space14 = H5Dget_space (dset14);
    ndims14 = H5Sget_simple_extent_dims (space14, dims14, NULL);
    rdata14 = (char **) malloc (dims14[0] * sizeof (char *));
    space15 = H5Dget_space (dset15);
    ndims15 = H5Sget_simple_extent_dims (space15, dims15, NULL);
    rdata15 = (char **) malloc (dims15[0] * sizeof (char *));
    space16 = H5Dget_space (dset16);
    ndims16 = H5Sget_simple_extent_dims (space16, dims16, NULL);
    rdata16 = (char **) malloc (dims16[0] * sizeof (char *));
    space17 = H5Dget_space (dset17);
    ndims17 = H5Sget_simple_extent_dims (space17, dims17, NULL);
    rdata17 = (char **) malloc (dims17[0] * sizeof (char *));
    space18 = H5Dget_space (dset18);
    ndims18 = H5Sget_simple_extent_dims (space18, dims18, NULL);
    rdata18 = (char **) malloc (dims18[0] * sizeof (char *));
    space19 = H5Dget_space (dset19);
    ndims19 = H5Sget_simple_extent_dims (space19, dims19, NULL);
    rdata19 = (char **) malloc (dims19[0] * sizeof (char *));
    space20 = H5Dget_space (dset20);
    ndims20 = H5Sget_simple_extent_dims (space20, dims20, NULL);
    rdata20 = (char **) malloc (dims20[0] * sizeof (char *));
    space21 = H5Dget_space (dset21);
    ndims21 = H5Sget_simple_extent_dims (space21, dims21, NULL);
    rdata21 = (char **) malloc (dims21[0] * sizeof (char *));
    space22 = H5Dget_space (dset22);
    ndims22 = H5Sget_simple_extent_dims (space22, dims22, NULL);
    rdata22 = (char **) malloc (dims22[0] * sizeof (char *));
    space23 = H5Dget_space (dset23);
    ndims23 = H5Sget_simple_extent_dims (space23, dims23, NULL);
    rdata23 = (char **) malloc (dims23[0] * sizeof (char *));
    space24 = H5Dget_space (dset24);
    ndims24 = H5Sget_simple_extent_dims (space24, dims24, NULL);
    rdata24 = (char **) malloc (dims24[0] * sizeof (char *));
    space25 = H5Dget_space (dset25);
    ndims25 = H5Sget_simple_extent_dims (space25, dims25, NULL);
    rdata25 = (char **) malloc (dims25[0] * sizeof (char *));
    space26 = H5Dget_space (dset26);
    ndims26 = H5Sget_simple_extent_dims (space26, dims26, NULL);
    rdata26 = (char **) malloc (dims26[0] * sizeof (char *));
    space27 = H5Dget_space (dset27);
    ndims27 = H5Sget_simple_extent_dims (space27, dims27, NULL);
    rdata27 = (char **) malloc (dims27[0] * sizeof (char *));
    space28 = H5Dget_space (dset28);
    ndims28 = H5Sget_simple_extent_dims (space28, dims28, NULL);
    rdata28 = (char **) malloc (dims28[0] * sizeof (char *));
    space29 = H5Dget_space (dset29);
    ndims29 = H5Sget_simple_extent_dims (space29, dims29, NULL);
    rdata29 = (char **) malloc (dims29[0] * sizeof (char *));
    space30 = H5Dget_space (dset30);
    ndims30 = H5Sget_simple_extent_dims (space30, dims30, NULL);
    rdata30 = (char **) malloc (dims30[0] * sizeof (char *));
    space31 = H5Dget_space (dset31);
    ndims31 = H5Sget_simple_extent_dims (space31, dims31, NULL);
    rdata31 = (char **) malloc (dims31[0] * sizeof (char *));
    space32 = H5Dget_space (dset32);
    ndims32 = H5Sget_simple_extent_dims (space32, dims32, NULL);                                                                                                                                      
    rdata32 = (char **) malloc (dims32[0] * sizeof (char *));                                                                                                                                         
    space33 = H5Dget_space (dset33);                                                                                                                                                                  
    ndims33 = H5Sget_simple_extent_dims (space33, dims33, NULL);                                                                                                                                      
    rdata33 = (char **) malloc (dims33[0] * sizeof (char *));
    space34 = H5Dget_space (dset34);
    ndims34 = H5Sget_simple_extent_dims (space34, dims34, NULL);
    rdata34 = (char **) malloc (dims34[0] * sizeof (char *));
    space35 = H5Dget_space (dset35);
    ndims35 = H5Sget_simple_extent_dims (space35, dims35, NULL);
    rdata35 = (char **) malloc (dims35[0] * sizeof (char *));
    space36 = H5Dget_space (dset36);
    ndims36 = H5Sget_simple_extent_dims (space36, dims36, NULL);
    rdata36 = (char **) malloc (dims36[0] * sizeof (char *));


    //space37 = H5Dget_space (dset37);
    //ndims37 = H5Sget_simple_extent_dims (space37, dims37, NULL);
    //rdata37 = (float **) malloc (dims37[0][0] * sizeof (float *));


    /*
     * Create the memory datatype.
     */
    memtype = H5Tcopy (H5T_C_S1);
    status = H5Tset_size (memtype, H5T_VARIABLE);
    memtype2 = H5Tcopy (H5T_C_S1);
    status = H5Tset_size (memtype2, H5T_VARIABLE);
    memtype3 = H5Tcopy (H5T_C_S1);
    status = H5Tset_size (memtype3, H5T_VARIABLE);
    memtype4 = H5Tcopy (H5T_C_S1);
    status = H5Tset_size (memtype4, H5T_VARIABLE);
    memtype5 = H5Tcopy (H5T_C_S1);
    status = H5Tset_size (memtype5, H5T_VARIABLE);
    memtype6 = H5Tcopy (H5T_C_S1);
    status = H5Tset_size (memtype6, H5T_VARIABLE);
    memtype7 = H5Tcopy (H5T_C_S1);
    status = H5Tset_size (memtype7, H5T_VARIABLE);
    memtype8 = H5Tcopy (H5T_C_S1);
    status = H5Tset_size (memtype8, H5T_VARIABLE);
    memtype9 = H5Tcopy (H5T_C_S1);
    status = H5Tset_size (memtype9, H5T_VARIABLE);
    memtype10 = H5Tcopy (H5T_C_S1);
    status = H5Tset_size (memtype10, H5T_VARIABLE);
    memtype11 = H5Tcopy (H5T_C_S1);
    status = H5Tset_size (memtype11, H5T_VARIABLE);
    memtype12 = H5Tcopy (H5T_C_S1);
    status = H5Tset_size (memtype12, H5T_VARIABLE);
    memtype13 = H5Tcopy (H5T_C_S1);
    status = H5Tset_size (memtype13, H5T_VARIABLE);
    memtype14 = H5Tcopy (H5T_C_S1);
    status = H5Tset_size (memtype14, H5T_VARIABLE);
    memtype15 = H5Tcopy (H5T_C_S1);
    status = H5Tset_size (memtype15, H5T_VARIABLE);
    memtype16 = H5Tcopy (H5T_C_S1);
    status = H5Tset_size (memtype16, H5T_VARIABLE);
    memtype17 = H5Tcopy (H5T_C_S1);
    status = H5Tset_size (memtype17, H5T_VARIABLE);
    memtype18 = H5Tcopy (H5T_C_S1);
    status = H5Tset_size (memtype18, H5T_VARIABLE);
    memtype19 = H5Tcopy (H5T_C_S1);
    status = H5Tset_size (memtype19, H5T_VARIABLE);
    memtype20 = H5Tcopy (H5T_C_S1);
    status = H5Tset_size (memtype20, H5T_VARIABLE);
    memtype21 = H5Tcopy (H5T_C_S1);                                                                                                                                                                   
    status = H5Tset_size (memtype21, H5T_VARIABLE);
    memtype22 = H5Tcopy (H5T_C_S1);
    status = H5Tset_size (memtype22, H5T_VARIABLE);
    memtype23 = H5Tcopy (H5T_C_S1);                                                                                                                                                                   
    status = H5Tset_size (memtype23, H5T_VARIABLE);                                                                                                                                                   
    memtype24 = H5Tcopy (H5T_C_S1);                                                                                                                                                                   
    status = H5Tset_size (memtype24, H5T_VARIABLE);                                                                                                                                                   
    memtype25 = H5Tcopy (H5T_C_S1);                                                                                                                                                                   
    status = H5Tset_size (memtype25, H5T_VARIABLE);                                                                                                                                                   
    memtype26 = H5Tcopy (H5T_C_S1);
    status = H5Tset_size (memtype26, H5T_VARIABLE);
    memtype27 = H5Tcopy (H5T_C_S1);                                                                                                                                                                   
    status = H5Tset_size (memtype27, H5T_VARIABLE);                                                                                                                                                   
    memtype28 = H5Tcopy (H5T_C_S1);                                                                                                                                                                   
    status = H5Tset_size (memtype28, H5T_VARIABLE);                                                                                                                                                   
    memtype29 = H5Tcopy (H5T_C_S1);                                                                                                                                                                   
    status = H5Tset_size (memtype29, H5T_VARIABLE);                                                                                                                                                   
    memtype30 = H5Tcopy (H5T_C_S1);
    status = H5Tset_size (memtype30, H5T_VARIABLE);
    memtype31 = H5Tcopy (H5T_C_S1);                                                                                                                                                                   
    status = H5Tset_size (memtype31, H5T_VARIABLE);                                                                                                                                                   
    memtype32 = H5Tcopy (H5T_C_S1);                                                                                                                                                                   
    status = H5Tset_size (memtype32, H5T_VARIABLE);                                                                                                                                                   
    memtype33 = H5Tcopy (H5T_C_S1);                                                                                                                                                                   
    status = H5Tset_size (memtype33, H5T_VARIABLE);
    memtype34 = H5Tcopy (H5T_C_S1);
    status = H5Tset_size (memtype34, H5T_VARIABLE);
    memtype35 = H5Tcopy (H5T_C_S1);
    status = H5Tset_size (memtype35, H5T_VARIABLE);
    memtype36 = H5Tcopy (H5T_C_S1);
    status = H5Tset_size (memtype36, H5T_VARIABLE);


    //memtype37 = H5Tcopy (H5T_C_S1);
    //status = H5Tset_size (memtype37, H5T_VARIABLE);

    /*
     * Read the data.
     */
    status = H5Dread (dset, memtype, H5S_ALL, H5S_ALL, H5P_DEFAULT, rdata);
    status = H5Dread (dset2, memtype2, H5S_ALL, H5S_ALL, H5P_DEFAULT, rdata2);
    status = H5Dread (dset3, memtype3, H5S_ALL, H5S_ALL, H5P_DEFAULT, rdata3);
    status = H5Dread (dset4, memtype4, H5S_ALL, H5S_ALL, H5P_DEFAULT, rdata4);
    status = H5Dread (dset5, memtype5, H5S_ALL, H5S_ALL, H5P_DEFAULT, rdata5);
    status = H5Dread (dset6, memtype6, H5S_ALL, H5S_ALL, H5P_DEFAULT, rdata6);
    status = H5Dread (dset7, memtype7, H5S_ALL, H5S_ALL, H5P_DEFAULT, rdata7);
    status = H5Dread (dset8, memtype8, H5S_ALL, H5S_ALL, H5P_DEFAULT, rdata8);
    status = H5Dread (dset9, memtype9, H5S_ALL, H5S_ALL, H5P_DEFAULT, rdata9);
    status = H5Dread (dset10, memtype10, H5S_ALL, H5S_ALL, H5P_DEFAULT, rdata10);
    status = H5Dread (dset11, memtype11, H5S_ALL, H5S_ALL, H5P_DEFAULT, rdata11);
    status = H5Dread (dset12, memtype12, H5S_ALL, H5S_ALL, H5P_DEFAULT, rdata12);
    status = H5Dread (dset13, memtype13, H5S_ALL, H5S_ALL, H5P_DEFAULT, rdata13);
    status = H5Dread (dset14, memtype14, H5S_ALL, H5S_ALL, H5P_DEFAULT, rdata14);
    status = H5Dread (dset15, memtype15, H5S_ALL, H5S_ALL, H5P_DEFAULT, rdata15);
    status = H5Dread (dset16, memtype16, H5S_ALL, H5S_ALL, H5P_DEFAULT, rdata16);
    status = H5Dread (dset17, memtype17, H5S_ALL, H5S_ALL, H5P_DEFAULT, rdata17);
    status = H5Dread (dset18, memtype18, H5S_ALL, H5S_ALL, H5P_DEFAULT, rdata18);
    status = H5Dread (dset19, memtype19, H5S_ALL, H5S_ALL, H5P_DEFAULT, rdata19);
    status = H5Dread (dset20, memtype20, H5S_ALL, H5S_ALL, H5P_DEFAULT, rdata20);
    status = H5Dread (dset21, memtype21, H5S_ALL, H5S_ALL, H5P_DEFAULT, rdata21);
    status = H5Dread (dset22, memtype22, H5S_ALL, H5S_ALL, H5P_DEFAULT, rdata22);
    status = H5Dread (dset23, memtype23, H5S_ALL, H5S_ALL, H5P_DEFAULT, rdata23);
    status = H5Dread (dset24, memtype24, H5S_ALL, H5S_ALL, H5P_DEFAULT, rdata24);
    status = H5Dread (dset25, memtype25, H5S_ALL, H5S_ALL, H5P_DEFAULT, rdata25);
    status = H5Dread (dset26, memtype26, H5S_ALL, H5S_ALL, H5P_DEFAULT, rdata26);
    status = H5Dread (dset27, memtype27, H5S_ALL, H5S_ALL, H5P_DEFAULT, rdata27);
    status = H5Dread (dset28, memtype28, H5S_ALL, H5S_ALL, H5P_DEFAULT, rdata28);
    status = H5Dread (dset29, memtype29, H5S_ALL, H5S_ALL, H5P_DEFAULT, rdata29);
    status = H5Dread (dset30, memtype30, H5S_ALL, H5S_ALL, H5P_DEFAULT, rdata30);
    status = H5Dread (dset31, memtype31, H5S_ALL, H5S_ALL, H5P_DEFAULT, rdata31);
    status = H5Dread (dset32, memtype32, H5S_ALL, H5S_ALL, H5P_DEFAULT, rdata32);
    status = H5Dread (dset33, memtype33, H5S_ALL, H5S_ALL, H5P_DEFAULT, rdata33);
    status = H5Dread (dset34, memtype34, H5S_ALL, H5S_ALL, H5P_DEFAULT, rdata34);
    status = H5Dread (dset35, memtype35, H5S_ALL, H5S_ALL, H5P_DEFAULT, rdata35);
    status = H5Dread (dset36, memtype36, H5S_ALL, H5S_ALL, H5P_DEFAULT, rdata36);

    //status = H5Dread (dset37, memtype37, H5S_ALL, H5S_ALL, H5P_DEFAULT, rdata37);


    /*
     * Output the data to the screen.
     */
    for (i=0; i<dims[0]; i++)
        printf ("%s[%d]: %s\n", DATASET, i, rdata[i]);
    for (i2=0; i2<dims2[0]; i2++)
        printf ("%s[%d]: %s\n", DATASET2, i2, rdata2[i2]);
    for (i3=0; i3<dims3[0]; i3++)
        printf ("%s[%d]: %s\n", DATASET3, i3, rdata3[i3]);
    for (i4=0; i4<dims4[0]; i4++)
        printf ("%s[%d]: %s\n", DATASET4, i4, rdata4[i4]);
    for (i5=0; i5<dims5[0]; i5++)
        printf ("%s[%d]: %s\n", DATASET5, i5, rdata5[i5]);
    for (i6=0; i6<dims6[0]; i6++)
        printf ("%s[%d]: %s\n", DATASET6, i6, rdata6[i6]);
    for (i7=0; i7<dims7[0]; i7++)
        printf ("%s[%d]: %s\n", DATASET7, i7, rdata7[i7]);
    for (i8=0; i8<dims8[0]; i8++)
        printf ("%s[%d]: %s\n", DATASET8, i8, rdata8[i8]);
    for (i9=0; i9<dims9[0]; i9++)
        printf ("%s[%d]: %s\n", DATASET9, i9, rdata9[i9]);
    for (i10=0; i10<dims10[0]; i10++)
        printf ("%s[%d]: %s\n", DATASET10, i10, rdata10[i10]);
    for (i11=0; i11<dims11[0]; i11++)
        printf ("%s[%d]: %s\n", DATASET11, i11, rdata11[i11]);
    for (i12=0; i12<dims12[0]; i12++)
        printf ("%s[%d]: %s\n", DATASET12, i12, rdata12[i12]);
    for (i13=0; i13<dims13[0]; i13++)
        printf ("%s[%d]: %s\n", DATASET13, i13, rdata13[i13]);
    for (i14=0; i14<dims14[0]; i14++)
        printf ("%s[%d]: %s\n", DATASET14, i14, rdata14[i14]);
    for (i15=0; i15<dims15[0]; i15++)
        printf ("%s[%d]: %s\n", DATASET15, i15, rdata15[i15]);
    for (i16=0; i16<dims16[0]; i16++)
        printf ("%s[%d]: %s\n", DATASET16, i16, rdata16[i16]);
    for (i17=0; i17<dims17[0]; i17++)
        printf ("%s[%d]: %s\n", DATASET17, i17, rdata17[i17]);
    for (i18=0; i18<dims18[0]; i18++)
        printf ("%s[%d]: %s\n", DATASET18, i18, rdata18[i18]);
    for (i19=0; i19<dims19[0]; i19++)
        printf ("%s[%d]: %s\n", DATASET19, i19, rdata19[i19]);
    for (i20=0; i20<dims20[0]; i20++)
        printf ("%s[%d]: %s\n", DATASET20, i20, rdata20[i20]);
    for (i21=0; i21<dims21[0]; i21++)                                                                                                                                                                 
        printf ("%s[%d]: %s\n", DATASET21, i21, rdata21[i21]);                                                                                                                                        
    for (i22=0; i22<dims22[0]; i22++)
        printf ("%s[%d]: %s\n", DATASET22, i22, rdata22[i22]);
    for (i23=0; i23<dims23[0]; i23++)                                                                                                                                                                 
        printf ("%s[%d]: %s\n", DATASET23, i23, rdata23[i23]);                                                                                                                                        
    for (i24=0; i24<dims24[0]; i24++)                                                                                                                                                                 
        printf ("%s[%d]: %s\n", DATASET24, i24, rdata24[i24]);                                                                                                                                        
    for (i25=0; i25<dims25[0]; i25++)                                                                                                                                                                 
        printf ("%s[%d]: %s\n", DATASET25, i25, rdata25[i25]);                                                                                                                                        
    for (i26=0; i26<dims26[0]; i26++)
        printf ("%s[%d]: %s\n", DATASET26, i26, rdata26[i26]);
    for (i27=0; i27<dims27[0]; i27++)                                                                                                                                                                 
        printf ("%s[%d]: %s\n", DATASET27, i27, rdata27[i27]);                                                                                                                                        
    for (i28=0; i28<dims28[0]; i28++)                                                                                                                                                                 
        printf ("%s[%d]: %s\n", DATASET28, i28, rdata28[i28]);                                                                                                                                        
    for (i29=0; i29<dims29[0]; i29++)                                                                                                                                                                 
        printf ("%s[%d]: %s\n", DATASET29, i29, rdata29[i29]);
    for (i30=0; i30<dims30[0]; i30++)                                                                                                                                                                 
        printf ("%s[%d]: %s\n", DATASET30, i30, rdata30[i30]);
    for (i31=0; i31<dims31[0]; i31++)
        printf ("%s[%d]: %s\n", DATASET31, i31, rdata31[i31]);
    for (i32=0; i32<dims32[0]; i32++)                                                                                                                                                                 
        printf ("%s[%d]: %s\n", DATASET32, i32, rdata32[i32]);                                                                                                                                        
    for (i33=0; i33<dims33[0]; i33++)
        printf ("%s[%d]: %s\n", DATASET33, i33, rdata33[i33]);
    for (i34=0; i34<dims34[0]; i34++)
        printf ("%s[%d]: %s\n", DATASET34, i34, rdata34[i34]);
    for (i35=0; i35<dims35[0]; i35++)
        printf ("%s[%d]: %s\n", DATASET35, i35, rdata35[i35]);
    for (i36=0; i36<dims36[0]; i36++)                                                                                                                                                                 
        printf ("%s[%d]: %s\n", DATASET36, i36, rdata36[i36]);                                                                                                                                        
    //for (i37=0; i37<dims37[0]; i37++)
    //    printf ("%s[%d]: %s\n", DATASET37, i37, rdata37[i37]);


/*    printf ("\n");
    for (myi=0; myi<NX2; myi++)
    {
        for (myj=0; myj<NY2; myj++)
        {
         printf (" %s[%d][%d]: %f", DATASET37, myi, myj, rdataf2[myi][myj]);
        }
        printf ("\n");
   }
*/





    /*
     * Close and release resources.  Note that H5Dvlen_reclaim works
     * for variable-length strings as well as variable-length arrays.
     * Also note that we must still free the array of pointers stored
     * in rdata, as H5Tvlen_reclaim only frees the data these point to.
     */
    status = H5Dvlen_reclaim (memtype, space, H5P_DEFAULT, rdata);
    status = H5Dvlen_reclaim (memtype2, space2, H5P_DEFAULT, rdata2);
    status = H5Dvlen_reclaim (memtype3, space3, H5P_DEFAULT, rdata3);
    status = H5Dvlen_reclaim (memtype4, space4, H5P_DEFAULT, rdata4);
    status = H5Dvlen_reclaim (memtype5, space5, H5P_DEFAULT, rdata5);
    status = H5Dvlen_reclaim (memtype6, space6, H5P_DEFAULT, rdata6);
    status = H5Dvlen_reclaim (memtype7, space7, H5P_DEFAULT, rdata7);
    status = H5Dvlen_reclaim (memtype8, space8, H5P_DEFAULT, rdata8);
    status = H5Dvlen_reclaim (memtype9, space9, H5P_DEFAULT, rdata9);
    status = H5Dvlen_reclaim (memtype10, space10, H5P_DEFAULT, rdata10);
    status = H5Dvlen_reclaim (memtype11, space11, H5P_DEFAULT, rdata11);
    status = H5Dvlen_reclaim (memtype12, space12, H5P_DEFAULT, rdata12);
    status = H5Dvlen_reclaim (memtype13, space13, H5P_DEFAULT, rdata13);
    status = H5Dvlen_reclaim (memtype14, space14, H5P_DEFAULT, rdata14);
    status = H5Dvlen_reclaim (memtype15, space15, H5P_DEFAULT, rdata15);
    status = H5Dvlen_reclaim (memtype16, space16, H5P_DEFAULT, rdata16);
    status = H5Dvlen_reclaim (memtype17, space17, H5P_DEFAULT, rdata17);
    status = H5Dvlen_reclaim (memtype18, space18, H5P_DEFAULT, rdata18);
    status = H5Dvlen_reclaim (memtype19, space19, H5P_DEFAULT, rdata19);
    status = H5Dvlen_reclaim (memtype20, space20, H5P_DEFAULT, rdata20);
    status = H5Dvlen_reclaim (memtype21, space21, H5P_DEFAULT, rdata21);
    status = H5Dvlen_reclaim (memtype22, space22, H5P_DEFAULT, rdata22);
    status = H5Dvlen_reclaim (memtype23, space23, H5P_DEFAULT, rdata23);
    status = H5Dvlen_reclaim (memtype24, space24, H5P_DEFAULT, rdata24);
    status = H5Dvlen_reclaim (memtype25, space25, H5P_DEFAULT, rdata25);
    status = H5Dvlen_reclaim (memtype26, space26, H5P_DEFAULT, rdata26);
    status = H5Dvlen_reclaim (memtype27, space27, H5P_DEFAULT, rdata27);
    status = H5Dvlen_reclaim (memtype28, space28, H5P_DEFAULT, rdata28);
    status = H5Dvlen_reclaim (memtype29, space29, H5P_DEFAULT, rdata29);
    status = H5Dvlen_reclaim (memtype30, space30, H5P_DEFAULT, rdata30);
    status = H5Dvlen_reclaim (memtype31, space31, H5P_DEFAULT, rdata31);
    status = H5Dvlen_reclaim (memtype32, space32, H5P_DEFAULT, rdata32);
    status = H5Dvlen_reclaim (memtype33, space33, H5P_DEFAULT, rdata33);
    status = H5Dvlen_reclaim (memtype34, space34, H5P_DEFAULT, rdata34);
    status = H5Dvlen_reclaim (memtype35, space35, H5P_DEFAULT, rdata35);
    status = H5Dvlen_reclaim (memtype36, space36, H5P_DEFAULT, rdata36);

    //status = H5Dvlen_reclaim (memtype37, space37, H5P_DEFAULT, rdata37);



    free (rdata);
    free (rdata2);
    free (rdata3);
    free (rdata4);
    free (rdata5);
    free (rdata6);
    free (rdata7);
    free (rdata8);
    free (rdata9);
    free (rdata10);
    free (rdata11);
    free (rdata12);
    free (rdata13);
    free (rdata14);
    free (rdata15);
    free (rdata16);
    free (rdata17);
    free (rdata18);
    free (rdata19);
    free (rdata20);
    free (rdata21);
    free (rdata22);
    free (rdata23);
    free (rdata24);
    free (rdata25);
    free (rdata26);
    free (rdata27);
    free (rdata28);
    free (rdata29);
    free (rdata30);
    free (rdata31);
    free (rdata32);
    free (rdata33);
    free (rdata34);
    free (rdata35);
    free (rdata36);
    //free (rdata37);
    //free (rdataf2);
   // status = H5Sclose(space37);
   // status = H5Dclose(dset37);


    status = H5Dclose (dset);
    status = H5Dclose (dset2);
    status = H5Dclose (dset3);
    status = H5Dclose (dset4);
    status = H5Dclose (dset5);
    status = H5Dclose (dset6);
    status = H5Dclose (dset7);
    status = H5Dclose (dset8);
    status = H5Dclose (dset9);
    status = H5Dclose (dset10);
    status = H5Dclose (dset11);
    status = H5Dclose (dset12);
    status = H5Dclose (dset13);
    status = H5Dclose (dset14);
    status = H5Dclose (dset15);
    status = H5Dclose (dset16);
    status = H5Dclose (dset17);
    status = H5Dclose (dset18);
    status = H5Dclose (dset19);
    status = H5Dclose (dset20);
    status = H5Dclose (dset21);
    status = H5Dclose (dset22);
    status = H5Dclose (dset23);
    status = H5Dclose (dset24);
    status = H5Dclose (dset25);
    status = H5Dclose (dset26);
    status = H5Dclose (dset27);
    status = H5Dclose (dset28);
    status = H5Dclose (dset29);
    status = H5Dclose (dset30);
    status = H5Dclose (dset31);
    status = H5Dclose (dset32);
    status = H5Dclose (dset33);
    status = H5Dclose (dset34);
    status = H5Dclose (dset35);
    status = H5Dclose (dset36);
    //status = H5Dclose(dset37); 


    status = H5Sclose (space);
    status = H5Sclose (space2);
    status = H5Sclose (space3);
    status = H5Sclose (space4);
    status = H5Sclose (space5);
    status = H5Sclose (space6);
    status = H5Sclose (space7);
    status = H5Sclose (space8);
    status = H5Sclose (space9);
    status = H5Sclose (space10);
    status = H5Sclose (space11);
    status = H5Sclose (space12);
    status = H5Sclose (space13);
    status = H5Sclose (space14);
    status = H5Sclose (space15);
    status = H5Sclose (space16);
    status = H5Sclose (space17);
    status = H5Sclose (space18);
    status = H5Sclose (space19);
    status = H5Sclose (space20);
    status = H5Sclose (space21);
    status = H5Sclose (space22);
    status = H5Sclose (space23);
    status = H5Sclose (space24);
    status = H5Sclose (space25);
    status = H5Sclose (space26);
    status = H5Sclose (space27);
    status = H5Sclose (space28);
    status = H5Sclose (space29);
    status = H5Sclose (space30);
    status = H5Sclose (space31);
    status = H5Sclose (space32);
    status = H5Sclose (space33);
    status = H5Sclose (space34);
    status = H5Sclose (space35);
    status = H5Sclose (space36);
  //  status = H5Sclose(space37);
    



    status = H5Tclose (filetype);
    status = H5Tclose (filetype2);
    status = H5Tclose (filetype3);
    status = H5Tclose (filetype4);
    status = H5Tclose (filetype5);
    status = H5Tclose (filetype6);
    status = H5Tclose (filetype7);
    status = H5Tclose (filetype8);
    status = H5Tclose (filetype9);
    status = H5Tclose (filetype10);
    status = H5Tclose (filetype11);
    status = H5Tclose (filetype12);
    status = H5Tclose (filetype13);
    status = H5Tclose (filetype14);
    status = H5Tclose (filetype15);
    status = H5Tclose (filetype16);
    status = H5Tclose (filetype17);
    status = H5Tclose (filetype18);
    status = H5Tclose (filetype19);
    status = H5Tclose (filetype20);
    status = H5Tclose (filetype21);
    status = H5Tclose (filetype22);
    status = H5Tclose (filetype23);
    status = H5Tclose (filetype24);
    status = H5Tclose (filetype25);
    status = H5Tclose (filetype26);
    status = H5Tclose (filetype27);
    status = H5Tclose (filetype28);
    status = H5Tclose (filetype29);
    status = H5Tclose (filetype30);
    status = H5Tclose (filetype31);
    status = H5Tclose (filetype32);
    status = H5Tclose (filetype33);
    status = H5Tclose (filetype34);
    status = H5Tclose (filetype35);
    status = H5Tclose (filetype36);
//    status = H5Tclose (filetype37);


    status = H5Tclose (memtype);
    status = H5Tclose (memtype2);
    status = H5Tclose (memtype3);
    status = H5Tclose (memtype4);
    status = H5Tclose (memtype5);
    status = H5Tclose (memtype6);
    status = H5Tclose (memtype7);
    status = H5Tclose (memtype8);
    status = H5Tclose (memtype9);
    status = H5Tclose (memtype10);
    status = H5Tclose (memtype11);
    status = H5Tclose (memtype12);
    status = H5Tclose (memtype13);
    status = H5Tclose (memtype14);
    status = H5Tclose (memtype15);
    status = H5Tclose (memtype16);
    status = H5Tclose (memtype17);
    status = H5Tclose (memtype18);
    status = H5Tclose (memtype19);
    status = H5Tclose (memtype20);
    status = H5Tclose (memtype21);
    status = H5Tclose (memtype22);
    status = H5Tclose (memtype23);
    status = H5Tclose (memtype24);
    status = H5Tclose (memtype25);
    status = H5Tclose (memtype26);
    status = H5Tclose (memtype27);
    status = H5Tclose (memtype28);
    status = H5Tclose (memtype29);
    status = H5Tclose (memtype30);
    status = H5Tclose (memtype31);
    status = H5Tclose (memtype32);
    status = H5Tclose (memtype33);
    status = H5Tclose (memtype34);
    status = H5Tclose (memtype35);
    status = H5Tclose (memtype36);
   // status = H5Tclose (memtype37);

    status = H5Pclose (gcpl1);
    status = H5Pclose (gcpl2);
    status = H5Pclose (gcpl3);
    status = H5Pclose (gcpl4);
    status = H5Pclose (gcpl5);
    status = H5Pclose (gcpl6);
    status = H5Pclose (gcpl7);
    status = H5Pclose (gcpl8);
    //status = H5Pclose (gcpl9);
    status = H5Gclose (group1);
    status = H5Gclose (group2);
    status = H5Gclose (group3);
    status = H5Gclose (group4);
    status = H5Gclose (group5);
    status = H5Gclose (group6);
    status = H5Gclose (group7);
    status = H5Gclose (group8);
    //status = H5Gclose (group9);
    status = H5Fclose (file);

    return 0;
   
}



/************************************************************

  Operator function for H5Ovisit.  This function prints the
  name and type of the object passed to it.

 ************************************************************/
herr_t op_func (hid_t loc_id, const char *name, const H5O_info_t *info,
            void *operator_data)
{ 
    printf ("/");               /* Print root group in object path */

    /*
     * Check if the current object is the root group, and if not print
     * the full path name and type.
     */
    if (name[0] == '.')         /* Root group, do not print '.' */
        printf ("  (Group)\n");
    else
        switch (info->type) {
            case H5O_TYPE_GROUP:
                printf ("%s  (Group)\n", name);
                break;
            case H5O_TYPE_DATASET:
                printf ("%s  (Dataset)\n", name);
                break;
            case H5O_TYPE_NAMED_DATATYPE:
                printf ("%s  (Datatype)\n", name);
                break;
            default:
                printf ("%s  (Unknown)\n", name);
        }

    return 0;
}