-------------------------------------------------------------------------------- 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 --------------------------------------------------------------------------------