/data1/mhar/model/cetimb_1/results/dummy /data1/mhar/model/cetimb_1/results/HD209458b.12.SP /data1/mhar/model/cetimb_1/input_data/full_solar.dat /data1/mhar/model/cetimb_1/input_data/full_solar_xsect.dat t .... sw1 - H3+ cooling switch. t .... sw2 - Tidal locking switch. t .... sw3 - Non-LTE cooling switch. f .... sw4 - Previous simulation input switch. f .... sw5 - Globally average the startup t .... sw6 - Vertical diffusion switch t .... sw7 - Ion and neutral chemistry switch t .... sw8 - H cooling f .... sw9 - Tidal Forcing t .... sw10 - H recombination cooling f .... sw11 - H2 cooling t .... sw12 - H+ cooling 12 .... Number of days to run 1 .... Season (1=equinox; 2=northern summer; 3=southern summer) 1 .... X ray multiplication factor (0 = no X rays, 1 = current Sun) 0.685 .... Planet mass (Jupiters) 1.01 .... Stellar mass (Suns) 0.d0 .... Solar Zenith Angle (1-d model) 0.04707 .... Distance from the host star (AU) 1.32 .... Radius of the planet (Jupiters) 0.04707 .... Semi-major axis 0 .... Obliquity 7.0d3 .... Lower boundary altitude (km) 1.d0 .... Lower boundary pressure (Pa) 0.d0 .... Upper boundary pressure (Pa) 700 .... Temperature of base (if T profile in input, set NEGATIVE) 1.628d-1 f .... Starting He bottom level; allow chemistry to update -1.00E+00 f .... Starting H bottom level (-1 equilib.); allow chem to update -1.00E+00 f .... Starting H2 bottom level(-1 equilib.); allow chem to update 1.d-10 t .... Starting H2+ bottom level; allow chem to update 1.d-10 t .... Starting H+ bottom level; allow chem to update 1.d-10 t .... Starting He+ bottom level; allow chem to update 1.d-20 t .... Starting H3+ bottom level; allow chem to update 1.d-10 t .... Starting e bottom level; allow chem to update 1.d-10 t .... Starting HeH+ bottom level; allow chem to update 3.58d-6 t .... Starting Na bottom level; allow chem to update 1.d-10 t .... Starting Na+ bottom level; allow chem to update 5.94d-5 t .... Starting Si bottom level; allow chem to update 1.d-10 t .... Starting Si+ bottom level; allow chem to update 1.d-10 t .... Starting Si++ bottom level; allow chem to update 1.4d-3 t .... Starting O bottom level; allow chem to update 1.d-10 t .... Starting O+ bottom level; allow chem to update 6.1d-5 t .... Starting C bottom level; allow chem to update 1.d-10 t .... Starting C+ bottom level; allow chem to update 1.d-10 t .... Starting Si+++ bottom level; allow chem to update 1.d-10 t .... Starting CH ; allow chemistry to update 1.d-10 t .... Starting CH+ ; allow chemistry to update 1.d-10 t .... Starting CH2 ; allow chemistry to update 1.d-10 t .... Starting CH2+ ; allow chemistry to update 1.d-10 t .... Starting CH3 ; allow chemistry to update 1.d-10 t .... Starting CH3+ ; allow chemistry to update 1.d-10 t .... Starting CH4 ; allow chemistry to update 1.d-10 t .... Starting CH4+ ; allow chemistry to update 1.d-10 t .... Starting CH5+ ; allow chemistry to update 1.d-10 t .... Starting CO ; allow chemistry to update 1.d-10 t .... Starting CO+ ; allow chemistry to update 1.d-10 t .... Starting D ; allow chemistry to update 1.d-10 t .... Starting D+ ; allow chemistry to update 1.d-10 t .... Starting H2D+ ; allow chemistry to update 1.d-10 t .... Starting H2O ; allow chemistry to update 1.d-10 t .... Starting H2O+ ; allow chemistry to update 1.d-10 t .... Starting H3O+ ; allow chemistry to update 1.d-10 t .... Starting HCO+ ; allow chemistry to update 1.d-10 t .... Starting HD ; allow chemistry to update 1.d-10 t .... Starting HD+ ; allow chemistry to update 1.d-10 t .... Starting HN2+ ; allow chemistry to update 1.d-10 t .... Starting N ; allow chemistry to update 1.d-10 t .... Starting N+ ; allow chemistry to update 1.d-10 t .... Starting N2 ; allow chemistry to update 1.d-10 t .... Starting N2+ ; allow chemistry to update 1.d-10 t .... Starting NH ; allow chemistry to update 1.d-10 t .... Starting NH+ ; allow chemistry to update 1.d-10 t .... Starting NH2 ; allow chemistry to update 1.d-10 t .... Starting NH2+ ; allow chemistry to update 1.d-10 t .... Starting O2 ; allow chemistry to update 1.d-10 t .... Starting O2+ ; allow chemistry to update 1.d-10 t .... Starting OH ; allow chemistry to update 1.d-10 t .... Starting OH+ ; allow chemistry to update 0.0 f .... Starting bottom level gas fraction; allow chem to update 0.0 f .... Starting bottom level gas fraction; allow chem to update 0.0 f .... Starting bottom level gas fraction; allow chem to update 0.0 f .... Starting bottom level gas fraction; allow chem to update 0.0 f .... Starting bottom level gas fraction; allow chem to update 0.0 f .... Starting bottom level gas fraction; allow chem to update 0.0 f .... Starting bottom level gas fraction; allow chem to update 0.0 f .... Starting bottom level gas fraction; allow chem to update 0.0 f .... Starting bottom level gas fraction; allow chem to update 0.0 f .... Starting bottom level gas fraction; allow chem to update 0.0 f .... Starting bottom level gas fraction; allow chem to update 0.0 f .... Starting bottom level gas fraction; allow chem to update 10.d3 .... Vertical altitude spacing at bottom level (m) 1.014 .... Vertical grid stretch factor 2.0d4 .... Eddy diffusion coefficient 1.0 .... Time step (in seconds) 84.59397 .... Day length in hours 3.524749 .... Period (in Earth days) 720 .... Horizontal smoothing frequency for composition (time step frequency, 0 for none) 720 .... Horizontal smoothing frequency for thermal energy (CvT) 720 .... Horizontal smoothing frequency for winds 10 .... Vertical smoothing frequency 1 .... Output frequency (days) 1 .... Chemistry call frequency 10 .... Radiative heating/cooling call frequency 0.5 .... Photoelectron Neutral heating efficiency (fraction of 1) 0.5 0.00 .... Molecular line heating efficiency 5 .... advection: 16-CN; 15-Spherical Wilson; 5-Van-Leer 0 .... plotting: 0-off;1-plotfiles no tranfer;2-plotfiles+transfer 1 .... verbose debugging - more diagnostic output when running 3000 .... Verbose debug output frequency 8000000. .... maximum vertical velocity cap 2 .... Chemistry integrator. 1-Old CETIM; 2-ASAD Impact; 3-ASAD SVODE 3 .... Upper boundary 0-set_to_lev_below;1-Conserve; 2-Jean's flux; 3-Extrap f .... f-Diffuse all gases as majors; t-Treat some as minors if applicable 0 .... Dummy placeholder until new switcch needed 0 .... Dummy placeholder until new switcch needed 0 .... Dummy placeholder until new switcch needed 0 .... Dummy placeholder until new switcch needed *************************** EXOTIM INPUT FILE *************************** - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Switches/other input are: - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - NOTE: Please read the descriptions below before running the model!!! sw1 .... H3+ cooling switch. (On/Off) True = Include infrared cooling by H3+ ions. False = Generate a model excluding infrared cooling. sw2 .... Tidal locking switch. (On/Off) True = The planet is rotationally synchronised with the parent star. Likely to be the case within 0.1 AU from the host for a Jupiter-type EGP. Dummy daylength value should still be specified to control the run duration. False = Asynchronous rotation rate is based on the user-specified daylength. sw3 .... Non-LTE cooling switch. (On/Off) True = H3+ infrared emissions are damped by a non-lte factor that is based on detailed balance calculations performed by Steve Miller. False = LTE emission rates are used for H3+ emissions. This is not normarly appropriate for the upper thermosphere. sw4 .... Previous simulation input switch. (On/Off) True = Use one of the previous runs as a startup atmosphere. The top line in the instream above is the input file name. False = Standard average atmosphere is set up with a default temperature (temperature of base), winds (=0) and composition. The user specifies the mole fractions for He and H in relation to H2 at the lower boundary and an initial composition is then calculated by assuming hydrostatic and diffusive equilibria. The top line in the instream is ignored. sw5 .... Globally average the startup atmosphere. (On/Off) True = Globally average the startup atmosphere input. This switch should be false if a start up atmosphere is used instead of a previous run. sw6 .... Vertical diffusion switch sw7 .... Ion and neutral chemistry switch sw8 .... Hydrogen Lyman alpha cooling sw9 .... Stellar tidal forcing NOTE: Horizontal smoothing frequency is set as times per day for energy, winds and composition. Note that if day lenght is, say, halved, smoothing is doubled for tidally locked targets. This affects the results at close-in orbits where the model is hot. SOL_ACT.. Solar activity level: 1 = low ; 2 = high SEASON... 1 = equinox 2 = northern hemisphere solstice 3 = southern hemisphere solstice NN_DAYS.. Number of days to run code