Skip to content

Commit c2d6be7

Browse files
author
Rory Barnes
committed
Added new tests for backward integrations with STELLAR and AtmEsc.
Fixed missing iModuleBit assignments in POISE.
1 parent f76811e commit c2d6be7

File tree

8 files changed

+356
-0
lines changed

8 files changed

+356
-0
lines changed

‎src/poise.c‎

Lines changed: 4 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -4064,6 +4064,7 @@ void InitializeOutputPoise(OUTPUT *output, fnWriteOutput fnWrite[]) {
40644064
output[OUT_TEMPLANDLAT].dNeg = 1;
40654065
output[OUT_TEMPLANDLAT].iNum = 1;
40664066
output[OUT_TEMPLANDLAT].bGrid = 1;
4067+
output[OUT_TEMPMAXWATER].iModuleBit = POISE;
40674068
fnWrite[OUT_TEMPLANDLAT] = &WriteTempLandLat;
40684069

40694070
fvFormattedString(&output[OUT_TEMPWATERLAT].cName, "TempWaterLat");
@@ -4075,6 +4076,7 @@ void InitializeOutputPoise(OUTPUT *output, fnWriteOutput fnWrite[]) {
40754076
output[OUT_TEMPWATERLAT].dNeg = 1;
40764077
output[OUT_TEMPWATERLAT].iNum = 1;
40774078
output[OUT_TEMPWATERLAT].bGrid = 1;
4079+
output[OUT_TEMPMAXWATER].iModuleBit = POISE;
40784080
fnWrite[OUT_TEMPWATERLAT] = &WriteTempWaterLat;
40794081

40804082
fvFormattedString(&output[OUT_LATITUDE].cName, "Latitude");
@@ -4100,6 +4102,7 @@ void InitializeOutputPoise(OUTPUT *output, fnWriteOutput fnWrite[]) {
41004102
output[OUT_ALBEDOLANDLAT].bNeg = 0;
41014103
output[OUT_ALBEDOLANDLAT].iNum = 1;
41024104
output[OUT_ALBEDOLANDLAT].bGrid = 1;
4105+
output[OUT_TEMPMAXWATER].iModuleBit = POISE;
41034106
fnWrite[OUT_ALBEDOLANDLAT] = &WriteAlbedoLandLat;
41044107

41054108
fvFormattedString(&output[OUT_ALBEDOWATERLAT].cName, "AlbedoWaterLat");
@@ -4108,6 +4111,7 @@ void InitializeOutputPoise(OUTPUT *output, fnWriteOutput fnWrite[]) {
41084111
output[OUT_ALBEDOWATERLAT].bNeg = 0;
41094112
output[OUT_ALBEDOWATERLAT].iNum = 1;
41104113
output[OUT_ALBEDOWATERLAT].bGrid = 1;
4114+
output[OUT_TEMPMAXWATER].iModuleBit = POISE;
41114115
fnWrite[OUT_ALBEDOWATERLAT] = &WriteAlbedoWaterLat;
41124116

41134117
fvFormattedString(&output[OUT_ANNUALINSOL].cName, "AnnInsol");

‎tests/AtmescStellar/Backward/a.in‎

Lines changed: 15 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,15 @@
1+
# Star
2+
sName a # Body's name
3+
saModules stellar # Modules to apply, exact spelling required
4+
5+
# Physical Parameters
6+
dAge 1e9
7+
dMass 1
8+
dRotPeriod -1
9+
10+
# Stellar Parameters
11+
sStellarModel baraffe
12+
sMagBrakingModel briemann
13+
14+
# Output
15+
saOutputOrder Time Age -RotPer -Luminosity -Radius Temperature -TotEn -TotAngMom

‎tests/AtmescStellar/Backward/b.in‎

Lines changed: 14 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,14 @@
1+
# This planet is experiencing Roche lobe overflow
2+
sName b
3+
saModules atmesc
4+
dXFrac 1.0
5+
dSurfWaterMass 0
6+
bHaltEnvelopeGone 0
7+
sPlanetRadiusModel lopez
8+
dAtmXAbsEffH 0.1
9+
dSemi -0.1
10+
dEcc 0.
11+
dEnvelopeMass -0.81818
12+
dMass -1.81818
13+
bAtmEscAuto 1
14+
saOutputOrder Time -EnvelopeMass Radius -Mass DeltaTime

‎tests/AtmescStellar/Backward/test_Backward.py‎

Lines changed: 180 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,180 @@
1+
from benchmark import Benchmark, benchmark
2+
import astropy.units as u
3+
4+
@benchmark(
5+
{
6+
"log.initial.system.Age": {"value": 3.155760e+16, "unit": u.sec},
7+
"log.initial.system.Time": {"value": -0.000000, "unit": u.sec},
8+
"log.initial.system.TotAngMom": {"value": 5.058770e+42, "unit": (u.kg * u.m ** 2) / u.sec},
9+
"log.initial.system.TotEnergy": {"value": -2.463825e+41, "unit": u.erg},
10+
"log.initial.system.PotEnergy": {"value": -2.465658e+41, "unit": u.Joule},
11+
"log.initial.system.KinEnergy": {"value": 1.833857e+38, "unit": u.Joule},
12+
"log.initial.system.DeltaTime": {"value": 0.000000, "unit": u.sec},
13+
"log.initial.a.Mass": {"value": 1.988416e+30, "unit": u.kg},
14+
"log.initial.a.Radius": {"value": 100.680639, "unit": u.Rearth},
15+
"log.initial.a.RadGyra": {"value": 0.290831},
16+
"log.initial.a.RotAngMom": {"value": 5.043469e+42, "unit": (u.kg * u.m ** 2) / u.sec},
17+
"log.initial.a.RotVel": {"value": 4.669855e+04, "unit": u.m / u.sec},
18+
"log.initial.a.BodyType": {"value": 0.000000},
19+
"log.initial.a.RotRate": {"value": 7.272205e-05, "unit": 1 / u.sec},
20+
"log.initial.a.RotPer": {"value": 1.000000, "unit": u.day},
21+
"log.initial.a.Density": {"value": 1792.696544, "unit": u.kg / u.m ** 3},
22+
"log.initial.a.HZLimitDryRunaway": {"value": 1.175525e+11, "unit": u.m},
23+
"log.initial.a.HZLimRecVenus": {"value": 9.760054e+10, "unit": u.m},
24+
"log.initial.a.HZLimRunaway": {"value": 1.279108e+11, "unit": u.m},
25+
"log.initial.a.HZLimMoistGreenhouse": {"value": 1.291407e+11, "unit": u.m},
26+
"log.initial.a.HZLimMaxGreenhouse": {"value": 2.207286e+11, "unit": u.m},
27+
"log.initial.a.HZLimEarlyMars": {"value": 2.408255e+11, "unit": u.m},
28+
"log.initial.a.Instellation": {"value": -1.000000, "unit": u.kg / u.sec ** 3},
29+
"log.initial.a.CriticalSemiMajorAxis": {"value": -1.000000, "unit": u.m},
30+
"log.initial.a.LXUVTot": {"value": 1.697389e+22, "unit": u.kg / u.sec ** 3},
31+
"log.initial.a.LostEnergy": {"value": 5.562685e-309, "unit": u.Joule},
32+
"log.initial.a.LostAngMom": {"value": 5.562685e-309, "unit": (u.kg * u.m ** 2) / u.sec},
33+
"log.initial.a.EscapeVelocity": {"value": 6.429130e+05, "unit": u.m / u.sec},
34+
"log.initial.a.Luminosity": {"value": 0.749501, "unit": u.LSUN},
35+
"log.initial.a.LXUVStellar": {"value": 1.697389e+22, "unit": u.W},
36+
"log.initial.a.Temperature": {"value": 5591.770652, "unit": u.K},
37+
"log.initial.a.LXUVFrac": {"value": 5.888437e-05},
38+
"log.initial.a.RossbyNumber": {"value": 0.063834},
39+
"log.initial.a.DRotPerDtStellar": {"value": 4.087843e-11},
40+
"log.initial.a.WindTorque": {"value": 2.390187e+27},
41+
"log.initial.b.Mass": {"value": 1.818180, "unit": u.Mearth},
42+
"log.initial.b.Radius": {"value": 6.987673e+07, "unit": u.m},
43+
"log.initial.b.RadGyra": {"value": 0.500000},
44+
"log.initial.b.BodyType": {"value": 0.000000},
45+
"log.initial.b.Density": {"value": 7.597731, "unit": u.kg / u.m ** 3},
46+
"log.initial.b.HZLimitDryRunaway": {"value": 1.175525e+11, "unit": u.m},
47+
"log.initial.b.HZLimRecVenus": {"value": 9.760054e+10, "unit": u.m},
48+
"log.initial.b.HZLimRunaway": {"value": 1.279108e+11, "unit": u.m},
49+
"log.initial.b.HZLimMoistGreenhouse": {"value": 1.291407e+11, "unit": u.m},
50+
"log.initial.b.HZLimMaxGreenhouse": {"value": 2.207286e+11, "unit": u.m},
51+
"log.initial.b.HZLimEarlyMars": {"value": 2.408255e+11, "unit": u.m},
52+
"log.initial.b.Instellation": {"value": 1.024993e+05, "unit": u.kg / u.sec ** 3},
53+
"log.initial.b.MeanMotion": {"value": 6.296061e-06, "unit": 1 / u.sec},
54+
"log.initial.b.OrbPeriod": {"value": 9.979550e+05, "unit": u.sec},
55+
"log.initial.b.SemiMajorAxis": {"value": 1.495979e+10, "unit": u.m},
56+
"log.initial.b.LXUVTot": {"value": -1.000000, "unit": u.kg / u.sec ** 3},
57+
"log.initial.b.EscapeVelocity": {"value": 4554.450923, "unit": u.m / u.sec},
58+
"log.initial.b.SurfWaterMass": {"value": 0.000000, "unit": u.kg},
59+
"log.initial.b.EnvelopeMass": {"value": 0.818180, "unit": u.Mearth},
60+
"log.initial.b.OxygenMass": {"value": 0.000000, "unit": u.kg},
61+
"log.initial.b.RGLimit": {"value": 1.228298e+11, "unit": u.m},
62+
"log.initial.b.XO": {"value": 0.000000},
63+
"log.initial.b.EtaO": {"value": 0.000000},
64+
"log.initial.b.PlanetRadius": {"value": 6.987673e+07, "unit": u.m},
65+
"log.initial.b.OxygenMantleMass": {"value": 0.000000, "unit": u.kg},
66+
"log.initial.b.RadXUV": {"value": -1.000000, "unit": u.m},
67+
"log.initial.b.RadSolid": {"value": -1.000000, "unit": u.m},
68+
"log.initial.b.PresXUV": {"value": 5.000000, "unit": (null)},
69+
"log.initial.b.ScaleHeight": {"value": -1.000000, "unit": u.m},
70+
"log.initial.b.ThermTemp": {"value": 400.000000, "unit": u.K},
71+
"log.initial.b.AtmGasConst": {"value": 4124.000000, "unit": (null)},
72+
"log.initial.b.PresSurf": {"value": -1.000000, "unit": u.Pa},
73+
"log.initial.b.DEnvMassDt": {"value": -1.965548e+09, "unit": u.kg / u.sec},
74+
"log.initial.b.FXUV": {"value": 6.035605, "unit": u.W / u.m ** 2},
75+
"log.initial.b.AtmXAbsEffH2O": {"value": 0.300000},
76+
"log.initial.b.RocheRadius": {"value": 1.826583e+08, "unit": u.m},
77+
"log.initial.b.BondiRadius": {"value": 7.395371e+07, "unit": u.m},
78+
"log.initial.b.HEscapeRegime": {"value": 3.000000},
79+
"log.initial.b.RRCriticalFlux": {"value": 0.062660, "unit": u.W / u.m ** 2},
80+
"log.initial.b.CrossoverMass": {"value": 0.000000, "unit": u.kg},
81+
"log.initial.b.WaterEscapeRegime": {"value": 8.000000},
82+
"log.initial.b.FXUVCRITDRAG": {"value": 0.000670, "unit": u.W / u.m ** 2},
83+
"log.initial.b.HREFFLUX": {"value": 5.787364e+19, "unit": 1 / u.m ** 2 / u.sec},
84+
"log.initial.b.XO2": {"value": 0.000000},
85+
"log.initial.b.XH2O": {"value": 0.000000},
86+
"log.initial.b.HDiffFlux": {"value": 2.874041e+15, "unit": 1 / u.m ** 2 / u.sec},
87+
"log.initial.b.HRefODragMod": {"value": 1.000000},
88+
"log.initial.b.KTide": {"value": 0.454161},
89+
"log.initial.b.CumulativeXUVFlux": {"value": 0.000000, "unit": u.kg / u.sec ** 3},
90+
"log.initial.b.RunawayGreenhouseFlux": {"value": 1520.422432, "unit": u.kg / u.sec ** 3},
91+
"log.initial.b.RGDuration": {"value": 0.00000e+00, "unit": u.yr},
92+
"log.final.system.Age": {"value": 1.577880e+14, "unit": u.sec, "rtol": 1e-4},
93+
"log.final.system.Time": {"value": -3.139981e+16, "unit": u.sec, "rtol": 1e-4},
94+
"log.final.system.TotAngMom": {"value": 5.072047e+42, "unit": (u.kg * u.m ** 2) / u.sec, "rtol": 1e-4},
95+
"log.final.system.TotEnergy": {"value": -2.464150e+41, "unit": u.erg, "rtol": 1e-4},
96+
"log.final.system.PotEnergy": {"value": -1.587322e+41, "unit": u.Joule, "rtol": 1e-4},
97+
"log.final.system.KinEnergy": {"value": 5.860165e+40, "unit": u.Joule, "rtol": 1e-4},
98+
"log.final.a.Mass": {"value": 1.988416e+30, "unit": u.kg, "rtol": 1e-4},
99+
"log.final.a.Radius": {"value": 156.391718, "unit": u.Rearth, "rtol": 1e-4},
100+
"log.final.a.RadGyra": {"value": 0.449294, "rtol": 1e-4},
101+
"log.final.a.RotAngMom": {"value": 2.163512e+44, "unit": (u.kg * u.m ** 2) / u.sec, "rtol": 1e-4},
102+
"log.final.a.RotVel": {"value": 5.403630e+05, "unit": u.m / u.sec, "rtol": 1e-4},
103+
"log.final.a.BodyType": {"value": 0.000000, "rtol": 1e-4},
104+
"log.final.a.RotRate": {"value": 0.000542, "unit": 1 / u.sec, "rtol": 1e-4},
105+
"log.final.a.RotPer": {"value": 0.134241, "unit": u.day, "rtol": 1e-4},
106+
"log.final.a.Density": {"value": 478.303342, "unit": u.kg / u.m ** 3, "rtol": 1e-4},
107+
"log.final.a.HZLimitDryRunaway": {"value": 1.084460e+11, "unit": u.m, "rtol": 1e-4},
108+
"log.final.a.HZLimRecVenus": {"value": 9.446714e+10, "unit": u.m, "rtol": 1e-4},
109+
"log.final.a.HZLimRunaway": {"value": 1.254168e+11, "unit": u.m, "rtol": 1e-4},
110+
"log.final.a.HZLimMoistGreenhouse": {"value": 1.249951e+11, "unit": u.m, "rtol": 1e-4},
111+
"log.final.a.HZLimMaxGreenhouse": {"value": 2.269129e+11, "unit": u.m, "rtol": 1e-4},
112+
"log.final.a.HZLimEarlyMars": {"value": 2.475151e+11, "unit": u.m, "rtol": 1e-4},
113+
"log.final.a.Instellation": {"value": -1.000000, "unit": u.kg / u.sec ** 3, "rtol": 1e-4},
114+
"log.final.a.CriticalSemiMajorAxis": {"value": -1.000000, "unit": u.m, "rtol": 1e-4},
115+
"log.final.a.LXUVTot": {"value": 2.453266e+23, "unit": u.kg / u.sec ** 3, "rtol": 1e-4},
116+
"log.final.a.LostEnergy": {"value": -1.462843e+41, "unit": u.Joule, "rtol": 1e-4},
117+
"log.final.a.LostAngMom": {"value": -2.113232e+44, "unit": (u.kg * u.m ** 2) / u.sec, "rtol": 1e-4},
118+
"log.final.a.EscapeVelocity": {"value": 5.158439e+05, "unit": u.m / u.sec, "rtol": 1e-4},
119+
"log.final.a.Luminosity": {"value": 0.637875, "unit": u.LSUN, "rtol": 1e-4},
120+
"log.final.a.LXUVStellar": {"value": 2.453266e+23, "unit": u.W, "rtol": 1e-4},
121+
"log.final.a.Temperature": {"value": 4309.979229, "unit": u.K, "rtol": 1e-4},
122+
"log.final.a.LXUVFrac": {"value": 0.001000, "rtol": 1e-4},
123+
"log.final.a.RossbyNumber": {"value": 0.003889, "rtol": 1e-4},
124+
"log.final.a.DRotPerDtStellar": {"value": -4.682593e-11, "rtol": 1e-4},
125+
"log.final.a.WindTorque": {"value": 1.342986e+28, "rtol": 1e-4},
126+
"log.final.b.Mass": {"value": 5.236487, "unit": u.Mearth, "rtol": 1e-4},
127+
"log.final.b.Radius": {"value": 1.902090e+08, "unit": u.m, "rtol": 1e-4},
128+
"log.final.b.RadGyra": {"value": 0.500000, "rtol": 1e-4},
129+
"log.final.b.BodyType": {"value": 0.000000, "rtol": 1e-4},
130+
"log.final.b.Density": {"value": 1.084905, "unit": u.kg / u.m ** 3, "rtol": 1e-4},
131+
"log.final.b.HZLimitDryRunaway": {"value": 1.084460e+11, "unit": u.m, "rtol": 1e-4},
132+
"log.final.b.HZLimRecVenus": {"value": 9.446714e+10, "unit": u.m, "rtol": 1e-4},
133+
"log.final.b.HZLimRunaway": {"value": 1.254168e+11, "unit": u.m, "rtol": 1e-4},
134+
"log.final.b.HZLimMoistGreenhouse": {"value": 1.249951e+11, "unit": u.m, "rtol": 1e-4},
135+
"log.final.b.HZLimMaxGreenhouse": {"value": 2.269129e+11, "unit": u.m, "rtol": 1e-4},
136+
"log.final.b.HZLimEarlyMars": {"value": 2.475151e+11, "unit": u.m, "rtol": 1e-4},
137+
"log.final.b.Instellation": {"value": 8.723364e+04, "unit": u.kg / u.sec ** 3, "rtol": 1e-4},
138+
"log.final.b.MeanMotion": {"value": 6.296093e-06, "unit": 1 / u.sec, "rtol": 1e-4},
139+
"log.final.b.OrbPeriod": {"value": 9.979499e+05, "unit": u.sec, "rtol": 1e-4},
140+
"log.final.b.SemiMajorAxis": {"value": 1.495979e+10, "unit": u.m, "rtol": 1e-4},
141+
"log.final.b.LXUVTot": {"value": -1.000000, "unit": u.kg / u.sec ** 3, "rtol": 1e-4},
142+
"log.final.b.EscapeVelocity": {"value": 4684.770900, "unit": u.m / u.sec, "rtol": 1e-4},
143+
"log.final.b.SurfWaterMass": {"value": 0.000000, "unit": u.kg, "rtol": 1e-4},
144+
"log.final.b.EnvelopeMass": {"value": 4.236487, "unit": u.Mearth, "rtol": 1e-4},
145+
"log.final.b.OxygenMass": {"value": 0.000000, "unit": u.kg, "rtol": 1e-4},
146+
"log.final.b.RGLimit": {"value": 1.174025e+11, "unit": u.m, "rtol": 1e-4},
147+
"log.final.b.XO": {"value": 0.000000, "rtol": 1e-4},
148+
"log.final.b.EtaO": {"value": 0.000000, "rtol": 1e-4},
149+
"log.final.b.PlanetRadius": {"value": 1.902090e+08, "unit": u.m, "rtol": 1e-4},
150+
"log.final.b.OxygenMantleMass": {"value": 0.000000, "unit": u.kg, "rtol": 1e-4},
151+
"log.final.b.RadXUV": {"value": -1.000000, "unit": u.m, "rtol": 1e-4},
152+
"log.final.b.RadSolid": {"value": -1.000000, "unit": u.m, "rtol": 1e-4},
153+
"log.final.b.PresXUV": {"value": 5.000000, "unit": (null), "rtol": 1e-4},
154+
"log.final.b.ScaleHeight": {"value": -1.000000, "unit": u.m, "rtol": 1e-4},
155+
"log.final.b.ThermTemp": {"value": 400.000000, "unit": u.K, "rtol": 1e-4},
156+
"log.final.b.AtmGasConst": {"value": 4124.000000, "unit": (null), "rtol": 1e-4},
157+
"log.final.b.PresSurf": {"value": -1.000000, "unit": u.Pa, "rtol": 1e-4},
158+
"log.final.b.DEnvMassDt": {"value": 3.419958e+09, "unit": u.kg / u.sec, "rtol": 1e-4},
159+
"log.final.b.FXUV": {"value": 87.233642, "unit": u.W / u.m ** 2, "rtol": 1e-4},
160+
"log.final.b.AtmXAbsEffH2O": {"value": 0.300000, "rtol": 1e-4},
161+
"log.final.b.RocheRadius": {"value": 2.598805e+08, "unit": u.m, "rtol": 1e-4},
162+
"log.final.b.BondiRadius": {"value": 2.217192e+08, "unit": u.m, "rtol": 1e-4},
163+
"log.final.b.HEscapeRegime": {"value": 6.000000, "rtol": 1e-4},
164+
"log.final.b.RRCriticalFlux": {"value": 0.001249, "unit": u.W / u.m ** 2, "rtol": 1e-4},
165+
"log.final.b.CrossoverMass": {"value": 0.000000, "unit": u.kg, "rtol": 1e-4},
166+
"log.final.b.WaterEscapeRegime": {"value": 8.000000, "rtol": 1e-4},
167+
"log.final.b.FXUVCRITDRAG": {"value": 0.000275, "unit": u.W / u.m ** 2, "rtol": 1e-4},
168+
"log.final.b.HREFFLUX": {"value": 3.590455e+21, "unit": 1 / u.m ** 2 / u.sec, "rtol": 1e-4},
169+
"log.final.b.XO2": {"value": 0.000000, "rtol": 1e-4},
170+
"log.final.b.XH2O": {"value": 0.000000, "rtol": 1e-4},
171+
"log.final.b.HDiffFlux": {"value": 1.117118e+15, "unit": 1 / u.m ** 2 / u.sec, "rtol": 1e-4},
172+
"log.final.b.HRefODragMod": {"value": 1.000000, "rtol": 1e-4},
173+
"log.final.b.KTide": {"value": 0.100000, "rtol": 1e-4},
174+
"log.final.b.CumulativeXUVFlux": {"value": 8.269339e+17, "unit": u.kg / u.sec ** 3, "rtol": 1e-4},
175+
"log.final.b.RunawayGreenhouseFlux": {"value": 1416.380837, "unit": u.kg / u.sec ** 3, "rtol": 1e-4},
176+
"log.final.b.RGDuration": {"value": 0.00000e+00, "unit": u.yr, "rtol": 1e-4},
177+
}
178+
)
179+
class Test_Backward(Benchmark):
180+
pass

‎tests/AtmescStellar/Backward/vpl.in‎

Lines changed: 25 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,25 @@
1+
#
2+
sSystemName star # System Name
3+
iVerbose 5 # Verbosity level
4+
bOverwrite 1 # Allow file overwrites?
5+
6+
# List of "body files" that contain body-specific parameters
7+
saBodyFiles a.in b.in # The first star
8+
9+
# Input/Output Units
10+
sUnitMass solar # Options: gram, kg, Earth, Neptune, Jupiter, solar
11+
sUnitLength aU # Options: cm, m, km, Earth, Jupiter, solar, AU
12+
sUnitTime YEARS # Options: sec, day, year, Myr, Gyr
13+
sUnitAngle d # Options: deg, rad
14+
15+
# Input/Output
16+
bDoLog 1 # Write a log file?
17+
iDigits 6 # Maximum number of digits to right of decimal
18+
dMinValue 1e-10 # Minimum value of eccentricity/obliquity
19+
20+
# Evolution Parameters
21+
bDoBackward 1 # Perform a forward evolution?
22+
bVarDt 1 # Use variable timestepping?
23+
dEta 0.01 # Coefficient for variable timestepping
24+
dStopTime 9.95e8 # Stop time for evolution
25+
dOutputTime 1e6 # Output timesteps (assuming in body files)

‎tests/Stellar/BaraffeBackward/a.in‎

Lines changed: 15 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,15 @@
1+
# Star
2+
sName a # Body's name
3+
saModules stellar # Modules to apply, exact spelling required
4+
5+
# Physical Parameters
6+
dAge 1e9
7+
dMass 1
8+
dRotPeriod -1
9+
10+
# Stellar Parameters
11+
sStellarModel baraffe
12+
sMagBrakingModel briemann
13+
14+
# Output
15+
saOutputOrder Time Age -RotPer -Luminosity -Radius Temperature -TotEn -TotAngMom

0 commit comments

Comments
 (0)