Two non-dimensional parameters are added: "dissipation_number" and "gruneisen" under the Solver component. One can use the original incompressible solver by setting "gruneisen=0". The code will treat this as "gruneisen=infinity". Setting non-zero value to "gruneisen" will switch to compressible solver. One can use the TALA solver for incompressible case by setting "gruneisen" to a non-zero value while setting "dissipation_number=0". This is useful when debugging the compressible solver. Two implementations are available: one by Wei Leng (U. Colorado) and one by Eh Tan (CIG). Leng's version uses the original conjugate gradient method for the Uzawa iteration and moves the contribution of compressibility to the RHS, similar to the method of Ita and King, JGR, 1994. Tan's version uses the bi-conjugate gradient stablized method for the Uzawa iteration, similar to the method of Tan and Gurnis, JGR, 2007. Both versions agree very well. In the benchmark case, 33x33x33 nodes per cap, Di/gamma=1.0, Ra=1.0, delta function of load at the mid mantle, the peak velocity differs by only 0.007%. Leng's version is enabled by default. Edit function solve_Ahat_p_fhat() in lib/Stokes_flow_Incomp.c to switch to Tan's version.
To reference or cite the objects present in the Software Heritage archive, permalinks based on SoftWare Heritage persistent IDentifiers (SWHIDs) must be used instead of copying and pasting the url from the address bar of the browser (as there is no guarantee the current URI scheme will remain the same over time).
Select below a type of object currently browsed in order to display its associated SWHID and permalink.
Finished the compressible Stokes solver for TALA.
CitComS.py NEWS -- history of user-visible changes. 2006-12-04 Please send CitComS.py bug reports to firstname.lastname@example.org. Version 2.2.1 * One example, Cookbook7, was missing in v2.2.0. Version 2.2.0 * A new feature to track particles in the flow (code donated by Allen McNamara and Shijie Zhong). * A new feature to use the tracers to compute composition field and compute thermo-chemical convection (code donated by Allen McNamara and Shijie Zhong). A new cookcook (Cookbook7) in the manual introduces how to use this feature. * New feature to use the tracers to define the low viscosity wedges/channels in the subduction zones (code donated by Vlad Manea). * Fixes of sevearal bugs: issue 72, 73, 75, 77 and 81. See the Geodynamics Roundup issue tracker (http://geodynamics.org/roundup) and the file ChangeLog for details. Important bug fixes include: - 'make install' broken for non-pyred CitcomS. (Issue 72) - 'autocombine.py' did not work in Ubuntu 6.10. (Issue 73) - Velocity b.c not scaled correctly when age becomes negative. When using file_vbcs=on and the age becomes negative (passing the present time), the read-in velocity BC was not scaled properly. Usually the boundary velocity became ~1000x smaller than it should be. (Issue 77) - Errors in the manual. The manual failed to document the unit of the vbc_file is in cm/yr. The manual contained a small error in the description of the header of the coor_file. (Issue 75 and 81) Version 2.1.0 * The filenames of post-processed and combined data are changed slightly. For example, "model.cap0.100" is renamed to "model.cap00.100", "model.cap6.100" is renamed to "model.cap06.100", and "model.cap11.100" is unchanged. * A new input parameter "datadir" specifies the path of output files. The "datafile" input parameter (cannot contain path) specifies the filename prefix of the output files. In v2.0.x and before, "datafile" is used to specify the path and filename prefix of output files. This change is not backward compatible. Users need to change their old input file to run it. * Improved build procedure. The CitComS executable embeds python interpreter by default and is statically linked with MPI library, which improves the portability of the code. * New input format (INI-styled .cfg) * No need to specify launcher.nodes anymore, it is calculated from the mesher properties. * Support of binary format output (using HDF5 library). * Additional data output (including: geoid, stress, pressure, and horizontal average) can be enabled/disabled by changing the input parameters. * Scripts (using GMT commands) to plot 2-D cross sections of temperature field. * Added/removed a few input parameters (all in solver.const component) for computing geoid. Removed "layerd" and "wdensity" and added "radius", "density_above" and "density_below". * Fixes for several bugs: issue 30, 35, 42, 44 and 55. See the Geodynamics Roundup issue tracker (http://geodynamics.org/roundup) and the file ChangeLog for details. Important bug fixes include: - Inconsistent times in C and Pyre. When 'fixed_timestep' is set and the grid resolution is insufficient, the Pyre part of CitComS could give the wrong model time. (Issue 35) - Accidental deletion of surf files by batchsurf.py. batchsurf.py is a script to combine the surf files. If the current working directory happens to be the output directory of CitcomS, the script would delete the original surf files. (Issue 55) Version 2.0.2 * Fixes for several minor bugs: issues 7, 8, 9, 10, 11, 12, and 13. Most of these are installation-related. See the Geodynamics Roundup issue tracker (http://geodynamics.org/roundup) and the file ChangeLog for details. Version 2.0.1 * CitComS.py is now distributed under the GNU General Public License. See the file COPYING for copying conditions. * CitComS.py now uses the GNU Build System (Autoconf, Automake, and Libtool). See the file INSTALL for building and installation instructions. Local variables: version-control: never mode: indented-text end:
Computing file changes ...