afrt_phs.cpp

Go to the documentation of this file.

/**************************************************************************
*
* NAME: afrt_phs.cpp
*
* DESCRIPTION: Source file to run Ahmad Frazier Top Of Atmosphere calculation
*
* REFERENCES:
*
* REVISION HISTORY:
*   DATE:       PR#         AUTHOR         Description
*   --------    ------      --------       -----------------
*   05-02-2018              S. Anderson
*
* NOTES (MISCELLANEOUS) SECTION:
* none
*
**************************************************************************/
 
#include <genutils.h>
#include <AfrtOptions.h>
#include <AfrtProcess.h>
#include <afrt_input.h>
 
//-----------------------------------------------------------------------------
//
// Main Function
//
//-----------------------------------------------------------------------------
 
int main(int argc, char* argv[])
{
    int status = 0;
 
    clo_optionList_t* list = clo_createList();
    clo_setEnablePositionOptions(1);
 
    //char softwareVersion[200];
    //sprintf(softwareVersion, "%d.%d.%d-r%d", L3MAPGEN_VERSION_MAJOR, L3MAPGEN_VERSION_MINOR,
    //        L3MAPGEN_VERSION_PATCH_LEVEL, SVN_REVISION);
    afrt_init_options(list, "2.0");
    if(argc == 1) {
        clo_printUsage(list);
        exit(1);
    }
    afrt_read_options(list, argc, argv);
 
    if(clo_getBool(list, "verbose"))
        want_verbose = 1;
    else
        want_verbose = 0;
 
    AfrtProcess* generator = new AfrtProcessPhs();
    generator->setHistory(afrt_get_history(argc, argv));
    
    status = generator->initialize();
 
    if ( status != RT_SUCCESS ) {
        return(status);
    }
 
    status = generator->process();
 
    if ( status != RT_SUCCESS ) {
        return(status);
    }
 
    string ofile = afrt_get_option(OFILE_PHS);
    NcFile* nc_output = generator->open_nc4(ofile);
    if (nc_output == 0L) {
        std::cerr << "Failed to open NetCDF output file" << std::endl;
        status = RT_FAIL;
    } else {
        status = generator->write_nc4(nc_output);
        delete nc_output;
    }
    delete generator;
 
    if ( status != RT_SUCCESS ) {
        std::cerr << "PHS NetCDF tables NOT created" << std::endl;
    } else {
        std::cerr << "\nPHS NetCDF tables created. \n" << std::endl;
    }
 
   return(status);
}