Search NASA⌕ Search

SEARCH · Search NASA

Results for “Tidal Modeling”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
NASA NTRS2001

Orbital, Rotational and Climatic Interactions: Energy Dissipation and Angular Momentum Exchange in the Earth-Moon System

A numerical ocean tide model has been developed and tested using highly accurate TOPEX/Poseidon (T/P) tidal solutions. The hydrodynamic model is based on time stepping a finite difference approximation to the non-linear shallow water equations. Two novel features of our implementation are a rigorous treatment of self attraction and loading (SAL), and a physically based parameterization for internal tide (IT) radiation drag. The model was run for a range of grid resolutions, and with variations in model parameters and bathymetry. For a rational treatment of SAL and IT drag, the model run at high resolution (1/12 degree) fits the T/P solutions to within 5 cm RMS in the open ocean. Both the rigorous SAL treatment and the IT drag parameterization are required to obtain solutions of this quality. The sensitivity of the solution to perturbations in bathymetry suggest that the fit to T/P is probably now limited by errors in this critical input. Since the model is not constrained by any data, we can test the effect of dropping sea-level to match estimated bathymetry from the last glacial maximum (LGM). Our results suggest that the 100 m drop in sea-level in the LGM would have significantly increased tidal amplitudes in the North Atlantic, and increased overall tidal dissipation by about 40%. However, details in tidal solutions for the past 20 ka are sensitive to the assumed stratification. IT drag accounts for a significant fraction of dissipation, especially in the LGM when large areas of present day shallow sea were exposed, and this parameter is poorly constrained at present.

Egbert, Gary D.↗
NASA NTRS1986

Simulation of the thermospheric tides by use of the NCAR thermospheric general circulation model

Numerical calculations of the thermospheric tidal winds and temperatures at equinox are presented. The calculations were made using the National Center for Atmospheric Research Thermospheric General Circulation Model (TGCM), which includes the effects of viscosity, conductivity, diffusion, ion drag, winds, and temperature gradients. The thermospheric diurnal and semidiurnal tides are excited in situ by solar heating and by ion-neutral momentum coupling. The semidiurnal tidal calculations also include the effects of upward propagating waves generated by heating in the lower atmosphere. This semidiurnal propagating component is modeled by use of the classical tidal perturbations as lower boundary conditions. The model is tuned by adjusting the propagating tidal forcing term until calculated semidiurnal wind and temperature fields best approximate incoherent scatter observations. The tidal TGCM results are consistent with previous theoretical work and successfully reproduce high altitude temperature and meridional velocity data, but they give significantly lower magnitudes for velocities and temperatures near 160 km than are seen by observations.

Fensen, C. G.↗
NASA NTRS2017

Shallow Lunar Seismic Activity and the Current Stress State of the Moon

A vast, global network of more than 3200 lobate thrust fault scarps has been revealed in high resolution Lunar Reconnaissance Orbiter Camera (LROC) images. The fault scarps very young, less than 50 Ma, based on their small scale and crisp appearance, crosscutting relations with small-diameter impact craters, and rates of infilling of associated small, shallow graben and may be actively forming today. The population of young thrust fault scarps provides a window into the recent stress state of the Moon and offers insight into the origin of global lunar stresses. The distribution of orientations of the fault scarps is non-random, inconsistent with isotropic stresses from late-stage global contraction as the sole source of stress Modeling shows that tidal stresses contribute significantly to the current stress state of the lunar crust. Tidal stresses (orbital recession and diurnal tides) superimposed on stresses from global contraction result in non-isotropic compressional stress and thrust faults consistent with lobate scarp orientations. Stresses due to orbital recession do not change with orbital position, thus it is with the addition of diurnal stresses that peak stresses are reached. At apogee, diurnal and recession stresses are most compressive near the tidal axis, while at perigee they are most compressive 90 degrees away from the tidal axis. Coseismic slip events on currently active thrust faults are expected to be triggered when peak stresses are reached. Analysis of the timing of the 28 the shallow moonquakes recorded by the Apollo seismic network shows that 19 indeed occur when the Moon is closer to apogee, while only 9 shallow events occur when the Moon is closer to perigee. Here we show the results of relocating the shallow moonquake using an algorithm designed for sparse networks to better constrain their epicentral locations in order to compare them with stress models. The model for the current stress state of the Moon is refined by investigating the contribution of polar wander.

siesmology↗
NASA NTRS2019

KELT-22Ab: A Massive, Short-Period Hot Jupiter Transiting a Near-solar Twin

We present the discovery of KELT-22Ab, a hot Jupiter from the KELT-South survey. KELT-22Ab transits the moderately bright (V11.1) Sun-like G2V star TYC 7518-468-1. The planet has an orbital period of P = 1.3866529 0.0000027 days, a radius of = -R 1.285+ R P 0.071 J 0.12 , and a relatively large mass of = - M 3.47+ M P 0.14 J 0.15 . The star has = - + R 1.099 0.046 R 0.079 , = - + M 1.092 0.041 M 0.045 , = - T 5767+ eff 49 50 K, = - log g 4.393+0.0600.039 (cgs), and [m/H]=+ -0.259+0.0830.085; thus other than its slightly super-solar metallicity, it appears to be a near-solar twin. Surprisingly, KELT-22A exhibits kinematics and a Galactic orbit that are somewhat atypical for thin-disk stars. Nevertheless, the star is rotating rapidly for its estimated age, and shows evidence of chromospheric activity. Imaging reveals a slightly fainter companion to KELT-22A that is likely bound, with a projected separation of 6″ (∼1400 au). In addition to the orbital motion caused by the transiting planet, we detect a possible linear trend in the radial velocity of KELT-22A, suggesting the presence of another relatively nearby body that is perhaps non-stellar. KELT-22Ab is highly irradiated (as a consequence of the small semimajor axis of a R = 4.97), and is mildly inflated. At such small separations, tidal forces become significant. The configuration of this system is optimal for measuring the rate of tidal dissipation within the host star. Our models predict that, due to tidal forces, the semimajor axis is decreasing rapidly, and KELT-22Ab is predicted to spiral into the star within the next Gyr.

Knicole Colon↗
NASA NTRS2017

Shallow Lunar Seismic Activity and the Current Stress State of the Moon

A vast, global network of more than 3200 lobate thrust fault scarps has been revealed in high resolution Lunar Reconnaissance Orbiter Camera (LROC) images. The fault scarps are very young, less than 50 Ma, based on their small scale and crisp appearance, crosscutting relations with small-diameter impact craters, and rates of infilling of associated small, shallow graben and may be actively forming today. The population of young thrust fault scarps provides a window into the recent stress state of the Moon and offers insight into the origin of global lunar stresses. The distribution of orientations of the fault scarps is non-random, inconsistent with isotropic stresses from late-stage global contraction as the sole source of stress. Modeling shows that tidal stresses contribute significantly to the current stress state of the lunar crust. Tidal stresses (orbital recession and diurnal tides) superimposed on stresses from global contraction result in non-isotropic compressional stress and may produce thrust faults consistent with lobate scarp orientations. At any particular point on the lunar surface, peak compressive stress will be reached at a certain time in the diurnal cycle. Coseismic slip events on currently active thrust faults are expected to be triggered when peak stresses are reached. Analysis of the timing of the 28 the shallow moonquakes recorded by the Apollo seismic network shows that 19 indeed occur when the Moon is closer to apogee, while only 9 shallow events occur when the Moon is closer to perigee. Here we report efforts to refine the model for the current stress state of the Moon by investigating the contribution of polar wander. Progress on relocating the epicentral locations of the shallow moonquakes using an algorithm designed for sparse networks is also reported.

structure↗
DOE OSTI2025

StreamGen: Connecting Populations of Streams and Shells to Their Host Galaxies

In this work, we study how the abundance and dynamics of populations of disrupting satellite galaxies change systematically as a function of host galaxy properties. We apply a theoretical model of the phase-mixing process to classify intact satellite galaxies and stellar streamlike and shell-like debris in ∼1500 Milky Way–mass systems generated by a semi-analytic galaxy formation code, SatGen. In particular, we test the effect of host galaxy halo mass, disk mass, ratio of disk scale height to length, and stellar feedback model on disrupting satellite populations. We find that the counts of tidal debris are consistent across all host galaxy models, within a given host mass range, and that all models can have streamlike debris on low-energy orbits, consistent with that observed around the Milky Way. However, we find a preference for streamlike debris on lower-energy orbits in models with a thicker (lower-density) host disk or on higher-energy orbits in models with a more massive host disk. Importantly, we observe significant halo-to-halo variance across all models. These results highlight the importance of simulating and observing large samples of Milky Way–mass galaxies and accounting for variations in host properties when using disrupting satellites in studies of near-field cosmology.

dark matter↗
NASA NTRS1973

Evolution of satellite resonances by tidal dissipation.

Analysis of a realistic model shows how satellites' gravitational interaction can halt their differential tidal evolution when resonant commensurabilities of their orbital periods are reached. The success of this study lends support to the hypothesis that orbit-orbit resonances among satellites in the solar system, including the Titan-Hyperion case, did evolve as a result of tidal energy dissipation. Consideration of the time scale for this evolution process, possible now that the capture mechanism has been revealed, can offer more sophisticated constraints on the tidal dissipation function, Q, and on past orbital conditions.

Greenberg, R.↗
DOE OSTI2016

ROTATING STARS AND THE FORMATION OF BIPOLAR PLANETARY NEBULAE. II. TIDAL SPIN-UP

We present new binary stellar evolution models that include the effects of tidal forces, rotation, and magnetic torques with the goal of testing planetary nebulae (PNs) shaping via binary interaction. We explore whether tidal interaction with a companion can spin-up the asymptotic giant brach (AGB) envelope. To do so, we have selected binary systems with main-sequence masses of 2.5 M {sub ⊙} and 0.8 M {sub ⊙} and evolve them allowing initial separations of 5, 6, 7, and 8 au. The binary stellar evolution models have been computed all the way to the PNs formation phase or until Roche lobe overflow (RLOF) is reached, whatever happens first. We show that with initial separations of 7 and 8 au, the binary avoids entering into RLOF, and the AGB star reaches moderate rotational velocities at the surface (∼3.5 and ∼2 km s{sup −1}, respectively) during the inter-pulse phases, but after the thermal pulses it drops to a final rotational velocity of only ∼0.03 km s{sup −1}. For the closest binary separations explored, 5 and 6 au, the AGB star reaches rotational velocities of ∼6 and ∼4 km s{sup −1}, respectively, when the RLOF is initiated. We conclude that the detached binary models that avoid entering the RLOF phase during the AGB will not shape bipolar PNs, since the acquired angular momentum is lost via the wind during the last two thermal pulses. This study rules out tidal spin-up in non-contact binaries as a sufficient condition to form bipolar PNs.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Stress Models of the Annual Hydrospheric, Atmospheric, Thermal, and Tidal Loading Cycles on California Faults: Perturbation of Background Stress and Changes in Seismicity

Stresses in the lithosphere arise from multiple natural loading sources that include both surface and body forces. The largest surface loads include near-surface water storage, snow and ice, atmosphere pressure, ocean loading, and temperature changes. The solid Earth also deforms from celestial body interactions and variations in Earth's rotation. We model the seasonal stress changes in California from 2006 through 2014 for seven different loading sources with annual periods to produce an aggregate stressing history for faults in the study area. Our modeling shows that the annual water loading, atmosphere, temperature, and Earth pole tides are the largest loading sources and should each be evaluated to fully describe seasonal stress changes. In California we find that the hydrological loads are the largest source of seasonal stresses. We explore the seasonal stresses with respect to the background principal stress orientation constrained with regional focal mechanisms and analyze the modulation of seismicity. Our results do not suggest a resolvable seasonal variation for the ambient stress orientation in the shallow crust. When projecting the seasonal stresses into the background stress orientation we find that the timing of microseismicity modestly increases from an ~8 kPa seasonal mean-normal-stress perturbation. The results suggest that faults in California are optimally oriented with the background stress field and respond to subsurface pressure changes, possibly due to processes we have not considered in this study. At any time a population of faults are near failure as evident from earthquakes triggered by these slight seasonal stress perturbations.

58 GEOSCIENCES↗
DOE OSTI2015

A HYDRODYNAMICAL SOLUTION FOR THE ''TWIN-TAILED'' COLLIDING GALAXY CLUSTER ''EL GORDO''

The distinctive cometary X-ray morphology of the recently discovered massive galaxy cluster ''El Gordo'' (ACT-CT J0102–4915; z = 0.87) indicates that an unusually high-speed collision is ongoing between two massive galaxy clusters. A bright X-ray ''bullet'' leads a ''twin-tailed'' wake, with the Sunyaev-Zel'dovich (SZ) centroid at the end of the northern tail. We show how the physical properties of this system can be determined using our FLASH-based, N-body/hydrodynamic model, constrained by detailed X-ray, SZ, and Hubble lensing and dynamical data. The X-ray morphology and the location of the two dark matter components and the SZ peak are accurately described by a simple binary collision viewed about 480 million years after the first core passage. We derive an impact parameter of ≅300 kpc, and a relative initial infall velocity of ≅2250 km s{sup –1} when separated by the sum of the two virial radii assuming an initial total mass of 2.15 × 10{sup 15} M {sub ☉} and a mass ratio of 1.9. Our model demonstrates that tidally stretched gas accounts for the northern X-ray tail along the collision axis between the mass peaks, and that the southern tail lies off axis, comprising compressed and shock heated gas generated as the less massive component plunges through the main cluster. The challenge for ΛCDM will be to find out if this physically extreme event can be plausibly accommodated when combined with the similarly massive, high-infall-velocity case of the Bullet cluster and other such cases being uncovered in new SZ based surveys.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

capella (α Aurigae) revisited: New binary orbit, physical properties, and evolutionary state

Knowledge of the chemical composition and absolute masses of Capella are key to understanding the evolutionary state of this benchmark binary system comprising two giant stars. Previous efforts, including our own 2009 study, have largely failed to reach an acceptable agreement between the observations and current stellar evolution models, preventing us from assessing the status of the primary. Here we report a revision of the physical properties of the components incorporating recently published high-precision radial velocity measurements, and a new detailed chemical analysis providing abundances for more than 20 elements in both stars. We obtain highly precise (∼0.3%) masses of 2.5687 ± 0.0074 M{sub ⊙} and 2.4828 ± 0.0067 M{sub ⊙}, radii of 11.98 ± 0.57 R{sub ⊙} and 8.83 ± 0.33 R{sub ⊙}, effective temperatures of 4970 ± 50 and 5730 ± 60 K, and independently measured luminosities based on the orbital parallax (78.7 ± 4.2 L{sub ⊙} and 72.7 ± 3.6 L{sub ⊙}). We find an excellent match to stellar evolution models at the measured composition of [Fe/H]=−0.04±0.06. Three different sets of models place the primary star firmly at the end of the core helium-burning phase (clump), while the secondary is known to be evolving rapidly across the Hertzprung gap. The measured lithium abundance, the C/N ratio, and the {sup 12}C/{sup 13}C isotopic carbon abundance ratio, which change rapidly in the giant phase, are broadly in agreement with expectations from models. Predictions from tidal theory for the spin rates, spin–orbit alignment, and other properties do not fare as well, requiring a 40-fold increase in the efficiency of the dissipation mechanisms in order to match the observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

The dynamical tides of spinning Newtonian stars

ABSTRACT We carefully develop the framework required to model the dynamical tidal response of a spinning neutron star in an inspiralling binary system, in the context of Newtonian gravity, making sure to include all relevant details and connections to the existing literature. The tidal perturbation is decomposed in terms of the normal oscillation modes, used to derive an expression for the effective Love number which is valid for any rotation rate. In contrast to previous work on the problem, our analysis highlights subtle issues relating to the orthogonality condition required for the mode-sum representation of the dynamical tide and shows how the prograde and retrograde modes combine to provide the overall tidal response. Utilizing a slow-rotation expansion, we show that the dynamical tide (the effective Love number) is corrected at first order in rotation, whereas in the case of the static tide (the static Love number) the rotational corrections do not enter until second order.

Pnigouras, P. (ORCID:0000000318959431)↗
DOE OSTI2019

The hidden giant: discovery of an enormous Galactic dwarf satellite in Gaia DR2

ABSTRACT We report the discovery of a Milky Way satellite in the constellation of Antlia. The Antlia 2 dwarf galaxy is located behind the Galactic disc at a latitude of b ∼ 11° and spans 1.26°, which corresponds to ∼2.9 kpc at its distance of 130 kpc. While similar in spatial extent to the Large Magellanic Cloud, Antlia 2 is orders of magnitude fainter at MV = −9 mag, making it by far the lowest surface brightness system known (at ∼31.9 mag arcsec−2), ∼100 times more diffuse than the so-called ultra diffuse galaxies. The satellite was identified using a combination of astrometry, photometry, and variability data from Gaia Data Release 2, and its nature confirmed with deep archival DECam imaging, which revealed a conspicuous BHB signal. We have also obtained follow-up spectroscopy using AAOmega on the AAT, identifying 159 member stars, and we used them to measure the dwarf’s systemic velocity, 290.9 ± 0.5 km s−1, its velocity dispersion, 5.7 ± 1.1 km s−1, and mean metallicity, [Fe/H] = −1.4. From these properties we conclude that Antlia 2 inhabits one of the least dense dark matter (DM) haloes probed to date. Dynamical modelling and tidal-disruption simulations suggest that a combination of a cored DM profile and strong tidal stripping may explain the observed properties of this satellite. The origin of this core may be consistent with aggressive feedback, or may even require alternatives to cold dark matter (such as ultra-light bosons).

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

From Existing and New Nuclear and Astrophysical Constraints to Stringent Limits on the Equation of State of Neutron-Rich Dense Matter

Through continuous progress in nuclear theory and experiment and an increasing number of neutron-star (NS) observations, a multitude of information about the equation of state (EOS) for matter at extreme densities is available. To constrain the EOS across its entire density range, this information needs to be combined consistently. However, the impact and model dependency of individual observations vary. Given their growing number, assessing the various methods is crucial to compare the respective effects on the EOS and discover potential biases. For this purpose, we present a broad compendium of different constraints and apply them individually to a large set of EOS candidates within a Bayesian framework. Specifically, we explore different ways of how chiral effective field theory and perturbative quantum chromodynamics can be used to place a likelihood on EOS candidates. We also investigate the impact of nuclear experimental constraints, as well as different radio and x-ray observations of NS masses and radii. This is augmented by reanalyses of the existing data from binary neutron star coalescences, in particular of GW170817, with improved models for the tidal waveform and kilonova light curves, which we also utilize to construct a tight upper limit of 2.39 M ⊙ on the TOV mass based on GW170817’s remnant. Our diverse set of constraints is eventually combined to obtain stringent limits on NS properties. We organize the combination in a way to distinguish between constraints where the systematic uncertainties are deemed small and those that rely on less conservative assumptions. For the former, we find the radius of the canonical 1.4 M ⊙ neutron star to be R 1.4 = 12.2 6 − 0.91 + 0.80 km and the TOV mass at M TOV = 2.2 5 − 0.22 + 0.42 M ⊙ (95% credibility). Including all the presented constraints yields R 1.4 = 12.2 0 − 0.48 + 0.50 km and M TOV = 2.3 0 − 0.20 + 0.07 M ⊙ . When comparing these limits to individual data points, we find that the quoted radius of HESS J1731-347 displays noticeable tension with other constraints. Constraining microphysical properties of the EOS proves more challenging. For instance, the symmetry energy slope is restricted to L sym = 48 − 25 + 21 MeV , where this constraint is mainly dominated by our reanalysis of the PREX-II and CREX experiment. Published by the American Physical Society 2025

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2016

A HUBBLE SPACE TELESCOPE STUDY OF THE ENIGMATIC MILKY WAY HALO GLOBULAR CLUSTER CRATER

We analyze the resolved stellar populations of the faint stellar system, Crater, based on deep optical imaging taken with the Advanced Camera for Surveys on the Hubble Space Telescope. Crater’s color–magnitude diagram (CMD) extends ∼4 mag below the oldest main-sequence (MS) turnoff. Structurally, we find that Crater has a half-light radius of ∼20 pc and no evidence for tidal distortions. We model Crater’s CMD as a simple stellar population (SSP) and alternatively by solving for its full star formation history. In both cases, Crater is well described by an SSP with an age of ∼7.5 Gyr, a metallicity of [M/H] ∼ −1.65, a total stellar mass of M{sub ⋆}∼1e4 M{sub ⊙}, and a luminosity of M{sub V}∼−5.3, located at a distance of d ∼ 145 kpc, with modest uncertainties due to differences in the underlying stellar evolution models. We argue that the sparse sampling of stars above the turnoff and subgiant branch are likely to be 1.0–1.4 M{sub ⊙} blue stragglers and their evolved descendants, as opposed to intermediate-age MS stars. We find that Crater is an unusually young cluster given its location in the Galaxy’s outer halo. We discuss scenarios for Crater’s origin, including the possibility of being stripped from the SMC or the accretion from lower-mass dwarfs such as Leo I or Carina. Despite uncertainty over its progenitor system, Crater appears to have been incorporated into the Galaxy more recently than z ∼ 1 (8 Gyr ago), providing an important new constraint on the accretion history of the Galaxy.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Constraints on the Physical Properties of GW190814 through Simulations Based on DECam Follow-up Observations by the Dark Energy Survey

On 2019 August 14, the LIGO and Virgo Collaborations detected gravitational waves from a black hole and a 2.6 solar mass compact object, possibly the first neutron star–black hole merger. In search of an optical counterpart, the Dark Energy Survey (DES) obtained deep imaging of the entire 90% confidence level localization area with Blanco/DECam 0, 1, 2, 3, 6, and 16 nights after the merger. Objects with varying brightness were detected by the DES Pipeline, and we systematically reduced the candidate counterparts through catalog matching, light-curve properties, host-galaxy photometric redshifts, Southern Astrophysical Research spectroscopic follow-up observations, and machine-learning-based photometric classification. All candidates were rejected as counterparts to the merger. To quantify the sensitivity of our search, we applied our selection criteria to full light-curve simulations of supernovae and kilonovae as they would appear in the DECam observations. Because the source class of the merger was uncertain, we utilized an agnostic, three-component kilonova model based on tidally disrupted neutron star (NS) ejecta properties to quantify our detection efficiency of a counterpart if the merger included an NS. We find that, if a kilonova occurred during this merger, configurations where the ejected matter is greater than 0.07 solar masses, has lanthanide abundance less than 10 -8.56 , and has a velocity between 0.18c and 0.21c are disfavored at the 2σ level. Furthermore, we estimate that our background reduction methods are capable of associating gravitational wave signals with a detected electromagnetic counterpart at the 4σ level in 95% of future follow-up observations.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
NASA NTRS1991

Libration in the earth's rotation

External luni-solar torque exerted on the difference (B-A) of the earth's two equatorial principal moments of inertia gives rise to two types of librational motions in the earth's rotation: the semidiurnal libration in spin and the prograde diurnal libration in polar motion. Formulas for the librations considering a realistic earth model and their tidal decompositions are derived and evaluated. The spin libration has a maximum peal-to-peak amplitude of 0.90 milliarcseconds, that of the polar libration is 0.06 milliarcseconds. Implications concerning their detectability and role in the tidal variation of earth rotation are discussed.

Chao, B. F.↗
NASA NTRS1993

Mars Global Reference Atmospheric Model (Mars-GRAM): Release No. 2 - Overview and applications

The Mars Global Reference Atmospheric Model (Mars-GRAM), a science and engineering model for empirically parameterizing the temperature, pressure, density, and wind structure of the Martian atmosphere, is described with particular attention to the model's newest version, Mars-GRAM, Release No. 2 and to the improvements incorporated into the Release No. 2 model as compared with the Release No. 1 version. These improvements include (1) an addition of a new capability to simulate local-scale Martian dust storms and the growth and decay of these storms; (2) an addition of the Zurek and Haberle (1988) wave perturbation model, for simulating tidal perturbation effects; and (3) a new modular version of Mars-GRAM, for incorporation as a subroutine into other codes.

James, B.↗