calculate_radiative_loss_rate(goests, force_download=False, download_dir=None)¶
Calculates radiative loss rate from GOES observations.
This function calculates the radiative loss rate as a function of time of solar soft X-ray-emitting plasma across all wavelengths given a LightCurve object containing GOES data. The radiative loss rate is determined from the GOES isothermal temperature and volume emission measure as a function of time, as calculated by
calculate_temperature_em(). See docstring of that function for more details. If the LightCurve object does not contain the temperatures and emission measures, but only contain the GOES fluxes, then the temperature and emission measures are calculated using calculate_temperature_em(). The unit of the resulting radiative loss rates is W. Once the radiative loss rates have been found, they are returned as part of a new LightCurve object also containing the metadata, GOES fluxes and corresponding temperatures and emission measures of the input LightCurve object.
XRSTimeSeries) – TimeSeries object containing GOES data. The units of these data MUST be W/m^2 (flux), MK (temperature) and cm^-3 (emission measure). If LightCurve object does not contain temperature and emission measure values, they are calculated from the flux values using calculate_temperature_em().
bool) – If True, the GOES radiative loss data file is downloaded even if already locally stored. It is important to do this if a new version of the file has been generated due to a new CHIANTI version being released or the launch of new GOES satellites. Default=False
str) – The directory to download the GOES radiative loss data file to. Default=SunPy default download directory
XRSTimeSeries) – Contains same metadata and data as input LightCurve with the following additional data columns:ts_new.to_dataframe().temperature - Array of temperature values [MK]ts_new.to_dataframe().em - Array of volume emission measure values [cm**-3]ts_new.to_dataframe().rad_loss_rate - radiative loss rate of the coronal soft X-ray-emitting plasma across all wavelengths [W]
The GOES radiative loss rates are calculated using a csv file containing a table of radiative loss rate per unit emission measure at various temperatures. The appropriate values are then found via interpolation. This table was generated using CHIANTI atomic physics database employing the methods of 1. Coronal abundances, a default density of 10**10 cm**-3, and ionization equilibrium of 2 were used.
Cox, D.P., Tucker, W.H. 1969, ApJ, 157, 1157, DOI: 10.1086/150144
Mazzotta, P., Mazzitelli, G., Colafrancesco, S., & Vittorio, N. 1998, A&AS, 133, 339, DOI: 10.1051/aas:1998330
>>> import sunpy.timeseries as ts >>> from sunpy.instr.goes import calculate_radiative_loss_rate >>> from sunpy.data.sample import GOES_XRS_TIMESERIES >>> goests = ts.TimeSeries(GOES_XRS_TIMESERIES) >>> goests.to_dataframe()[0:10] xrsa xrsb 2011-06-06 23:59:59.961999893 1.000000e-09 1.887100e-07 2011-06-07 00:00:02.008999944 1.000000e-09 1.834600e-07 2011-06-07 00:00:04.058999896 1.000000e-09 1.860900e-07 2011-06-07 00:00:06.104999900 1.000000e-09 1.808400e-07 2011-06-07 00:00:08.151999950 1.000000e-09 1.860900e-07 2011-06-07 00:00:10.201999903 1.000000e-09 1.808400e-07 2011-06-07 00:00:12.248999953 1.000000e-09 1.860900e-07 2011-06-07 00:00:14.298999906 1.000000e-09 1.834600e-07 2011-06-07 00:00:16.344999909 1.000000e-09 1.808400e-07 2011-06-07 00:00:18.391999960 1.000000e-09 1.834600e-07 >>> goests_new = calculate_radiative_loss_rate(goests) >>> goests_new.to_dataframe()[0:10] xrsa xrsb temperature em rad_loss_rate 2011-06-06 23:59:59.961999893 1.000000e-09 1.887100e-07 3.503510 2.190626e+48 1.781001e+19 2011-06-07 00:00:02.008999944 1.000000e-09 1.834600e-07 3.534262 2.055847e+48 1.660031e+19 2011-06-07 00:00:04.058999896 1.000000e-09 1.860900e-07 3.518700 2.122771e+48 1.719931e+19 2011-06-07 00:00:06.104999900 1.000000e-09 1.808400e-07 3.550100 1.990333e+48 1.601718e+19 2011-06-07 00:00:08.151999950 1.000000e-09 1.860900e-07 3.518700 2.122771e+48 1.719931e+19 2011-06-07 00:00:10.201999903 1.000000e-09 1.808400e-07 3.550100 1.990333e+48 1.601718e+19 2011-06-07 00:00:12.248999953 1.000000e-09 1.860900e-07 3.518700 2.122771e+48 1.719931e+19 2011-06-07 00:00:14.298999906 1.000000e-09 1.834600e-07 3.534262 2.055847e+48 1.660031e+19 2011-06-07 00:00:16.344999909 1.000000e-09 1.808400e-07 3.550100 1.990333e+48 1.601718e+19 2011-06-07 00:00:18.391999960 1.000000e-09 1.834600e-07 3.534262 2.055847e+48 1.660031e+19