started working on ring polymer contraction, create force_calc.f90
fixed barostat
misc - rewrite parts of energy estimator, try to fix pressure estimators
start working on ring polymer contraction (new forces routine)