#!/usr/bin/env python

import os 

cwd 		= os.getcwd()
eV2kJmol	= 96.48530749925794

# Pt(111) ############################################

Easym	= -.14209779E+03
Eads    = float( open( os.path.join( cwd, 'Pt111', 'OSZICAR' ) ).readlines()[-1].split()[4] ) 
E	= ( Easym - Eads + 0.041 ) * eV2kJmol

Zc	= 0.1582674527817457 * 22.2841654530509672

print " Pt(111)  %.1f kJ/mol  |  %.3f Ang " % ( E, Zc )

#######################################################
# Pt(211) ############################################

Easym	= -.11485880E+03
Eads    = float( open( os.path.join( cwd, 'Pt211', 'OSZICAR' ) ).readlines()[-1].split()[4] ) 
E	= ( Easym - Eads + 0.022 ) * eV2kJmol

Zc	= ( 0.5299351894442889-0.4048174595071856 ) * 21.8420385199999991

print " Pt(211)  %.1f kJ/mol  |  %.3f Ang " % ( E, Zc )

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