*/
public static void main(String[] args) throws SimuException, GaiaException, IOException, StarSystemSimuException {
// Use DMSSimuDU437 to generate planetary system
DmsSimuDU437 dmsSimuDU437 = new DmsSimuDU437();
// Exoplanet orbital parameters
double eccentricity = 0.27;
double inclination = 0.838; // (rad)
double nodeAngle = 0.4; // position angle of the line of nodes (rad)
double omega2 = 0.1; // argument of periastron (rad)
double period = 30.26; // (days)
//double semiMajorAxis = 0.13; // semi-maj. axis (A.U.) of the secondary (a2)
double timePeriastron = 1.0; // (days)
double massPlanet = 1.0; // (Jupiter masses)
double bondAlbedo = 0.5; // Albedo; shouldn't matter for astrometry
// Star data
//long idRoot = 1223345453L;
double magMv = 10.17; // Absolute magnitude. 51 Peg: +4.7. GOG does nothing for G < 12.0
double vMinusI = 1.52; // Intrinsec Mean (V-I) colour. 51 Peg: +0.8
double absV = 0.0; // Interstellar absortion in the V band
double ra = 343.31972458; // RA ICRS (deg). 51 Peg: 344.3665854
double dec = -14.263700556; // DEC ICRS (deg). 51 Peg: 20.76883222
//double parallax = 65.10; // (mas)
//double muRa = 208.07; // Proper motion RA (mas/yr)
//double muDec = 60.96; // Proper motion DEC (mas/yr)
double parallax = 213.3; // (mas). 51 Peg: 64.07 mas
double distance = 1000 / parallax;// (pc)
double muRa = 0.0; // Proper motion RA (mas/yr)
double muDec = 0.0; // Proper motion DEC (mas/yr)
double vRad = -1.7; // Radial velocity (km/s)
double feH = 0.05; // [Fe/H]
double alphaE = 0.0; // Alpha elements enhancement
double mass = 0.334; // (Msun?)
double magBol = 12.1; // (Compute from bolometric corrections in Allen)
int pop = 6; // (Besancon model. 6: 5-7 Gyr)
String spType = "M4 V"; // String defining the spectral type
//double teff = 5500; // (G2)
//double logg = 4.4; // log g (cgs)
//double radius = 2.0; // (Rsun?)
// Generate solar system with star and planet
ExoplanetSystem system = new ExoplanetSystem("GJ 876");
Star star = system.createStar(spType, magMv, distance, vMinusI, absV, ra, dec, parallax, muRa, muDec, vRad, mass, magBol, pop, feH, alphaE);
OrbitalParams orbParam = dmsSimuDU437.generateOrbitalParams(
star, massPlanet / Math.sin(inclination), period, timePeriastron, eccentricity, omega2, nodeAngle, inclination);
// StellarSource source = dmsSimuDU437.generateSystem( star, orbitalParams );
// Copied for adjustment from DmsSimuDU437.java