/***********************************************************************
/
/  GRID CLASS (IDENTIFY THE MULTI-SPECIES FIELDS)
/
/  written by: Chia-Jung Hsu, Patch generated by naunet
/  date:       2021
/
************************************************************************/

#include "preincludes.h"
#include "macros_and_parameters.h"
#include "typedefs.h"
#include "global_data.h"
#include "Fluxes.h"
#include "GridList.h"
#include "ExternalBoundary.h"
#include "Grid.h"

/* function prototypes */

int FindField(int f, int farray[], int n);

#ifdef USE_NAUNET

int grid::IdentifyNaunetSpeciesFields(
    {{ arguments | wordwrap(80) | replace("\n", "\n    ") }})
{
    {{ initial | wordwrap(80) | replace("\n", "\n    ") }};

    // Find Fields for the Naunet species
    {% for f in findfield -%}
        {{ f }}
    {% endfor -%}

    // check fields exist
    {% for c in check -%}
        {{ c }};
    {% endfor -%}

    return SUCCESS;
}

#endif
