bc9a15c | Eh Tan | 03 June 2009, 21:08:01 UTC | Update SVTInlet to inherit from BaseSVTInlet | 03 June 2009, 21:08:01 UTC |
56e2c06 | Eh Tan | 03 June 2009, 21:07:04 UTC | adding files to makefile | 03 June 2009, 21:07:04 UTC |
b3dfd01 | Eh Tan | 02 June 2009, 22:58:45 UTC | update to current coupler | 02 June 2009, 22:58:45 UTC |
7fd3100 | Eh Tan | 02 June 2009, 22:58:24 UTC | reducing the accuracy of solvers, otherwise, they will converge to the wrong result | 02 June 2009, 22:58:24 UTC |
748c305 | Eh Tan | 02 June 2009, 22:56:46 UTC | Added parameters 'amending_outflow' and 'exchange_pressure' to help the convergence of esolver When 'amending_outflow' is set to true, the imposed velocity BC will be amended slightly to be divergence-free. (The divergence, e.g. outflow, is caused by the combination of solver inaccuracy and interpolation inaccuracy.) When 'exchange_pressure' is set to true, the initial pressure (at element level) of the embedded solver is taken from the pressure solution of the containing solver. | 02 June 2009, 22:56:46 UTC |
179463a | Eh Tan | 02 June 2009, 22:09:25 UTC | Boundary for global mesh, not sure whether it works with interpolator or not | 02 June 2009, 22:09:25 UTC |
2ef0c3d | Eh Tan | 02 June 2009, 21:58:40 UTC | init'd E->output.seismic and E->output.coord_bin | 02 June 2009, 21:58:40 UTC |
101123f | Eh Tan | 02 June 2009, 21:58:23 UTC | refining cookbook10 | 02 June 2009, 21:58:23 UTC |
e1af994 | Leif Strand | 28 May 2009, 01:50:40 UTC | CUDA-ized e_assemble_del2_u(), and ran it successfully -- and slowly! -- under the device emulator. Added a new 'configure' option: "--with-cuda". | 28 May 2009, 01:50:40 UTC |
500e3f5 | Leif Strand | 27 May 2009, 02:32:03 UTC | Parallelized strip_bcs_from_residual(). | 27 May 2009, 02:32:03 UTC |
511d3f3 | Leif Strand | 27 May 2009, 01:28:55 UTC | Sketched data-parallel version of e_assemble_del2_u(). | 27 May 2009, 01:28:55 UTC |
060b73c | Thorsten Becker | 26 May 2009, 15:21:42 UTC | Added blob_bc_persist parameter (default = off) to allow assigning blob-like temperature anomalies to BCs. | 26 May 2009, 15:21:42 UTC |
5a6accf | Thorsten Becker | 26 May 2009, 15:03:36 UTC | Fixed typo in blob temperature assignment. Changed conversion of coordinates, now assignment works. | 26 May 2009, 15:03:36 UTC |
accc9b4 | Eh Tan | 21 May 2009, 23:16:27 UTC | add changes in r14999 | 21 May 2009, 23:16:27 UTC |
f99c960 | Eh Tan | 15 May 2009, 22:32:48 UTC | ignoring the degree-0 and degree-1 components of geoid | 15 May 2009, 22:32:48 UTC |
1670e55 | Eh Tan | 10 May 2009, 23:18:05 UTC | add more chaanges in v3.1.0 | 10 May 2009, 23:18:05 UTC |
a622619 | Eh Tan | 10 May 2009, 16:37:31 UTC | Change version number for v3.1.0 alpha release | 10 May 2009, 16:37:31 UTC |
658d80a | Eh Tan | 10 May 2009, 16:35:07 UTC | Summary of changes in this release | 10 May 2009, 16:35:07 UTC |
573f5d0 | Eh Tan | 01 May 2009, 22:19:35 UTC | Adding cookbook for synthetic seismogram, mostly copied from cookbook7 | 01 May 2009, 22:19:35 UTC |
7a399ef | Eh Tan | 01 May 2009, 22:18:00 UTC | Depth is measured from r=1.0, not outer radius, because in coupled model, the outer radius can be < 1.0. | 01 May 2009, 22:18:00 UTC |
c9c08bb | Eh Tan | 01 May 2009, 22:17:27 UTC | Printing a big warning when plate velocity is imposed but remove_angular_momentum or remove_rigid_rotation is on. | 01 May 2009, 22:17:27 UTC |
24392fd | Eh Tan | 01 May 2009, 22:15:05 UTC | Removing angular momentum in pseudo free surface case | 01 May 2009, 22:15:05 UTC |
5ab9dda | Eh Tan | 01 May 2009, 22:14:29 UTC | ignore anything after ';' | 01 May 2009, 22:14:29 UTC |
b215c86 | Eh Tan | 01 May 2009, 22:13:26 UTC | explicitly add Python.h to requirement | 01 May 2009, 22:13:26 UTC |
05da830 | Leif Strand | 28 April 2009, 22:53:38 UTC | For my own sanity, simplified code assuming the following: E->sphere.caps_per_proc == 1 E->parallel.nproc == 1 E->mesh.levmax == 0 E->mesh.nsd == 3 | 28 April 2009, 22:53:38 UTC |
2064543 | Leif Strand | 27 April 2009, 20:15:10 UTC | Moved malloc/free calls outside of kernel. (The amount of memory allocated appears to be invariant.) | 27 April 2009, 20:15:10 UTC |
fa871b7 | Leif Strand | 27 April 2009, 18:39:15 UTC | Consolidated calls to malloc/free. | 27 April 2009, 18:39:15 UTC |
2d995dd | Thorsten Becker | 24 April 2009, 18:37:49 UTC | Added flag to allow using the more precise method of computing strain-rates throughout the globe, and not just close to the poles. precise_strain_rate=on will select this, default is "odd", i.e. poles only as before. | 24 April 2009, 18:37:49 UTC |
0eb5c37 | Leif Strand | 23 April 2009, 02:32:01 UTC | Don't call through function pointers -- otherwise, 'nvcc' crashes! | 23 April 2009, 02:32:01 UTC |
185800d | Leif Strand | 23 April 2009, 02:22:19 UTC | Proposed CUDA kernel for conjugate gradient solver. | 23 April 2009, 02:22:19 UTC |
a6bc514 | Leif Strand | 22 April 2009, 00:59:20 UTC | Avoid core dump in record() when E->fp hasn't been set yet. | 22 April 2009, 00:59:20 UTC |
e45eab4 | Leif Strand | 21 April 2009, 22:29:26 UTC | The new Python.h header check reorders this script somewhat. CIT_CHECK_PYTHON_HEADER must follow AC_PROG_CC, and precede CIT_PYTHON_EGG_SETUP. | 21 April 2009, 22:29:26 UTC |
c46afd8 | Leif Strand | 21 April 2009, 19:33:59 UTC | Issue45: Check for Python.h. | 21 April 2009, 19:33:59 UTC |
756b97b | Thorsten Becker | 17 April 2009, 03:37:28 UTC | Added viscosity options 9 and 10 for TDPEPV which are identical to 3 and 8 but temperature is not limited to [0;1]. This is cumbersome, but perhaps needed for backward compatibility. | 17 April 2009, 03:37:28 UTC |
9ea8a44 | Eh Tan | 16 April 2009, 18:24:28 UTC | Move parameter "mantle_temp" from CitcomS.solver.param to CitcomS.solver.ic. Renamed E->control.lith_age_mantle_temp to E->control.mantle_temp. | 16 April 2009, 18:24:28 UTC |
0591133 | Eh Tan | 16 April 2009, 18:23:06 UTC | update cookbook8 with the manual | 16 April 2009, 18:23:06 UTC |
487eeea | Eh Tan | 16 April 2009, 18:22:22 UTC | Add -DUSE_HDF5 to CPPFLAGS in module/Exchanger | 16 April 2009, 18:22:22 UTC |
075b734 | Eh Tan | 16 April 2009, 18:21:50 UTC | change the viscosity section to be the same as that in cookbook3. | 16 April 2009, 18:21:50 UTC |
2c35d25 | Eh Tan | 16 April 2009, 18:20:08 UTC | Removed domain extents, since they are determined by coor_file. | 16 April 2009, 18:20:08 UTC |
a7389af | Eh Tan | 16 April 2009, 18:07:53 UTC | Removed instructions to 'make install', 'make clean' and 'make distclean'. Added dependencies. | 16 April 2009, 18:07:53 UTC |
a2a39ef | Eh Tan | 09 April 2009, 20:53:58 UTC | Fixed a bug in cgrad solver introduced in r13270. | 09 April 2009, 20:53:58 UTC |
7dbc76c | Thorsten Becker | 09 April 2009, 20:39:24 UTC | Left debugging output in, now cleaned up | 09 April 2009, 20:39:24 UTC |
d5b89ca | Thorsten Becker | 09 April 2009, 20:37:59 UTC | Fixed file output issue which caused vtk log file I/O to bomb during Citcom finalize. | 09 April 2009, 20:37:59 UTC |
9d6cfb4 | Eh Tan | 08 April 2009, 23:39:49 UTC | Sync'ing C and Pyre and adding comments | 08 April 2009, 23:39:49 UTC |
971856b | Eh Tan | 08 April 2009, 23:38:51 UTC | Added a function citcom_finalize() to close files and MPI before exit. Return 0 when the code is finished normally. Rearrang post_processing and stokes_flow_only. | 08 April 2009, 23:38:51 UTC |
6e42502 | Eh Tan | 08 April 2009, 23:34:43 UTC | Clean up. | 08 April 2009, 23:34:43 UTC |
811838a | Eh Tan | 08 April 2009, 23:31:11 UTC | Removed multigrid-el option of stokes solver, which is broken for a long time and not used by anyone. | 08 April 2009, 23:31:11 UTC |
f482a84 | Eh Tan | 08 April 2009, 23:27:51 UTC | Moved where to call read_rayleigh_from_file() earlier. I hope this doesn't break the ray_control stuff. | 08 April 2009, 23:27:51 UTC |
5ab6970 | Thorsten Becker | 08 April 2009, 20:24:59 UTC | Fixed bug in vtop assignment. | 08 April 2009, 20:24:59 UTC |
c6a6fd7 | Eh Tan | 07 April 2009, 22:50:02 UTC | Revert the unintended modification of examples/Makefile.am | 07 April 2009, 22:50:02 UTC |
c5e1db6 | Eh Tan | 07 April 2009, 22:45:32 UTC | Change the viscosity option in cookbook3 and cookbook4 | 07 April 2009, 22:45:32 UTC |
ce1c4ce | Eh Tan | 07 April 2009, 22:44:24 UTC | Change the mesh size in cookbook2 | 07 April 2009, 22:44:24 UTC |
02125c8 | Eh Tan | 04 April 2009, 00:59:32 UTC | Renamed Seismic_model.c to Mineral_physics_models.c | 04 April 2009, 00:59:32 UTC |
1ff1990 | Eh Tan | 04 April 2009, 00:57:47 UTC | Fixed a few bugs in seismic output * fixed errors in the coefficient table * fixed typos * normalized drho by reference density profile * more digits in the prem radius table * added two specfem flags in PREM calculation. * "dv" output for debugging purpose, disabled. | 04 April 2009, 00:57:47 UTC |
795f249 | Eh Tan | 04 April 2009, 00:52:15 UTC | Inconsisent parameters with the manual | 04 April 2009, 00:52:15 UTC |
9de0f51 | Thorsten Becker | 28 March 2009, 02:29:00 UTC | Modified the way the compositional init with ggrd tracer can now select either a range of layers (>0) or a specific layer (<0). | 28 March 2009, 02:29:00 UTC |
dc9672a | Eh Tan | 20 March 2009, 23:36:56 UTC | limiting the number of multigrid cycles in solve_del2_u(). A new input parameter solver.vsolver.max_mg_cycles (default to 50) controls the max. number of multigrid cycles when solving the velocity. Without this parameter, when the stiffness matrix is near singular, the the multigrid solver may become an infinite loop. | 20 March 2009, 23:36:56 UTC |
25a7493 | Eh Tan | 18 March 2009, 19:39:54 UTC | Remove mgunitx etc from pyre input. This restores the behavior in v3.0 and earlier version. In pyre version, the mesh size is always specified by nodex etc. In C version, the mesh size is specified by nodex if Solver=cgrad, and by mgunitx and levels if Solver=multigrid. | 18 March 2009, 19:39:54 UTC |
04c4125 | Eh Tan | 17 March 2009, 23:42:01 UTC | Rollback the unintended change in velo output. | 17 March 2009, 23:42:01 UTC |
6e50c43 | Eh Tan | 17 March 2009, 00:18:09 UTC | Check return values of sscanf/fscanf. This fixs issue265. | 17 March 2009, 00:18:09 UTC |
213f078 | Eh Tan | 17 March 2009, 00:17:04 UTC | Removed obsolete function read_previous_field() | 17 March 2009, 00:17:04 UTC |
c1fbc1d | Eh Tan | 17 March 2009, 00:16:11 UTC | remove angular momentum from stokes solution by setting solver.vsolver.remove_angular_mementum=1 | 17 March 2009, 00:16:11 UTC |
e8525af | Eh Tan | 16 March 2009, 23:17:30 UTC | Converting temperature/composition fields to seismic velocities, in a format that specfem3d portal and v4.1 can understand. * solver.param.mineral_physics_mode: default to 3, the model of Tramper, Vacher and Vlaar's PEPI 2001. * output_optional=seismic will write coordinates and seismic velocities in binary files for specfem comsumption. * domain bounds file "datafile.domain": binary file with nproc * 10 doubles. The 10 doubles are (rmin, rmax) and four (theta, phi) pairs for the four bottom corner nodes. This file is written by rank-0 processor only. | 16 March 2009, 23:17:30 UTC |
ebff88c | Eh Tan | 16 March 2009, 22:52:08 UTC | add test cases for the stokes solver under bousinessq approximation and truncated anelastic liquid approximation | 16 March 2009, 22:52:08 UTC |
e1f9320 | Eh Tan | 16 March 2009, 22:50:06 UTC | seperate initialization and parameter input | 16 March 2009, 22:50:06 UTC |
ac6cb95 | Eh Tan | 16 March 2009, 22:44:42 UTC | change the default of mgunit | 16 March 2009, 22:44:42 UTC |
8c586f0 | Eh Tan | 16 March 2009, 22:43:37 UTC | Make the initialization step in pyre conforms to C version. This fixs a bug in r13270. | 16 March 2009, 22:43:37 UTC |
b655ef9 | Eh Tan | 10 March 2009, 00:58:34 UTC | Disable HC related stuff if USE_GGRD is not set | 10 March 2009, 00:58:34 UTC |
aaec6e9 | Thorsten Becker | 04 March 2009, 23:45:21 UTC | Changed ic flavors from 100 to 99 as per Eh's request | 04 March 2009, 23:45:21 UTC |
c7d458e | Thorsten Becker | 04 March 2009, 23:00:56 UTC | Changed ggrd flavor init for tracers with checkpoint override from 2 to 100 as per Eh's request. | 04 March 2009, 23:00:56 UTC |
0085e63 | Thorsten Becker | 28 February 2009, 21:23:21 UTC | Added ic_method_for_flavors=2, which works like the ggrd based tracer flavor init method ic_method_for_flavors=1 but overrides restart settings for tracers. This allows a compositional restart of a purely thermal run. | 28 February 2009, 21:23:21 UTC |
a9ab99d | Eh Tan | 20 February 2009, 00:09:53 UTC | Partially back out r19268,r9272,r11215, since the AVM stuff is redesigned. | 20 February 2009, 00:09:53 UTC |
66ad9f1 | Eh Tan | 19 February 2009, 23:46:52 UTC | Partially back out r11218, since the AVM stuff is redesigned. | 19 February 2009, 23:46:52 UTC |
c8a7398 | Eh Tan | 19 February 2009, 23:41:16 UTC | Partially back out r11221, since the AVM stuff is redesigned. | 19 February 2009, 23:41:16 UTC |
edc858e | Eh Tan | 19 February 2009, 23:35:07 UTC | Back out r11279, since the AVM stuff is redesigned. | 19 February 2009, 23:35:07 UTC |
f26ee7a | Leif Strand | 13 February 2009, 03:46:04 UTC | Fixed minor bug exposed by C-to-C++ translation (r14045). | 13 February 2009, 03:46:04 UTC |
8e778fa | Eh Tan | 06 February 2009, 19:51:18 UTC | Emit the warning message from a different place to get rid of static variables. Also, the macros TRUE and FALSE are undefined unless configured with GMT. | 06 February 2009, 19:51:18 UTC |
2dbb83a | Eh Tan | 22 December 2008, 18:48:50 UTC | Fixed a bug in viscosity rheol=7. Bug reported by Laura Alisic. | 22 December 2008, 18:48:50 UTC |
e6414a7 | Eh Tan | 19 November 2008, 06:49:04 UTC | Fixed two bugs in lg_pow(a, n), which computes a^n. One bug is when n=0, it returned a, instead of 1 Another bug is a^(n-1) was returned | 19 November 2008, 06:49:04 UTC |
00f2357 | Eh Tan | 18 November 2008, 23:56:37 UTC | One more fix to zero geoid arrays | 18 November 2008, 23:56:37 UTC |
d11440a | Eh Tan | 18 November 2008, 22:58:29 UTC | Zero'ed geoid arrays. | 18 November 2008, 22:58:29 UTC |
cd0f4b8 | Thorsten Becker | 18 November 2008, 06:55:46 UTC | Patched up Ggrd_handling for velocity grids close to pole. Addded flag to suppress check of incompressibility and pressure convergence. Partially, this is because for kinematic BCs pressure will only be constrained up to a constant, partially because mixed density/plate flow models show very poor convergence, still to be checked. | 18 November 2008, 06:55:46 UTC |
b1d50c6 | Thorsten Becker | 11 November 2008, 02:55:07 UTC | blob parameters were required even if tic_method != 2, changed such that the parameters were only read in if tic_method == 2 | 11 November 2008, 02:55:07 UTC |
54a6331 | Thorsten Becker | 11 November 2008, 02:48:36 UTC | Minor fix to allow compilation with USE_GGRD | 11 November 2008, 02:48:36 UTC |
3d888ba | Eh Tan | 07 November 2008, 23:35:06 UTC | forget to check in module/setProperties.c | 07 November 2008, 23:35:06 UTC |
bcddb74 | Eh Tan | 07 November 2008, 23:32:20 UTC | * Reserved 'tic_method=100' for user-defined initial temperature. * Read all parameters related to initial temperature regardless which tic_method is used. This will make the life easier when adding new tic_method. * Fixed a bug in 'tic_method=0' for regional model. The bug causes that the sinosoidal temperature perturbation is applied at 0th processor only. * Fixed a bug in 'tic_method=1' for regional model. 'tic_method=1' should generate a top TBL according to the input parameter 'half_space_age'. However, the code effectively multiply 25x to the age. Cookbook6 is using this tic_method, so its half_space_age needs to become 2500 to get the original temperature. * In 'tic_method=2', compute distance in Cartesian coordinate, instead of spherical coordinate. * Refactoring codes for temperature initial conditions. * Merged regional_construct_tic_from_input() and full_construct_tic_from_input() to construct_tic_from_input(). | 07 November 2008, 23:32:20 UTC |
ca13d76 | Eh Tan | 07 November 2008, 23:32:00 UTC | Reserved 'rheol=100' for user-defined viscosity law. | 07 November 2008, 23:32:00 UTC |
a0d3b96 | Eh Tan | 07 November 2008, 23:31:29 UTC | Added parameter 'use_cbf_topo' to pyre binding, default to off. | 07 November 2008, 23:31:29 UTC |
8ad1a60 | Eh Tan | 07 November 2008, 23:14:19 UTC | * Added back parameters mgunitx/mgunity/mgunitz. This partly reverts r13256. * Under pyre, uses mgunitx/mgunity/mgunitz and levels to compute nodex/nodey/nodez in multigrid solver. * Merged regional_global_derived_values() and full_global_derived_values() to global_derived_values(). | 07 November 2008, 23:14:19 UTC |
70e5d7d | Eh Tan | 05 November 2008, 00:19:05 UTC | Minor fix on node counting | 05 November 2008, 00:19:05 UTC |
00db2ff | Eh Tan | 04 November 2008, 22:50:02 UTC | Added a missing import | 04 November 2008, 22:50:02 UTC |
3c39d9c | Eh Tan | 04 November 2008, 21:02:08 UTC | zslice.py will slice opt files if present | 04 November 2008, 21:02:08 UTC |
629713a | Eh Tan | 04 November 2008, 21:01:26 UTC | Removed input parameters "mgunitx", "mgunity" and "mgunitz". Their values are inferred from other parameters. Also, the code checks that "levels" > 1 when multigrid solver is used. | 04 November 2008, 21:01:26 UTC |
dc54393 | Eh Tan | 04 November 2008, 21:00:01 UTC | Fixed a bug in outputting comp_el | 04 November 2008, 21:00:01 UTC |
4b191bb | Thorsten Becker | 01 November 2008, 18:37:20 UTC | Fixed two typos, velocity norm function was called where pressure norm should have been called. This leads to core dumps, all seems to be working now. | 01 November 2008, 18:37:20 UTC |
f8e52fc | Eh Tan | 31 October 2008, 20:18:13 UTC | Reporting correct nno and nel (number of nodes and elements) | 31 October 2008, 20:18:13 UTC |
ef7299c | Eh Tan | 29 October 2008, 23:22:10 UTC | (for r13196) Fixing convergence criterion for Stokes solver. * New functions global_v_norm2(), global_p_norm2(), global_div_norm2() compute the L2 norm of velocity, pressure, and divergence. * The L2 norm of the vector f is defined as: \sum_{all elements} ( \int_{element} (f*f dA) ) / A where A is the volume of the domain. * The convergence criterion is controlled by parameter "accuracy" under "[CitcomS.solver.vsolver]". The iteration stops if (norm(div(rho*V))/norm(V) < accuracy) or ( (norm(dV)/norm(V) < accuracy) and (norm(dP)/norm(P) < accuracy) for the last two iterations ) * The default value of "accuracy" becomes 1e-4. * The input parameters "tole_compressibility" and "relative_err_accuracy" under "[CitcomS.solver.vsolver]" are gone. * The mass matrix E->MASS becomes double precision (was single precision). | 29 October 2008, 23:22:10 UTC |
2ba507e | Eh Tan | 29 October 2008, 23:17:11 UTC | Fixed two bugs in vtk velocity output | 29 October 2008, 23:17:11 UTC |
f52d984 | Thorsten Becker | 23 October 2008, 23:35:27 UTC | Initial test implementation of netcdf grd based assignment of designated Euler vectors based on a code grd. | 23 October 2008, 23:35:27 UTC |
f2bb425 | Eh Tan | 04 September 2008, 22:46:39 UTC | multi-component chemical viscosity | 04 September 2008, 22:46:39 UTC |