Script for calculating interfacial tension

Dear All,

I would like to share a script for calculating interfacial tension. The script header describes how it works:

###################################################################

#

# INTERFACIAL TENSION

#

# Author:     Jason DeJoannis

# Date:     Sept. 23, 2012

#

# Calculates the average interfacial tension from an atomistic

# or Mesocite trajectory using the stress tensor stored in

# each frame (Stress->Eij). The assumption is that there

# are two interfaces and the box size is fixed.

#

#    gamma      =  0.5 L < P_normal - P_lateral >

#    L         = box length in direction normal to surfaces

#    P_normal      = pressure in normal direction

#    P_lateral     = pressure in lateral direction

#

########## USER DEFINED ###########################################

# Surface normal direction X, Y or Z

my \$surfnormal = "Z";

# Trajectory with stress tensor

my \$doc = \$Documents{"water slab.xtd"};

# Begin with this frame number

my \$startFrame = 1;

###################################################################

ift.bmp