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

/*
modification history
--------------------
01a,1Feb02         Original
*/

/* 
*    Usage: testCreateRecord
*    This program
*        
*/

#include <stdlib.h>
#include <string.h>
#include <stdio.h>
#include <madrec.h>
#include <cedar.h>

int
main (argc, argv)
     int     argc;
     char    *argv[];
{
    Madrec *madrecp;    /* Output File - type 0-4 */
    int iotype=0, stat=0, record=0, i=0;
    int lprol=0, jpar=0, mpar=0, nrow=0, krec=0, kinst=0, kindat=0, ibyr=0, 
        ibmo=0, ibdy=0, ibh=0, ibm=0, ibs=0, ibcs=0, ieyr=0, iemo=0, iedy=0, 
        ieh=0, iem=0, ies=0, iecs=0;
    double kp[8] = {2.0, 3.3, 4.0, 3.0, 8.0, 7.7, 6.3, 5.0};
    double ap[8] = {10.0, 15.0, 16.0, 15.0, 24.0, 21.0, 18.0, 15.0};
    char text[161] = "";


    /* Get the input parameters */
/*
    if (argc == 4) {
        infile = argv[1];
        outfile = argv[2];
        iotype = atoi(argv[3]);
        if (iotype < 0 || iotype > 4) {
            fprintf(stderr, "%s\n", err3);
            return(1);
        } else {
            iotype = 20 + iotype;
        }
    } else {
        fprintf(stderr, "%s\n", err2);
        return(1);
    }
*/
  

    /* Create a madrec object for the output file */
    if ((madrecp = madrecCreate()) == (Madrec *) NULL) {
        fprintf(stderr, "create madrecw: %s\n", madrecGetError(madrecp));
        return(1);
    }

    /* Connect the output madrec object to a madrigal file */
    iotype = 20;
    stat = madrecOpen(madrecp, iotype, "madout");
    fprintf(stderr, "open madrecw: %s\n", madrecGetError(madrecp));
    
    /* first, add a catalog record - lenght must be 80*n */
    /* create some text */
    text[160] = '\0';
    for (i=0; i<160; i++)
        text[i] = ' ';
    strcpy(text, "This is a catalog record line 0");
    text[strlen("This is a catalog record line 0")] = ' ';
    strcpy(text + 80, "This is a catalog record line 1");
    text[80 + strlen("This is a catalog record line 1")] = ' ';
    madrecp->recordp = cedarCreateCatalogRecord(31, 30007,
                                2001, 8, 20,
                                0, 0, 0, 0,
                                2001, 8, 21,
                                23, 59, 59, 99,
				text);
				
    /* append some more text to the catalog record */
    for (i=0; i<160; i++)
        text[i] = ' ';
    strcpy(text, "This is a catalog record line 2");
    text[strlen("This is a catalog record line 2")] = ' ';
    strcpy(text + 80, "This is a catalog record line 3");
    text[80 + strlen("This is a catalog record line 3")] = ' ';
    stat = cedarAppendCatalogRecord(&(madrecp->recordp), text);
				
    stat = madrecPutNextRec(madrecp);
    fprintf(stderr, "putNextRec madrecw: %s\n", madrecGetError(madrecp));
    free(madrecp->recordp);
    
    /* next, add a header record - lenght must be 80*n */
    /* create some text */
    for (i=0; i<160; i++)
        text[i] = ' ';
    strcpy(text, "This is a header record line 0");
    text[strlen("This is a header record line 0")] = ' ';
    strcpy(text + 80, "This is a header record line 1");
    text[80 + strlen("This is a header record line 1")] = ' ';
    madrecp->recordp = cedarCreateHeaderRecord(31, 30007,
                                2001, 8, 20,
                                0, 0, 0, 0,
                                2001, 8, 21,
                                23, 59, 59, 99,
				2,2,
				text);
				
    /* append some more text to the header record */
    for (i=0; i<160; i++)
        text[i] = ' ';
    strcpy(text, "This is a header record line 2");
    text[strlen("This is a header record line 2")] = ' ';
    strcpy(text + 80, "This is a header record line 3");
    text[80 + strlen("This is a header record line 3")] = ' ';
    stat = cedarAppendHeaderRecord(&(madrecp->recordp), text);
				
    stat = madrecPutNextRec(madrecp);
    fprintf(stderr, "putNextRec madrecw: %s\n", madrecGetError(madrecp));
	
    lprol = 16;
    jpar = 3;
    mpar = 2;
    nrow = 8;
    krec = 1002;
    kinst = 210;
    kindat = 30007;
    ibyr = 2001;
    ibmo = 8;
    ibdy = 20;
    ibh = 0;
    ibm = 0;
    ibs = 0;
    ibcs = 0;
    ieyr = 2001;
    iemo = 8;
    iedy = 20;
    ieh = 23;
    iem = 59;
    ies = 59;
    iecs = 59;
    for (record=0; record<5; record++) {

        ibdy++;
        iedy++;

	/* Create a Cedar record in the madrec object */
        if (madrecp->recordp != (Int16 *)NULL) {
            free(madrecp->recordp);
        }
	madrecp->recordp = cedarCreateRecord(lprol, jpar, mpar, nrow, krec,
					     kinst, kindat, ibyr, ibmo, ibdy, 
					     ibh, ibm, ibs, ibcs, ieyr,
					     iemo, iedy, ieh, iem, ies,
					     iecs);
    
	/* Set 1d parameters */
	stat = cedarSet1dParm(madrecp->recordp, 340, 12.0, 0);
	stat = cedarSet1dParm(madrecp->recordp, 354, 1.5e-20, 1);
	stat = cedarSet1dParm(madrecp->recordp, 356, 1.2e-20, 2);
    
	/* Set 2d parms */
	stat = cedarSet2dParm(madrecp->recordp, 310, kp, 0);
	stat = cedarSet2dParm(madrecp->recordp, 335, ap, 1);
    
	/* cedarPrintRecord(madrecp->recordp); */
    
	stat = madrecPutNextRec(madrecp);
	fprintf(stderr, "putNextRec madrecw: %s\n", madrecGetError(madrecp));

    }

    stat = madrecClose(madrecp);      
    fprintf(stderr, "close madrecw: %s\n", madrecGetError(madrecp));

    madrecDestroy(madrecp);

    return(0);
}