| 36f5b32 | Leif Strand | 29 April 2009, 20:47:56 UTC | Merged r14275:14351 from trunk. | 29 April 2009, 20:47:56 UTC |
| 5ceef5d | Leif Strand | 22 April 2009, 19:42:53 UTC | Merged r14108:14275 from trunk. | 22 April 2009, 19:42:53 UTC |
| 5759f63 | Leif Strand | 09 April 2009, 03:05:05 UTC | For each function, place a unique prototype in a header file which corresponds to the source file in which the function is defined. Declare functions as "static" wherever possible. Removed "obsolete" source files. | 09 April 2009, 03:05:05 UTC |
| 13f829b | Leif Strand | 27 February 2009, 21:34:47 UTC | Merged r14046:14108 from trunk. | 27 February 2009, 21:34:47 UTC |
| 1a1c833 | Leif Strand | 13 February 2009, 03:39:35 UTC | New, unauthorized, experimental branch. Ported C code to C++, updating old K&R style functions, adding function prototypes, adding casts, etc. | 13 February 2009, 03:39:35 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 |
| fed8acd | Eh Tan | 04 September 2008, 20:18:53 UTC | Fixed issue 158: incorrect E->control.fi_max with coor=1 | 04 September 2008, 20:18:53 UTC |
| e2213ae | Eh Tan | 04 September 2008, 00:20:50 UTC | Some basic VTK parallel XML output * output in vts/pvts Structured VTK XML format * only these fields are available: coordinate, temperature, velocity and viscosity * coordinate and velocity are in Cartesian coordinates Some notes for future tasks: * adding other node- or element-based fields (pressure, stress, comp_nd, comp_el) is easy * tracer output should be in vtp format * surf, botm, and horiz_avg will be in vts format * geoid will still be as an ascii table | 04 September 2008, 00:20:50 UTC |
| 46cfb7a | Eh Tan | 04 September 2008, 00:12:24 UTC | forgot to checkin this change in r12799 | 04 September 2008, 00:12:24 UTC |
| 6ef0cc0 | Eh Tan | 04 September 2008, 00:11:36 UTC | pyre input for "self_gravitation" | 04 September 2008, 00:11:36 UTC |
| 9eff47c | Eh Tan | 04 September 2008, 00:10:51 UTC | self-gravitation geoid stuff * compute geoid due to internal buoyancy at the bottom * get rid of stress arrays * seperating topo effect on surface geoid into two part: surface topo and CMB topo. | 04 September 2008, 00:10:51 UTC |
| 867f7be | Eh Tan | 24 July 2008, 18:42:45 UTC | E->viscosity.zbase_layer is not init'd in pyre version. This fixes issue 155. | 24 July 2008, 18:42:45 UTC |
| 972d8eb | Thorsten Becker | 24 July 2008, 18:21:37 UTC | Both Eh and I were right and now the z_layer check should be finally fixed. | 24 July 2008, 18:21:37 UTC |
| ccc8ab9 | Thorsten Becker | 24 July 2008, 18:13:22 UTC | Fixed typo that provided assignment of zbase_layer thanks to Eh. | 24 July 2008, 18:13:22 UTC |
| cc2e6d9 | Thorsten Becker | 24 July 2008, 01:00:16 UTC | Added rheology option eight. The other merges are surprising. | 24 July 2008, 01:00:16 UTC |
| 55333fb | Thorsten Becker | 30 June 2008, 15:49:56 UTC | Moved the stress computation into stress output subroutine, for consistency with the style of the Output...c files. | 30 June 2008, 15:49:56 UTC |
| 871fd84 | Thorsten Becker | 30 June 2008, 02:12:56 UTC | Made CBF topography method and call to E->solver.parallel_communication_routs_s contingent on a "use_cbf_topo" parameter which is for now, by default, off, until the potential remaining memory bug (?) is tracked down. Made sure the stress tensor is computed before stress output, as CBF topo does not require the computation of the stress tensor. Added a z_layer input flag for zbase_layer other than the four layers used for the control of phase boundary depth. (Default is backward compatible.) | 30 June 2008, 02:12:56 UTC |
| eb573f1 | Thorsten Becker | 27 June 2008, 00:53:30 UTC | Check in of proper version for elliptical net rotation correction. | 27 June 2008, 00:53:30 UTC |
| 9f7396f | Thorsten Becker | 26 June 2008, 23:56:32 UTC | Reordered loop in get_buoyancy to speed up array access. Fixed omission of one elliptical correction, added a flag. | 26 June 2008, 23:56:32 UTC |
| fa90e33 | Thorsten Becker | 26 June 2008, 23:04:06 UTC | Replaced several trig function evaluations in Size_does_matter.c witht he precomputed SinCos arrays for speedup. Experimental ellipticity implementation has been improved somewhat, and made accessible only when -DALLOW_ELLIPTICAL is set. This was necessary to not slow the spherical version of the code down with if statements. | 26 June 2008, 23:04:06 UTC |
| cc26d94 | Thorsten Becker | 25 June 2008, 01:13:32 UTC | Modified net rotation computation for ellipticity. Minor additions (prototypes) to make the code compile smoothly under PGI. | 25 June 2008, 01:13:32 UTC |
| 6ada82f | Eh Tan | 24 June 2008, 21:18:37 UTC | Modified routs_s for lines. Added some debugging output and updated comments. | 24 June 2008, 21:18:37 UTC |
| a0ca2a8 | Eh Tan | 24 June 2008, 21:11:18 UTC | Only top/bottom processors need to communicate | 24 June 2008, 21:11:18 UTC |
| 30e6081 | Thorsten Becker | 21 June 2008, 22:28:51 UTC | Assignment of velocity boundary condition only called once at the top multigrid level during boundary condition assignment. Experimental implementation of mixed velocity/free slip boundary condition based on netcdf grids, in Ggrd_handling. | 21 June 2008, 22:28:51 UTC |
| e6aaa4d | Eh Tan | 20 June 2008, 23:01:28 UTC | Enable CBF topography | 20 June 2008, 23:01:28 UTC |
| 9cab369 | Eh Tan | 20 June 2008, 22:51:56 UTC | Added back get_CBF_topo(), from Shijie Zhong, with some minor modifications. Minor modifications: * some double arrays -> float arrays * calling exchange_snode_d() -> calling *_exchange_snode_f() | 20 June 2008, 22:51:56 UTC |
| 1f78588 | Eh Tan | 20 June 2008, 22:50:51 UTC | Resurrected full_parallel_communication_routs_s(), which is needed by CBF topo. Renamed exchange_snode_f() to full_exchange_snode_f() and regional_exchange_snode_f(). | 20 June 2008, 22:50:51 UTC |
| b0228c2 | Eh Tan | 20 June 2008, 22:42:30 UTC | Added get_global_1d_shape_fn_L() from Shijie Zhong | 20 June 2008, 22:42:30 UTC |
| d98cafb | Eh Tan | 20 June 2008, 22:41:57 UTC | Compute E->L and E->Lx 1d shape functions (1d is a misnomer, they are 2d actually). | 20 June 2008, 22:41:57 UTC |
| 2e17e1d | Eh Tan | 20 June 2008, 19:08:03 UTC | The function declaration was put at a wrong place. | 20 June 2008, 19:08:03 UTC |
| bf432bb | Eh Tan | 19 June 2008, 22:20:07 UTC | Change back the order of input arguments to match the manual. | 19 June 2008, 22:20:07 UTC |
| fff43af | Eh Tan | 19 June 2008, 21:43:25 UTC | The comment in rheol=4 doesn't match the code and is removed. Added option rheol=7, with a complete temperature- and pressure-dependent Arrhenius law. This fixed issue152. | 19 June 2008, 21:43:25 UTC |
| c3a51e5 | Eh Tan | 19 June 2008, 21:42:27 UTC | Added debugging (but disabled) function for spherical harmonics expansion. | 19 June 2008, 21:42:27 UTC |
| ab74087 | Eh Tan | 19 June 2008, 01:05:12 UTC | Send geoid arrays to vertical columns of processors. After this fix, the geoid output is the same for nprocz=1 and nprocz>1. | 19 June 2008, 01:05:12 UTC |
| 08ab04b | Eh Tan | 19 June 2008, 01:04:06 UTC | Compute E->surf_det (top surface area, ie. surface jacobian, of an element) on a unit sphere. Use E->surf_det, instead of the surface jacobian returned by get_global_1d_shape_fn(), to expand spherical harmonics. The surface jacobian returned by get_global_1d_shape_fn() is defined at the top radius of the current mesh. When nprocz>1, the top radius is not the same as the outer radius of the sphere for processors. As a result, the spherical harmonics expansion gave incorrect result when nprocz>1 in the previous revisions. | 19 June 2008, 01:04:06 UTC |
| d3025a8 | Thorsten Becker | 18 June 2008, 22:35:04 UTC | Speed up of computation by precomputing constants. | 18 June 2008, 22:35:04 UTC |
| 2d8894e | Eh Tan | 17 June 2008, 21:23:45 UTC | Compute geoid with self gravitation. Code contributed by Shijie Zhong. | 17 June 2008, 21:23:45 UTC |
| de569a3 | Thorsten Becker | 17 June 2008, 20:31:30 UTC | Added a few more comments and made a few loops more efficient. | 17 June 2008, 20:31:30 UTC |
| 67ed449 | Thorsten Becker | 14 June 2008, 00:00:58 UTC | Added a few comments on Rayleigh number netcdf control. | 14 June 2008, 00:00:58 UTC |
| 9735d41 | Eh Tan | 13 June 2008, 23:36:53 UTC | update NEWS item | 13 June 2008, 23:36:53 UTC |
| d5dc9d3 | Eh Tan | 13 June 2008, 23:35:38 UTC | Replaced MPI_COMM_WORLD by E->parallel.world. MPI_COMM_WORLD is the global communicator. E->parallel.world is the communicator for the current solver. These two are different in a coupled run. Before this fix, full-regional coupling will dead lock. | 13 June 2008, 23:35:38 UTC |
| adc7607 | Eh Tan | 13 June 2008, 23:34:08 UTC | Added some comments | 13 June 2008, 23:34:08 UTC |
| dbd7942 | Eh Tan | 13 June 2008, 23:33:29 UTC | To be consistent with C version, disable checkpoint point output at 0th step if restarting. | 13 June 2008, 23:33:29 UTC |
| 35878a7 | Eh Tan | 13 June 2008, 23:32:21 UTC | Correctly counting the global number of equations. This will slightly affect the results of the stokes solver. This avoid double counting the overlapping nodes between processors and the global norm of velocity vector. | 13 June 2008, 23:32:21 UTC |
| e5d4428 | Eh Tan | 13 June 2008, 23:25:41 UTC | Using firewall channel, since error channel doesn't stop the code. | 13 June 2008, 23:25:41 UTC |
| 94e1195 | Eh Tan | 13 June 2008, 23:24:32 UTC | Decalred variables as double to avoid float-double conversion later in the loop. | 13 June 2008, 23:24:32 UTC |
| b118cc2 | Eh Tan | 13 June 2008, 19:14:12 UTC | update the help message | 13 June 2008, 19:14:12 UTC |
| d2eccfa | Thorsten Becker | 12 June 2008, 21:22:41 UTC | Added tentative implementation of a rotating, elliptical Earth. Made restart for constant viscosity cases possible. Added alternative geoid computation place holder for self-gravitating mantle, to be filled out by Shijie's old code. | 12 June 2008, 21:22:41 UTC |
| 6f05c7a | Thorsten Becker | 30 May 2008, 18:42:57 UTC | Fixed typo in memory allocation procedure that led to core dump. | 30 May 2008, 18:42:57 UTC |
| de1b3e0 | Michael Gurnis | 12 May 2008, 21:22:42 UTC | fixed typo | 12 May 2008, 21:22:42 UTC |
| 1fbd9ec | Michael Gurnis | 12 May 2008, 21:16:55 UTC | Just a place holder for the next release | 12 May 2008, 21:16:55 UTC |
| 9a0fcce | Michael Gurnis | 05 May 2008, 19:17:31 UTC | Bug Fix | 05 May 2008, 19:17:31 UTC |
| 8ffc5f0 | Eh Tan | 01 May 2008, 22:59:14 UTC | Add python binding to update top surface temperatures | 01 May 2008, 22:59:14 UTC |
| e7216e6 | Michael Gurnis | 01 May 2008, 22:56:12 UTC | More additions to allow top surface temperatures to be read in a each time step. Should read in every time step for C only version. Additional changes for Pyre version to follow; will test more once these changes are made. Small changes to viscosity structures when reading in materials for each element. More documentation to follow. | 01 May 2008, 22:56:12 UTC |
| e137764 | Eh Tan | 01 May 2008, 00:51:25 UTC | Inputing tbcs_file & temperature_boundary_file | 01 May 2008, 00:51:25 UTC |
| 7499aea | Michael Gurnis | 30 April 2008, 23:11:18 UTC | Start to add new features for reading surface temperature for all top surface node points Added new variables to the CONTROL structure vbcs_file & temperature_boundary_fi le | 30 April 2008, 23:11:18 UTC |
| 1d4d90a | Thorsten Becker | 06 April 2008, 21:58:45 UTC | GMT/Netcdf grd input can now deal with velocity boundary conditions, material dependence, and local Rayleigh number in surface layers. | 06 April 2008, 21:58:45 UTC |
| 1866b6f | Eh Tan | 05 March 2008, 21:15:02 UTC | Close trace_log files when finalizing | 05 March 2008, 21:15:02 UTC |
| adc1f2d | Eh Tan | 05 March 2008, 21:07:02 UTC | Clean up unused arrays and functions | 05 March 2008, 21:07:02 UTC |
| b33dc3c | Eh Tan | 05 March 2008, 21:06:11 UTC | Write exclusive tracer time every 20 steps | 05 March 2008, 21:06:11 UTC |
| 331291e | Eh Tan | 05 March 2008, 21:04:21 UTC | Revert r11328 | 05 March 2008, 21:04:21 UTC |
| 1709362 | Eh Tan | 05 March 2008, 20:50:31 UTC | Write exclusive tracer time every 20 steps | 05 March 2008, 20:50:31 UTC |
| 37cc271 | Eh Tan | 28 February 2008, 01:03:50 UTC | Removed artificial tolerance of empty elements. This greatly decreases the minimal number of tracers per element required for a given mesh. This fixed issue140. | 28 February 2008, 01:03:50 UTC |
| cc3e43d | Eh Tan | 28 February 2008, 01:00:54 UTC | Minor tweak | 28 February 2008, 01:00:54 UTC |
| e15117c | Thorsten Becker | 27 February 2008, 07:42:44 UTC | Built surface velocity BC support via Netcdf grd files back in. This needs more testing, but it would be good to have the framework in place before new additions are made. Also, age control is not yet implemented, but will be soon. | 27 February 2008, 07:42:44 UTC |
| 332d337 | Eh Tan | 21 February 2008, 01:25:34 UTC | Removed unnecessary call to construct_mat_group() | 21 February 2008, 01:25:34 UTC |
| c42f6fd | Eh Tan | 21 February 2008, 01:03:43 UTC | Finished interpolating fields onto tracers | 21 February 2008, 01:03:43 UTC |
| 222c941 | Eh Tan | 21 February 2008, 01:02:41 UTC | Renamed fix_phi() to fix_angle() | 21 February 2008, 01:02:41 UTC |
| 2bd8f42 | Eh Tan | 21 February 2008, 01:01:09 UTC | Fixed issue139. Bug for tracers near South pole. | 21 February 2008, 01:01:09 UTC |
| c0e2815 | Eh Tan | 21 February 2008, 01:00:21 UTC | Get new tracers in post processing stage | 21 February 2008, 01:00:21 UTC |
| 346ba4e | Eh Tan | 21 February 2008, 00:56:29 UTC | Read multiple extraq | 21 February 2008, 00:56:29 UTC |
| 5cc21b8 | Eh Tan | 21 February 2008, 00:55:41 UTC | Write averaged composition in horiz_avg output. The HDF5 part of output is not implemented... | 21 February 2008, 00:55:41 UTC |
| 892dd0d | Eh Tan | 21 February 2008, 00:54:25 UTC | Removed horizontal average from output_interpolated_fields() | 21 February 2008, 00:54:25 UTC |
| 89dc1c4 | Eh Tan | 21 February 2008, 00:53:35 UTC | Compute horizontal average for composition fields | 21 February 2008, 00:53:35 UTC |
| a9e6ba7 | Thorsten Becker | 17 February 2008, 00:14:06 UTC | Modified Citcom.c such that E->mat gets assigned based on element depth (as was default before), even if mat_control is not selected. Added (back in?) ggrd netcdf type VIP prefactor material control via ggrd_mat_control flags (as documented in Instructions.c) Renamed some of the ggrd stuff to make it more modular. | 17 February 2008, 00:14:06 UTC |
| 797ffec | Eh Tan | 09 February 2008, 00:14:35 UTC | Interpolate fields in post_processing() * The parameter "itracer_interpolate_fields" is the same as "model_type" in the AVM code. * If itracer_interpolate_fields is 1, 2, or 3, it interpolates temperature and compositions. The horizonal average is not removed yet. | 09 February 2008, 00:14:35 UTC |
| 0262f7b | Eh Tan | 08 February 2008, 23:58:57 UTC | Refactored the functions. * get shape functions on tracer location * interpolate data using acquired shape functions | 08 February 2008, 23:58:57 UTC |
