47 static double density_initial,
66 density_initial =
dense.xMassDensity;
78 if( !
wind.lgStatic() )
159 double angle_wrt_los=-1. , Border_init=-1.;
184 angle_wrt_los *=
PI2 / 360.;
187 Bpar_init = Border_init*cos(angle_wrt_los);
188 Btan_init = Border_init*sin(angle_wrt_los);
205 " This value of gamma (%.3e) is impossible. Must have gamma = 0 or > 1.\n Sorry.\n",
220 strcpy(
optimize.chVarFmt[
optimize.nparm],
"MAGNETIC FIELD TANGLED= %f LOG GAMMA= %f" );
229 strcpy(
optimize.chVarFmt[
optimize.nparm],
"MAGNETIC FIELD ORDERED= %f LOG ANGLE RADIANS= %f" );
#define DEBUG_ENTRY(funcname)
bool nMatch(const char *chKey) const
double getNumberCheckAlwaysLog(const char *chDesc)
double getNumberCheck(const char *chDesc)
double getNumberDefault(const char *chDesc, double fdef)
static double Btangl_here
void Magnetic_evaluate(void)
static double Btangl_init
static bool lgBinitialized
void Magnetic_reinit(void)
void ParseMagnet(Parser &p)