MtkFileGridFieldToDataType.c

Go to the documentation of this file.

/*===========================================================================
=                                                                           =
=                        MtkFileGridFieldToDataType                         =
=                                                                           =
=============================================================================

                         Jet Propulsion Laboratory
                                   MISR
                               MISR Toolkit

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

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

#include "MisrFileQuery.h"
#include "MisrUtil.h"
#include "MisrError.h"
#include <stdlib.h>
#include <mfhdf.h>
#include <HdfEosDef.h>

MTKt_status MtkFileGridFieldToDataType(
  const char *filename, 
  const char *gridname, 
  const char *fieldname, 
  MTKt_DataType *datatype  )
{
  MTKt_status status;       /* Return status */

  status = MtkFileGridFieldToDataTypeNC(filename, gridname, fieldname, datatype); // try netCDF
  if (status != MTK_NETCDF_OPEN_FAILED) return status;

  return MtkFileGridFieldToDataTypeHDF(filename, gridname, fieldname, datatype); // try HDF
}

MTKt_status MtkFileGridFieldToDataTypeNC(
  const char *filename, 
  const char *gridname, 
  const char *fieldname, 
  MTKt_DataType *datatype  )
{
  MTKt_status status;
  MTKt_status status_code;
  int ncid = 0;

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

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

  /* Read data type for the field. */
  status = MtkFileGridFieldToDataTypeNcid(ncid, gridname, fieldname, datatype);
  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 MtkFileGridFieldToDataTypeHDF(
  const char *filename, 
  const char *gridname, 
  const char *fieldname, 
  MTKt_DataType *datatype  )
{
  MTKt_status status_code;      /* Return status of this function */
  intn status;          /* HDF-EOS return status */
  int32 fid = FAIL;     /* HDF-EOS file identifier */

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

  fid = GDopen((char *)filename, DFACC_READ);
  if (fid == FAIL) MTK_ERR_CODE_JUMP(MTK_HDFEOS_GDOPEN_FAILED);

  /* Read data type for the field. */
  status = MtkFileGridFieldToDataTypeFid(fid, gridname, fieldname, datatype);
  MTK_ERR_COND_JUMP(status);

  status = GDclose(fid);
  if (status == FAIL) MTK_ERR_CODE_JUMP(MTK_HDFEOS_GDCLOSE_FAILED);

  return MTK_SUCCESS;

ERROR_HANDLE:
  if (fid != FAIL) GDclose(fid);
  return status_code;
}

MTKt_status MtkFileGridFieldToDataTypeFid(
  int32 fid,            
  const char *gridname, 
  const char *fieldname, 
  MTKt_DataType *datatype  )
{
  MTKt_status status_code;      /* Return status of this function */
  intn status;          /* HDF-EOS return status */
  int32 gid = FAIL;     /* HDF-EOS grid identifier */
  int32 hdf_datatype;       /* HDF-EOS datatype number */
  int32 rank;           /* Not used */
  int32 dims[10];       /* Not used */
  char dimlist[80];     /* Not used */
  char *basefield = NULL;   /* Base fieldname */
  int nextradims;               /* Number of extra dimensions */
  int *extradims = NULL;    /* Extra dimension list */

  if (gridname == NULL ||
      fieldname == NULL || datatype == NULL)
    MTK_ERR_CODE_JUMP(MTK_NULLPTR);

  gid = GDattach(fid, (char *)gridname);
  if (gid == FAIL) MTK_ERR_CODE_JUMP(MTK_HDFEOS_GDATTACH_FAILED);

  status = MtkParseFieldname(fieldname, &basefield, &nextradims, &extradims);
  MTK_ERR_COND_JUMP(status);

  status = GDfieldinfo(gid, basefield, &rank, dims, &hdf_datatype, dimlist);
  if (status == FAIL) MTK_ERR_CODE_JUMP(MTK_HDFEOS_GDFIELDINFO_FAILED);

  status = GDdetach(gid);
  if (status == FAIL) MTK_ERR_CODE_JUMP(MTK_HDFEOS_GDDETACH_FAILED);

  status = MtkHdfToMtkDataTypeConvert(hdf_datatype, datatype);
  MTK_ERR_COND_JUMP(status);

  free(basefield);
  free(extradims);
  return MTK_SUCCESS;

ERROR_HANDLE:
  if (gid != FAIL) GDdetach(gid);
  if (basefield != NULL) free(basefield);
  if (extradims != NULL) free(extradims);
  return status_code;
}

MTKt_status MtkFileGridFieldToDataTypeNcid(
  int ncid,            
  const char *gridname, 
  const char *fieldname, 
  MTKt_DataType *datatype  )
{
  MTKt_status status_code;      /* Return status of this function */
  intn status;          /* HDF-EOS return status */
  char *basefield = NULL;   /* Base fieldname */
  int nextradims;               /* Number of extra dimensions */
  int *extradims = NULL;    /* Extra dimension list */

  if (gridname == NULL ||
      fieldname == NULL || datatype == NULL)
    MTK_ERR_CODE_JUMP(MTK_NULLPTR);

  status = MtkParseFieldname(fieldname, &basefield, &nextradims, &extradims);
  MTK_ERR_COND_JUMP(status);

  int group_id;
  {
    int nc_status = nc_inq_grp_ncid(ncid, gridname, &group_id);
    if (nc_status != NC_NOERR) MTK_ERR_CODE_JUMP(MTK_NETCDF_READ_FAILED);
  }

  MTKt_ncvarid var;
  status = MtkNCVarId(group_id, basefield, &var);
  MTK_ERR_COND_JUMP(status);

  nc_type nc_datatype;
  {
    int nc_status = nc_inq_vartype(var.gid, var.varid, &nc_datatype);
    if (nc_status != NC_NOERR) MTK_ERR_CODE_JUMP(MTK_NETCDF_READ_FAILED);
  }

  status = MtkNcToMtkDataTypeConvert(nc_datatype, datatype);
  MTK_ERR_COND_JUMP(status);

  free(basefield);
  free(extradims);
  return MTK_SUCCESS;

ERROR_HANDLE:
  if (basefield != NULL) free(basefield);
  if (extradims != NULL) free(extradims);
  return status_code;
}