afrt_nc4.cpp

Go to the documentation of this file.

/**************************************************************************
*
* NAME: afrt_nc4.cpp
*
* DESCRIPTION: Source file to run Ahmad Frazier Radiative Transfer 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.1");
    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;
 
    AfrtProcessNC4* generator = new AfrtProcessNC4();
    generator->setHistory(afrt_get_history(argc, argv));
 
    status = generator->initialize();
 
    if ( status != RT_SUCCESS ) {
        std::cerr << "ERROR!!! Initialization failure." << std::endl;
        return(status);
    }
 
#ifdef WITH_MPI
    status = generator->process_mpi();
#else
    status = generator->process();
#endif
 
    if ( status != RT_SUCCESS ) {
        std::cerr << "ERROR!!! Process failure." << std::endl;
        return(status);
    }
 
    delete generator;
 
    if ( status != RT_SUCCESS ) {
        std::cerr << "WARNING!!! RT tables NOT created" << std::endl;
    } else {
        std::cerr << "RT2 NetCDF4 tables created: " << std::endl;
    }
 
   return(status);
}