sim_surface93.h

/*

  Copyright 2013, Vasudevan Venkateshwaran, Garde group @ RPI
  
  This program is free software: you can redistribute it and/or modify
  it under the terms of the GNU General Public License as published by
  the Free Software Foundation, either version 3 of the License, or
  (at your option) any later version.
  
  This program is distributed in the hope that it will be useful,
  but WITHOUT ANY WARRANTY; without even the implied warranty of
  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  GNU General Public License for more details.
  
  You should have received a copy of the GNU General Public License
  along with this program.  If not, see <http://www.gnu.org/licenses/>.

  sim_harmonicpin.h

  Header file containing the class which represents a uniform field applied 
  along the interface normal.

*/

#include <iostream>
#include <fstream>
#include "grid.h"
#include "constants.h"


#ifndef _surface93_h
#define _surface93_h

#include "sim_field.h"

class surface93: public Field
{
private:
    double zref;
    double sigma;    // LJ sigma of head groups (nm)   
    double epsilon;  // LJ epsilon of head groups (kJ/mol)
    double rho;      // Surface density of head groups (#/nm^2)
    double R_0;
    double rho_l;
    double d_skin;
    double sigmasq;

public:
    surface93();
    virtual ~surface93();
    void init ( int argc, char **argv, std::ofstream& fp );
    void calc_forces ( double* h, double *f, t_Grid* grid);
      
};

#endif /* _surface93_h */