GyotoDynamicalDisk3D.h
Go to the documentation of this file.
/*
Copyright 2011-2016, 2018 Frederic Vincent, Thibaut Paumard
This file is part of Gyoto.
Gyoto is free software: you can redistribute it and/or modify
it under the terms of the GNU General Public License as published by
the Free Software Foundation, either version 3 of the License, or
(at your option) any later version.
Gyoto is distributed in the hope that it will be useful,
but WITHOUT ANY WARRANTY; without even the implied warranty of
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
GNU General Public License for more details.
You should have received a copy of the GNU General Public License
along with Gyoto. If not, see <http://www.gnu.org/licenses/>.
*/
#ifndef __GyotoDynamicalDisk3D_H_
#define __GyotoDynamicalDisk3D_H_
#include <iostream>
#include <fstream>
#include <iomanip>
#include <cstring>
namespace Gyoto{
namespace Astrobj { class DynamicalDisk3D; }
}
#include <GyotoDisk3D.h>
#include <GyotoBlackBodySpectrum.h>
class Gyoto::Astrobj::DynamicalDisk3D : public Astrobj::Disk3D {
friend class Gyoto::SmartPointer<Gyoto::Astrobj::DynamicalDisk3D>;
protected:
SmartPointer<Spectrum::BlackBody> spectrumBB_;
private:
bool temperature_;
char* dirname_;
double tinit_;
double dt_;
int nb_times_;
double PLindex_;
bool novel_;
double floortemperature_;
double ** emission_array_;
double ** absorption_array_;
double ** velocity_array_;
// Constructors - Destructor
// -------------------------
public:
GYOTO_OBJECT;
GYOTO_OBJECT_THREAD_SAFETY;
DynamicalDisk3D();
DynamicalDisk3D(const DynamicalDisk3D& ) ;
virtual DynamicalDisk3D* clone () const;
virtual ~DynamicalDisk3D() ;
// Accessors
// ---------
public:
using Generic::metric;
void metric(SmartPointer<Metric::Generic> gg);
void file(std::string const &f);
std::string file() const;
void tinit(double t);
double tinit()const;
void dt(double t);
double dt()const;
void PLindex(double t);
double PLindex()const;
void floorTemperature(double t);
double floorTemperature()const;
void temperature(bool t);
bool temperature() const;
void withVelocity(bool t);
bool withVelocity() const;
// Stuff
// -----
double emission1date(double nu_em, double dsem,
state_t const &c_ph, double const c_obj[8]) const;
using Disk3D::emission;
virtual double emission(double nu_em, double dsem,
state_t const &c_ph, double const c_obj[8]=NULL) const;
double transmission1date(double nu_em, double dsem,
state_t const &c_ph, double const c_obj[8]) const;
double transmission(double nu_em, double dsem,
state_t const &cph, double const co[8]) const;
void getVelocity(double const pos[4], double vel[4]);
double const * getVelocity() const;
protected:
void copyQuantities(int iq) ;
};
#endif
GyotoBlackBodySpectrum.h
I_nu(nu, T) = cst_*2*h*nu^3/c^2/(exp(h*nu/k*T)-1.);.
GyotoDisk3D.h
A geometrically thick, optically thin disk.
GYOTO_OBJECT
#define GYOTO_OBJECT
Declare class::properties and class::getProperties().
Definition GyotoObject.h:84
Gyoto::Astrobj::Disk3D
Geometrically thick disk read from FITS file.
Definition GyotoDisk3D.h:56
Gyoto::Astrobj::DynamicalDisk3D
Geometrically thick optically thin disk read from a set of FITS files.
Definition GyotoDynamicalDisk3D.h:58
Gyoto::Astrobj::DynamicalDisk3D::velocity_array_
double ** velocity_array_
velocity(r, z, phi)
Definition GyotoDynamicalDisk3D.h:93
Gyoto::Astrobj::DynamicalDisk3D::clone
virtual DynamicalDisk3D * clone() const
Cloner.
Gyoto::Astrobj::DynamicalDisk3D::emission1date
double emission1date(double nu_em, double dsem, state_t const &c_ph, double const c_obj[8]) const
Compute emission at one grid date.
Gyoto::Astrobj::DynamicalDisk3D::spectrumBB_
SmartPointer< Spectrum::BlackBody > spectrumBB_
emission law
Definition GyotoDynamicalDisk3D.h:61
Gyoto::Astrobj::DynamicalDisk3D::metric
void metric(SmartPointer< Metric::Generic > gg)
Set the Metric gg_.
Gyoto::Astrobj::DynamicalDisk3D::absorption_array_
double ** absorption_array_
Definition GyotoDynamicalDisk3D.h:85
Gyoto::Astrobj::DynamicalDisk3D::tinit_
double tinit_
Time of the first FITS file.
Definition GyotoDynamicalDisk3D.h:66
Gyoto::Astrobj::DynamicalDisk3D::getVelocity
void getVelocity(double const pos[4], double vel[4])
Get fluid 4-velocity at point.
Gyoto::Astrobj::DynamicalDisk3D::dirname_
char * dirname_
FITS files directory.
Definition GyotoDynamicalDisk3D.h:65
Gyoto::Astrobj::DynamicalDisk3D::transmission1date
double transmission1date(double nu_em, double dsem, state_t const &c_ph, double const c_obj[8]) const
Compute transmission at one grid date.
Gyoto::Astrobj::DynamicalDisk3D::PLindex_
double PLindex_
power law index such that density_elec(E) ∝ E-p
Definition GyotoDynamicalDisk3D.h:69
Gyoto::Astrobj::DynamicalDisk3D::emission
virtual double emission(double nu_em, double dsem, state_t const &c_ph, double const c_obj[8]=NULL) const
Interpolate emission between grid dates.
Gyoto::Astrobj::DynamicalDisk3D::novel_
bool novel_
put to true if velocity of emitting particle is not provided
Definition GyotoDynamicalDisk3D.h:70
Gyoto::Astrobj::DynamicalDisk3D::transmission
double transmission(double nu_em, double dsem, state_t const &cph, double const co[8]) const
Interpolate transmission between grid dates.
Gyoto::Astrobj::DynamicalDisk3D::floortemperature_
double floortemperature_
if non-zero, emission and absorption are 0 for temperatures below this floor, emission=blackbody and ...
Definition GyotoDynamicalDisk3D.h:71
Gyoto::Astrobj::DynamicalDisk3D::emission_array_
double ** emission_array_
Definition GyotoDynamicalDisk3D.h:78
Gyoto::Astrobj::DynamicalDisk3D::copyQuantities
void copyQuantities(int iq)
Set underlying Disk3D pointers to a specific date slice.
Gyoto::Astrobj::DynamicalDisk3D::dt_
double dt_
Time increment between two FITS (assumed constant).
Definition GyotoDynamicalDisk3D.h:67
Gyoto::Astrobj::DynamicalDisk3D::nb_times_
int nb_times_
Number of times.
Definition GyotoDynamicalDisk3D.h:68
Gyoto::Astrobj::DynamicalDisk3D::DynamicalDisk3D
DynamicalDisk3D()
Standard constructor.
Gyoto::Astrobj::DynamicalDisk3D::temperature_
bool temperature_
1 if temperature is given in fits data file, 0 if emission coef is directly given
Definition GyotoDynamicalDisk3D.h:64
Gyoto::Astrobj::Generic::metric
virtual SmartPointer< Metric::Generic > metric() const
Get the Metric gg_.
Gyoto::Astrobj::Generic::emission
virtual double emission(double nu_em, double dsem, state_t const &coord_ph, double const coord_obj[8]=NULL) const
Specific intensity Iν.
Gyoto::SmartPointer
Pointers performing reference counting.
Definition GyotoSmartPointer.h:135
Gyoto::Astrobj
Access to astronomical objects.
Definition GyotoAstrobj.h:49
Gyoto
Namespace for the Gyoto library.
Definition GyotoAstrobj.h:44