# Test Coulomb potential                                                              
                                                                                           
units real                                                                                 
dimension 3                                                                                
boundary p p p                                                                             
atom_style full                                                                            	
timestep 1
                                                                                                                                          
pair_style lj/cut/coul/long 4                                                           
pair_modify mix arithmetic                                                                 
bond_style harmonic                                                                        
angle_style harmonic                                                                       
kspace_style ewald 0.0001                                                                  
                                                                                             
read_data Test.data                           
neighbor 0.2 bin                                                                                 
neigh_modify delay 0 every 5 check yes 

group Na type 1                            
group Cl type 2 

compute 1 Na group/group Cl pair yes kspace yes

thermo_style custom step c_1 
thermo 1   
run 1