6d90070 | Eh Tan | 30 October 2007, 22:29:44 UTC | Minor modification | 30 October 2007, 22:29:44 UTC |
1f562ce | Eh Tan | 30 October 2007, 21:54:14 UTC | Added the C program project_geoid to the makefile | 30 October 2007, 21:54:14 UTC |
1ea69a0 | Eh Tan | 30 October 2007, 21:53:50 UTC | A post-processing program to project geoid coefficents onto a regular (longitude, latitude) mesh | 30 October 2007, 21:53:50 UTC |
4d7f4a3 | Eh Tan | 30 October 2007, 21:50:52 UTC | Fixed a bug in dimensionalizing density. Provided the formula of geoid calculation in the comments. Rearranged the order of functions. | 30 October 2007, 21:50:52 UTC |
f40dafa | Eh Tan | 30 October 2007, 21:49:58 UTC | Compute d(rho)/dr/rho from rho(r) | 30 October 2007, 21:49:58 UTC |
8102449 | Eh Tan | 22 October 2007, 20:57:45 UTC | Updated file ChangeLog to r8170 | 22 October 2007, 20:57:45 UTC |
26c3d4a | Eh Tan | 22 October 2007, 20:26:44 UTC | Updated NEWS | 22 October 2007, 20:26:44 UTC |
a05de44 | Eh Tan | 22 October 2007, 19:46:26 UTC | Minor update | 22 October 2007, 19:46:26 UTC |
7823da3 | Eh Tan | 22 October 2007, 19:44:28 UTC | Simpilify INSTALL | 22 October 2007, 19:44:28 UTC |
4832d3d | Eh Tan | 22 October 2007, 19:03:41 UTC | Eliminated unnecessary computation | 22 October 2007, 19:03:41 UTC |
e8a9d98 | Eh Tan | 22 October 2007, 19:03:00 UTC | Variable material properties in latent heat correction on adiabatic cooling | 22 October 2007, 19:03:00 UTC |
5b2e55b | Eh Tan | 22 October 2007, 19:02:29 UTC | Added some comments | 22 October 2007, 19:02:29 UTC |
a27d3c6 | Eh Tan | 19 October 2007, 19:24:45 UTC | Renamed control.width* -> control.inv_width* for clarity | 19 October 2007, 19:24:45 UTC |
88a2e23 | Eh Tan | 18 October 2007, 18:53:49 UTC | Renamed coupler inventory: dimensional -> si_unit and transformational -> cartesian | 18 October 2007, 18:53:49 UTC |
f614f83 | Eh Tan | 17 October 2007, 22:44:38 UTC | Minor change in the examples | 17 October 2007, 22:44:38 UTC |
25c9b53 | Eh Tan | 17 October 2007, 20:56:25 UTC | Added comments on where C code is replaced by python code in Citcom.c | 17 October 2007, 20:56:25 UTC |
8d2bbbe | Eh Tan | 17 October 2007, 20:55:17 UTC | Put the rest of r8111 back | 17 October 2007, 20:55:17 UTC |
9676e09 | Eh Tan | 17 October 2007, 20:33:21 UTC | Fixed input diffusivity parameter in the examples, from r8111. | 17 October 2007, 20:33:21 UTC |
69e4ad8 | Eh Tan | 17 October 2007, 20:32:54 UTC | Added back minor modification from r8111 | 17 October 2007, 20:32:54 UTC |
7732890 | Eh Tan | 17 October 2007, 18:40:42 UTC | Fixed a typo in the example | 17 October 2007, 18:40:42 UTC |
6a8e4b5 | Eh Tan | 17 October 2007, 18:36:29 UTC | Fixed a bug in remove_rigid_rot() | 17 October 2007, 18:36:29 UTC |
93471fe | Eh Tan | 17 October 2007, 18:35:03 UTC | Roll back r8111, reading velo files for init. temperature stays as tic_method=-1. Other changes in r8111 will be put back later. | 17 October 2007, 18:35:03 UTC |
1602d87 | Thorsten Becker | 13 October 2007, 19:07:53 UTC | Fixed rigid rotation code bug. Added restart=2 option which will not use checkpoints, but tic_method=-1 style temperature input (formally realized via tic_method=-1, but this should be cleaner) fixed input diffusivity parameter in full example file. | 13 October 2007, 19:07:53 UTC |
334bd72 | Eh Tan | 12 October 2007, 20:24:17 UTC | Turn off monitor_max_T in the test | 12 October 2007, 20:24:17 UTC |
fdf7f72 | Eh Tan | 12 October 2007, 20:23:53 UTC | Boundary class is for regional model only, add a checking for it. | 12 October 2007, 20:23:53 UTC |
17b3276 | Eh Tan | 11 October 2007, 01:23:50 UTC | Cookbook 9: nested solver coupling + coor_refine + tsolver tuning | 11 October 2007, 01:23:50 UTC |
0f0b81a | Thorsten Becker | 11 October 2007, 00:03:40 UTC | Wiser now, commented in the E->fp output | 11 October 2007, 00:03:40 UTC |
4af480b | Thorsten Becker | 10 October 2007, 23:40:38 UTC | Added coor=3 option for radial node spacing a la CitcomCU. Streamlined ggrd temperature init. Added file flush to heat flow output. | 10 October 2007, 23:40:38 UTC |
9c14928 | Eh Tan | 10 October 2007, 21:00:13 UTC | Add an assertion | 10 October 2007, 21:00:13 UTC |
3c5bee2 | Eh Tan | 10 October 2007, 20:11:00 UTC | Remove rigid body rotation from the velocity solution from global model. A input parameter 'solver.vsolver.remove_rigid_rotation', default to on, indicate whether to remove the rigid body rotation. For free-slip model, this parameter should be on. For model with imposed plate velocity, it is advised to turn off this parameter. | 10 October 2007, 20:11:00 UTC |
ee814bc | Eh Tan | 05 October 2007, 20:02:02 UTC | Renamed excludeTop -> exclude_top, and excludeBottom -> exclude_bottom. (No other CitcomS parameter uses CamelCase.) Reworked the interaction of initial temperature and coupling. Created a new parameter exchange_initial_temperature. | 05 October 2007, 20:02:02 UTC |
b0f818a | Eh Tan | 05 October 2007, 20:00:32 UTC | Strengthen the default convergence criterion for the 3rd iteration in compressible cg solver | 05 October 2007, 20:00:32 UTC |
72ad845 | Eh Tan | 05 October 2007, 19:58:54 UTC | Disable a debugging output | 05 October 2007, 19:58:54 UTC |
4bcf0ad | Eh Tan | 05 October 2007, 19:58:31 UTC | Add licence text to the multi-couplers | 05 October 2007, 19:58:31 UTC |
820d139 | Eh Tan | 05 October 2007, 19:57:48 UTC | Turn off filter_temp by default since it's not in the energy equation | 05 October 2007, 19:57:48 UTC |
4e87245 | Thorsten Becker | 04 October 2007, 19:35:04 UTC | Fxied material assignment bug in PDEPV (0 vs. 1 array issue...) | 04 October 2007, 19:35:04 UTC |
cc73a80 | Eh Tan | 01 October 2007, 20:38:22 UTC | Updated the cookbook scripts | 01 October 2007, 20:38:22 UTC |
f4f2b13 | Eh Tan | 01 October 2007, 19:26:40 UTC | Added Cookbook8 for compressible model + geoid + checkpoint + vsolver tuning | 01 October 2007, 19:26:40 UTC |
ae31e95 | Thorsten Becker | 28 September 2007, 23:36:15 UTC | Fixed restart I/O bug in Output_gzdir. | 28 September 2007, 23:36:15 UTC |
e429e26 | Eh Tan | 27 September 2007, 19:33:57 UTC | Added function prototype and fixed a few errors | 27 September 2007, 19:33:57 UTC |
a79e4a9 | Eh Tan | 27 September 2007, 19:07:57 UTC | Shrink the size of checkpoint files | 27 September 2007, 19:07:57 UTC |
02acf94 | Eh Tan | 27 September 2007, 19:07:33 UTC | Enable coor_refine for regional model, only radial coordinate is refined | 27 September 2007, 19:07:33 UTC |
e29d3c2 | Eh Tan | 25 September 2007, 19:40:08 UTC | Fixed a bug in vtkio | 25 September 2007, 19:40:08 UTC |
f392355 | Eh Tan | 25 September 2007, 01:31:30 UTC | Introducing choice of reference state: When solver.param.reference_state=0, the reference state is read from a file (solver.param.refstate_file) When solver.param.reference_state=1 (default), the reference state is calculated using Adams-Williamson equation of state. | 25 September 2007, 01:31:30 UTC |
4170440 | Eh Tan | 24 September 2007, 21:29:24 UTC | Detect the presence of libz and enable ascii-gz at configure time | 24 September 2007, 21:29:24 UTC |
95c453f | Eh Tan | 24 September 2007, 21:28:26 UTC | Comment out varialbe thermal conductivity. Removed unused struct members. | 24 September 2007, 21:28:26 UTC |
e4cd242 | Eh Tan | 24 September 2007, 21:28:05 UTC | Multiply rho*cp to DT/dt | 24 September 2007, 21:28:05 UTC |
6a4ff8e | Eh Tan | 19 September 2007, 20:34:37 UTC | Don't write to log file, since only rank-0 proc can access it | 19 September 2007, 20:34:37 UTC |
37b7dd0 | Eh Tan | 19 September 2007, 20:34:22 UTC | Set default value for tracer/composition flags earlier | 19 September 2007, 20:34:22 UTC |
a3387ca | Eh Tan | 19 September 2007, 20:33:42 UTC | Shrink the size of arrays | 19 September 2007, 20:33:42 UTC |
3488a3f | Eh Tan | 19 September 2007, 20:31:11 UTC | Shrink the size of arrays | 19 September 2007, 20:31:11 UTC |
37e37f8 | Eh Tan | 19 September 2007, 20:30:44 UTC | Fixed uninit'd heat flux and geoid arrays bugs, as reported by valgrind | 19 September 2007, 20:30:44 UTC |
dec899b | Thorsten Becker | 19 September 2007, 00:28:00 UTC | Minor bug fix in Output_gzidir, output_heating function didn't work. | 19 September 2007, 00:28:00 UTC |
c89da81 | Eh Tan | 18 September 2007, 23:08:06 UTC | Fixed a memory bug in checkpoint | 18 September 2007, 23:08:06 UTC |
995fff8 | Eh Tan | 18 September 2007, 23:07:48 UTC | Setup composition for regional model | 18 September 2007, 23:07:48 UTC |
f23ffbf | Leif Strand | 18 September 2007, 01:52:50 UTC | Fix for Issue129: "CitcomS configured without exchanger still links with -lExchanger". Also, use CIT_HEADER_MPI so that configuring --with-exchanger using MPICH2 works (see Issue14 and Issue57). | 18 September 2007, 01:52:50 UTC |
a2fd2a4 | Eh Tan | 17 September 2007, 20:30:30 UTC | Generalized init. composition from neighbors for multi-component composition | 17 September 2007, 20:30:30 UTC |
c10e388 | Eh Tan | 17 September 2007, 19:07:29 UTC | Using neighboring elements to determine the initial composition, if only a few elements are empty. | 17 September 2007, 19:07:29 UTC |
e1af673 | Eh Tan | 15 September 2007, 00:17:24 UTC | Revert accidental commit in part of r7957 | 15 September 2007, 00:17:24 UTC |
94cd779 | Eh Tan | 14 September 2007, 23:39:59 UTC | Rearranged computation of latent heat and the test condition of phase change | 14 September 2007, 23:39:59 UTC |
fedce4a | Eh Tan | 14 September 2007, 23:11:46 UTC | Added the contribution of latent heat to the diffusion term | 14 September 2007, 23:11:46 UTC |
ffe1a6f | Eh Tan | 14 September 2007, 23:11:27 UTC | Replace abs() by fabs() | 14 September 2007, 23:11:27 UTC |
d4bc483 | Eh Tan | 14 September 2007, 23:10:59 UTC | Consevert energy (sum of rho*cp*T) instead of temperature (sum of T) in the Lenardic filter, also fixed a long-standing bug (abs -> fabs) | 14 September 2007, 23:10:59 UTC |
d13bdd5 | Eh Tan | 14 September 2007, 23:09:14 UTC | Changed the magic numbers (0.0 and 1.0) to const variables, reindent the code. | 14 September 2007, 23:09:14 UTC |
9189bab | Eh Tan | 13 September 2007, 22:36:57 UTC | Comment out adiabatic temperature profile for now, since it is not used in the code. | 13 September 2007, 22:36:57 UTC |
fcd95c8 | Eh Tan | 12 September 2007, 19:24:56 UTC | Fixed a problem when the input string is empty | 12 September 2007, 19:24:56 UTC |
249d790 | Eh Tan | 12 September 2007, 19:24:29 UTC | Added heating file to optional output | 12 September 2007, 19:24:29 UTC |
cb39951 | Eh Tan | 12 September 2007, 19:24:15 UTC | Keep the load vector intact in the bicg solver | 12 September 2007, 19:24:15 UTC |
0b02b96 | Eh Tan | 12 September 2007, 19:23:56 UTC | Using drand48() instead of rand() to improve the quality of random tracers. Added some debugging output if the tracer count is low. | 12 September 2007, 19:23:56 UTC |
f8219a2 | Eh Tan | 12 September 2007, 19:22:05 UTC | Add new input boolean parameter 'itracer_warnings'. Default to on. When it is off, the exit condition of tracer is relaxed. | 12 September 2007, 19:22:05 UTC |
6b934e2 | Eh Tan | 12 September 2007, 19:19:50 UTC | Removed the test for F=0 load vector bug. The bug was fixed in r7940. | 12 September 2007, 19:19:50 UTC |
d6e4a43 | Eh Tan | 07 September 2007, 22:53:22 UTC | Fixed tsolver. Fixed a bug in TMass and another in convert element # to nz. | 07 September 2007, 22:53:22 UTC |
0d53949 | Eh Tan | 05 September 2007, 22:35:30 UTC | Fixed a few memory problems. Shrink the size of several arrays related to Stokes eqn. Moved jacobi() to Obsolete.c | 05 September 2007, 22:35:30 UTC |
7273549 | Eh Tan | 05 September 2007, 22:21:12 UTC | Forgot to update surfacee_tem | 05 September 2007, 22:21:12 UTC |
c77dfd0 | Eh Tan | 05 September 2007, 22:20:46 UTC | Set the min of 2nd strain rate invariant to 1e-16 to prevent infinite viscosity when SDEPV=on | 05 September 2007, 22:20:46 UTC |
681786d | Eh Tan | 05 September 2007, 22:20:08 UTC | Init'd timestep variables | 05 September 2007, 22:20:08 UTC |
e3ccb79 | Thorsten Becker | 04 September 2007, 23:09:14 UTC | Added test for F=0 load vector bug, moved open_qfiles to after set_elapsed time. | 04 September 2007, 23:09:14 UTC |
e3f6845 | Eh Tan | 29 August 2007, 23:50:51 UTC | Forgot to update the makefile | 29 August 2007, 23:50:51 UTC |
4b3831f | Eh Tan | 29 August 2007, 23:39:17 UTC | Removed unused velocity bc files | 29 August 2007, 23:39:17 UTC |
f0d9ebe | Eh Tan | 29 August 2007, 23:18:52 UTC | Restructing the repository layout for CitcomS documentation | 29 August 2007, 23:18:52 UTC |
e06f544 | Eh Tan | 29 August 2007, 20:55:12 UTC | Forgot to update r7902 | 29 August 2007, 20:55:12 UTC |
3711a41 | Eh Tan | 29 August 2007, 20:48:05 UTC | Replaced output_open() with output_open_mode() Append the log and time files if restarting | 29 August 2007, 20:48:05 UTC |
a981bca | Eh Tan | 29 August 2007, 20:47:28 UTC | Minor fix | 29 August 2007, 20:47:28 UTC |
d1f0b86 | Thorsten Becker | 29 August 2007, 01:46:41 UTC | Missed to merge one modification after the conflict. | 29 August 2007, 01:46:41 UTC |
13492bf | Thorsten Becker | 29 August 2007, 00:51:20 UTC | When checking in my changes, I had to resolve a conflict for lib/Instructions.c by hand, which I hope I did properly. Here are my changes: - renamed CONTOL structure members ORTHO and ORTHOZ to CITCOM_ORTHO and CITCOM_ORTHOZ Those were never used and conflicted with definitions in GMT gmt.h - Added a higher frequency heat flow output option like so: write_q_files=1 # option to write heat flux to files qt.dat and qb.dat # at intervals smaller than storage_spacing (0) - Added the capability to read in initial temperatures from netcdf grd files, if -USE_GGRD is used. Lot of options, like so: # # read initial temperature conditions from grd files (default values in parentheses) # tic_method=4 # read initial temperature from netcdf GRD files (off) ggrd_tinit_scale_with_prem=off # scale the temperature with PREM densities (off) ggrd_tinit_scale=1.0 # scaling factor to apply to read in scalars f (1.0) ggrd_tinit_offset=-0.5 # offset, T = f * scale + offset + tm (0.0) # where tm is the mean between top and bottom TBC values # if the bottom is flux, will use 1 for bottom TBC value ggrd_tinit_gfile="../../data/tomography/s20a_smean_new_age/t" # prefix of grd files, will # try to read t.1.grd, t.2.grd ... t.n.grd # where n is the number of layers in the depth file ggrd_tinit_dfile="../../data/tomography/s20a_smean_new_age/tdepth.dat" # file with layer depths in km from bottom up ggrd_tinit_prem_file="../../progs/src/hc-svn/prem/prem.dat" # PREM data file ggrd_tinit_override_tbc=on # override temperature boundary conditions (off) | 29 August 2007, 00:51:20 UTC |
3ba2a04 | Eh Tan | 28 August 2007, 01:10:18 UTC | Defaulting datadir_old to the current directory | 28 August 2007, 01:10:18 UTC |
220393d | Eh Tan | 28 August 2007, 01:09:14 UTC | Added various heating term to the energy equation. A new input parameter: solver.surfaceT for non-dimensional surface temperature. Removed old (and unused) input parameter: solver.const.surf_temp. | 28 August 2007, 01:09:14 UTC |
557acb3 | Eh Tan | 27 August 2007, 19:34:44 UTC | Fixed a index bug in stokes solver | 27 August 2007, 19:34:44 UTC |
8347039 | Eh Tan | 27 August 2007, 19:34:26 UTC | Add pyre input for new parameters introduced by TWB | 27 August 2007, 19:34:26 UTC |
1a41a1d | Eh Tan | 27 August 2007, 19:34:08 UTC | Clean up the code and declared function prototypes to avoid compiler warnings | 27 August 2007, 19:34:08 UTC |
12591d6 | Eh Tan | 23 August 2007, 20:20:24 UTC | Compute total visc and adi heating | 23 August 2007, 20:20:24 UTC |
fce4d3d | Eh Tan | 23 August 2007, 20:19:33 UTC | Minor refactoring | 23 August 2007, 20:19:33 UTC |
4bc0ea7 | Eh Tan | 23 August 2007, 20:18:52 UTC | Compute phase change wrt reference density and gravity profiles. | 23 August 2007, 20:18:52 UTC |
d0098c3 | Eh Tan | 23 August 2007, 20:18:18 UTC | Avoid reading unused input parameters | 23 August 2007, 20:18:18 UTC |
1268a26 | Eh Tan | 23 August 2007, 20:17:33 UTC | assert solver.tsolver.monitor_max_T is off in coupled run and remove debugging output | 23 August 2007, 20:17:33 UTC |
ae9b86a | Eh Tan | 23 August 2007, 20:16:25 UTC | Renamed E->control.filter_temperature to E->advection.filter_temperature | 23 August 2007, 20:16:25 UTC |
8233668 | Eh Tan | 23 August 2007, 20:16:10 UTC | Added new input parameter solver.tsolver.monitor_max_T When it is on (default), if the max temperature changes too much between timestep, the temperature field is restored at the tsolver is called using half of the timestep size. | 23 August 2007, 20:16:10 UTC |
447f924 | Eh Tan | 23 August 2007, 20:15:50 UTC | Removed unused var. in struct advection | 23 August 2007, 20:15:50 UTC |
80a070b | Eh Tan | 23 August 2007, 20:14:52 UTC | A new way to compute TMass | 23 August 2007, 20:14:52 UTC |
6d31d51 | Eh Tan | 23 August 2007, 20:14:04 UTC | Minor refactoring: add comments, rename variables | 23 August 2007, 20:14:04 UTC |