MacroscopicProblem.C

#include <stdio.h>
#include <iostream>
#include "modena.h"

using namespace std;

int
main(int argc, char *argv[])
{
    double T = 270;
    double Tend = 290.0;
    
    // Instantiate index set
    //modena_index_set_t *indexSet = modena_index_set_new("species");


    // Instantiate a model
    modena_model_t *model = modena_model_new("SurfaceTension[A=AIR,B=THF]");   //muss das FunctionModule genau so heißen??
    if(modena_error_occurred())
    {
        return modena_error();
    }

    // Allocate memory and fetch arg positions
    modena_inputs_t *inputs = modena_inputs_new(model);           //How many inputs and outputs is defined in the function module!! 
    modena_outputs_t *outputs = modena_outputs_new(model);

    
    size_t Tpos = modena_model_inputs_argPos(model, "T");

    modena_model_argPos_check(model);

    
    while(T < Tend)
    {
      // Set input vector
      modena_inputs_set(inputs, Tpos, T);

      // Call the model
     int ret = modena_model_call(model, inputs, outputs);

     // Terminate, if requested
     if(modena_error_occurred())
            {
                modena_inputs_destroy(inputs);
                modena_outputs_destroy(outputs);
                modena_model_destroy(model);

                return modena_error();
            }

     // Fetch result
     double rho = modena_outputs_get(outputs, 0);
     
     cout << "T = " << T;
     
     
     T = T + 10.0;
    }    
    
       
    modena_inputs_destroy(inputs);
    modena_outputs_destroy(outputs);
    modena_model_destroy(model);

    return 0;
}