int DBObjData::Read(const char *fileName) {
DBInt ret;
FILE *file;
if (strncmp(CMfileExtension(fileName), "nc", 2) == 0) {
Type(DBTypeGridContinuous); // TODO: Limiting to Continuous grid
ret = DBImportNetCDF(this, fileName);
}
else if (strncmp(CMfileExtension(fileName), "gz", 2) == 0) {
char pCommand[strlen(fileName) + 16];
sprintf(pCommand, "gunzip -c %s", fileName);
if ((file = popen(pCommand, "r")) == (FILE *) NULL) {
CMmsgPrint(CMmsgSysError, "File (%s) Opening Error in: %s %d", fileName, __FILE__, __LINE__);
return (DBFault);
}
ret = Read(file);
pclose(file);
}
else {
if ((file = fopen(fileName, "r")) == (FILE *) NULL) {
CMmsgPrint(CMmsgAppError, "File (%s) Opening Error in: %s %d", fileName, __FILE__, __LINE__);
return (DBFault);
}
ret = Read(file);
fclose(file);
}
FileName(fileName);
return (ret);
}
int DBObjData::Read (const char *fileName)
{
FILE *file;
if (strncmp(CMfileExtension(fileName),"nc",2) == 0)
{
Type (DBTypeGridContinuous); // TODO: Limiting to Continuous grid
if (DBImportNetCDF (this,fileName) == DBFault) return (DBFault);
}
else
{
if ((file = fopen (fileName,"r")) == (FILE *) NULL)
{
CMmsgPrint (CMmsgAppError, "File (%s) Opening Error in: %s %d",fileName, __FILE__,__LINE__);
return (DBFault);
}
if (Read (file) == DBFault) return (DBFault);
}
FileName (fileName);
return (DBSuccess);
}
int main (int argc,char *argv [])
{
int argPos, argNum = argc, ret, verbose = false;
char *title = (char *) NULL, *subject = (char *) NULL;
char *domain = (char *) NULL, *version = (char *) NULL;
DBInt shadeSet = DBFault;
DBObjData *grdData;
for (argPos = 1;argPos < argNum; )
{
if (CMargTest (argv [argPos],"-t","--title"))
{
if ((argNum = CMargShiftLeft (argPos,argv,argNum)) <= argPos)
{ CMmsgPrint (CMmsgUsrError,"Missing title!"); return (CMfailed); }
title = argv [argPos];
if ((argNum = CMargShiftLeft (argPos,argv,argNum)) <= argPos) break;
continue;
}
if (CMargTest (argv [argPos],"-u","--subject"))
{
if ((argNum = CMargShiftLeft (argPos,argv,argNum)) <= argPos)
{ CMmsgPrint (CMmsgUsrError,"Missing subject!"); return (CMfailed); }
subject = argv [argPos];
if ((argNum = CMargShiftLeft (argPos,argv,argNum)) <= argPos) break;
continue;
}
if (CMargTest (argv [argPos],"-d","--domain"))
{
if ((argNum = CMargShiftLeft (argPos,argv,argNum)) <= argPos)
{ CMmsgPrint (CMmsgUsrError,"Missing domain!"); return (CMfailed); }
domain = argv [argPos];
if ((argNum = CMargShiftLeft (argPos,argv,argNum)) <= argPos) break;
continue;
}
if (CMargTest (argv [argPos],"-v","--version"))
{
if ((argNum = CMargShiftLeft (argPos,argv,argNum)) <= argPos)
{ CMmsgPrint (CMmsgUsrError,"Missing version!"); return (CMfailed); }
version = argv [argPos];
if ((argNum = CMargShiftLeft (argPos,argv,argNum)) <= argPos) break;
continue;
}
if (CMargTest (argv [argPos],"-s","--shadeset"))
{
int shadeCodes [] = { DBDataFlagDispModeContStandard,
DBDataFlagDispModeContGreyScale,
DBDataFlagDispModeContBlueScale,
DBDataFlagDispModeContBlueRed,
DBDataFlagDispModeContElevation };
const char *shadeSets [] = { "standard","grey","blue","blue-to-red","elevation", (char *) NULL };
if ((argNum = CMargShiftLeft (argPos,argv,argNum)) <= argPos)
{ CMmsgPrint (CMmsgUsrError,"Missing shadeset!"); return (CMfailed); }
if ((shadeSet = CMoptLookup (shadeSets,argv [argPos],true)) == DBFault)
{ CMmsgPrint (CMmsgUsrError,"Invalid shadeset!"); return (CMfailed); }
shadeSet = shadeCodes [shadeSet];
if ((argNum = CMargShiftLeft (argPos,argv,argNum)) <= argPos) break;
continue;
}
if (CMargTest (argv [argPos],"-V","--verbose"))
{
verbose = true;
if ((argNum = CMargShiftLeft (argPos,argv,argNum)) <= argPos) break;
continue;
}
if (CMargTest (argv [argPos],"-h","--help"))
{
CMmsgPrint (CMmsgInfo,"%s [options] <input netcdf> <output data>",CMfileName(argv[0]));
CMmsgPrint (CMmsgInfo," -t,--title [dataset title]");
CMmsgPrint (CMmsgInfo," -u,--subject [subject]");
CMmsgPrint (CMmsgInfo," -d,--domain [domain]");
CMmsgPrint (CMmsgInfo," -v,--version [version]");
CMmsgPrint (CMmsgInfo," -s,--shadeset [standard|grey|blue|blue-to-red|elevation]");
CMmsgPrint (CMmsgInfo," -V,--verbose");
CMmsgPrint (CMmsgInfo," -h,--help");
return (DBSuccess);
}
if ((argv [argPos][0] == '-') && (strlen (argv [argPos]) > 1))
{ CMmsgPrint (CMmsgUsrError,"Unknown option: %s!",argv [argPos]); return (CMfailed); }
argPos++;
}
if (argNum < 2) { CMmsgPrint (CMmsgUsrError,"Missing input file"); return (CMfailed); }
if (argNum > 3) { CMmsgPrint (CMmsgUsrError,"Extra arguments!"); return (CMfailed); }
grdData = new DBObjData ("Noname",DBTypeGridContinuous);
ret = DBImportNetCDF (grdData,argv[1]);
if (title != (char *) NULL) grdData->Name (title);
if (subject != (char *) NULL) grdData->Document (DBDocSubject,subject);
if (domain != (char *) NULL) grdData->Document (DBDocGeoDomain,domain);
if (version != (char *) NULL) grdData->Document (DBDocVersion,version);
if (shadeSet != DBFault) { grdData->Flags (DBDataFlagDispModeContShadeSets,DBClear); grdData->Flags (shadeSet,DBSet); }
ret = (argNum > 2) && (strcmp (argv [2],"-") != 0) ? grdData->Write (argv [2]) : grdData->Write (stdout);
delete grdData;
return (ret);
}
