Newer
Older
Dave Whipp
committed
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
--------------------------------------------------------------------------------
DOUAR-WSMP v0.1 - 23/03/2011
--------------------------------------------------------------------------------
Changes listed with revision number in the douar_newiso branch (r###) in the
following format:
--------------------------------------------------------------------------------
r###:
<Description of the change or additional feature>
- <Files affected, file2, file3...>
--------------------------------------------------------------------------------
Major changes since 2009:
--------------------------------------------------------------------------------
r175-183,195-196,198,216,218-223,225,227-233,236,244,247,249-255,258-261:
Added new define_isostasy_bc.f90 routine to allow velocity boundary
conditions to adjust to the crustal subsidence/uplift.
- define_cloud.f90, define_isostasy_bc.f90, DOUAR.f90,
interpolate_ov_on_osolve.f90, isostasy.f90, Makefile,
module_definitions.f90, octree_interpolate_six.f90 (new@251),
read_input_file.f90, write_global_output.f90
r234:
Added modifications suggested by J. Braun to be able to vertically compress
the FE mesh by factor vex
- define_cloud.f90, define_surface.f90, erosion.f90, make_matrix.f90,
make_pressure.f90, module_definitions.f90, move_cloud.f90,
move_surface.f90, read_input_file.f90, write_global_output.f90
r257,262,267,271-274,292-293,296:
Calculate and apply isostasy after surfaces are advected
- build_system_wsmp.f90, DOUAR.f90, define_cloud.f90, define_ov.f90,
define_surface.f90, interpolate_ov_on_osolve.f90, isostasy.f90,
make_cut.f90, make_matrix.f90, module_definitions.f90, move_surface.f90,
write_global_output.f90
r297-302:
Added ability of points to be removed when density becomes too high in a
given surface (surface level + 4)
- DOUAR.f90, refine_surface.f90
r306,310-311,315 (make_cut.f90 rewrite),317 (pressure_cut.f90 rewrite),327,330:
Added ability to choose to use a volumetric majority or minority to define
elemental material properties, in addition to using divFEM
- make_cut.f90, module_definitions.f90, pressure_cut.f90, read_input_file.f90
r331-333,336:
Revised vrm.f90 to have consistent formatting and 7 Drucker-Prager options.
Some are known to be bad, but available for comparison with older versions
of the code.
- vrm.f90
--------------------------------------------------------------------------------
Minor changes:
--------------------------------------------------------------------------------
r164:
e2d is now output as sqrt(J2d), not J2d (r164)
- do_leaf_measurements.f90
r165-173,184-188,192,194,197,199,256,263:
Elemental pressure, smoothed pressure, viscosity, pre-isostasy velocities,
isostasy-only velocities, is_plastic and countnode now output
- define_ov.f90, DOUAR.f90, input_of_outputs.txt, isostasy.f90,
module_definitions.f90, post.f90, write_global_output.f90
r189-190:
Gave post a Makefile like the rest of the code
- Makefile, VTK/Makefile
r212:
Fixed output formatting to be more flexible
- post.f90,
r235,238-242,328:
Added surface type(s)
12: rectangular emboss with specified slope
13: Simple plateau geometry with variable width
- create_surfaces.f90
r265:
Modified MPI usage to call mpi_ssend() rather than mpi_send()
- solve_with_pwgsmp.f90, solve_with_pwssmp.f90
r270:
Implemented fix from Anshul to make the temperature calculation work
properly
- solve_with_pwgsmp.f90
r305,309,314,321,323,325-327,329,441:
Added storage/output of elemental dilatation rate
- define_cloud.f90, define_ov.f90, do_leaf_measurements.f90, DoRuRe.f90,
DOUAR.f90, input_of_outputs.txt, module_definitions.f90, post.f90,
write_global_output.f90
r307,311,316,318,326-327:
Added ability to use an initial strain rate to place plastic materials on
yield
- DOUAR.f90, make_matrix.f90, module_definitions.f90, read_input_file.f90,
vrm.f90
r307, 311-312, 327:
Added option to specify an elemental bulk viscosity rather than penalty
value
- make_matrix.f90, make_pressure.f90, module_definitions.f90,
read_input_file.f90
r319:
Revised nonlinear iteration convergence criterion to determine convergence
has been reached when the maximum vector difference over the maximum vector
length from the previous iteration is below the given tolerance value
- compute_convergence_criterion.f90
r436:
Removed output of 'done visu.vtk' from post.f90 since it is not part of the
standard VTK file format
- post.f90
r439, 441-443:
Added input option for choosing to use a 2D form of the strain rate tensor
invariants. Assumes plane strain in the x-z plane.
- do_leaf_measurements.f90, improve_osolve.f90, make_matrix.f90,
module_definitions.f90, module_invariants.f90, read_input_file.f90
--------------------------------------------------------------------------------