input.h

Go to the documentation of this file.

#ifndef _INPUT_H
#define _INPUT_H

#include <iostream>
#include <string>
#include <cstring>

#include "constants.h"

using namespace std;

class Input {
 public:
  Input() { } 
  Input(char** argv, int argc) :
    Ekfact(atof(argv[2])),
    numr((atoi(argv[3])%2 == 1) ? atoi(argv[3]): atoi(argv[3])+1),
    numz((atoi(argv[4])%2 == 1) ? atoi(argv[4]): atoi(argv[4])+1),
    D0(atof(argv[5])*1.e28/kpc/kpc*Myr),
    zt(atof(argv[6])),
    zmax(atof(argv[7])),
    delta(atof(argv[8])),
    index_radial(atof(argv[9])),
#ifdef HAVE_DS
      mx(atof(argv[10])),
      taudec(atof(argv[11])),
#else
    ab_C(atof(argv[10])),
    ab_N(atof(argv[11])),
#endif
#ifdef REAC
    vAlfven(atof(argv[12])*km*Myr/kpc),
#endif
    v0(atof(argv[13])*km*Myr/kpc),
    dvdz(atof(argv[14])*km*Myr/kpc),
    etaT(atof(argv[15])),
#ifdef HAVE_DS
    sigmav(atof(argv[16])),
    dmmode(atoi(argv[17])),
    dmprof(DMprofile(atoi(argv[18]))),
#endif
    filename(argv[1])
      { 

        if (!strcmp(argv[argc-1], "const")) set_profile = Constant;
        else if (!strcmp(argv[argc-1], "exp")) set_profile = Exp;
        else if (!strcmp(argv[argc-1], "expr")) set_profile = Expr;
        else if (!strcmp(argv[argc-1], "radial")) set_profile = Radial;
        else if (!strcmp(argv[argc-1], "exp_radial")) set_profile = ExpRadial;
        else {
          cerr << "What DPerpType did you choose?" << endl;
          exit(-2);
        }
      }
  Input(char* run, double& Ekfact_, int& Nr, int& Nz, double& D0_, double& zt_, double& zmax_, double& delta_, double& index_radial_, double& ab_C_, double& ab_N_, double& vA_, double& vC_, double& dvdz_conv_, double etaT_, char* profile) :
    Ekfact(Ekfact_),
    numr((Nr%2 == 1) ? Nr: Nr+1),
    numz((Nz%2 == 1) ? Nz: Nz+1),
    D0(D0_*1.e28/kpc/kpc*Myr),
    zt(zt_),
    zmax(zmax_),
    delta(delta_),
    index_radial(index_radial_),
#ifndef HAVE_DS
    ab_C(ab_C_),
    ab_N(ab_N_),
#endif
#ifdef REAC
    vAlfven(vA_*km*Myr/kpc),
#endif
    v0(vC_*km*Myr/kpc),
    dvdz(dvdz_conv_*km*Myr/kpc),
    etaT(etaT_),
    filename(run)
      { 

        int ipos = filename.find(' ');
        if (ipos < string::npos) filename.erase(ipos);

        if (!strcmp(profile, "const")) set_profile = Constant;
        else if (!strcmp(profile, "exp")) set_profile = Exp;
        else if (!strcmp(profile, "expr")) set_profile = Expr;
        else if (!strcmp(profile, "radial")) set_profile = Radial;
        else if (!strcmp(profile, "exp_radial")) set_profile = ExpRadial;
        else {
          cerr << "What DPerpType did you choose?" << endl;
          exit(-2);
        }
      }
    virtual ~Input() { }
    double Ekfact; 
    int numr; 
    int numz; 
    double D0; 
    double zt; 
    double zmax; 
    double delta; 
    double index_radial; 
    double vAlfven; 
    double ab_C; 
    double ab_N; 
    double v0; 
    double dvdz; 
    double etaT; 
    DPerpType set_profile; 
    string filename; 
    double mx; 
    double taudec; 
    double sigmav; 
    int dmmode; 
    DMprofile dmprof; 
};

#endif