Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “variational analysis”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,135 records · Page 63

Streamlining Energy Sprawl: Assessment of Geothermal Impacts on Public Lands

Renewable energy power generation, including geothermal, solar photovoltaics (PV), and wind, can have significant land use impacts depending on the technology and the size of the facility. At the same time, innovations across these renewable technologies are working to improve plant efficiency per acre of disturbance. Through a survey of satellite imagery and other sources, this analysis reevaluates the impact of geothermal power generation and examines the variations in land use for different geothermal technologies (binary, flash, and dry steam). Geothermal operations on Bureau of Land Management (BLM) sites serve as the focus for this analysis - i.e., only active sites with federal mineral ownership for which the BLM receives royalties from power production. The results show significant variation in direct land disturbance across the operations sites and across these three geothermal technologies. Although there is variation across plants within the same geothermal technology category, dry steam plants (The Geysers) were shown to have the lowest impact, in acres per megawatt capacity and acres per gigawatt-hour generation, followed by flash steam and binary plants. Overall, direct land impacts relative to capacity and generation have declined 16-27% for facilities constructed since 2000 as compared to those constructed prior, even though most new facilities are binary plants (which traditionally had higher land disturbance). Additionally, capacity factors based on actual net production data were found to average only 50% across all operations sites, suggesting that previous estimates based on nameplate capacity only may lead to an undercounting of generation-based land disturbance metrics. Should the geothermal industry capture additional market share, improving land disturbance metrics through accurate net production and capacity data and further minimizing impacts will be increasingly important as competition over land use intensifies.

GEOTHERMAL ENERGY↗

Hybrid Mobile Communication Networks for Planetary Exploration

A paper discusses the continuing work of the Mobile Exploration System Project, which has been performing studies toward the design of hybrid communication networks for future exploratory missions to remote planets. A typical network could include stationary radio transceivers on a remote planet, mobile radio transceivers carried by humans and robots on the planet, terrestrial units connected via the Internet to an interplanetary communication system, and radio relay transceivers aboard spacecraft in orbit about the planet. Prior studies have included tests on prototypes of these networks deployed in Arctic and desert regions chosen to approximate environmental conditions on Mars. Starting from the findings of the prior studies, the paper discusses methods of analysis, design, and testing of the hybrid communication networks. It identifies key radio-frequency (RF) and network engineering issues. Notable among these issues is the study of wireless LAN throughput loss due to repeater use, RF signal strength, and network latency variations. Another major issue is that of using RF-link analysis to ensure adequate link margin in the face of statistical variations in signal strengths.

Alena, Richard↗

Improved Measurements of Galaxy Star Formation Stochasticity from the Intrinsic Scatter of Burst Indicators

Abstract Measurements of short-timescale star formation variations (i.e., “burstiness” or star formation stochasticity) are integral to our understanding of star formation feedback mechanisms and the assembly of stellar populations in galaxies. We expand upon the work of Broussard et al. by introducing a new analysis of galaxy star formation burstiness that accounts for variations in the Q sg = E B − V stars / E B − V gas distribution, a major confounding factor. We use Balmer decrements from the MOSFIRE Deep Evolution Field (MOSDEF) survey to measure Q sg , which we use to construct mock catalogs from the Santa Cruz Semi-Analytic Models and Mufasa cosmological hydrodynamical simulation based on 3D-HST, Fiber Multi-Object Spectrograph (FMOS)-COSMOS, and MOSDEF galaxies with H α detections. The results of the mock catalogs are compared against observations using the burst indicator η = log 10 ( SFR H α / SFR NUV ) , with the standard deviation of the η distribution indicating burstiness. We find decent agreement between mock and observed η distribution shapes; however, the FMOS-COSMOS and MOSDEF mocks show a systematically low median and scatter in η in comparison to the observations. This work also presents the novel approach of analytically deriving the relationship between the intrinsic scatter in η , scatter added by measurement uncertainties, and observed scatter, resulting in an intrinsic burstiness measurement of 0.06–0.16 dex.

79 ASTRONOMY AND ASTROPHYSICS↗

Monitoring Global Precipitation Using Satellite Observations: Status and Future

The current status of monitoring global precipitation amounts and patterns is described using data sets from the Global Precipitation Climatology Project (GPCP) of the World Climate Research Program (WCRP) and from recent research satellites, especially the Tropical Rainfall Measuring Mission (TRMM). The GPCP monthly (and pentad) data set is a 23-year, globally complete precipitation analysis that is used to explore global and regional variations and trends. The data set is a blend of data mainly from low-orbit microwave satellites and geosynchronous infrared satellites, with additional input from satellite sounder data, Outgoing Longwave Radiation (OLR) data and raingauges. The monthly GPCP data set shows no significant global trend in precipitation over the twenty years, unlike the positive trend in global surface temperatures over the past century. Regional trends are also analyzed. A trend pattern that is a combination of both El Nino and La Nina precipitation features is evident in the 23-year data set. This pattern is related to an increase with time in the number of combined months of El Nino and La Nina during the 23-year period. This apparent trend may be a short-term variation, but also might be related to the increase with time of extreme precipitation events reported elsewhere. Patterns of precipitation variation related to ENSO and other phenomena are shown with clear signals extending from the Tropics into middle and high latitudes of both hemispheres. Also shown, as an example of higher time resolution data is the GPCP daily analysis, which is available for the last six years. A second focus of the talk is on TRMM precipitation data and how these newer data sets incorporating information from the first space-borne meteorological radar compare with the established GPCP data sets.

Adler, Robert F.↗

Range‐separated density functional theory using multiresolution analysis and quantum computing

Quantum computers are expected to outperform classical computers for specific problems in quantum chemistry. Such calculations remain expensive, but costs can be lowered through the partition of the molecular system. In the present study, partition was achieved with range-separated density functional theory (RS-DFT). The use of RS-DFT reduces both the basis set size and the active space size dependence of the ground state energy in comparison with the use of wave function theory (WFT) alone. The utilization of pair natural orbitals (PNOs) in place of canonical molecular orbitals (MOs) results in more compact qubit Hamiltonians. To test this strategy, a basis-set independent framework, known as multiresolution analysis (MRA), was employed to generate PNOs. Tests were conducted with the variational quantum eigensolver for a number of molecules. The results show that the proposed approach reduces the number of qubits needed to reach a target energy accuracy.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A Search for Vector Magnetic Field Variations Associated with the M-Class Flares of 1991 June 10 in AR 6659

A careful analysis of a 6-hour time sequence of vector magnetograms of AR 6659, observed on 1991 June 10 with the MSFC vector magnetograph, has revealed only minor changes in the vector magnetic field azimuths in the vicinity of two M-class flares, and the association of these changes with the flares is not unambiguous. In this paper we present our analysis of the data which includes comparison of vector magnetograms prior to and during the flares, calculation of distributions of the rms variation of the azimuth at each pixel in the field of view of the active region, and examination of the variation with time of the azimuths at every pixel covered by the main flare emissions as observed with the H-alpha telescope coaligned with the vector magnetograph. We also present results of an analysis of evolutionary changes in the azimuth over the field of view of the active region.

Hagyard, Mona J.↗

A Search for Vector Magnetic Field Variations Associated with the M-Class Flares of 1991 June 10 in AR 6659

A careful analysis of a 6-hour time sequence of vector magnetograms of AR 6659, observed on 1991 June 10 with the MSFC vector magnetograph, has revealed only minor changes in the vector magnetic field azimuths in the vicinity of two M-class flares, and the association of these changes with the flares is not unambiguous. In this paper we present our analysis of the data which includes comparison of vector magnetograms prior to and during the flares, calculation of distributions of the rms variation of the azimuth at each pixel in the field of view of the active region, and examination of the variation with time of the azimuths at every pixel covered by the main flare emissions as observed with the H-alpha telescope coaligned with the vector magnetograph. We also present results of an analysis of evolutionary changes in the azimuth over the field of view of the active region.

Hagyard, Mona J.↗

Development of new baseline models for U.S. medium office buildings based on commercial buildings energy consumption survey data

Building energy estimation for the building sector under various scenarios are needed for building energy regulation and policy making. This often starts with representative baselines (either empirical baseline or modeled baseline). Commercial Buildings Energy Consumption Survey (CBECS) data is a widely used empirical baseline for U.S. commercial buildings, but none of the existing baseline model are developed to represent the CBECS data. This paper aims to develop new baseline models for the U.S. medium office buildings, which can produce modeled baselines consistent with the CBECS data. Here, we introduced the methodology to create baseline models and the criteria to evaluate the performance of baseline models. The methodology consists of three phases: (1) identification of model inputs, (2) model calibration, and (3) model validation with uncertainty analysis. The evaluation index is the coefficient of variation of the root-mean-square deviation (CV(RMSD)) of site energy use intensities (EUIs) between the modeled baseline and empirical baseline. Then 30 new baseline models for two vintages (pre- and post-1980) and 15 climate zones were created. The evaluation shows that the CV(RMSD) is lower than 0.05 for the modeled baselines produced by the new baseline models. As a comparison, the CV(RMSD) is higher than 0.1 for the existing modeled baselines generated by DOE Commercial Reference Building Models. Further analysis shows that the new baseline models are able to capture the uncertainties of the representative features of existing buildings.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Dynamic and Regression Modeling of Ocean Variability in the Tide-Gauge Record at Seasonal and Longer Periods

Comparison of monthly mean tide-gauge time series to corresponding model time series based on a static inverted barometer (IB) for pressure-driven fluctuations and a ocean general circulation model (OM) reveals that the combined model successfully reproduces seasonal and interannual changes in relative sea level at many stations. Removal of the OM and IB from the tide-gauge record produces residual time series with a mean global variance reduction of 53%. The OM is mis-scaled for certain regions, and 68% of the residual time series contain a significant seasonal variability after removal of the OM and IB from the tide-gauge data. Including OM admittance parameters and seasonal coefficients in a regression model for each station, with IB also removed, produces residual time series with mean global variance reduction of 71%. Examination of the regional improvement in variance caused by scaling the OM, including seasonal terms, or both, indicates weakness in the model at predicting sea-level variation for constricted ocean regions. The model is particularly effective at reproducing sea-level variation for stations in North America, Europe, and Japan. The RMS residual for many stations in these areas is 25-35 mm. The production of "cleaner" tide-gauge time series, with oceanographic variability removed, is important for future analysis of nonsecular and regionally differing sea-level variations. Understanding the ocean model's strengths and weaknesses will allow for future improvements of the model.

Hill, Emma M.↗

Insight into ideal shear strength of Ni-based dilute alloys using first-principles calculations and correlational analysis

Here the present work examines the effect of alloying elements (denoted X) on the ideal shear strength for 26 dilute Ni-based alloys, Ni11X, as determined by first-principles calculations of pure alias shear deformations. The variations in ideal shear strength are quantitatively explored with correlational analysis techniques, showing the importance of atomic properties such as size and electronegativity. The shear moduli of the alloys are affirmed to show a strong linear relationship with their ideal shear strengths, while the shear moduli of the individual alloying elements were not indicative of alloy shear strength. Through combination with available ideal shear strength data on Mg alloys, a potential application of the Ni alloy data is demonstrated in the search for a set of atomic features suitable for machine learning applications to mechanical properties. As another illustration, the calculated Ni ideal shear strengths play a key role in a predictive multiscale framework for deformation behavior of single crystal alloys at large strains, as shown by simulated stress–strain curves.

36 MATERIALS SCIENCE↗

Application of a Navier-Stokes Solver to the Analysis of Multielement Airfoils and Wings Using Multizonal Grid Techniques

A computational study was performed to determine the predictive capability of a Reynolds averaged Navier-Stokes code (CFL3D) for two-dimensional and three-dimensional multielement high-lift systems. Three configurations were analyzed: a three-element airfoil, a wing with a full span flap and a wing with a partial span flap. In order to accurately model these complex geometries, two different multizonal structured grid techniques were employed. For the airfoil and full span wing configurations, a chimera or overset grid technique was used. The results of the airfoil analysis illustrated that although the absolute values of lift were somewhat in error, the code was able to predict reasonably well the variation with Reynolds number and flap position. The full span flap analysis demonstrated good agreement with experimental surface pressure data over the wing and flap. Multiblock patched grids were used to model the partial span flap wing. A modification to an existing patched- grid algorithm was required to analyze the configuration as modeled. Comparisons with experimental data were very good, indicating the applicability of the patched-grid technique to analyses of these complex geometries.

Jones, Kenneth M.↗

Two Successive Martian Years on the Orbit: Similarities and Differences of CO2 Seasonal Cycle from HEND/ODYSSEY Data

The three years of Mars Odyssey successful work on the martian orbit provide a lot of new information about peculiarities of long term variations of CO2 seasonal cycle. To start such analysis we have used observations of neutron albedo of Mars obtained by High Energy Neutron detector (HEND) mounted onboard Mars Odyssey spacecraft. The high latitude northern and southern regions of Mars are affected by global redistribution of atmospheric CO2 which resulted in 25% of atmospheric mass condensed on martian surface of these regions during winter period of time. The seasonal deposit is formed starting from 60N/60S latitudes and achieve its maximal thickness about 1 m at latitudes close to martian poles. Changes of CO2 deposit thickness is the reason for significant variations of neutron flux above martian poles from summer to winter seasons because CO2 frost effectively hides upper water rich surface layers from the orbit observations in neutrons and gamma-rays. This effect was used to estimate column density of CO2 deposit at different latitudes on North and South of Mars and reconstruct multidimensional model of CO2 deposit showing how snow depth varies as function of latitude, longitude and time. In this presentation we tried to make a next step in our study of martian seasonal CO2 cycle and look for similarities and differences between two successive martian years.

Litvak, M. L.↗

Ectomycorrhizal effects on decomposition are highly dependent on fungal traits, climate, and litter properties: A model-based assessment. Dataset.

To simulate the effects of mycorrhizal fungi on soil organic matter cycling, we incorporated mycorrhizal processes into the Carbon, Organisms, Rhizosphere, and Protection in the Soil Environment (CORPSE) model to develop a new soil model Myco-CORPSE. The new model was calibrated and evaluated against soil measurements taken at temperate forests in New Hampshire (NH) and Georgia (GA). A series of scenario analysis were also conducted to explore the conditions under which ectomycorrhizal (ECM) N acquisition processes can induce different soil C accumulation in ECM systems compared to arbuscular (AM) systems.In this data package, we included:-The Python codes of the standard Myco-CORPSE model we developed: "Standard Myco_CORPSE python codes.zip". The main program is the "gradient_sim.py" which calculates the bulk soil microbes and CN content along a user defined gradient of clay, soil temperature, soil moisture and mycorrhizal dominance, and relies on two subprograms "CORPSE_deriv.py" and "CORPSE_integrate.py". "CORPSE_deriv.py" calculated the changes in all simulated soil stock within every time step and "CORPSE_integrate.py" integrate the changes in all simulated soil stock within simulated time period. The program "Plot.py" is used to plot the major outputs produced by the main program "gradient_sim.py".-The modified Python codes of Myco-CORPSE models with site-level environmental inputs (in NH and GA) used to conduct simulations in NH and GA sites: "NH_GA model simulations.zip". -The Python codes used to evaluate the Myco-CORPSE simulation outputs in NH and GA sites against site-level measurements: "Plot NH_GA simulation against measurements.zip". It includes both the evaluation Python code, the model outputs on NH and GA sites, and the measured soil properties in both sites.-The modified Python codes of Myco-CORPSE models "Scenario analysis_model simulations.zip" that is used to conduct scenario analysis of how different litter properties, mycorrhizal fungal traits, climate, and seasonal variation in temperature and vegetation phenology impact the mycorrhizal effects on soil CN properties. The sub file folder "Scenario analysis_litter traits" contains the codes for scenario analysis of different litter properties; The sub file folder "Scenario analysis_ECM types" contains the codes for scenario analysis of different ECM fungal traits; The sub file folder "Scenario analysis_climate&seasonality" contains the codes for scenario analysis of different climate and seasonalities;-"Scenario analysis_model results and plotting codes.zip" contains all the output files from the the scenario analysis of Myco-CORPSE model as described above and the plotting codes used to the generate the heatmaps and scatterplots shown in the manuscript "Ectomycorrhizal effects on decomposition are highly dependent on fungal traits, climate, and litter properties: A model-based assessment"The majority of the model outputs did not have specific geographic information or temporal coverage because the analysis we conducted are mainly hypothetical model simulations. We only provided geographic description, coordinates and temporal coverage for those soil measurements which we used for model evaluations (included in the "Plot NH_GA simulation against measurements.zip").

54 ENVIRONMENTAL SCIENCES↗

Advances and Limitations of Atmospheric Boundary Layer Observations with GPS Occultation over Southeast Pacific Ocean

The typical atmospheric boundary layer (ABL) over the southeast (SE) Pacific Ocean is featured with a strong temperature inversion and a sharp moisture gradient across the ABL top. The strong moisture and temperature gradients result in a sharp refractivity gradient that can be precisely detected by the Global Positioning System (GPS) radio occultation (RO) measurements. In this paper, the Constellation Observing System for Meteorology, Ionosphere & Climate (COSMIC) GPS RO soundings, radiosondes and the high-resolution ECMWF analysis over the SE Pacific are analyzed. COSMIC RO is able to detect a wide range of ABL height variations (1-2 kilometer) as observed from the radiosondes. However, the ECMWF analysis systematically underestimates the ABL heights. The sharp refractivity gradient at the ABL top frequently exceeds the critical refraction (e.g., −157 N-unit per kilometer) and becomes the so-called ducting condition, which results in a systematic RO refractivity bias (or called N-bias) inside the ABL. Simulation study based on radiosonde profiles reveals the magnitudes of the N-biases are vertical resolution dependent. The N-bias is also the primary cause of the systematically smaller refractivity gradient (rarely exceeding −110 N-unit per kilometer) at the ABL top from RO measurement. However, the N-bias seems not affect the ABL height detection. Instead, the very large RO bending angle and the sharp refractivity gradient due to ducting allow reliable detection of the ABL height from GPS RO. The seasonal mean climatology of ABL heights derived from a nine-month composite of COSMIC RO soundings over the SE Pacific reveals significant differences from the ECMWF analysis. Both show an increase of ABL height from the shallow stratocumulus near the coast to a much higher trade wind inversion further off the coast. However, COSMIC RO shows an overall deeper ABL and reveals different locations of the minimum and maximum ABL heights as compared to the ECMWF analysis. At low latitudes, despite the decreasing number of COSMIC RO soundings and the lower percentage of soundings that penetrate into the lowest 500-m above the mean-sea-level, there are small sampling errors in the mean ABL height climatology. The difference of ABL height climatology between COSMIC RO and ECMWF analysis over SE Pacific is significant and requires further studies.

COSMIC RO↗

HST Multicolor (255-1042 nm) Photometry of Saturn's Main Rings: Radial Profiles, Phase and Opening Angle Variations, and Regional Spectra - 1

The main rings of Saturn were observed with the Planetary Camera of the WFPC2 instrument on the Hubble Space Telescope (HST) from September 1996 to August 2000 as the'ring opening angle to Earth and Sun increased from 4 deg to 24 deg, with a spread of phase angles between 0.3 deg and 6 deg at each opening angle. The rings were routinely observed in the five HST wideband UBVRI filters (F336W, F439W, F555W, F675W, and F814W) and occasionally in the F255W, F785LP, and F1042M filters. The emphasis in this series of papers will be on radial color (implying compositional) variations. In this first paper we describe the analysis technique and calibration procedure, note revisions in a previously published Voyager ring color data analysis, and present new results based on over 100 HST images. In the 300-600 nm spectral range where the rings are red, the 555nm/336nm ratio increases by about 14% as the phase angle increases from 0.3 deg to 6 deg. This effect, never reported previously for the rings, is significantly larger than the phase reddening which characterizes other icy objects, primarily because of the redness of the rings. However, there is no discernible tendency for color to vary with ring opening angle at a given phase angle, and there is no phase variation of color where the spectrum is flat. We infer from this combination of facts that multiple intraparticle scattering, either in a regolith or between facets of an unusually rough surface, is important in these geometries, but that multiple interparticle scattering in a vertically extended layer is not. Voyager color ratios at a phase angle of 14 deg are compatible with this trend, but calibration uncertainties prevent their use in quantitative modeling. Overall ring-average spectra are compatible with those of earlier work within calibration uncertainties, but ring spectra vary noticeably with region. We refine and subdivide the regions previously defined by others. The variation seen between radial profiles of ratios between different wavelength suggests the presence of multiple compositional components with different radial distributions. We present new radial profiles of far UV color ratio (F336W/F255W) showing, substantial global variations having a different radial structure than seen between 555 and 336nm. We also find evidence for absorption in the 850nm spectral range (a feature previously only weakly indicated in ring-averaged spectra) primarily through its radial variation; it is located primarily in the C ring, where the particles are known to have lower albedo, and is consistent with "interplanetary pollution" of the rings.

Cuzzi, Jeffrey N.↗

Dynamic model of the Earth's upper atmosphere

An initial modification to the MSF/J70 Thermospheric Model, in which the variations due to sudden geomagnetic disturbances upon the Earth's upper atmospheric density structure were modeled is presented. This dynamic model of the geomagnetic variation included is an improved version of one which SAO developed from the analysis of the ESRO 4 mass spectrometer data that was incorporated in the Jacchia 1977 model. The variation with geomagnetic local time as well as with geomagnetic latitude are included, and also the effects due to disturbance of the temperature profiles in the region of energy deposition.

Slowey, J. W.↗

Diurnal variation of marine stratocumulus over San Nicolas Island during the FIRE IFO

Preliminary analysis was made of data collected at San Nicolas Island during the Intensive Field Observation (IFO) phase of the First International Satellite Cloud Climatology Program's Regional Experiment (FIRE). Of particular interest was an examination of a distinct diurnal variation in the cloud properties, despite an apparent absence of diurnal forcing from the surface. Direct or indirect radiative modulation of such clouds, as proposed by Fravalo at el. (1981) and Turton and Nicholls (1987) indeed seems likely. Preliminary observational evidence for diurnal change in the marine stratocumulus adjacent to San Nicolas Island is presented. A comparison is then made between the observed behavior and predictions from theoretical models of the interactive effect of radiation on boundary layer clouds.

Davies, R.↗

High charge state carbon and oxygen ions in Earth's equatorial quasi-trapping region

Observations of energetic (1.5 - 300 keV/e) medium-to-high charge state (+3 less than or equal to Q less than or equal to +7) solar wind origin C and O ions made in the quasi-trapping region (QTR) of Earth's magnetosphere are compared to ion trajectories calculated in model equatorial magnetospheric magnetic and electric fields. These comparisons indicate that solar wind ions entering the QTR on the nightside as an energetic component of the plasma sheet exit the region on the dayside, experiencing little or no charge exchange on the way. Measurements made by the CHarge Energy Mass (CHEM) ion spectrometer on board the Active Magnetospheric Particle Tracer Explorer/Charge Composition Explorer (AMPTE/CCE) spacecraft at 7 less than L less than 9 from September 1984 to January 1989 are the source of the new results contained herein: quantitative long-term determination of number densities, average energies, energy spectra, local time distributions, and their variation with geomagnetic disturbance level as indexed by Kp. Solar wind primaries (ions with charge states unchanged) and their secondaries (ions with generally lower charge states produced from primaries in the magnetosphere via charge exchange)are observed throughout the QTR and have distinctly different local time variations that persist over the entire 4-year analysis interval. During Kp larger than or equal to 3 deg intervals, primary ion (e.g., O(+6)) densities exhibit a pronounced predawn maximum with average energy minimum and a broad near-local-noon density minimum with average energy maximum. Secondary ion (e.g., O(+5)) densities do not have an identifiable predawn peak, rather they have a broad dayside maximum peaked in local morning and a nightside minimum. During Kp less than or equal to 2(-) intervals, primary ion density peaks are less intense, broader in local time extent, and centered near midnight, while secondary ion density local time variations diminish. The long-time-interval baseline helps to refine and extend previous observations; for example, we show that ionospheric contribution to O(+3)) is negligible. Through comparison with model ion trajectories, we interpret the lack of pronounced secondary ion density peaks colocated with the primary density peaks to indicate that: (1) negligible charge exchange occurs at L greater than 7, that is, solar wind secondaries are produced at L less than 7, and (2) solar wind secondaries do not form a significant portion of the plasma sheet population injected into the QTR. We conclude that little of the energetic solar wind secondary ion population is recirculated through the magnetosphere.

Christon, S. P.↗