GEO_get_ephatt_inputs.c

Go to the documentation of this file.

#include <errno.h>
#include <stdio.h>
#include "PGS_IO.h"
#include "smfio.h"
#include "GEO_earth.h"
#include "GEO_product.h"
#include "PGS_MODIS_35251.h"
 
/* Used internally by the EPH portion of the SDP Toolkit library. */
extern int byteswap(char *, int);
 
PGSt_SMF_status GEO_get_ephatt_inputs(
    PGSt_PC_Logical file_logical,
    PGSt_integer    file_version,
    char        universal_references[][PGSd_UR_FIELD_SIZE]
){
    PGSt_IO_Gen_FileHandle  *handle;
    PGSt_ephemHeader        fileHeader;
    PGSt_SMF_status     retval = PGS_S_SUCCESS;
    char            filefunc[] = __FILE__ ", GEO_get_ephatt_inputs";
    char            msgbuf[256];
 
    if(universal_references == NULL)
    {
    modsmf(MODIS_E_BAD_INPUT_ARG, "universal_references is null",
        filefunc);
    
    return MODIS_E_BAD_INPUT_ARG;
    }
 
    if(PGS_IO_Gen_Open(file_logical, PGSd_IO_Gen_Read, &handle, file_version)
    != PGS_S_SUCCESS)
    {
    sprintf(msgbuf, "PGS_IO_Gen_Open(%ld, %ld)",
        (long)file_logical, (long)file_version);
    modsmf(MODIS_E_GEO, msgbuf, filefunc);
    
    return MODIS_E_GEO;
    }
 
    errno = 0;
    if(fread(&fileHeader, sizeof fileHeader, 1, handle) != 1)
    {
    if(errno)
        PGS_SMF_SetUNIXMsg(errno, "during fread(fileHeader)", filefunc);
    else
        modsmf(MODIS_E_GEO, "fread(fileHeader)", filefunc);
 
    retval = PGS_E_UNIX;
    }
    else
    {
    /* Simplistic byte-order check, borrowed from PGS_EPH_GetEphMet() */
    if(fileHeader.nURs > 100)
        byteswap((char*)&fileHeader.nURs, sizeof fileHeader.nURs);
 
    if(fileHeader.nURs > MAX_EA_INPUTS)
        fileHeader.nURs = MAX_EA_INPUTS;
 
    if(fread(universal_references, PGSd_UR_FIELD_SIZE, fileHeader.nURs,
        handle) != fileHeader.nURs)
    {
        if(errno)
        PGS_SMF_SetUNIXMsg(errno, "during fread(universal_references",
            filefunc);
        else
        modsmf(MODIS_E_GEO, "fread(universal_references)", filefunc);
 
        retval = PGS_E_UNIX;
    }
    else if(fileHeader.nURs < MAX_EA_INPUTS)
        /*Terminate the array with an empty string. */
        universal_references[fileHeader.nURs][0] = '\0';
    }
 
    if(PGS_IO_Gen_Close(handle) != PGS_S_SUCCESS)
    modsmf(MODIS_E_GEO, "PGE_IO_Gen_Close()", filefunc);
 
    return retval;
}