surface tension calculation part is revised