https://github.com/geodynamics/citcoms
Revision bcf06ab870d4cfd4a7c8594146ed51e41b23d5f9 authored by Eh Tan on 09 August 2007, 22:57 UTC, committed by Eh Tan on 09 August 2007, 22:57 UTC
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.

1 parent 91bcb85
Raw File
Tip revision: bcf06ab870d4cfd4a7c8594146ed51e41b23d5f9 authored by Eh Tan on 09 August 2007, 22:57 UTC
Finished the compressible Stokes solver for TALA.
Tip revision: bcf06ab
NEWS
CitComS.py NEWS -- history of user-visible changes.  2006-12-04

Please send CitComS.py bug reports to cig-mc@geodynamics.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:
back to top