MtkFileToPath.c

Go to the documentation of this file.

/*===========================================================================
=                                                                           =
=                              MtkFileToPath                                =
=                                                                           =
=============================================================================

                         Jet Propulsion Laboratory
                                   MISR
                               MISR Toolkit

            Copyright 2005, California Institute of Technology.
                           ALL RIGHTS RESERVED.
                 U.S. Government Sponsorship acknowledged.

============================================================================*/

#include "MisrFileQuery.h"
#include "MisrError.h"
#include <hdf.h>
#include <mfhdf.h>

MTKt_status MtkFileToPath( const char *filename, 
               int *path              )
{
  MTKt_status status;       /* Return status */

  status = MtkFileToPathNC(filename, path); // try netCDF
  if (status != MTK_NETCDF_OPEN_FAILED) return status;

  return MtkFileToPathHDF(filename, path); // try HDF
}

MTKt_status MtkFileToPathNC( const char *filename, 
               int *path              )
{
  MTKt_status status_code;      /* Return code of this function */
  MTKt_status status;           /* Return status */

  if (filename == NULL) return MTK_NULLPTR;

  /* Open file */
  int ncid = 0;
  {
    int nc_status = nc_open(filename, NC_NOWRITE, &ncid);
    if (nc_status != NC_NOERR) MTK_ERR_CODE_JUMP(MTK_NETCDF_OPEN_FAILED);
  }

  /* Read grid attribute */
  status = MtkFileToPathNcid(ncid, path);
  MTK_ERR_COND_JUMP(status);

  /* Close file */
  {
    int nc_status = nc_close(ncid);
    if (nc_status != NC_NOERR) MTK_ERR_CODE_JUMP(MTK_NETCDF_CLOSE_FAILED);
  }
  ncid = 0;

  return MTK_SUCCESS;

 ERROR_HANDLE:
  if (ncid != 0) nc_close(ncid);
  return status_code;
}

MTKt_status MtkFileToPathHDF( const char *filename, 
               int *path              )
{
  MTKt_status status_code;      /* Return status of this function */
  int32 sid = FAIL;     /* HDF SD file identifier */
  MTKt_status status;       /* Return status */
  intn hdfstatus;       /* HDF return status */

  if (filename == NULL)
    MTK_ERR_CODE_JUMP(MTK_NULLPTR);

  sid = SDstart(filename, DFACC_READ);
  if (sid == FAIL)
    MTK_ERR_CODE_JUMP(MTK_HDF_SDSTART_FAILED);

  /* Read path number. */
  status = MtkFileToPathFid(sid, path);
  MTK_ERR_COND_JUMP(status);

  hdfstatus = SDend(sid);
  if (hdfstatus == FAIL)
    MTK_ERR_CODE_JUMP(MTK_HDF_SDEND_FAILED);
  sid = FAIL;

  return MTK_SUCCESS;

ERROR_HANDLE:
  if (sid != FAIL) {
    hdfstatus = SDend(sid);
  }
  return status_code;
}

MTKt_status MtkFileToPathFid( int32 sid,            
                  int *path              )
{
  MTKt_status status_code;      /* Return status of this function */
  intn status;          /* HDF return status */
  int32 attr_index;     /* HDF attribute index */
  int32 path_tmp;       /* Temp path number */

  if (path == NULL)
    MTK_ERR_CODE_JUMP(MTK_NULLPTR);

  attr_index = SDfindattr(sid, "Path_number");
  if (attr_index == FAIL)
    MTK_ERR_CODE_JUMP(MTK_HDF_SDFINDATTR_FAILED);

  status = SDreadattr(sid, attr_index, &path_tmp);
  if (status == FAIL)
    MTK_ERR_CODE_JUMP(MTK_HDF_SDREADATTR_FAILED);

  *path = path_tmp;

  return MTK_SUCCESS;

ERROR_HANDLE:
  return status_code;
}

MTKt_status MtkFileToPathNcid( int ncid,            
                  int *path              )
{
  MTKt_status status_code;      /* Return status of this function */

  if (path == NULL)
    MTK_ERR_CODE_JUMP(MTK_NULLPTR);

  int path_tmp;
  {
    int nc_status = nc_get_att(ncid, NC_GLOBAL, "Path_number", &path_tmp);
    if (nc_status != NC_NOERR) MTK_ERR_CODE_JUMP(MTK_NETCDF_READ_FAILED);
  }

  *path = path_tmp;

  return MTK_SUCCESS;

ERROR_HANDLE:
  return status_code;
}