Citcom

Last updated
CitCom
Original author(s) Louis Moresi (Monash U.)
Initial releaseearly 1990s
Stable release
2.2 / 27 March 2007;16 years ago (2007-03-27)
Written in C
License GNU General Public License
Website geodynamics.org/cig/software/citcomcu/   OOjs UI icon edit-ltr-progressive.svg

CitCom (for California Institute of Technology Convection in the Mantle) is a finite element code designed to solve thermal convection problems relevant to Earth's mantle released under the GNU General Public License. Written in C, the code's latest version, CitComS, runs on a variety of parallel processing computers, including shared and distributed memory platforms.

Contents

History

CitCom was originally written in the early 1990s by Louis Moresi (Monash U.). Although the code for three-dimensional problems was incorporated from its inception, early versions of the software only solved for time-dependent convection problems within two-dimensional Cartesian domains. Moresi's original code turned out to be incredibly modular and easily extensible. Consequently, the fundamental finite element infrastructure which Louis wrote is still in place and forms the basis for much of the code contained in the present release.

In the mid-1990s Moresi wrote versions of the code that solved the equations within three-dimensional Cartesian domains. Then Shijie Zhong (U. of Colorado, Boulder) successfully parallelized CitCom using message passing routines on a limited release Intel supercomputer. Zhong then created a spherical version of the code which he named CitComS. Lijie Han (Planetary Science Institute) then created a regional version of CitComS as well as an alternate version of message passing for an arbitrarily large number of processors. Clint Conrad (Johns Hopkins) created the first Beowulf implementations of the code, then Conrad and Eh Tan (Computational Infrastructure for Geodynamics) re-coded the message passing of the fully spherical version so that problems run on arbitrarily large numbers of processors could also be solved. A plethora of different versions of CitCom exist both on computers at the California Institute of Technology and around the world.

Consequently, by 2002, there were so many different versions of the code that some rationalization was in order. The software was migrated into a version control system and Eh Tan and Eun-seo Choi (Caltech) created a version of CitComS that generates either a fully spherical or regional model, CitcomSFull and CitcomSRegional respectively. CitComS was released to the community through the former GeoFramework project as version 1.0 and 1.1.

By 2004, in order to increase the functionality of CitComS, the developers began to reengineer the code into an object-oriented environment specifically so it could work with a Python-based modeling framework called Pyre. This release of the software, now named CitComS.py, is essentially the product of those reengineering efforts. Eh Tan was the principal developer of CitComS.py, with considerable help from Eun-seo Choi and Michael Aivazis (Caltech).

CitComS is one component of a larger collection of software encompassed by the former GeoFramework project, a collaboration between the Center for Advanced Computing Research (CACR) [1] and the Seismological Laboratory, [2] both at Caltech, and the Victorian Partnership for Advanced Computing [3] in Australia. The GeoFramework project developed a suite of tools to model multi-scale deformation for Earth science problems. This effort was motivated by the need to understand interactions between the long-term evolution of plate tectonics and shorter term processes such as the evolution of faults during and between earthquakes. During 2005 and 2006 much of the remaining software developed by GeoFramework was released under a GPL license and made available from Computational Infrastructure for Geodynamics (CIG). [4]

The second major release of CitComS (2.0) incorporated the software framework Pyre, free surface modeling methods, and stress boundary conditions on the top and bottom surfaces. In the summer of 2005, as part of the 2.0.1 release, CIG replaced the old build procedure with the GNU Build System. A subsequent release, version 2.0.2, could compile and run on 64-bit systems.

The third major release of CitComS (2.1) incorporated new features and functionality, the most important being the use of HDF5 (a parallel version of the Hierarchical Data Format). The HDF5 format allows you to deal with the massive data output created for production runs. This version accepted .cfg files on input, which are easier to create and read.

Other improvements included the incorporation of geoid calculations that had been left out of earlier releases, as well as new scripts to allow results to be visualized with MayaVi2 [5] in addition to Generic Mapping Tools (GMT) [6] and OpenDX. [7] Instructions were provided on using this version as a preinstalled package on some of the NSF TeraGrid sites.

The latest release of CitComS (2.2, 3/27/07) incorporates the ability of tracing particles in the flow. The tracer code was developed by Allen McNamara and Shijie Zhong in 2004 and donated to CIG in early 2007. The tracer code has a wide range of applications in the mantle convection. It can be used in tracing the trajectory of passive particles, in delineating the top boundary of subducted slabs to define the low viscosity wedges, or in tracking the evolution of the chemical composition field.

Related Research Articles

<span class="mw-page-title-main">Plate tectonics</span> Movement of Earths lithosphere

Plate tectonics is the scientific theory that Earth's lithosphere comprises a number of large tectonic plates which have been slowly moving since about 3.4 billion years ago. The model builds on the concept of continental drift, an idea developed during the first decades of the 20th century. Plate tectonics came to be accepted by geoscientists after seafloor spreading was validated in the mid-to-late 1960s.

<span class="mw-page-title-main">Wolfram Mathematica</span> Computational software program

Wolfram Mathematica is a software system with built-in libraries for several areas of technical computing that allow machine learning, statistics, symbolic computation, data manipulation, network analysis, time series analysis, NLP, optimization, plotting functions and various types of data, implementation of algorithms, creation of user interfaces, and interfacing with programs written in other programming languages. It was conceived by Stephen Wolfram, and is developed by Wolfram Research of Champaign, Illinois. The Wolfram Language is the programming language used in Mathematica. Mathematica 1.0 was released on June 23, 1988 in Champaign, Illinois and Santa Clara, California.

<span class="mw-page-title-main">Geoid</span> Ocean shape without winds and tides

The geoid is the shape that the ocean surface would take under the influence of the gravity of Earth, including gravitational attraction and Earth's rotation, if other influences such as winds and tides were absent. This surface is extended through the continents. According to Gauss, who first described it, it is the "mathematical figure of the Earth", a smooth but irregular surface whose shape results from the uneven distribution of mass within and on the surface of Earth. It can be known only through extensive gravitational measurements and calculations. Despite being an important concept for almost 200 years in the history of geodesy and geophysics, it has been defined to high precision only since advances in satellite geodesy in the late 20th century.

Parallel Virtual Machine (PVM) is a software tool for parallel networking of computers. It is designed to allow a network of heterogeneous Unix and/or Windows machines to be used as a single distributed parallel processor. Thus large computational problems can be solved more cost effectively by using the aggregate power and memory of many computers. The software is very portable; the source code, available free through netlib, has been compiled on everything from laptops to Crays.

Louis-Noël Moresi is a Professor of Computational Mathematics & Geophysics at The Australian National University. He has deeply influenced the understanding of the Geophysics community through his own research as well as providing software for the community to use.

SALOME is a multi-platform open source (LGPL-2.1-or-later) scientific computing environment, allowing the realization of industrial studies of physics simulations.

<span class="mw-page-title-main">COIN-OR</span>

Computational Infrastructure for Operations Research (COIN-OR), is a project that aims to "create for mathematical software what the open literature is for mathematical theory." The open literature provides the operations research (OR) community with a peer-review process and an archive. Papers in operations research journals on mathematical theory often contain supporting numerical results from computational studies. The software implementations, models, and data used to produce the numerical results are typically not published. The status quo impeded researchers needing to reproduce computational results, make fair comparisons, and extend the state of the art.

<span class="mw-page-title-main">Model photosphere</span>

The photosphere denotes those solar or stellar surface layers from which optical radiation escapes. These stellar outer layers can be modeled by different computer programs. Often, calculated models are used, together with other programs, to calculate synthetic spectra for stars. For example, in varying the assumed abundance of a chemical element, and comparing the synthetic spectra to observed ones, the abundance of that element in that particular star can be determined. As computers have evolved, the complexity of the models has deepened, becoming more realistic in including more physical data and excluding more of the simplifying assumptions. This evolution of the models has also made them applicable to different kinds of stars.

The Computational Infrastructure for Geodynamics (CIG) is a community-driven organization that advances Earth science by developing and disseminating software for geophysics and related fields. It is a National Science Foundation-sponsored collaborative effort to improve geodynamic modelling and develop, support, and disseminate open-source software for the geodynamics research and higher education communities.

MADNESS is a high-level software environment for the solution of integral and differential equations in many dimensions using adaptive and fast harmonic analysis methods with guaranteed precision based on multiresolution analysis and separated representations .

<span class="mw-page-title-main">KIVA (software)</span>

KIVA is a family of Fortran-based Computational Fluid Dynamics software developed by Los Alamos National Laboratory (LANL). The software predicts complex fuel and air flows as well as ignition, combustion, and pollutant-formation processes in engines. The KIVA models have been used to understand combustion chemistry processes, such as auto-ignition of fuels, and to optimize diesel engines for high efficiency and low emissions. General Motors has used KIVA in the development of direct-injection, stratified charge gasoline engines as well as the fast burn, homogeneous-charge gasoline engine. Cummins reduced development time and cost by 10%–15% using KIVA to develop its high-efficiency 2007 ISB 6.7-L diesel engine that was able to meet 2010 emission standards in 2007. At the same time, the company realized a more robust design and improved fuel economy while meeting all environmental and customer constraints.

<span class="mw-page-title-main">HeuristicLab</span> Software environment

HeuristicLab is a software environment for heuristic and evolutionary algorithms, developed by members of the Heuristic and Evolutionary Algorithm Laboratory (HEAL) at the University of Applied Sciences Upper Austria, in Hagenberg im Mühlkreis. HeuristicLab has a strong focus on providing a graphical user interface so that users are not required to have comprehensive programming skills to adjust and extend the algorithms for a particular problem. In HeuristicLab algorithms are represented as operator graphs and changing or rearranging operators can be done by drag-and-drop without actually writing code. The software thereby tries to shift algorithm development capability from the software engineer to the user and practitioner. Developers can still extend the functionality on code level and can use HeuristicLab's plug-in mechanism that allows them to integrate custom algorithms, solution representations or optimization problems.

<span class="mw-page-title-main">Earth's internal heat budget</span> Accounting of the energy flows at and below the planets crust

Earth's internal heat budget is fundamental to the thermal history of the Earth. The flow of heat from Earth's interior to the surface is estimated at 47±2 terawatts (TW) and comes from two main sources in roughly equal amounts: the radiogenic heat produced by the radioactive decay of isotopes in the mantle and crust, and the primordial heat left over from the formation of Earth.

<span class="mw-page-title-main">GPlates</span> Open-source application software for interactive plate-tectonic reconstructions

GPlates is open-source application software offering a novel combination of interactive plate-tectonic reconstructions, geographic information system (GIS) functionality and raster data visualisation.

ILNumerics is a mathematical class library for Common Language Infrastructure (CLI) developers and a domain specific language (DSL) for the implementation of numerical algorithms on the .NET platform. While algebra systems with graphical user interfaces focus on prototyping of algorithms, implementation of such algorithms into distribution-ready applications is done using development environments and general purpose programming languages (GPL). ILNumerics is an extension to Visual Studio and aims at supporting the creation of technical applications based on .NET.

The Flatiron Institute is an American internal research division of the Simons Foundation, launched in 2016. It comprises five centers for computational science: the Center for Computational Astrophysics (CCA); the Center for Computational Biology (CCB); the Center for Computational Quantum Physics (CCQ); the Center for Computational Mathematics (CCM); and the Center for Computational Neuroscience (CCN). It also has a group called the Scientific Computing Core (SCC). The institute takes its name from the Flatiron District in New York City where it's based.

<span class="mw-page-title-main">Numerical modeling (geology)</span> Technique to solve geological problems by computational simulation

In geology, numerical modeling is a widely applied technique to tackle complex geological problems by computational simulation of geological scenarios.

This glossary of computer science is a list of definitions of terms and concepts used in computer science, its sub-disciplines, and related fields, including terms relevant to software, data science, and computer programming.

<span class="mw-page-title-main">Qiskit</span>

Qiskit is an open-source software development kit (SDK) for working with quantum computers at the level of circuits, pulses, and algorithms. It provides tools for creating and manipulating quantum programs and running them on prototype quantum devices on IBM Quantum Platform or on simulators on a local computer. It follows the circuit model for universal quantum computation, and can be used for any quantum hardware that follows this model.

References

  1. "Center for Advanced Computing Research (CACR)". Archived from the original on 2008-12-23. Retrieved 2020-07-22.
  2. Seismological Laboratory
  3. Victorian Partnership for Advanced Computing
  4. Computational Infrastructure for Geodynamics (CIG)
  5. MayaVi2
  6. "Generic Mapping Tools (GMT)". Archived from the original on 2013-08-13. Retrieved 2022-03-05.
  7. "OpenDX". Archived from the original on 2006-12-05. Retrieved 2019-07-03.