doing NPT simulations varying p