/*  $Id: testParmArray.c 3304 2011-01-17 15:25:59Z brideout $ */

/*
modification history
--------------------
00a,22Apr00         original
*/

/* 
* USAGE: testParmArray cedarFile fileType
*    e.f.: testParmArray mil991102g.001 3
*    This program is a simple example illustrating the use of the madrec
*    library. It is hardwired to open and process
*    $MADROOT/experiments/1998/mlh/20jan98/mil980120g.002. 
* 
*    Supported formats:
*
*       Open Cedar file for sequential reading:
*         1 - Determine file type automatically
*        10 - Madrigal
*        11 - Blocked Binary
*        12 - Cbf
*        13 - Unblocked Binary
*        14 - Ascii
*
*       Open Cedar file for random reading using memory image of file; 
*         3 - Determine file type automatically
*        30 - Madrigal file
*        31 - Blocked Binary file
*        32 - Cbf file
*        33 - Unblocked Binary file
*        34 - Ascii file
*/

#include <stdlib.h>
#include <string.h>
#include <sys/types.h>
#include <sys/stat.h>
#include <fcntl.h>
#include <stdio.h>
#include <cedarIO.h>
#include <madrec.h>
#include <cedar.h>

int
main (argc, argv)
     int     argc;
     char    *argv[];
{
    Madrec *madrecp;
    char madfile[128];
    int iotype, i, j, k, nlines, nrecs, ind;

    double *parp;
    Ffspec fspec;
    int pparms[7] = {34, 120, 132, 142, 550, 10, 160};

    int fparms[4] = {34, 120, 142, 110};
    double fmin[4] = { 2.0, 200.0,  80.0, 300.0};
    double fmax[4] = { 2.5, 600.0,  89.0, 400.0};


    /* Get the name of the madrigal file */
    if (argc != 3) {
        cedarGetMadroot(madfile);
        strcat(madfile, "/experiments/1998/mlh/20jan98/mil980120g.002");
        iotype = 1;
    } else {
        strcpy(madfile, argv[1]);
        iotype = atoi(argv[2]);
    }

    /* Create a madrec object */
    madrecp = madrecCreate();
    printf("create: %s\n", madrecGetError(madrecp));

    /* Connect the madrec object to a madrigal file */
    madrecOpen(madrecp, iotype, madfile);
    printf("open: %s\n", madrecGetError(madrecp));

    /* Specify cedarGetParmArray arguments */
    fspec.nparms = 7;
    fspec.pparms = pparms;  
    fspec.nfilters = 4;
    fspec.fparms = fparms;
    fspec.fmin = fmin;
    fspec.fmax = fmax;


    /* Read the file */
    nrecs = 999999;
    for (j=0; j<nrecs; j++) {
        ind = madrecGetNextRec(madrecp);
        if (ind != 0) {
            if (ind == -1)
                break;
            printf("Error reading cedar file\n");
            return(0);
        }
        parp = cedarGetParmArray(madrecp->recordp, &fspec, &nlines);
        for (i=0; i<nlines; i++) {
            for (k=0; k<fspec.nparms; k++) {
                printf("%12.5e", parp[i+nlines*k]);
            }
            printf("\n");
        }
        free(parp);
    }

    /* Print the record */
    /*
    cedarPrintRecord(madrecp->recordp);
    printf("printRecord: %s\n", madrecGetError(madrecp));
    */

    /* Disconnect the madrigal file */
    madrecClose(madrecp);
    printf("close: %s\n", madrecGetError(madrecp));

    /* Destroy the madrec object */
    madrecDestroy(madrecp);

    return(0);
}