/home/kapf/tycho_docu/start_file_reader.c File Reference

#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <math.h>
#include "variables_global.h"
#include "prototypes.h"

Go to the source code of this file.

Functions
int	start_file_reader (char filename[])

---

Function Documentation

int start_file_reader

(

char

filename[]

)

The read routine for the parameterfile

Definition at line 19 of file start_file_reader.c.

References advection, boundary::back, bound, courant, dens_ic, dimension, boundary::down, endtime_sim, filetype, boundary::front, Gamma1, gasconstant, gravity_on_off, inflow_density, inflow_temperature, inflow_velocity, intial_velocity_file, boundary::left, make_ics, marker_density, marker_ic, minute_in_sec, number_of_threads, obst_ic, obstacle_density, obstacle_heat_conductivity, obstacle_temperature, obstalce_absorption_coefficient, output_dir, output_frequency, restart_frequency, boundary::right, sound_frequency, sound_ic, sound_pressure_level, sound_reflexion, specific_heat_capacity_gas, specific_heat_capacity_obstacle, starting_flow, strat_const_atmos, temp_ic, boundary::up, vel_ic, viscosity_on_off, wind_ic, wind_on_off, wind_speed, with_obstacles, with_one_pulse, with_sound_emitter, x, xmax, y, ymax, z, and zmax.

Referenced by main().

                                       {

    char data[55][200];
    int l;
    FILE *fd;
    char buffer[200];

    fd = fopen(filename, "r");
    if (fd == NULL) {
        printf("The Parameter File you specified does not exist.\n");
        exit(10);
    }

    l = 0;
    while (feof(fd) == 0) {
        fgets(data[l], 200, fd);
        if ((data[l][0] == '#') || (data[l][0] == '\n')) continue;
        data[l][strlen(data[l]) - 1] = '\0';
        l++;
    }
    fclose(fd);



    filetype = atoi(data[0]);

    make_ics = atoi(data[1]);

    sprintf(dens_ic, "%s", data[2]);

    sprintf(temp_ic, "%s", data[3]);

    sprintf(vel_ic, "%s", data[4]);

    sprintf(wind_ic, "%s", data[5]);

    sprintf(obst_ic, "%s", data[6]);
    
    sprintf(sound_ic, "%s", data[7]);

    sprintf(marker_ic, "%s", data[8]);

    intial_velocity_file = atoi(data[9]);

    sprintf(output_dir, "%s", data[10]);

    strat_const_atmos = atoi(data[11]);

    x = atoi(data[12]);
    y = atoi(data[13]);
    z = atoi(data[14]);

    dimension = atoi(data[15]);

    xmax = atof(data[16]);
    ymax = atof(data[17]);
    zmax = atof(data[18]);

    bound.down = atoi(data[19]);
    bound.up = atoi(data[20]);
    bound.left = atoi(data[21]);
    bound.right = atoi(data[22]);
    bound.front = atoi(data[23]);
    bound.back = atoi(data[24]);

    inflow_velocity = atof(data[25]);
    inflow_density = atof(data[26]);
    inflow_temperature = atof(data[27]);
    starting_flow = atoi(data[28]);

    endtime_sim = atof(data[29]);
    endtime_sim = endtime_sim*minute_in_sec;

    output_frequency = atof(data[30]);

    restart_frequency = atof(data[31]);

    courant = atof(data[32]);

    gravity_on_off = atoi(data[33]);

    gasconstant = atof(data[34]);

    Gamma1 = atof(data[35]);

    wind_on_off = atoi(data[36]);

    wind_speed = atof(data[37]);

    number_of_threads = atoi(data[38]);

    with_obstacles = atoi(data[39]);

    obstacle_density = atof(data[40]);

    obstacle_temperature = atof(data[41]);
    obstacle_temperature = obstacle_density * gasconstant * obstacle_temperature;

    obstacle_heat_conductivity = atof(data[42]);

    specific_heat_capacity_gas = atof(data[43]);

    specific_heat_capacity_obstacle = atof(data[44]);

    advection = atoi(data[45]);

    marker_density = atof(data[46]);

    viscosity_on_off = atoi(data[47]);
    
    with_sound_emitter = atoi(data[48]);

    sound_pressure_level = atof(data[49]);

    sound_pressure_level = pow(10, (sound_pressure_level / 20))*2.0E-5 * sqrt(2.0);
    
    with_one_pulse = atoi(data[50]);

    sound_frequency = atof(data[51]);
    
    sound_reflexion = atof(data[52]);

    obstalce_absorption_coefficient = atof(data[53]);
    obstalce_absorption_coefficient = pow(10,(obstalce_absorption_coefficient/20));
    obstalce_absorption_coefficient = 1.0/obstalce_absorption_coefficient;
    
    //printf("%g\n", obstalce_absorption_coefficient);

    //for debug
    //printf("%i\n", l);

    if (l != 54) {
        printf("Simulation stops, because of too few data in the start file.\n");
        exit(99);
    }

    return 0;
}