* Edit 'Makefile', modify the line 'F90 = ...' for your F90 compiler. Modify
  'CC = ...' for C compiler if necessary.
* Run "make" to build an optimized binary.
* Run "make debug=1" to build a debugging binary.
* The binary is called 'dynearthsol2'.


* Edit 'resolution' in 'parameters.h' to change mesh size.
* Edit 'strain_correction' in 'parameters.h' to disable De Micheli's correction.
* Edit 'maxvbcval' in 'parameters.h' to change the speed of boundaries.
* Edit 'apply_bcs()' in 'field.f90' to modify the boundary condition.
* Edit 'crate_plstrain()' in 'field.f90' to modify the initial weak zone.
* Edit 'create_temperature()' in 'field.f90' to modify the initial temperature.
* Edit 'rheol_type' in 'rheol.f90' to modify rheology law.
* Edit 'markerspacing' in 'marker.f90' to modify the average distance between
  initial location of markers.
* Edit material properties in 'matprop.f90'.
* Edit 'output_frequency' and 'nloops' in 'dynearthsol2.f90' to modify time 
* Edit 'restart' in 'dynearthsol2.f90' to restart from a previous a run.
* Edit 'NSTRESS' in 'constants.h' to choose 2D or 3D plain strain formulation.
  Also edit 'nstress' in 'dyn2vtk.py' accordingly.
* The output files '*.0' and '*.rs' will be saved in the current directory.

* Run "python dyn2vtk.py" to convert the binary output to VTK files.
* Some of the simulation outputs might be disabled. Edit dyn2vtk.py and
  'output()' in 'dynearthsol2.f90' to disable/enable them.
* Plot the VTK files with Paraview or LLNL's Visit program.


This program is free software: you can redistribute it and/or modify
it under the terms of the MIT / X Windows System license (see the
file LICENSE for the full text).

These files: triangle.c, triangle.h, qsort.f90 are distributed by their
own license(s).

by Eh Tan, Oct. 2012