Model Storage: Preservation and Reuse¶
Storing data is often tricky. You may be wondering what output to save, what output to throw away. You may also have questions surround what formats to use to store your data.
PICASO’s recommended model output is xarray
. There are several xarray
tutorials throughout the picaso documentation. The most thorough is in the data uniformity tutorial.
Here you will learn how to use a nifty picaso function that will store your the minimum data needed to reproduce your 1D spectrum. The one caveat is that we will not store opacity data. However, you can always get your reference opacity data using the auto citation tools.
For this tutorial we assume you already know the basics of running 1D models. Note that we have still not enabled this functionality for 3D and 4D runs. Please contact the developers if you are interested in this functionality.
[1]:
import picaso.justdoit as jdi
import picaso.justplotit as jpi
u = jdi.u #astropy units
Storing your run with xarray¶
Here we show a simple example of how to use the output_xarray
function for preservation
[2]:
opa = jdi.opannection(wave_range=[0.3,1])
pl = jdi.inputs()#calculation='brown')
pl.gravity(radius=1, radius_unit=u.Unit('R_jup'),
mass=1, mass_unit=u.Unit('M_jup'))
pl.atmosphere(filename=jdi.jupiter_pt(), sep='\s+')
pl.phase_angle(0)
pl.clouds(filename=jdi.jupiter_cld(), sep='\s+')
pl.star(opa, 5000,0,4, radius=1, radius_unit=u.Unit("R_sun"), semi_major=1, semi_major_unit=1*u.AU)
#MUST USE full output=True for this functionality
df= pl.spectrum(opa, calculation='reflected', full_output=True)
[3]:
preserve = jdi.output_xarray(df,pl,savefile='/data/picaso_dbs/model.nc')
[4]:
preserve
[4]:
<xarray.Dataset> Size: 579kB Dimensions: (pressure: 61, wavelength: 12039, pressure_layer: 60, wavenumber_layer: 196) Coordinates: * pressure (pressure) float64 488B 1e-06 1.41e-06 2e-06 ... 708.0 1e+03 * wavelength (wavelength) float64 96kB 0.9999 0.9998 ... 0.3001 0.3 Dimensions without coordinates: pressure_layer, wavenumber_layer Data variables: (12/21) temperature (pressure) float64 488B 150.9 149.7 ... 1.525e+03 1.675e+03 albedo (wavelength) float64 96kB 0.2648 0.2092 ... 0.7193 0.7193 fpfs_reflected (wavelength) float64 96kB 6.048e-08 4.777e-08 ... 1.643e-07 e- (pressure) float64 488B 4.5e-38 4.5e-38 ... 2.15e-11 H2 (pressure) float64 488B 0.837 0.837 0.837 ... 0.836 0.836 H (pressure) float64 488B 4.5e-38 4.5e-38 ... 3.37e-06 ... ... NH3 (pressure) float64 488B 0.000137 0.000137 ... 0.000155 N2 (pressure) float64 488B 5.42e-17 1.58e-17 ... 2.73e-05 PH3 (pressure) float64 488B 4.5e-38 4.5e-38 ... 5.03e-07 opd (pressure_layer, wavenumber_layer) float64 94kB 0.0 ... 0.0 ssa (pressure_layer, wavenumber_layer) float64 94kB 0.0 ... 0.0 asy (pressure_layer, wavenumber_layer) float64 94kB 0.0 ... 0.0 Attributes: planet_params: {"mp": {"value": 1.8981245973360504e+30, "unit": "g"}, "... stellar_params: {"database": "ck04models", "steff": 5000, "feh": 0, "log... orbit_params: {"sma": {"value": 14959787070000.0, "unit": "cm"}}
- pressure: 61
- wavelength: 12039
- pressure_layer: 60
- wavenumber_layer: 196
- pressure(pressure)float641e-06 1.41e-06 ... 708.0 1e+03
- units :
- bar
array([1.0000e-06, 1.4100e-06, 2.0000e-06, 2.8200e-06, 3.9800e-06, 5.6200e-06, 7.9400e-06, 1.1220e-05, 1.5850e-05, 2.2390e-05, 3.1620e-05, 4.4670e-05, 6.3100e-05, 8.9120e-05, 1.2589e-04, 1.7783e-04, 2.5119e-04, 3.5481e-04, 5.0119e-04, 7.0795e-04, 1.0000e-03, 1.4125e-03, 1.9953e-03, 2.8184e-03, 3.9811e-03, 5.6234e-03, 7.9433e-03, 1.1220e-02, 1.5849e-02, 2.2387e-02, 3.1623e-02, 4.4668e-02, 6.3096e-02, 8.9125e-02, 1.2589e-01, 1.7783e-01, 2.5119e-01, 3.5481e-01, 5.0119e-01, 7.0795e-01, 1.0000e+00, 1.4125e+00, 1.9953e+00, 2.8184e+00, 3.9811e+00, 5.6234e+00, 7.9433e+00, 1.0000e+01, 1.5849e+01, 2.2387e+01, 3.1623e+01, 4.4668e+01, 6.3096e+01, 8.9125e+01, 1.2589e+02, 1.7783e+02, 2.5119e+02, 3.5481e+02, 5.0119e+02, 7.0795e+02, 1.0000e+03])
- wavelength(wavelength)float640.9999 0.9998 0.9997 ... 0.3001 0.3
- units :
- micron
array([0.999907, 0.999807, 0.999707, ..., 0.300084, 0.300054, 0.300024])
- temperature(pressure)float64150.9 149.7 ... 1.525e+03 1.675e+03
- units :
- Kelvin
array([ 150.87, 149.68, 148.4 , 147. , 145.46, 143.89, 141.99, 140.28, 139.18, 138.29, 136.86, 134.53, 130.78, 122.56, 121.14, 145.1 , 156.53, 103.12, 94.06, 92.71, 103.98, 111.22, 126.88, 152.23, 144.24, 135.54, 120.64, 107.33, 107.5 , 116.89, 124.44, 127.07, 127.52, 128.96, 132.48, 136.92, 145.5 , 156.96, 175.53, 196.07, 218.7 , 243.55, 270.86, 300.88, 333.96, 370.48, 410.84, 440.11, 504.97, 559.74, 620.39, 687.48, 761.64, 843.45, 933.51, 1032.37, 1140.53, 1258.45, 1386.6 , 1525.37, 1675.21])
- albedo(wavelength)float640.2648 0.2092 ... 0.7193 0.7193
- units :
- none
array([0.2648301 , 0.20915205, 0.21309161, ..., 0.71930807, 0.71930484, 0.71930034])
- fpfs_reflected(wavelength)float646.048e-08 4.777e-08 ... 1.643e-07
- units :
- erg/cm**2/s/cm/(erg/cm**2/s/cm)
array([6.04827334e-08, 4.77668052e-08, 4.86665334e-08, ..., 1.64277847e-07, 1.64277110e-07, 1.64276082e-07])
- e-(pressure)float644.5e-38 4.5e-38 ... 2.15e-11
- units :
- v/v
array([4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 1.03e-37, 2.28e-34, 8.66e-31, 2.85e-27, 4.07e-24, 4.21e-21, 2.03e-18, 2.47e-16, 7.49e-15, 9.87e-14, 7.56e-13, 4.46e-12, 2.15e-11])
- H2(pressure)float640.837 0.837 0.837 ... 0.836 0.836
- units :
- v/v
array([0.837, 0.837, 0.837, 0.837, 0.837, 0.837, 0.837, 0.837, 0.837, 0.837, 0.837, 0.837, 0.837, 0.837, 0.837, 0.837, 0.837, 0.837, 0.837, 0.837, 0.837, 0.837, 0.837, 0.837, 0.837, 0.837, 0.837, 0.837, 0.837, 0.837, 0.837, 0.837, 0.837, 0.837, 0.837, 0.837, 0.837, 0.837, 0.837, 0.837, 0.837, 0.836, 0.836, 0.836, 0.836, 0.836, 0.836, 0.836, 0.836, 0.836, 0.836, 0.836, 0.836, 0.836, 0.836, 0.836, 0.836, 0.836, 0.836, 0.836, 0.836])
- H(pressure)float644.5e-38 4.5e-38 ... 3.37e-06
- units :
- v/v
array([4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 2.96e-36, 1.40e-32, 2.74e-29, 2.45e-26, 1.55e-24, 2.69e-21, 3.77e-19, 3.21e-17, 1.74e-15, 6.27e-14, 1.56e-12, 2.76e-11, 3.60e-10, 3.55e-09, 2.72e-08, 1.65e-07, 8.18e-07, 3.37e-06])
- H+(pressure)float644.5e-38 4.5e-38 ... 1.88e-37
- units :
- v/v
array([4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 1.88e-37])
- H-(pressure)float644.5e-38 4.5e-38 ... 8.8e-14
- units :
- v/v
array([4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 8.21e-36, 7.36e-32, 6.79e-28, 2.35e-24, 1.67e-21, 2.45e-19, 1.33e-17, 3.61e-16, 6.61e-15, 8.80e-14])
- VO(pressure)float644.5e-38 4.5e-38 ... 4.3e-12
- units :
- v/v
array([4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 5.95e-37, 6.25e-34, 6.23e-30, 2.28e-26, 3.36e-23, 2.17e-20, 6.63e-18, 1.04e-15, 8.83e-14, 4.30e-12])
- TiO(pressure)float644.5e-38 4.5e-38 ... 2.49e-13
- units :
- v/v
array([4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 6.49e-37, 1.36e-33, 1.55e-29, 6.29e-26, 1.00e-22, 6.86e-20, 2.19e-17, 3.45e-15, 2.49e-13])
- CO2(pressure)float644.5e-38 4.5e-38 ... 3.85e-10
- units :
- v/v
array([4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.18e-37, 3.35e-33, 2.56e-30, 8.97e-28, 1.60e-25, 1.55e-23, 2.74e-22, 4.76e-20, 1.44e-18, 3.06e-17, 4.68e-16, 5.29e-15, 4.55e-14, 3.01e-13, 1.57e-12, 6.53e-12, 2.20e-11, 6.16e-11, 1.51e-10, 3.85e-10])
- He(pressure)float640.163 0.163 0.163 ... 0.162 0.162
- units :
- v/v
array([0.163, 0.163, 0.163, 0.163, 0.163, 0.163, 0.163, 0.163, 0.163, 0.163, 0.163, 0.163, 0.163, 0.163, 0.163, 0.163, 0.163, 0.163, 0.163, 0.163, 0.163, 0.163, 0.163, 0.163, 0.163, 0.163, 0.163, 0.163, 0.163, 0.163, 0.163, 0.163, 0.163, 0.163, 0.163, 0.163, 0.163, 0.163, 0.163, 0.163, 0.163, 0.163, 0.162, 0.162, 0.162, 0.162, 0.162, 0.162, 0.162, 0.162, 0.162, 0.162, 0.162, 0.162, 0.162, 0.162, 0.162, 0.162, 0.162, 0.162, 0.162])
- H2O(pressure)float646.9e-05 3.53e-05 ... 0.000887
- units :
- v/v
array([6.90e-05, 3.53e-05, 1.75e-05, 8.34e-06, 3.79e-06, 1.69e-06, 6.72e-07, 2.80e-07, 1.40e-07, 7.45e-08, 3.31e-08, 1.07e-08, 2.03e-09, 6.07e-11, 2.38e-11, 7.64e-08, 1.21e-06, 1.15e-15, 2.55e-18, 6.92e-19, 6.67e-16, 2.25e-14, 1.51e-11, 3.53e-08, 2.65e-09, 1.20e-10, 3.06e-13, 3.80e-16, 2.94e-16, 2.10e-14, 3.66e-13, 7.25e-13, 6.09e-13, 7.40e-13, 1.87e-12, 5.99e-12, 6.07e-11, 9.51e-10, 4.32e-08, 1.22e-06, 2.24e-05, 2.83e-04, 8.19e-04, 8.19e-04, 8.19e-04, 8.19e-04, 8.19e-04, 8.19e-04, 8.19e-04, 8.19e-04, 8.19e-04, 8.19e-04, 8.19e-04, 8.19e-04, 8.19e-04, 8.20e-04, 8.20e-04, 8.20e-04, 8.21e-04, 8.32e-04, 8.87e-04])
- CH4(pressure)float640.000466 0.000466 ... 0.00054
- units :
- v/v
array([0.000466, 0.000466, 0.000466, 0.000466, 0.000466, 0.000466, 0.000466, 0.000466, 0.000466, 0.000466, 0.000466, 0.000466, 0.000466, 0.000466, 0.000466, 0.000466, 0.000466, 0.000466, 0.000466, 0.000466, 0.000466, 0.000466, 0.000466, 0.000466, 0.000466, 0.000466, 0.000466, 0.000466, 0.000466, 0.000466, 0.000466, 0.000466, 0.000466, 0.000466, 0.000466, 0.000466, 0.000466, 0.000466, 0.000466, 0.000466, 0.000466, 0.000466, 0.000465, 0.000465, 0.000465, 0.000465, 0.000465, 0.000465, 0.000465, 0.000465, 0.000465, 0.000465, 0.000465, 0.000465, 0.000465, 0.000465, 0.000465, 0.000465, 0.000466, 0.000478, 0.00054 ])
- CO(pressure)float644.5e-38 4.5e-38 ... 1.2e-06
- units :
- v/v
array([4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 6.16e-36, 2.90e-32, 5.16e-29, 3.93e-26, 1.42e-23, 5.49e-22, 3.90e-19, 2.97e-17, 1.43e-15, 4.50e-14, 9.67e-13, 1.46e-11, 1.58e-10, 1.27e-09, 7.69e-09, 3.59e-08, 1.33e-07, 4.09e-07, 1.20e-06])
- NH3(pressure)float640.000137 0.000137 ... 0.000155
- units :
- v/v
array([1.37e-04, 1.37e-04, 1.37e-04, 1.37e-04, 1.37e-04, 1.37e-04, 1.37e-04, 1.37e-04, 1.37e-04, 1.37e-04, 1.37e-04, 1.37e-04, 1.37e-04, 1.37e-04, 1.37e-04, 1.37e-04, 1.37e-04, 8.94e-06, 2.08e-07, 8.35e-08, 4.25e-06, 2.64e-05, 1.37e-04, 1.37e-04, 1.37e-04, 1.37e-04, 6.07e-05, 1.14e-06, 8.47e-07, 9.22e-06, 4.05e-05, 5.24e-05, 4.23e-05, 4.27e-05, 5.33e-05, 8.14e-05, 1.31e-04, 1.37e-04, 1.37e-04, 1.37e-04, 1.37e-04, 1.37e-04, 1.37e-04, 1.37e-04, 1.37e-04, 1.37e-04, 1.37e-04, 1.37e-04, 1.37e-04, 1.37e-04, 1.37e-04, 1.37e-04, 1.35e-04, 1.33e-04, 1.31e-04, 1.24e-04, 1.20e-04, 1.23e-04, 1.32e-04, 1.43e-04, 1.55e-04])
- N2(pressure)float645.42e-17 1.58e-17 ... 2.73e-05
- units :
- v/v
array([5.42e-17, 1.58e-17, 4.37e-18, 1.14e-18, 2.74e-19, 6.35e-20, 1.23e-20, 2.55e-21, 7.23e-22, 2.26e-22, 5.24e-23, 7.22e-24, 4.12e-25, 1.14e-27, 2.17e-28, 1.15e-22, 1.01e-20, 4.32e-37, 4.51e-38, 4.51e-38, 1.37e-37, 1.94e-36, 3.82e-29, 1.26e-23, 1.51e-25, 8.02e-28, 6.85e-34, 4.51e-38, 4.51e-38, 8.60e-37, 4.92e-34, 3.30e-33, 9.88e-34, 7.08e-34, 5.76e-33, 1.32e-31, 3.02e-29, 6.06e-27, 3.24e-24, 8.30e-22, 1.09e-19, 8.01e-18, 3.54e-16, 1.01e-14, 2.00e-13, 2.84e-12, 3.02e-11, 1.26e-10, 1.61e-09, 8.49e-09, 3.66e-08, 1.31e-07, 3.89e-07, 9.83e-07, 2.11e-06, 3.72e-06, 5.96e-06, 9.44e-06, 1.47e-05, 2.09e-05, 2.73e-05])
- PH3(pressure)float644.5e-38 4.5e-38 ... 5.03e-07
- units :
- v/v
array([4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.55e-38, 5.45e-35, 5.34e-31, 1.60e-27, 2.05e-24, 1.25e-21, 6.18e-20, 6.85e-17, 7.03e-15, 4.51e-13, 1.89e-11, 5.41e-10, 1.08e-08, 1.33e-07, 4.94e-07, 5.33e-07, 5.28e-07, 5.20e-07, 5.11e-07, 5.03e-07])
- opd(pressure_layer, wavenumber_layer)float640.0 0.0 0.0 0.0 ... 0.0 0.0 0.0 0.0
- units :
- unitless
array([[0. , 0. , 0. , ..., 0. , 0. , 0. ], [0. , 0. , 0. , ..., 0. , 0. , 0. ], [0. , 0. , 0. , ..., 0. , 0. , 0. ], ..., [0.59295, 0.67848, 0.67467, ..., 0.47412, 0.47386, 0.47365], [0.29647, 0.33924, 0.33733, ..., 0.23706, 0.23693, 0.23683], [0. , 0. , 0. , ..., 0. , 0. , 0. ]])
- ssa(pressure_layer, wavenumber_layer)float640.0 0.0 0.0 0.0 ... 0.0 0.0 0.0 0.0
- units :
- unitless
array([[0. , 0. , 0. , ..., 0. , 0. , 0. ], [0. , 0. , 0. , ..., 0. , 0. , 0. ], [0. , 0. , 0. , ..., 0. , 0. , 0. ], ..., [0.948187 , 0.9617072, 0.818675 , ..., 0.5906784, 0.5935144, 0.5974549], [0.948187 , 0.9617072, 0.818675 , ..., 0.5906784, 0.5935144, 0.5974549], [0. , 0. , 0. , ..., 0. , 0. , 0. ]])
- asy(pressure_layer, wavenumber_layer)float640.0 0.0 0.0 0.0 ... 0.0 0.0 0.0 0.0
- units :
- unitless
array([[0. , 0. , 0. , ..., 0. , 0. , 0. ], [0. , 0. , 0. , ..., 0. , 0. , 0. ], [0. , 0. , 0. , ..., 0. , 0. , 0. ], ..., [0.3800323, 0.3667444, 0.4358756, ..., 0.8781489, 0.8735224, 0.8670378], [0.3800323, 0.3667444, 0.4358756, ..., 0.8781489, 0.8735224, 0.8670378], [0. , 0. , 0. , ..., 0. , 0. , 0. ]])
- pressurePandasIndex
PandasIndex(Index([ 1e-06, 1.41e-06, 2e-06, 2.82e-06, 3.98e-06, 5.62e-06, 7.94e-06, 1.122e-05, 1.585e-05, 2.239e-05, 3.162e-05, 4.467e-05, 6.31e-05, 8.912e-05, 0.00012589, 0.00017783, 0.00025119, 0.00035481, 0.00050119, 0.00070795, 0.001, 0.0014125, 0.0019953, 0.0028184, 0.0039811, 0.0056234, 0.0079433, 0.01122, 0.015849, 0.022387, 0.031623, 0.044668, 0.063096, 0.089125, 0.12589, 0.17783, 0.25119, 0.35481, 0.5011899999999999, 0.70795, 1.0, 1.4125, 1.9953, 2.8184, 3.9811, 5.6234, 7.9433, 10.0, 15.849, 22.387, 31.623, 44.668, 63.096, 89.125, 125.89, 177.83, 251.19, 354.81, 501.19, 707.95, 1000.0], dtype='float64', name='pressure'))
- wavelengthPandasIndex
PandasIndex(Index([ 0.999907199934528, 0.9998072142139076, 0.9997072384913596, 0.9996072727658833, 0.9995073170364811, 0.9994073713021508, 0.9993074355618947, 0.9992075098147134, 0.999107594059607, 0.9990076882955765, ... 0.30029414410706, 0.30026411619407106, 0.3002340912837236, 0.3002040693757168, 0.30017405046975054, 0.3001440345655245, 0.3001140216627392, 0.3000840117610942, 0.3000540048602893, 0.3000240009600245], dtype='float64', name='wavelength', length=12039))
- planet_params :
- {"mp": {"value": 1.8981245973360504e+30, "unit": "g"}, "rp": {"value": 7149200000.0, "unit": "cm"}, "p_reference": {"value": 1.0, "unit": "bar"}}
- stellar_params :
- {"database": "ck04models", "steff": 5000, "feh": 0, "logg": 4, "rs": {"value": 69570000000.0, "unit": "cm"}}
- orbit_params :
- {"sma": {"value": 14959787070000.0, "unit": "cm"}}
Reusing your run with xarray¶
We often revisit models time and time again. Maybe you want to change your wavelength range, or observing geometry (e.g. transit vs emission).
Here we show a simple example of how to use the input_xarray
function for reuse. In this simple example, we will take our previous input but instead of computing reflected light, add thermal emission and transmission
[5]:
opa = jdi.opannection(wave_range=[0.3,14])
ds = jdi.xr.load_dataset('/data/picaso_dbs/model.nc')
reuse = jdi.input_xarray(ds, opa)
new_model = reuse.spectrum(opa, calculation='reflected+thermal', full_output=True)
[6]:
new_output= jdi.output_xarray(new_model, reuse)
[7]:
new_output
[7]:
<xarray.Dataset> Size: 1MB Dimensions: (pressure: 61, wavelength: 28717, pressure_layer: 60, wavenumber_layer: 196) Coordinates: * pressure (pressure) float64 488B 1e-06 1.41e-06 2e-06 ... 708.0 1e+03 * wavelength (wavelength) float64 230kB 5.3 5.299 5.299 ... 0.3001 0.3 Dimensions without coordinates: pressure_layer, wavenumber_layer Data variables: (12/22) temperature (pressure) float64 488B 150.9 149.7 ... 1.525e+03 1.675e+03 flux_emission (wavelength) float64 230kB 1.838e+05 1.838e+05 ... 1.006e-08 fpfs_emission (wavelength) float64 230kB 1.774e-09 1.774e-09 ... 2.899e-24 albedo (wavelength) float64 230kB 0.09834 0.09837 ... 0.7193 0.7193 e- (pressure) float64 488B 4.5e-38 4.5e-38 ... 4.46e-12 2.15e-11 H2 (pressure) float64 488B 0.837 0.837 0.837 ... 0.836 0.836 ... ... NH3 (pressure) float64 488B 0.000137 0.000137 ... 0.000155 N2 (pressure) float64 488B 5.42e-17 1.58e-17 ... 2.73e-05 PH3 (pressure) float64 488B 4.5e-38 4.5e-38 ... 5.11e-07 5.03e-07 opd (pressure_layer, wavenumber_layer) float64 94kB 0.0 ... 0.0 ssa (pressure_layer, wavenumber_layer) float64 94kB 0.0 ... 0.0 asy (pressure_layer, wavenumber_layer) float64 94kB 0.0 ... 0.0 Attributes: planet_params: {"effective_temp": 19.572917134869886, "mp": {"value": 1... stellar_params: {"database": "ck04models", "steff": 5000, "feh": 0, "log...
- pressure: 61
- wavelength: 28717
- pressure_layer: 60
- wavenumber_layer: 196
- pressure(pressure)float641e-06 1.41e-06 ... 708.0 1e+03
- units :
- bar
array([1.0000e-06, 1.4100e-06, 2.0000e-06, 2.8200e-06, 3.9800e-06, 5.6200e-06, 7.9400e-06, 1.1220e-05, 1.5850e-05, 2.2390e-05, 3.1620e-05, 4.4670e-05, 6.3100e-05, 8.9120e-05, 1.2589e-04, 1.7783e-04, 2.5119e-04, 3.5481e-04, 5.0119e-04, 7.0795e-04, 1.0000e-03, 1.4125e-03, 1.9953e-03, 2.8184e-03, 3.9811e-03, 5.6234e-03, 7.9433e-03, 1.1220e-02, 1.5849e-02, 2.2387e-02, 3.1623e-02, 4.4668e-02, 6.3096e-02, 8.9125e-02, 1.2589e-01, 1.7783e-01, 2.5119e-01, 3.5481e-01, 5.0119e-01, 7.0795e-01, 1.0000e+00, 1.4125e+00, 1.9953e+00, 2.8184e+00, 3.9811e+00, 5.6234e+00, 7.9433e+00, 1.0000e+01, 1.5849e+01, 2.2387e+01, 3.1623e+01, 4.4668e+01, 6.3096e+01, 8.9125e+01, 1.2589e+02, 1.7783e+02, 2.5119e+02, 3.5481e+02, 5.0119e+02, 7.0795e+02, 1.0000e+03])
- wavelength(wavelength)float645.3 5.299 5.299 ... 0.3001 0.3
- units :
- micron
array([5.300002, 5.299472, 5.298942, ..., 0.300084, 0.300054, 0.300024])
- temperature(pressure)float64150.9 149.7 ... 1.525e+03 1.675e+03
- units :
- Kelvin
array([ 150.87, 149.68, 148.4 , 147. , 145.46, 143.89, 141.99, 140.28, 139.18, 138.29, 136.86, 134.53, 130.78, 122.56, 121.14, 145.1 , 156.53, 103.12, 94.06, 92.71, 103.98, 111.22, 126.88, 152.23, 144.24, 135.54, 120.64, 107.33, 107.5 , 116.89, 124.44, 127.07, 127.52, 128.96, 132.48, 136.92, 145.5 , 156.96, 175.53, 196.07, 218.7 , 243.55, 270.86, 300.88, 333.96, 370.48, 410.84, 440.11, 504.97, 559.74, 620.39, 687.48, 761.64, 843.45, 933.51, 1032.37, 1140.53, 1258.45, 1386.6 , 1525.37, 1675.21])
- flux_emission(wavelength)float641.838e+05 1.838e+05 ... 1.006e-08
- units :
- erg/cm**2/s/cm
array([1.83789348e+05, 1.83837129e+05, 1.83095852e+05, ..., 1.01901023e-08, 1.01285742e-08, 1.00634510e-08])
- fpfs_emission(wavelength)float641.774e-09 1.774e-09 ... 2.899e-24
- units :
- erg/cm**2/s/cm/(erg/cm**2/s/cm)
array([1.77410300e-09, 1.77404851e-09, 1.76638173e-09, ..., 2.94936760e-24, 2.92482714e-24, 2.89936270e-24])
- albedo(wavelength)float640.09834 0.09837 ... 0.7193 0.7193
- units :
- none
array([0.09834002, 0.09837396, 0.09821968, ..., 0.71930807, 0.71930484, 0.71930034])
- e-(pressure)float644.5e-38 4.5e-38 ... 2.15e-11
- units :
- v/v
array([4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 1.03e-37, 2.28e-34, 8.66e-31, 2.85e-27, 4.07e-24, 4.21e-21, 2.03e-18, 2.47e-16, 7.49e-15, 9.87e-14, 7.56e-13, 4.46e-12, 2.15e-11])
- H2(pressure)float640.837 0.837 0.837 ... 0.836 0.836
- units :
- v/v
array([0.837, 0.837, 0.837, 0.837, 0.837, 0.837, 0.837, 0.837, 0.837, 0.837, 0.837, 0.837, 0.837, 0.837, 0.837, 0.837, 0.837, 0.837, 0.837, 0.837, 0.837, 0.837, 0.837, 0.837, 0.837, 0.837, 0.837, 0.837, 0.837, 0.837, 0.837, 0.837, 0.837, 0.837, 0.837, 0.837, 0.837, 0.837, 0.837, 0.837, 0.837, 0.836, 0.836, 0.836, 0.836, 0.836, 0.836, 0.836, 0.836, 0.836, 0.836, 0.836, 0.836, 0.836, 0.836, 0.836, 0.836, 0.836, 0.836, 0.836, 0.836])
- H(pressure)float644.5e-38 4.5e-38 ... 3.37e-06
- units :
- v/v
array([4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 2.96e-36, 1.40e-32, 2.74e-29, 2.45e-26, 1.55e-24, 2.69e-21, 3.77e-19, 3.21e-17, 1.74e-15, 6.27e-14, 1.56e-12, 2.76e-11, 3.60e-10, 3.55e-09, 2.72e-08, 1.65e-07, 8.18e-07, 3.37e-06])
- H+(pressure)float644.5e-38 4.5e-38 ... 1.88e-37
- units :
- v/v
array([4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 1.88e-37])
- H-(pressure)float644.5e-38 4.5e-38 ... 8.8e-14
- units :
- v/v
array([4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 8.21e-36, 7.36e-32, 6.79e-28, 2.35e-24, 1.67e-21, 2.45e-19, 1.33e-17, 3.61e-16, 6.61e-15, 8.80e-14])
- VO(pressure)float644.5e-38 4.5e-38 ... 4.3e-12
- units :
- v/v
array([4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 5.95e-37, 6.25e-34, 6.23e-30, 2.28e-26, 3.36e-23, 2.17e-20, 6.63e-18, 1.04e-15, 8.83e-14, 4.30e-12])
- TiO(pressure)float644.5e-38 4.5e-38 ... 2.49e-13
- units :
- v/v
array([4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 6.49e-37, 1.36e-33, 1.55e-29, 6.29e-26, 1.00e-22, 6.86e-20, 2.19e-17, 3.45e-15, 2.49e-13])
- CO2(pressure)float644.5e-38 4.5e-38 ... 3.85e-10
- units :
- v/v
array([4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.18e-37, 3.35e-33, 2.56e-30, 8.97e-28, 1.60e-25, 1.55e-23, 2.74e-22, 4.76e-20, 1.44e-18, 3.06e-17, 4.68e-16, 5.29e-15, 4.55e-14, 3.01e-13, 1.57e-12, 6.53e-12, 2.20e-11, 6.16e-11, 1.51e-10, 3.85e-10])
- He(pressure)float640.163 0.163 0.163 ... 0.162 0.162
- units :
- v/v
array([0.163, 0.163, 0.163, 0.163, 0.163, 0.163, 0.163, 0.163, 0.163, 0.163, 0.163, 0.163, 0.163, 0.163, 0.163, 0.163, 0.163, 0.163, 0.163, 0.163, 0.163, 0.163, 0.163, 0.163, 0.163, 0.163, 0.163, 0.163, 0.163, 0.163, 0.163, 0.163, 0.163, 0.163, 0.163, 0.163, 0.163, 0.163, 0.163, 0.163, 0.163, 0.163, 0.162, 0.162, 0.162, 0.162, 0.162, 0.162, 0.162, 0.162, 0.162, 0.162, 0.162, 0.162, 0.162, 0.162, 0.162, 0.162, 0.162, 0.162, 0.162])
- H2O(pressure)float646.9e-05 3.53e-05 ... 0.000887
- units :
- v/v
array([6.90e-05, 3.53e-05, 1.75e-05, 8.34e-06, 3.79e-06, 1.69e-06, 6.72e-07, 2.80e-07, 1.40e-07, 7.45e-08, 3.31e-08, 1.07e-08, 2.03e-09, 6.07e-11, 2.38e-11, 7.64e-08, 1.21e-06, 1.15e-15, 2.55e-18, 6.92e-19, 6.67e-16, 2.25e-14, 1.51e-11, 3.53e-08, 2.65e-09, 1.20e-10, 3.06e-13, 3.80e-16, 2.94e-16, 2.10e-14, 3.66e-13, 7.25e-13, 6.09e-13, 7.40e-13, 1.87e-12, 5.99e-12, 6.07e-11, 9.51e-10, 4.32e-08, 1.22e-06, 2.24e-05, 2.83e-04, 8.19e-04, 8.19e-04, 8.19e-04, 8.19e-04, 8.19e-04, 8.19e-04, 8.19e-04, 8.19e-04, 8.19e-04, 8.19e-04, 8.19e-04, 8.19e-04, 8.19e-04, 8.20e-04, 8.20e-04, 8.20e-04, 8.21e-04, 8.32e-04, 8.87e-04])
- CH4(pressure)float640.000466 0.000466 ... 0.00054
- units :
- v/v
array([0.000466, 0.000466, 0.000466, 0.000466, 0.000466, 0.000466, 0.000466, 0.000466, 0.000466, 0.000466, 0.000466, 0.000466, 0.000466, 0.000466, 0.000466, 0.000466, 0.000466, 0.000466, 0.000466, 0.000466, 0.000466, 0.000466, 0.000466, 0.000466, 0.000466, 0.000466, 0.000466, 0.000466, 0.000466, 0.000466, 0.000466, 0.000466, 0.000466, 0.000466, 0.000466, 0.000466, 0.000466, 0.000466, 0.000466, 0.000466, 0.000466, 0.000466, 0.000465, 0.000465, 0.000465, 0.000465, 0.000465, 0.000465, 0.000465, 0.000465, 0.000465, 0.000465, 0.000465, 0.000465, 0.000465, 0.000465, 0.000465, 0.000465, 0.000466, 0.000478, 0.00054 ])
- CO(pressure)float644.5e-38 4.5e-38 ... 1.2e-06
- units :
- v/v
array([4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 6.16e-36, 2.90e-32, 5.16e-29, 3.93e-26, 1.42e-23, 5.49e-22, 3.90e-19, 2.97e-17, 1.43e-15, 4.50e-14, 9.67e-13, 1.46e-11, 1.58e-10, 1.27e-09, 7.69e-09, 3.59e-08, 1.33e-07, 4.09e-07, 1.20e-06])
- NH3(pressure)float640.000137 0.000137 ... 0.000155
- units :
- v/v
array([1.37e-04, 1.37e-04, 1.37e-04, 1.37e-04, 1.37e-04, 1.37e-04, 1.37e-04, 1.37e-04, 1.37e-04, 1.37e-04, 1.37e-04, 1.37e-04, 1.37e-04, 1.37e-04, 1.37e-04, 1.37e-04, 1.37e-04, 8.94e-06, 2.08e-07, 8.35e-08, 4.25e-06, 2.64e-05, 1.37e-04, 1.37e-04, 1.37e-04, 1.37e-04, 6.07e-05, 1.14e-06, 8.47e-07, 9.22e-06, 4.05e-05, 5.24e-05, 4.23e-05, 4.27e-05, 5.33e-05, 8.14e-05, 1.31e-04, 1.37e-04, 1.37e-04, 1.37e-04, 1.37e-04, 1.37e-04, 1.37e-04, 1.37e-04, 1.37e-04, 1.37e-04, 1.37e-04, 1.37e-04, 1.37e-04, 1.37e-04, 1.37e-04, 1.37e-04, 1.35e-04, 1.33e-04, 1.31e-04, 1.24e-04, 1.20e-04, 1.23e-04, 1.32e-04, 1.43e-04, 1.55e-04])
- N2(pressure)float645.42e-17 1.58e-17 ... 2.73e-05
- units :
- v/v
array([5.42e-17, 1.58e-17, 4.37e-18, 1.14e-18, 2.74e-19, 6.35e-20, 1.23e-20, 2.55e-21, 7.23e-22, 2.26e-22, 5.24e-23, 7.22e-24, 4.12e-25, 1.14e-27, 2.17e-28, 1.15e-22, 1.01e-20, 4.32e-37, 4.51e-38, 4.51e-38, 1.37e-37, 1.94e-36, 3.82e-29, 1.26e-23, 1.51e-25, 8.02e-28, 6.85e-34, 4.51e-38, 4.51e-38, 8.60e-37, 4.92e-34, 3.30e-33, 9.88e-34, 7.08e-34, 5.76e-33, 1.32e-31, 3.02e-29, 6.06e-27, 3.24e-24, 8.30e-22, 1.09e-19, 8.01e-18, 3.54e-16, 1.01e-14, 2.00e-13, 2.84e-12, 3.02e-11, 1.26e-10, 1.61e-09, 8.49e-09, 3.66e-08, 1.31e-07, 3.89e-07, 9.83e-07, 2.11e-06, 3.72e-06, 5.96e-06, 9.44e-06, 1.47e-05, 2.09e-05, 2.73e-05])
- PH3(pressure)float644.5e-38 4.5e-38 ... 5.03e-07
- units :
- v/v
array([4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.51e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.50e-38, 4.55e-38, 5.45e-35, 5.34e-31, 1.60e-27, 2.05e-24, 1.25e-21, 6.18e-20, 6.85e-17, 7.03e-15, 4.51e-13, 1.89e-11, 5.41e-10, 1.08e-08, 1.33e-07, 4.94e-07, 5.33e-07, 5.28e-07, 5.20e-07, 5.11e-07, 5.03e-07])
- opd(pressure_layer, wavenumber_layer)float640.0 0.0 0.0 0.0 ... 0.0 0.0 0.0 0.0
- units :
- unitless
array([[0. , 0. , 0. , ..., 0. , 0. , 0. ], [0. , 0. , 0. , ..., 0. , 0. , 0. ], [0. , 0. , 0. , ..., 0. , 0. , 0. ], ..., [0.59295, 0.67848, 0.67467, ..., 0.47412, 0.47386, 0.47365], [0.29647, 0.33924, 0.33733, ..., 0.23706, 0.23693, 0.23683], [0. , 0. , 0. , ..., 0. , 0. , 0. ]])
- ssa(pressure_layer, wavenumber_layer)float640.0 0.0 0.0 0.0 ... 0.0 0.0 0.0 0.0
- units :
- unitless
array([[0. , 0. , 0. , ..., 0. , 0. , 0. ], [0. , 0. , 0. , ..., 0. , 0. , 0. ], [0. , 0. , 0. , ..., 0. , 0. , 0. ], ..., [0.948187 , 0.9617072, 0.818675 , ..., 0.5906784, 0.5935144, 0.5974549], [0.948187 , 0.9617072, 0.818675 , ..., 0.5906784, 0.5935144, 0.5974549], [0. , 0. , 0. , ..., 0. , 0. , 0. ]])
- asy(pressure_layer, wavenumber_layer)float640.0 0.0 0.0 0.0 ... 0.0 0.0 0.0 0.0
- units :
- unitless
array([[0. , 0. , 0. , ..., 0. , 0. , 0. ], [0. , 0. , 0. , ..., 0. , 0. , 0. ], [0. , 0. , 0. , ..., 0. , 0. , 0. ], ..., [0.3800323, 0.3667444, 0.4358756, ..., 0.8781489, 0.8735224, 0.8670378], [0.3800323, 0.3667444, 0.4358756, ..., 0.8781489, 0.8735224, 0.8670378], [0. , 0. , 0. , ..., 0. , 0. , 0. ]])
- pressurePandasIndex
PandasIndex(Index([ 1e-06, 1.41e-06, 2e-06, 2.82e-06, 3.98e-06, 5.62e-06, 7.94e-06, 1.122e-05, 1.585e-05, 2.239e-05, 3.162e-05, 4.467e-05, 6.31e-05, 8.912e-05, 0.00012589, 0.00017783, 0.00025119, 0.00035481, 0.00050119, 0.00070795, 0.001, 0.0014125, 0.0019953, 0.0028184, 0.0039811, 0.0056234, 0.0079433, 0.01122, 0.015849, 0.022387, 0.031623, 0.044668, 0.063096, 0.089125, 0.12589, 0.17783, 0.25119, 0.35481, 0.5011899999999999, 0.70795, 1.0, 1.4125, 1.9953, 2.8184, 3.9811, 5.6234, 7.9433, 10.0, 15.849, 22.387, 31.623, 44.668, 63.096, 89.125, 125.89, 177.83, 251.19, 354.81, 501.19, 707.95, 1000.0], dtype='float64', name='pressure'))
- wavelengthPandasIndex
PandasIndex(Index([ 5.300001987536683, 5.29947201383708, 5.298942093132193, 5.298412225416728, 5.297882410685383, 5.297352648932866, 5.296822940153871, 5.296293284343108, 5.295763681495279, 5.295234131605083, ... 0.30029414410706, 0.30026411619407106, 0.3002340912837236, 0.3002040693757168, 0.30017405046975054, 0.3001440345655245, 0.3001140216627392, 0.3000840117610942, 0.3000540048602893, 0.3000240009600245], dtype='float64', name='wavelength', length=28717))
- planet_params :
- {"effective_temp": 19.572917134869886, "mp": {"value": 1.8981245973360504e+30, "unit": "g"}, "rp": {"value": 7149200000.0, "unit": "cm"}, "p_reference": {"value": 1.0, "unit": "bar"}}
- stellar_params :
- {"database": "ck04models", "steff": 5000, "feh": 0, "logg": 4, "rs": {"value": 69570000000.0, "unit": "cm"}}
Adding meta data¶
We have a tutorial that shows what some recommended meta data might be. output_xarray
has a field add_output
that allows you to add extra arguments to the xarray.Dataset.attrs
. However, it has to be in a very specific format.
You can see the basic format by running this function:
[8]:
jdi.standard_metadata()
[8]:
{'author': 'optional',
'contact': 'optional',
'code': 'optional',
'doi': 'optional',
'planet_params': {'rp': 'usually taken from picaso',
'mp': 'usually taken from picaso',
'mh': 'optional',
'cto': 'optional',
'heat_redis': 'optional',
'p_reference': 'usually taken from picaso',
'tint': 'optional'},
'stellar_params': {'logg': 'usually taken from picaso',
'feh': 'usually taken from picaso',
'steff': 'usually taken from picaso',
'rs': 'usually taken from picaso',
'ms': 'optional',
'database': 'usually taken from picaso'},
'orbit_params': {'sma': 'usually taken from picaso'}}
[9]:
opa = jdi.opannection(wave_range=[0.3,1])
pl = jdi.inputs()#calculation='brown')
pl.gravity(radius=1, radius_unit=u.Unit('R_jup'),
mass=1, mass_unit=u.Unit('M_jup'))
mh = 0
cto = 1
pl.atmosphere(filename=jdi.jupiter_pt(), sep='\s+')
pl.chemeq_visscher(cto, mh)
pl.phase_angle(0)
pl.clouds(filename=jdi.jupiter_cld(), sep='\s+')
pl.star(opa, 5000,0,4, radius=1, radius_unit=u.Unit("R_sun"), semi_major=1, semi_major_unit=1*u.AU)
#MUST USE full output=True for this functionality
df= pl.spectrum(opa, calculation='reflected', full_output=True)
preserve = jdi.output_xarray(df,pl,add_output={
'author':"Awesome Scientist",
'contact' : "awesomescientist@universe.com",
'code' : "picaso",
'planet_params':{'mh':mh, 'co':cto},
'cloud_params':{'fsed':3}
})
[10]:
preserve
[10]:
<xarray.Dataset> Size: 589kB Dimensions: (pressure: 61, wavelength: 12039, pressure_layer: 60, wavenumber_layer: 196) Coordinates: * pressure (pressure) float64 488B 1e-06 1.41e-06 2e-06 ... 708.0 1e+03 * wavelength (wavelength) float64 96kB 0.9999 0.9998 ... 0.3001 0.3 Dimensions without coordinates: pressure_layer, wavenumber_layer Data variables: (12/40) temperature (pressure) float64 488B 150.9 149.7 ... 1.525e+03 1.675e+03 albedo (wavelength) float64 96kB 0.2646 0.2089 ... 0.7193 0.7193 fpfs_reflected (wavelength) float64 96kB 6.043e-08 4.771e-08 ... 1.643e-07 e- (pressure) float64 488B 4.5e-38 4.5e-38 ... 2.278e-11 H2 (pressure) float64 488B 0.8359 0.8359 ... 0.8359 0.8358 H (pressure) float64 488B 4.5e-38 4.5e-38 ... 3.369e-06 ... ... MgH (pressure) float64 488B 4.5e-38 4.5e-38 ... 2.843e-11 OCS (pressure) float64 488B 4.5e-38 4.5e-38 ... 7.931e-12 graphite (pressure) float64 488B 4.5e-38 4.5e-38 ... 4.5e-38 4.5e-38 opd (pressure_layer, wavenumber_layer) float64 94kB 0.0 ... 0.0 ssa (pressure_layer, wavenumber_layer) float64 94kB 0.0 ... 0.0 asy (pressure_layer, wavenumber_layer) float64 94kB 0.0 ... 0.0 Attributes: author: Awesome Scientist code: picaso contact: awesomescientist@universe.com planet_params: {"mp": {"value": 1.8981245973360504e+30, "unit": "g"}, "... stellar_params: {"database": "ck04models", "steff": 5000, "feh": 0, "log... orbit_params: {"sma": {"value": 14959787070000.0, "unit": "cm"}} cloud_params: {'fsed': 3}
- pressure: 61
- wavelength: 12039
- pressure_layer: 60
- wavenumber_layer: 196
- pressure(pressure)float641e-06 1.41e-06 ... 708.0 1e+03
- units :
- bar
array([1.0000e-06, 1.4100e-06, 2.0000e-06, 2.8200e-06, 3.9800e-06, 5.6200e-06, 7.9400e-06, 1.1220e-05, 1.5850e-05, 2.2390e-05, 3.1620e-05, 4.4670e-05, 6.3100e-05, 8.9120e-05, 1.2589e-04, 1.7783e-04, 2.5119e-04, 3.5481e-04, 5.0119e-04, 7.0795e-04, 1.0000e-03, 1.4125e-03, 1.9953e-03, 2.8184e-03, 3.9811e-03, 5.6234e-03, 7.9433e-03, 1.1220e-02, 1.5849e-02, 2.2387e-02, 3.1623e-02, 4.4668e-02, 6.3096e-02, 8.9125e-02, 1.2589e-01, 1.7783e-01, 2.5119e-01, 3.5481e-01, 5.0119e-01, 7.0795e-01, 1.0000e+00, 1.4125e+00, 1.9953e+00, 2.8184e+00, 3.9811e+00, 5.6234e+00, 7.9433e+00, 1.0000e+01, 1.5849e+01, 2.2387e+01, 3.1623e+01, 4.4668e+01, 6.3096e+01, 8.9125e+01, 1.2589e+02, 1.7783e+02, 2.5119e+02, 3.5481e+02, 5.0119e+02, 7.0795e+02, 1.0000e+03])
- wavelength(wavelength)float640.9999 0.9998 0.9997 ... 0.3001 0.3
- units :
- micron
array([0.999907, 0.999807, 0.999707, ..., 0.300084, 0.300054, 0.300024])
- temperature(pressure)float64150.9 149.7 ... 1.525e+03 1.675e+03
- units :
- Kelvin
array([ 150.87, 149.68, 148.4 , 147. , 145.46, 143.89, 141.99, 140.28, 139.18, 138.29, 136.86, 134.53, 130.78, 122.56, 121.14, 145.1 , 156.53, 103.12, 94.06, 92.71, 103.98, 111.22, 126.88, 152.23, 144.24, 135.54, 120.64, 107.33, 107.5 , 116.89, 124.44, 127.07, 127.52, 128.96, 132.48, 136.92, 145.5 , 156.96, 175.53, 196.07, 218.7 , 243.55, 270.86, 300.88, 333.96, 370.48, 410.84, 440.11, 504.97, 559.74, 620.39, 687.48, 761.64, 843.45, 933.51, 1032.37, 1140.53, 1258.45, 1386.6 , 1525.37, 1675.21])
- albedo(wavelength)float640.2646 0.2089 ... 0.7193 0.7193
- units :
- none
array([0.26459774, 0.20888773, 0.21282954, ..., 0.71931084, 0.71930761, 0.71930311])
- fpfs_reflected(wavelength)float646.043e-08 4.771e-08 ... 1.643e-07
- units :
- erg/cm**2/s/cm/(erg/cm**2/s/cm)
array([6.04296671e-08, 4.77064382e-08, 4.86066824e-08, ..., 1.64278480e-07, 1.64277742e-07, 1.64276713e-07])
- e-(pressure)float644.5e-38 4.5e-38 ... 2.278e-11
- units :
- v/v
array([4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50060000e-38, 4.50060000e-38, 1.02883186e-37, 2.27414078e-34, 8.65636527e-31, 2.84538280e-27, 4.06788363e-24, 4.20503825e-21, 2.06732766e-18, 2.59831088e-16, 8.14244804e-15, 1.05826888e-13, 7.99085910e-13, 4.68631779e-12, 2.27829347e-11])
- H2(pressure)float640.8359 0.8359 ... 0.8359 0.8358
- units :
- v/v
array([0.83587277, 0.83587 , 0.83587 , 0.83587 , 0.83587 , 0.83587 , 0.83587 , 0.83587 , 0.83587 , 0.83587 , 0.83587 , 0.83587 , 0.83587 , 0.83587 , 0.83587 , 0.83587 , 0.83587108, 0.83587 , 0.83587 , 0.83587 , 0.83587 , 0.83587 , 0.83587 , 0.83587 , 0.83587 , 0.83587 , 0.83587 , 0.83587 , 0.83587 , 0.83587 , 0.83587 , 0.83587 , 0.83587 , 0.83587 , 0.83587 , 0.83587 , 0.83587 , 0.83587 , 0.83586433, 0.83584 , 0.83584 , 0.83584 , 0.83584 , 0.83611 , 0.83611 , 0.83611 , 0.83611 , 0.83611 , 0.83611 , 0.83611 , 0.83611 , 0.83611 , 0.83611196, 0.8361114 , 0.83610908, 0.83611421, 0.83611818, 0.83610398, 0.83609055, 0.83592727, 0.83576335])
- H(pressure)float644.5e-38 4.5e-38 ... 3.369e-06
- units :
- v/v
array([4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 2.95285732e-36, 1.39959350e-32, 2.73729402e-29, 2.45299869e-26, 1.54764519e-24, 2.68719363e-21, 3.76408590e-19, 3.20560345e-17, 1.73782033e-15, 6.26665283e-14, 1.55900523e-12, 2.76461654e-11, 3.60427397e-10, 3.55174222e-09, 2.71523945e-08, 1.65141593e-07, 8.17696496e-07, 3.36938878e-06])
- H+(pressure)float644.5e-38 4.5e-38 ... 1.805e-37
- units :
- v/v
array([4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50060000e-38, 4.50060000e-38, 4.50060000e-38, 4.50060000e-38, 4.50060000e-38, 4.50060000e-38, 4.50060000e-38, 4.50059557e-38, 4.50056950e-38, 4.50045220e-38, 4.50053139e-38, 4.50051779e-38, 4.50068765e-38, 4.50050000e-38, 1.80459872e-37])
- H-(pressure)float644.5e-38 4.5e-38 ... 9.306e-14
- units :
- v/v
array([4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50060000e-38, 4.50060000e-38, 4.50060000e-38, 4.50060000e-38, 4.50060000e-38, 8.21004331e-36, 7.35675799e-32, 6.78386272e-28, 2.38470543e-24, 1.74965033e-21, 2.66301810e-19, 1.42432626e-17, 3.81661388e-16, 6.94443618e-15, 9.30593118e-14])
- VO(pressure)float644.5e-38 4.5e-38 ... 3.492e-12
- units :
- v/v
array([4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50060000e-38, 4.50060000e-38, 4.50060000e-38, 4.50060000e-38, 4.50060000e-38, 5.94677373e-37, 6.24734876e-34, 6.22368618e-30, 2.28176685e-26, 3.36169853e-23, 2.16690836e-20, 6.60342890e-18, 1.01279908e-15, 7.94428453e-14, 3.49177473e-12])
- TiO(pressure)float644.5e-38 4.5e-38 ... 2.021e-13
- units :
- v/v
array([4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50060000e-38, 4.50060000e-38, 4.50060000e-38, 4.50060000e-38, 4.50060000e-38, 4.50060000e-38, 6.48684142e-37, 1.35853728e-33, 1.54878615e-29, 6.29302718e-26, 1.00251542e-22, 6.82354247e-20, 2.11726242e-17, 3.00346863e-15, 2.02143283e-13])
- CO2(pressure)float644.5e-38 4.5e-38 ... 2.746e-09
- units :
- v/v
array([4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.17042326e-37, 3.34592494e-33, 2.55370214e-30, 8.96269303e-28, 1.59776939e-25, 1.55234263e-23, 2.73513676e-22, 4.75774988e-20, 1.44290062e-18, 3.05928195e-17, 4.67120993e-16, 5.28918447e-15, 4.54724596e-14, 3.01183749e-13, 1.56997256e-12, 6.54597356e-12, 2.26263162e-11, 7.38654877e-11, 3.53995663e-10, 2.74603747e-09])
- He(pressure)float640.1624 0.1624 ... 0.1625 0.1625
- units :
- v/v
array([0.16244908, 0.16245 , 0.16245 , 0.16245 , 0.16245 , 0.16245 , 0.16245 , 0.16245 , 0.16245 , 0.16245 , 0.16245 , 0.16245 , 0.16245 , 0.16245 , 0.16245 , 0.16245 , 0.16245 , 0.16245 , 0.16245 , 0.16245 , 0.16245 , 0.16245 , 0.16245 , 0.16245 , 0.16245 , 0.16245 , 0.16245 , 0.16245 , 0.16245 , 0.16245 , 0.16245 , 0.16245 , 0.16245 , 0.16245 , 0.16245 , 0.16245 , 0.16245 , 0.16245 , 0.16245 , 0.16245 , 0.16245 , 0.16245 , 0.16245 , 0.16244 , 0.16244 , 0.16244 , 0.16244 , 0.16244 , 0.16244 , 0.16244 , 0.16244 , 0.16244 , 0.16244 , 0.16244 , 0.16244069, 0.16243968, 0.16243326, 0.16244052, 0.16243464, 0.16249103, 0.16252413])
- H2O(pressure)float646.895e-05 3.532e-05 ... 0.001348
- units :
- v/v
array([6.89517544e-05, 3.53183235e-05, 1.74440159e-05, 8.32376190e-06, 3.78098542e-06, 1.68526199e-06, 6.71857881e-07, 2.79855655e-07, 1.39894551e-07, 7.44319253e-08, 3.30556972e-08, 1.07133617e-08, 2.03457268e-09, 6.07506964e-11, 2.38279435e-11, 7.61693680e-08, 1.20514102e-06, 1.14635634e-15, 2.54058175e-18, 6.91657839e-19, 6.67517014e-16, 2.25411456e-14, 1.50764406e-11, 3.52528357e-08, 2.63985837e-09, 1.19805526e-10, 3.05737245e-13, 3.79556768e-16, 2.94263434e-16, 2.10031649e-14, 3.66369022e-13, 7.24273285e-13, 6.08681819e-13, 7.39950259e-13, 1.86899614e-12, 5.99569883e-12, 6.06099834e-11, 9.50531819e-10, 4.31682251e-08, 1.21717761e-06, 2.24066227e-05, 2.82636295e-04, 8.19160000e-04, 8.18920000e-04, 8.18920000e-04, 8.18920000e-04, 8.18610000e-04, 8.18610000e-04, 8.18610755e-04, 8.18618182e-04, 8.18619839e-04, 8.18623527e-04, 8.18634087e-04, 8.18699009e-04, 8.19009288e-04, 8.19558131e-04, 8.20068896e-04, 8.25211796e-04, 8.60493176e-04, 1.02865643e-03, 1.34769075e-03])
- CH4(pressure)float640.0004653 0.0004653 ... 0.001662
- units :
- v/v
array([0.00046533, 0.00046533, 0.00046533, 0.00046533, 0.00046533, 0.00046533, 0.00046533, 0.00046533, 0.00046533, 0.00046533, 0.00046533, 0.00046533, 0.00046533, 0.00046533, 0.00046533, 0.00046533, 0.00046533, 0.00046533, 0.00046533, 0.00046533, 0.00046533, 0.00046533, 0.00046533, 0.00046533, 0.00046533, 0.00046533, 0.00046533, 0.00046533, 0.00046533, 0.00046533, 0.00046533, 0.00046533, 0.00046533, 0.00046533, 0.00046533, 0.00046533, 0.00046533, 0.00046533, 0.00046533, 0.00046533, 0.00046533, 0.00046533, 0.00046533, 0.00046532, 0.00046532, 0.00046532, 0.00046532, 0.00046532, 0.00046532, 0.00046532, 0.00046532, 0.00046532, 0.00046532, 0.00046532, 0.00046532, 0.00046535, 0.00046584, 0.00047106, 0.000508 , 0.00073008, 0.00166184])
- CO(pressure)float644.5e-38 4.5e-38 ... 5.644e-06
- units :
- v/v
array([4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 6.16123462e-36, 2.89355093e-32, 5.14891747e-29, 3.93462915e-26, 1.41881488e-23, 5.49281572e-22, 3.89728262e-19, 2.96660647e-17, 1.42574163e-15, 4.49780867e-14, 9.66903095e-13, 1.45968794e-11, 1.58491229e-10, 1.27111690e-09, 7.70598919e-09, 3.66272269e-08, 1.51844422e-07, 7.74843859e-07, 5.64404558e-06])
- NH3(pressure)float640.0001372 0.0001372 ... 0.0001991
- units :
- v/v
array([1.37209077e-04, 1.37210000e-04, 1.37210000e-04, 1.37210000e-04, 1.37210000e-04, 1.37210000e-04, 1.37210000e-04, 1.37210000e-04, 1.37210000e-04, 1.37210000e-04, 1.37210000e-04, 1.37210000e-04, 1.37210000e-04, 1.37210000e-04, 1.37210000e-04, 1.37210000e-04, 1.37210000e-04, 8.91995522e-06, 2.07529906e-07, 8.34206327e-08, 4.24450076e-06, 2.64269045e-05, 1.37210000e-04, 1.37210000e-04, 1.37210000e-04, 1.37210000e-04, 6.06633702e-05, 1.13371313e-06, 8.46942691e-07, 9.21497557e-06, 4.05156502e-05, 5.23524469e-05, 4.23105627e-05, 4.26615347e-05, 5.32498759e-05, 8.13908818e-05, 1.30783777e-04, 1.37210000e-04, 1.37210000e-04, 1.37210000e-04, 1.37210000e-04, 1.37210000e-04, 1.37210000e-04, 1.37200000e-04, 1.37200000e-04, 1.37200000e-04, 1.37200000e-04, 1.37200000e-04, 1.37193107e-04, 1.37161304e-04, 1.37074673e-04, 1.36587407e-04, 1.35233654e-04, 1.33305013e-04, 1.37056405e-04, 1.51376093e-04, 1.70389349e-04, 1.83948436e-04, 1.91074761e-04, 1.94441564e-04, 1.99137556e-04])
- N2(pressure)float645.41e-17 1.583e-17 ... 4.505e-05
- units :
- v/v
array([5.41008642e-17, 1.58319221e-17, 4.36117631e-18, 1.13691545e-18, 2.73063327e-19, 6.35464068e-20, 1.22881236e-20, 2.55522047e-21, 7.20605797e-22, 2.25419502e-22, 5.23409326e-23, 7.22626665e-24, 4.13006331e-25, 1.14256734e-27, 2.17201078e-28, 1.14877922e-22, 1.00920708e-20, 4.30957431e-37, 4.50000000e-38, 4.50000000e-38, 1.37037798e-37, 1.94906994e-36, 3.80553594e-29, 1.25192889e-23, 1.50477556e-25, 8.01783664e-28, 6.82139520e-34, 4.50000000e-38, 4.50000000e-38, 8.59424105e-37, 4.93808130e-34, 3.29508183e-33, 9.90884633e-34, 7.09755760e-34, 5.72460165e-33, 1.32798764e-31, 3.01250381e-29, 6.05613276e-27, 3.23936312e-24, 8.27599640e-22, 1.09034958e-19, 7.99628171e-18, 3.54016036e-16, 1.01252645e-14, 1.99379876e-13, 2.83477668e-12, 3.01449665e-11, 1.26035127e-10, 1.61056671e-09, 8.48245050e-09, 3.65998213e-08, 1.30683112e-07, 3.89239618e-07, 9.82678361e-07, 2.32257599e-06, 5.51411223e-06, 1.19370839e-05, 2.10481718e-05, 3.06919379e-05, 3.87394198e-05, 4.50462199e-05])
- PH3(pressure)float644.5e-38 4.5e-38 ... 5.825e-07
- units :
- v/v
array([4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.54450666e-38, 5.44902680e-35, 5.33872349e-31, 1.59776826e-27, 2.05028762e-24, 1.24720887e-21, 6.18130182e-20, 6.84740027e-17, 7.02685645e-15, 4.50661656e-13, 1.89074369e-11, 5.41081168e-10, 1.08168751e-08, 1.36896598e-07, 5.61556101e-07, 5.36612588e-07, 5.34350231e-07, 5.39238424e-07, 5.54461583e-07, 5.82484946e-07])
- H2-(pressure)float644.5e-38 4.5e-38 ... 4.823e-17
- units :
- v/v
array([4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50060000e-38, 4.50060000e-38, 4.50060000e-38, 4.50060000e-38, 4.50060000e-38, 4.50060000e-38, 1.59021952e-37, 1.33124364e-34, 2.03312709e-30, 5.66505566e-27, 2.89946042e-24, 4.67266870e-22, 3.41968165e-20, 1.55617178e-18, 4.82288340e-17])
- H2+(pressure)float644.5e-38 4.5e-38 ... 8.177e-37
- units :
- v/v
array([4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50060000e-38, 4.50060000e-38, 4.50060000e-38, 4.50060000e-38, 4.50060000e-38, 4.50060000e-38, 4.50060000e-38, 4.50059557e-38, 4.50056950e-38, 4.50045220e-38, 4.50053139e-38, 4.50051779e-38, 4.50068765e-38, 4.50050000e-38, 8.17733576e-37])
- H3+(pressure)float644.5e-38 4.5e-38 ... 4.701e-27
- units :
- v/v
array([4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50060000e-38, 4.50060000e-38, 4.50060000e-38, 4.50060000e-38, 4.50060000e-38, 4.50060000e-38, 4.50060000e-38, 4.50059557e-38, 4.50056950e-38, 4.50045220e-38, 6.84028894e-38, 6.15017016e-36, 1.11151987e-32, 1.01590259e-29, 4.70084038e-27])
- H2S(pressure)float642.539e-05 2.54e-05 ... 2.792e-05
- units :
- v/v
array([2.53949077e-05, 2.53953128e-05, 2.53956311e-05, 2.53959439e-05, 2.52305711e-05, 2.48999340e-05, 2.43851396e-05, 2.13530828e-05, 1.19012327e-05, 6.51680168e-06, 3.03332980e-06, 1.06356522e-06, 2.30598081e-07, 9.47076074e-09, 3.94222221e-09, 5.03699229e-06, 2.01989097e-05, 4.64984597e-13, 1.84176135e-15, 5.52056767e-16, 2.57587867e-13, 5.89409212e-12, 1.97757803e-09, 2.03718177e-06, 1.92008144e-07, 1.14910950e-08, 5.16700864e-11, 1.21534669e-13, 9.33637633e-14, 4.18767381e-12, 5.29024435e-11, 9.43095445e-11, 7.79016253e-11, 8.96849832e-11, 1.99297220e-10, 5.49096546e-10, 4.24691079e-09, 4.87559063e-08, 1.45628142e-06, 1.86307568e-05, 2.53960000e-05, 2.53954357e-05, 2.53950000e-05, 2.53950000e-05, 2.53950000e-05, 2.53950000e-05, 2.53950000e-05, 2.53950000e-05, 2.53950000e-05, 2.53950000e-05, 2.53950000e-05, 2.53950000e-05, 2.53950000e-05, 2.53950000e-05, 2.53992471e-05, 2.53879218e-05, 2.60209405e-05, 2.72608958e-05, 2.72627367e-05, 2.72871294e-05, 2.79224083e-05])
- Fe(pressure)float644.5e-38 4.5e-38 ... 4.185e-09
- units :
- v/v
array([4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50060000e-38, 4.50060000e-38, 6.59204973e-37, 7.38200223e-33, 3.15897904e-29, 5.63642922e-26, 4.55093598e-23, 1.77546083e-20, 3.57501973e-18, 3.90398774e-16, 2.42427550e-14, 9.11291548e-13, 2.20925658e-11, 3.62087281e-10, 4.18494347e-09])
- FeH(pressure)float644.5e-38 4.5e-38 ... 1.401e-10
- units :
- v/v
array([4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50060000e-38, 4.50060000e-38, 6.37472707e-38, 1.70373232e-36, 2.99487818e-33, 1.17395404e-29, 1.96291142e-26, 1.50545084e-23, 5.68659113e-21, 1.11771467e-18, 1.20431284e-16, 7.60528583e-15, 3.00547873e-13, 7.81944881e-12, 1.40052400e-10])
- CrH(pressure)float644.5e-38 4.5e-38 ... 4.015e-10
- units :
- v/v
array([4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50060000e-38, 4.50060000e-38, 7.67446212e-38, 1.27431909e-35, 2.98632808e-32, 9.59808900e-29, 1.33991206e-25, 8.74022817e-23, 2.86286579e-20, 3.32621139e-18, 1.85576809e-16, 1.70448942e-14, 8.13964651e-13, 2.25609297e-11, 4.01476609e-10])
- Na(pressure)float644.5e-38 4.5e-38 ... 3.482e-06
- units :
- v/v
array([4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 1.45395686e-27, 2.58734072e-25, 2.73520219e-21, 1.20211951e-18, 2.74350333e-16, 3.45617548e-14, 2.54341177e-12, 1.14536284e-10, 3.31545031e-09, 6.40870685e-08, 8.53581121e-07, 3.36111870e-06, 3.36909743e-06, 3.40823731e-06, 3.48168606e-06])
- K(pressure)float644.5e-38 4.5e-38 ... 1.964e-07
- units :
- v/v
array([4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 5.22875233e-38, 2.85658091e-36, 5.13708416e-30, 7.14197530e-26, 4.25280317e-22, 1.27961649e-18, 2.09328034e-15, 1.77943987e-12, 4.11333527e-10, 1.42298758e-08, 7.68556506e-08, 1.41250192e-07, 1.55880673e-07, 1.72940566e-07, 1.96405172e-07])
- Rb(pressure)float644.5e-38 4.5e-38 ... 3.722e-10
- units :
- v/v
array([4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 1.43448026e-37, 2.53367597e-35, 3.33939062e-29, 3.72042080e-25, 1.81412795e-21, 3.30576523e-18, 9.66111022e-16, 4.92396332e-14, 1.49125307e-12, 2.34913070e-11, 1.36645555e-10, 2.63226230e-10, 2.89190435e-10, 3.22845627e-10, 3.72184987e-10])
- Cs(pressure)float644.5e-38 4.5e-38 ... 1.376e-11
- units :
- v/v
array([4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 9.94223405e-38, 1.14742559e-35, 1.65437213e-29, 1.95421291e-25, 2.61217865e-22, 6.36324351e-20, 4.91123329e-18, 2.86350078e-16, 1.06750994e-14, 2.11455830e-13, 1.82584257e-12, 5.50007764e-12, 7.12169710e-12, 9.42936159e-12, 1.37572325e-11])
- HCN(pressure)float644.5e-38 4.5e-38 ... 1.678e-08
- units :
- v/v
array([4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 5.82069558e-37, 2.79728957e-33, 5.61794653e-30, 5.61138807e-28, 2.06166226e-24, 4.43445520e-22, 5.24232472e-20, 3.59180461e-18, 1.49759408e-16, 3.97673270e-15, 7.34110056e-14, 1.00636340e-12, 1.00643883e-11, 7.19782524e-11, 4.02120450e-10, 2.29465873e-09, 1.67839046e-08])
- C2H2(pressure)float644.5e-38 4.5e-38 ... 1.848e-10
- units :
- v/v
array([4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50060000e-38, 4.50060000e-38, 9.76396877e-37, 4.37172557e-33, 8.38172009e-30, 7.40307914e-27, 3.23189424e-24, 7.37989696e-22, 9.24054670e-20, 6.66764471e-18, 2.89982942e-16, 8.06940777e-15, 1.67050318e-13, 4.17406379e-12, 1.84809932e-10])
- C2H4(pressure)float644.5e-38 4.5e-38 ... 9.272e-09
- units :
- v/v
array([4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 5.24480043e-37, 1.11151929e-33, 1.02954077e-30, 4.55525155e-28, 1.97767354e-26, 1.77724704e-23, 1.62730053e-21, 9.49180453e-20, 3.66705879e-18, 9.74346522e-17, 1.83475116e-15, 2.50860101e-14, 2.55678041e-13, 1.98948838e-12, 1.22955147e-11, 6.87580591e-11, 5.54032666e-10, 9.27188795e-09])
- C2H6(pressure)float641.247e-30 8.227e-31 ... 2.057e-08
- units :
- v/v
array([1.24675652e-30, 8.22708500e-31, 5.22181290e-31, 3.14739553e-31, 1.78322437e-31, 9.86772303e-32, 4.73817003e-32, 2.40705369e-32, 1.54334893e-32, 1.07163066e-32, 5.90461943e-33, 2.17588541e-33, 4.04960623e-34, 7.08619802e-36, 3.33226052e-36, 1.55872667e-31, 8.25836065e-30, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 7.31028414e-38, 6.34137978e-35, 1.98909245e-30, 1.12717688e-31, 3.36823623e-33, 2.54404273e-36, 4.50000000e-38, 4.50000000e-38, 6.10303421e-37, 1.87383435e-35, 6.95918407e-35, 8.66367921e-35, 1.72862663e-34, 8.78152139e-34, 6.05566105e-33, 1.81001024e-31, 9.48101133e-30, 1.93219814e-27, 2.14164334e-25, 1.37657877e-23, 5.45970492e-22, 1.43053955e-20, 2.61421553e-19, 3.50930348e-18, 3.60142496e-17, 2.89981860e-16, 1.06683054e-15, 1.14301757e-14, 5.63498106e-14, 2.41029749e-13, 9.05437590e-13, 3.02320619e-12, 9.05597895e-12, 2.44822920e-11, 6.02268468e-11, 1.35869726e-10, 2.87885562e-10, 6.42584941e-10, 2.37187351e-09, 2.05656528e-08])
- SiO(pressure)float644.5e-38 4.5e-38 ... 4.903e-08
- units :
- v/v
array([4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50060000e-38, 4.50060000e-38, 4.50060000e-38, 4.50060000e-38, 2.66939102e-35, 2.31154554e-31, 2.54756343e-27, 1.02125206e-23, 1.61398745e-20, 1.07979895e-17, 3.44229865e-15, 5.56032506e-13, 4.68715206e-11, 2.00888164e-09, 4.90300610e-08])
- MgH(pressure)float644.5e-38 4.5e-38 ... 2.843e-11
- units :
- v/v
array([4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50060000e-38, 4.50060000e-38, 4.50060000e-38, 4.50060000e-38, 4.50060000e-38, 1.82630806e-36, 5.43605288e-33, 4.95865552e-29, 1.73828847e-25, 2.56197350e-22, 1.68830641e-19, 5.36118714e-17, 8.62689826e-15, 6.73514905e-13, 2.84330739e-11])
- OCS(pressure)float644.5e-38 4.5e-38 ... 7.931e-12
- units :
- v/v
array([4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.32045632e-37, 5.04622115e-34, 2.64367947e-31, 7.05895446e-29, 2.27485180e-27, 1.18034881e-24, 7.41934105e-23, 3.03890514e-21, 8.40410565e-20, 1.62490605e-18, 2.25840625e-17, 2.29933116e-16, 1.75697091e-15, 1.05498246e-14, 5.13306247e-14, 2.09824673e-13, 1.06470819e-12, 7.93102334e-12])
- graphite(pressure)float644.5e-38 4.5e-38 ... 4.5e-38 4.5e-38
- units :
- v/v
array([4.5e-38, 4.5e-38, 4.5e-38, 4.5e-38, 4.5e-38, 4.5e-38, 4.5e-38, 4.5e-38, 4.5e-38, 4.5e-38, 4.5e-38, 4.5e-38, 4.5e-38, 4.5e-38, 4.5e-38, 4.5e-38, 4.5e-38, 4.5e-38, 4.5e-38, 4.5e-38, 4.5e-38, 4.5e-38, 4.5e-38, 4.5e-38, 4.5e-38, 4.5e-38, 4.5e-38, 4.5e-38, 4.5e-38, 4.5e-38, 4.5e-38, 4.5e-38, 4.5e-38, 4.5e-38, 4.5e-38, 4.5e-38, 4.5e-38, 4.5e-38, 4.5e-38, 4.5e-38, 4.5e-38, 4.5e-38, 4.5e-38, 4.5e-38, 4.5e-38, 4.5e-38, 4.5e-38, 4.5e-38, 4.5e-38, 4.5e-38, 4.5e-38, 4.5e-38, 4.5e-38, 4.5e-38, 4.5e-38, 4.5e-38, 4.5e-38, 4.5e-38, 4.5e-38, 4.5e-38, 4.5e-38])
- opd(pressure_layer, wavenumber_layer)float640.0 0.0 0.0 0.0 ... 0.0 0.0 0.0 0.0
- units :
- unitless
array([[0. , 0. , 0. , ..., 0. , 0. , 0. ], [0. , 0. , 0. , ..., 0. , 0. , 0. ], [0. , 0. , 0. , ..., 0. , 0. , 0. ], ..., [0.59295, 0.67848, 0.67467, ..., 0.47412, 0.47386, 0.47365], [0.29647, 0.33924, 0.33733, ..., 0.23706, 0.23693, 0.23683], [0. , 0. , 0. , ..., 0. , 0. , 0. ]])
- ssa(pressure_layer, wavenumber_layer)float640.0 0.0 0.0 0.0 ... 0.0 0.0 0.0 0.0
- units :
- unitless
array([[0. , 0. , 0. , ..., 0. , 0. , 0. ], [0. , 0. , 0. , ..., 0. , 0. , 0. ], [0. , 0. , 0. , ..., 0. , 0. , 0. ], ..., [0.948187 , 0.9617072, 0.818675 , ..., 0.5906784, 0.5935144, 0.5974549], [0.948187 , 0.9617072, 0.818675 , ..., 0.5906784, 0.5935144, 0.5974549], [0. , 0. , 0. , ..., 0. , 0. , 0. ]])
- asy(pressure_layer, wavenumber_layer)float640.0 0.0 0.0 0.0 ... 0.0 0.0 0.0 0.0
- units :
- unitless
array([[0. , 0. , 0. , ..., 0. , 0. , 0. ], [0. , 0. , 0. , ..., 0. , 0. , 0. ], [0. , 0. , 0. , ..., 0. , 0. , 0. ], ..., [0.3800323, 0.3667444, 0.4358756, ..., 0.8781489, 0.8735224, 0.8670378], [0.3800323, 0.3667444, 0.4358756, ..., 0.8781489, 0.8735224, 0.8670378], [0. , 0. , 0. , ..., 0. , 0. , 0. ]])
- pressurePandasIndex
PandasIndex(Index([ 1e-06, 1.41e-06, 2e-06, 2.82e-06, 3.98e-06, 5.62e-06, 7.94e-06, 1.122e-05, 1.585e-05, 2.239e-05, 3.162e-05, 4.467e-05, 6.31e-05, 8.912e-05, 0.00012589, 0.00017783, 0.00025119, 0.00035481, 0.00050119, 0.00070795, 0.001, 0.0014125, 0.0019953, 0.0028184, 0.0039811, 0.0056234, 0.0079433, 0.01122, 0.015849, 0.022387, 0.031623, 0.044668, 0.063096, 0.089125, 0.12589, 0.17783, 0.25119, 0.35481, 0.5011899999999999, 0.70795, 1.0, 1.4125, 1.9953, 2.8184, 3.9811, 5.6234, 7.9433, 10.0, 15.849, 22.387, 31.623, 44.668, 63.096, 89.125, 125.89, 177.83, 251.19, 354.81, 501.19, 707.95, 1000.0], dtype='float64', name='pressure'))
- wavelengthPandasIndex
PandasIndex(Index([ 0.999907199934528, 0.9998072142139076, 0.9997072384913596, 0.9996072727658833, 0.9995073170364811, 0.9994073713021508, 0.9993074355618947, 0.9992075098147134, 0.999107594059607, 0.9990076882955765, ... 0.30029414410706, 0.30026411619407106, 0.3002340912837236, 0.3002040693757168, 0.30017405046975054, 0.3001440345655245, 0.3001140216627392, 0.3000840117610942, 0.3000540048602893, 0.3000240009600245], dtype='float64', name='wavelength', length=12039))
- author :
- Awesome Scientist
- code :
- picaso
- contact :
- awesomescientist@universe.com
- planet_params :
- {"mp": {"value": 1.8981245973360504e+30, "unit": "g"}, "rp": {"value": 7149200000.0, "unit": "cm"}, "mh": 0, "co": 1, "p_reference": {"value": 1.0, "unit": "bar"}}
- stellar_params :
- {"database": "ck04models", "steff": 5000, "feh": 0, "logg": 4, "rs": {"value": 69570000000.0, "unit": "cm"}}
- orbit_params :
- {"sma": {"value": 14959787070000.0, "unit": "cm"}}
- cloud_params :
- {'fsed': 3}
[ ]: