add script to compute pressure a la paschalidis