/* $Id: StartStopFromData.c,v 1.1 2010/08/31 10:32:55 budnik Exp $ */ /** @file StartStopFromData.c * @brief */ #include #include #include #include #include #include #define MAX_FILE_NAME_LEN 120 #define TIMESHIFT 62167219200.0 /************************************* Global variables and structures **************************************/ long CDFDims, // Number of dimensions in a CDF file CDFVars, // Number of variables in a CDF file CDFDimSizes[CDF_MAX_DIMS], // Dimension Sizes in a CDF file CDFencoding, // Data encoding CDFmajority, // Variable majority CDFmaxRec, // max Record number CDFAttrs; // number of CDF Attributes // CDFdatatype; struct cdfvar { // CDF variable structure char name[CDF_VAR_NAME_LEN+1]; long num; //variable number long datatype; long numElem; // variable dimensionality long recVariance; // variable rec Variance long dimVariances[CDF_MAX_DIMS]; // # of data values in dimSizes } cdfVar; /*-------------------------------------------------------------------------- * Handles a CDF error. *--------------------------------------------------------------------------*/ void cdf_handle_error(CDFstatus status) { char message[CDF_STATUSTEXT_LEN+1]; CDFerror (status, message); /* Get the appropriate message */ fprintf (stderr, "%s\n", message); exit(1); } int main(int argc, char *argv[]) { char fileName[MAX_FILE_NAME_LEN]; long RecStart = 0, RecCount, RecInt = 1; long indices[1] = {0}, intervals[1] = {1}, counts[1] = {1}; long varN = -1; long datatype = 0; CDFstatus cstatus; // CDF status code CDFid id; double *value; int i; double epoch16[2]; long year, month, day, hour, min, sec, msec, mksec, nsec, psec; if (argc <= 1) exit(0); /* CDF input file name not specified */ else strcpy(fileName, argv[argc-1]); /* Get the input file name */ /*********************** Open CDF file *****************************/ if ((cstatus = CDFopen(fileName, &id)) != CDF_OK) cdf_handle_error(cstatus); /*********** treat all vars as zVars with eliminated false dimensionality and get Nrecs **********/ if ((cstatus = CDFlib(SELECT_, CDF_zMODE_, zMODEon2, NULL_)) != CDF_OK) cdf_handle_error (cstatus); while (datatype != CDF_EPOCH && datatype != CDF_EPOCH16){ varN++; if ((cstatus = CDFlib(SELECT_, zVAR_, varN, GET_, zVAR_DATATYPE_, &datatype, NULL_)) != CDF_OK) cdf_handle_error (cstatus); // if (datatype == CDF_EPOCH16) break; } if ((cstatus = CDFlib(SELECT_, zVAR_, varN, GET_, zVAR_MAXREC_, &CDFmaxRec, NULL_)) != CDF_OK) cdf_handle_error (cstatus); // fprintf(stdout,"Max Recs %d\n", CDFmaxRec+1); if (CDFmaxRec < 1) exit (0); value = (double *)malloc(sizeof(double)* (CDFmaxRec+1)*((datatype == CDF_EPOCH16)+1)); /******************* Get Epoch *************************/ if ((cstatus = CDFlib (SELECT_, zVAR_, varN, zVAR_RECNUMBER_, RecStart, zVAR_RECCOUNT_, CDFmaxRec+1, zVAR_RECINTERVAL_, RecInt, zVAR_DIMINDICES_, indices, zVAR_DIMCOUNTS_, counts, zVAR_DIMINTERVALS_, intervals, GET_, zVAR_HYPERDATA_, value, NULL_)) != CDF_OK) cdf_handle_error (cstatus); if (datatype == CDF_EPOCH16) { epoch16[0] = value[0]; epoch16[1] = value[1]; EPOCH16breakdown(epoch16, &year, &month, &day, &hour, &min, &sec, &msec, &mksec, &nsec, &psec); value[0] = computeEPOCH(year, month, day, hour, min, sec, msec); epoch16[0] = value[CDFmaxRec*2]; epoch16[1] = value[CDFmaxRec*2+1]; EPOCH16breakdown(epoch16, &year, &month, &day, &hour, &min, &sec, &msec, &mksec, &nsec, &psec); value[CDFmaxRec] = computeEPOCH(year, month, day, hour, min, sec, msec); } printf("%f %f\n", value[0]/1000.0-TIMESHIFT, value[CDFmaxRec]/1000.0-TIMESHIFT); if ((cstatus = CDFlib(CLOSE_, CDF_, NULL_)) != CDF_OK) cdf_handle_error (cstatus); free(value); }