diff --git a/.travis.yml b/.travis.yml
index 15f50f3..af6b2d3 100644
--- a/.travis.yml
+++ b/.travis.yml
@@ -2,8 +2,14 @@ language: bash
 before_install: sudo apt-get install python-tk
 
 language: python
-python: 2.7
-install: python2 setup -t
+python: 
+  - 2.7
+  - 3.5
+  - 3.9
+install:
+  - python --version
+  - pip install coverage
+  - python setup -t
 
 services:
   - xvfb
@@ -12,4 +18,9 @@ before_script:
   - "export DISPLAY=:99.0"
   - sleep 3 # give xvfb some time to start
 
-script: ./run
+script:
+  - python run
+  - coverage run -m unittest discover
+
+after_success:
+  - bash <(curl -s https://codecov.io/bash)
diff --git a/LICENSE b/LICENSE
index eaaff3b..16f2774 100644
--- a/LICENSE
+++ b/LICENSE
@@ -1,6 +1,6 @@
 MIT License
 
-Copyright (c) 2018 Hernan Chavez Thielemann
+Copyright (c) 2018-2021 Hernan Chavez Thielemann
 
 Permission is hereby granted, free of charge, to any person obtaining a copy
 of this software and associated documentation files (the "Software"), to deal
diff --git a/README.md b/README.md
index 5417746..478790d 100644
--- a/README.md
+++ b/README.md
@@ -1,6 +1,6 @@
- 
-python 2.7: [![Build Status](https://travis-ci.org/hernanchavezthielemann/GRO2LAM.svg?branch=27ene19)](https://travis-ci.org/hernanchavezthielemann/GRO2LAM)
+python 2.7 & 3.5-3.9: [![Build Status](https://travis-ci.com/hernanchavezthielemann/GRO2LAM.svg?branch=py3_compatibility)](https://travis-ci.com/hernanchavezthielemann/GRO2LAM)
 
+python 2.7 & 3.5-3.9: [![Build Status](https://travis-ci.org/hernanchavezthielemann/GRO2LAM.svg?branch=py3_compatibility)](https://travis-ci.org/hernanchavezthielemann/GRO2LAM)
 # GRO2LAM
    Gromacs to Lammps simulation converter
    
diff --git a/lib/grotolam.py b/lib/grotolam.py
index 9243f26..c84dd48 100644
--- a/lib/grotolam.py
+++ b/lib/grotolam.py
@@ -3,6 +3,7 @@
 from lib.misc.version import __version__
 from lib.misc.file import write_file
 from lib.gui.main_gui import launch_gui
+from lib.misc.display import show
 
 class Grotolam(object):
     def __init__(self, master=None):
@@ -27,14 +28,14 @@ def grotolam_launcher():
     ''' at this point this is just the name, quite deceptive ;) 
         maybe <run> can call directly the GUI in the future,
         but this handles the library'''
-    print '\n'+' '*20+__version__+'\n'
+    show( '\n' + ' '*20 + __version__ + '\n')
     
     launch_gui()
     
 def grotolam_launcher_test():
     
-    print '-'*21+'  TESTING  '+'-'*21+'\n\n'
-    print '\n'+' '*20+__version__+'\n'
+    show( '-'*21 + '  TESTING  ' + '-'*21 + '\n\n')
+    show( '\n' + ' '*20 + __version__ + '\n')
     launch_gui( True)
     
     #conv_data, ckbuttons, solva_tags,_quit_flag_ = launch_gui()
diff --git a/lib/gui/__init__.pyc b/lib/gui/__init__.pyc
new file mode 100644
index 0000000..6682dd6
Binary files /dev/null and b/lib/gui/__init__.pyc differ
diff --git a/lib/gui/conversion_gui.py b/lib/gui/conversion_gui.py
index 9683bf7..ab9fed8 100644
--- a/lib/gui/conversion_gui.py
+++ b/lib/gui/conversion_gui.py
@@ -7,13 +7,13 @@
 #///    Packages and globals definitions are here   ///
 #------------------------------------------------------
 
-from Tkinter import Frame, Label, StringVar, Button # , IntVar
+from lib.gui.tk_lib import Frame, Label, StringVar, Button # , IntVar
+from lib.gui.tk_lib import bottom_hline_deco, format_dec, Drop_Down_List
+from lib.gui.tk_lib import createmenubar, create_entry, get_entriesvalue
+from lib.gui.custom_row import File_Entry
 
-from tk_lib import bottom_hline_deco, format_dec, Drop_Down_List
-from tk_lib import createmenubar, create_entry, get_entriesvalue
-from custom_row import File_Entry
-
-from lib.misc.warn import print_dec_g, pop_wrg_1, pop_err_1
+from lib.misc.display import show_in_green, show
+from lib.misc.warn import pop_wrg_1, pop_err_1
 from lib.misc.file import check_file, check_in_file
 
 from lib.handling.gromacs import extract_gromacs_data, get_ffldfiles
@@ -59,7 +59,7 @@ def create_conversion_gui(self):
         _atomstyle_ = 'full'
         
         data_cont = self.master._convert_['setup']
-        if data_cont <> []:
+        if data_cont != []:
             _autoload_ = data_cont[0]
             eg_files = data_cont[1:-1]
             _atomstyle_ = data_cont[-1]
@@ -125,7 +125,7 @@ def create_conversion_gui(self):
         ######     END CONSTRUCTION SITE  -> PACK OUTSIDE #########
         
         if self.master.test:
-            print 'Seeing main gro2lam converter GUI'
+            show( 'Seeing main gro2lam converter GUI')
             self.after(2000, self.test_hook  )
 
     def atomstyle( self):
@@ -186,7 +186,7 @@ def autoloadstuff(self):
         #_autoload_ = self.objt_c[0].get()
         main_top_file = self.objt_c[2].get()
 
-        if self.autoload_buffer <> main_top_file:
+        if self.autoload_buffer != main_top_file:
         #    self.autoload_buffer = _autoload_
             self.autoload_buffer = main_top_file
             aux_cont, nonerr_flag = get_ffldfiles( main_top_file )
@@ -219,9 +219,9 @@ def load_top_file( self, *args ):#event=None,
         in order to avoid the the waving experienced with a first "''"
         in the entry an if is inside 
         '''
-        #print args
-        if self.objt_c[2].get() <> '':
-            #print ( self.objt_c[2].get())
+        #show( args
+        if self.objt_c[2].get() != '':
+            #show( ( self.objt_c[2].get())
             
             self.autol_b.configure( state = 'normal')
             for i in [2, 3, 4]:
@@ -252,7 +252,7 @@ def getdata_and_convert(self):
             
             data_cont = self.master._convert_['setup']
             root_folder = '/'.join(data_cont[1].split('/')[:-1]+[''])
-            print( 'Root folder: {}'.format( root_folder) )
+            show( 'Root folder: {}'.format( root_folder) )
             
             sim_data, _flags_ = extract_gromacs_data( data_cont[1:-1],
                                                      _autoload_)
@@ -279,7 +279,7 @@ def getdata_and_convert(self):
                     flag_done_ = False
                     
             if flag_done_:
-                print_dec_g( 'Data file generated as "data.gro2lam"' )
+                show_in_green( 'Data file generated as "data.gro2lam"' )
                 
                 self._convertdata_ = sim_data
                 self._convertdata_['filename'] = data_fnam
diff --git a/lib/gui/custom_row.py b/lib/gui/custom_row.py
index 5b438c6..737fdaf 100644
--- a/lib/gui/custom_row.py
+++ b/lib/gui/custom_row.py
@@ -2,9 +2,9 @@
 #    By Hernan Chavez Thielemann
 __author__ = 'Hernan Chavez Thielemann <hchavezthiele at gmail dot com>'
 
-from Tkinter import Entry, Button, Frame, StringVar
-from tk_lib import format_dec
-from popup import FilePopUp
+from lib.gui.tk_lib import Entry, Button, Frame, StringVar, format_dec
+from lib.misc.display import show
+from lib.gui.popup import FilePopUp
 
 
 class File_Entry(Frame):
@@ -85,8 +85,8 @@ def browsefile(self, entry, ext=None):
         '''Browse a file <button> action binder'''
         
         pop = FilePopUp( master = self)
-        if ext<>None and isinstance(ext, tuple):
-            #print '--- ', ext, type(ext)
+        if ext != None and isinstance(ext, tuple):
+            #show( '--- ', ext, type(ext)
             pop.filetypes['filetypes'] = ext #(("All","*"),) # 
         filepath = pop.getfilepath()
         
@@ -98,8 +98,8 @@ def browsefile(self, entry, ext=None):
             entry.xview_moveto(1)
         except:
             if filepath not in [ '', ()]:
-                print "Could not open File: ", filepath
-                print exc_info()[1]
+                show( "Could not open File: ", filepath)
+                show( exc_info()[1])
     
     def setter(self, value):
         
diff --git a/lib/gui/img/README.md b/lib/gui/img/README_img.md
similarity index 100%
rename from lib/gui/img/README.md
rename to lib/gui/img/README_img.md
diff --git a/lib/gui/main_gui.py b/lib/gui/main_gui.py
index 1a30811..4b70ec2 100644
--- a/lib/gui/main_gui.py
+++ b/lib/gui/main_gui.py
@@ -9,15 +9,12 @@
 from os.path import dirname, realpath
 from sys import exit
 
-from Tkinter import Tk, Frame, Label, TclError, PhotoImage
-
-from conversion_gui import Conversion
-from script_gui import Script_GUI
-from run_gui import Run_GUI
-
-from popup import AboutPopUp
-from tk_lib import createmenubar
-
+from lib.gui.tk_lib import Tk, Frame, Label, TclError, PhotoImage, createmenubar
+from lib.gui.conversion_gui import Conversion
+from lib.gui.script_gui import Script_GUI
+from lib.gui.run_gui import Run_GUI
+from lib.gui.popup import AboutPopUp
+from lib.misc.display import show
 from lib.misc.warn import wrg_3
 from lib.misc.file import run_command
 from lib.misc.version import __version__
@@ -74,7 +71,7 @@ def swapbody(self, _pbody_):# checked ok 16/09 -----------WF
             this, because it is just a "small" overlapping I gess
         '''
         
-        if self.prevailing_body <> _pbody_:
+        if self.prevailing_body != _pbody_:
             if self.body == None:
                 self.body = self.create_conversion_gui()
                 
@@ -82,15 +79,15 @@ def swapbody(self, _pbody_):# checked ok 16/09 -----------WF
                 self.body.destroy()
                     
                 if _pbody_==1:
-                    print 'Swapping to gro2lam converter GUI'
+                    show( 'Swapping to gro2lam converter GUI')
                     self.body = self.create_conversion_gui()
                     
                 elif _pbody_==2:
-                    print 'Swapping to input script generator GUI'
+                    show( 'Swapping to input script generator GUI')
                     self.body = self.create_script_gui()
                     
                 elif _pbody_==3:
-                    print 'Swapping to run script GUI'
+                    show( 'Swapping to run script GUI')
                     self.body = self.create_run_gui()
                     
                 else:
@@ -136,9 +133,9 @@ def launch_gui( started = False):
         Main GUI constructor
     '''
     
-    print wrg_3('Before you start, make sure there are no comments',
+    show( wrg_3('Before you start, make sure there are no comments',
                 '(;) in the middle of a line of the input GROMACS files.',
-                'Data after this symbol are not taken into account.')
+                'Data after this symbol are not taken into account.') )
     
     MasterWin = Tk()
     prompt = Gro2Lam_GUI( master= MasterWin, test = started)# xl_App
@@ -165,15 +162,16 @@ def launch_gui( started = False):
     ws = MasterWin.winfo_screenwidth() # width of the screen
     hs = MasterWin.winfo_screenheight() # height of the screen
     # calculate x and y coordinates for the Tk root window
-    x = (ws/6) - (w/2)
-    if x <100:
+    x = int( (ws/6.0) - (w/2.0))
+    if x < 100:
         x = 100
-    y = (hs/3) - (h/2)
-    if y< 40:
+    y = int( (hs/3.0) - (h/2.0))
+    if y < 40:
         y = 40
-        
+    
+    #show([ws, hs, w, h, x, y] )
     prompt.MAINVERTEX = [ws, hs, w, h, x, y]
-    #print MAINVERTEX
+    show( prompt.MAINVERTEX, v = 4)
     # set the dimensions of the screen 
     # and where it is placed
     MasterWin.geometry('{:d}x{:d}+{:d}+{:d}'.format( *prompt.MAINVERTEX[2:]))
@@ -188,15 +186,15 @@ def launch_gui( started = False):
 
 def showlicence():
     
-    print 'Opening licence file'
-    command = 'gedit ./lib/docs/COPYING'#
-    run_command(command)
+    show( 'Opening licence file')
+    command = 'gedit ./lib/docs/LICENSE'#
+    run_command( command)
 
 def launch_about( _master_window_):
     
-    print 'Launching about'
+    show( 'Launching about')
     
-    title_txt = ' '*17+'ABOUT GROTOLAM'
+    title_txt = ' '*10+'ABOUT GROTOLAM'
     
     pop = AboutPopUp(master = _master_window_,
                      title = title_txt,
@@ -205,7 +203,7 @@ def launch_about( _master_window_):
 
 def showuserman():
                      
-    print 'Opening readme file'
+    show( 'Opening readme file')
     command = 'gedit ./lib/docs/README.md'#
     run_command(command)
 
diff --git a/lib/gui/popup.py b/lib/gui/popup.py
index 9417bb4..0dcc5cb 100644
--- a/lib/gui/popup.py
+++ b/lib/gui/popup.py
@@ -2,12 +2,13 @@
 #    By Hernan Chavez Thielemann
 __author__ = 'Hernan Chavez Thielemann <hchavezthiele at gmail dot com>'
 
-from Tkinter import Frame, Toplevel, Label, Button, Checkbutton, Entry
-from Tkinter import SUNKEN, X, Y, PhotoImage, IntVar
-from tk_lib import create_entry, bottom_hline_deco, get_entriesvalue
-from tkFont import Font
-from webbrowser import open_new
-from lib.misc.version import __version__
+from lib.gui.tk_lib import Frame, Toplevel, Label, Button, Checkbutton, Entry
+from lib.gui.tk_lib import SUNKEN, X, Y, PhotoImage, IntVar, Font, Tk
+from lib.gui.tk_lib import create_entry, bottom_hline_deco, get_entriesvalue
+from lib.gui.tk_lib import link_row, polito_callback, email_me
+
+from lib.misc.version import __version__, __copyright__
+from lib.misc.display import show
 
 class Message_box(Frame):
     ''' It is a frame because at the end it is a message box launchpad'''
@@ -71,7 +72,7 @@ def filetypes(self):
         return self._options_
     @filetypes.setter
     def filetypes(self, **value):
-        print '-----------', value, type (value)
+        show( '-----------', value, type (value) )
         value += (("All","*"),)
         self._options_ = {'filetypes': ("All","*")}
         
@@ -97,7 +98,7 @@ class SaveAsPopUp():
     def __init__(self, master=None):
         self.master  = master
         self._options = {}
-        print 'Soon ...  but save what?'
+        show( 'Soon ...  but save what?')
         
     @property
     def options(self):
@@ -161,7 +162,7 @@ def set_self_vertex(self):
         
         x_pop = x_main + w_main - 30
         y_pop = y_main + h_main - h_pop
-        #print w_pop, h_pop, x_pop, y_pop
+        #show( w_pop, h_pop, x_pop, y_pop
         self._vertex_= [w_pop, h_pop, x_pop, y_pop]
 
     def create_content(self):
@@ -256,7 +257,7 @@ def __init__(self, master = None, **kwargs):
         self.set_self_vertex()
         
         self.wm_title(' '*5+title)
-        #print self._vertex_
+        #show( self._vertex_
         self.geometry('{:d}x{:d}+{:d}+{:d}'.format(*self._vertex_))
         
         self.ent_c = []
@@ -271,7 +272,7 @@ def __init__(self, master = None, **kwargs):
     def set_self_vertex(self ):
         
         ws, hs, w_main, h_main, x_main, y_main = self.master.MAINVERTEX
-        #print ws, hs, w_main, h_main, x_main, y_main
+        #show( ws, hs, w_main, h_main, x_main, y_main
         # calculate the greatness
         
         
@@ -282,7 +283,7 @@ def set_self_vertex(self ):
         
         x_pop = x_main + w_main - 30
         y_pop = y_main + h_main - self._height_
-        #print w_pop, h_pop, x_pop, y_pop
+        #show( w_pop, h_pop, x_pop, y_pop
         self._vertex_= [self._width_, self._height_, x_pop, y_pop]
 
     def create_content(self):
@@ -295,7 +296,7 @@ def create_content(self):
         # lets create space to do some stuff ...
         self.entr_maxlen = int(len(max(self._entries_, key=len))*2.5/3.0)
         row2fill = Frame(self)
-        #print self._entries_,'\n', self._defvals_,
+        #show( self._entries_,'\n', self._defvals_,
         #'\n',len(self._entries_),'\n', len(self._defvals_)
         
         aline=0
@@ -303,7 +304,7 @@ def create_content(self):
             if self._entries_[e] == '---':
                 aline+=1
                 bottom_hline_deco(row2fill)
-                #print self._entries_[e]
+                #show( self._entries_[e]
             else:
                 self.ent_c.append(create_entry(row2fill,
                                                self._entries_[e],
@@ -393,7 +394,7 @@ def __init__(self, master = None, **kwargs):
         self.set_self_vertex()
         
         self.wm_title(' '*5+title)
-        #print self._vertex_
+        #show( self._vertex_
         self.geometry('{:d}x{:d}+{:d}+{:d}'.format(*self._vertex_))
         
         self.ent_c = []
@@ -413,11 +414,11 @@ def set_self_vertex(self ):
         if  self._width_ == None:
             self._width_ = w_main - 40
         if self._height_ == None:
-            self._height_ = h_main*12/21
+            self._height_ = int( h_main*12/21.0)
         
-        x_pop = x_main + w_main - 60
-        y_pop = y_main + h_main - self._height_
-        #print w_pop, h_pop, x_pop, y_pop
+        x_pop = int( x_main + w_main - 60)
+        y_pop = int( y_main + h_main - self._height_)
+        #show( w_pop, h_pop, x_pop, y_pop
         self._vertex_= [self._width_, self._height_, x_pop, y_pop]
 
     def create_content(self):
@@ -450,7 +451,7 @@ def create_content(self):
                 self.add_new_line()
                 self.ch_init+=[0]
             else:
-                #print self.extra_checkb
+                #show( self.extra_checkb
                 gr, ids, k_xyz, rns, cks =self.extra_checkb
                 
                 for _e_ in range( len( gr)):
@@ -574,7 +575,7 @@ def checkbuttonstuff(self):
         # First detect wich one change, comparing the fresh value
         for v in range(ckbc_len):
             _ckvar_g_.append(self.ckb_c[v][0].get())
-            if self.ch_init[v]<>_ckvar_g_[v]:
+            if self.ch_init[v] != _ckvar_g_[v]:
                 auxi=v
         
         if auxi>=loc and auxi == ckbc_len-1:
@@ -607,7 +608,7 @@ def save_it(self, event=None):
         for g_e in range(len(self.ent_c))[:lo]:
             gr, ids, k_xyz, rns = get_entriesvalue( self.ent_c[g_e])
             
-            #print gr, ids, k_xyz, self.ch_init[g_e]
+            #show( gr, ids, k_xyz, self.ch_init[g_e]
             
             g_names.append(gr)
             g_aids.append(ids)
@@ -672,12 +673,12 @@ def set_self_vertex(self):
         
         ws, hs, w_main, h_main, x_main, y_main = self.master.MAINVERTEX
         # calculate the greatness
-        w_pop = w_main
-        h_pop = h_main*14/21
+        w_pop = int( w_main*6/7.0)
+        h_pop = int( h_main*2/3.0)
         
         x_pop = x_main + w_main + 30
         y_pop = y_main + h_main - h_pop - 140
-        #print w_pop, h_pop, x_pop, y_pop
+        #show( w_pop, h_pop, x_pop, y_pop
         self._vertex_= [w_pop, h_pop, x_pop, y_pop]
 
     def create_content(self):
@@ -685,92 +686,63 @@ def create_content(self):
         self.pop = Toplevel(self.master, bg = "white")
         self.pop.grab_set()# when you show the popup
         
-        self.pop.wm_title(' '*5+self._title_)
-        self.pop.geometry('{:d}x{:d}+{:d}+{:d}'.format(*self._vertex_))
+        self.pop.wm_title(' '*5 + self._title_)
+        self.pop.geometry('{:d}x{:d}+{:d}+{:d}'.format( *self._vertex_))
         
+        ###########################      Left   column      ###################
         leftcolumn = Frame(self.pop ,bg = "white")
         Label(leftcolumn, bg = "white").pack(side="top", fill="both", pady=30)
         Label(leftcolumn, bg = "white",
               image= self.im_logo2).pack(side="top", fill="both", padx=10) #side= LEFT,
-        leftcolumn.pack(side="left", fill='both', padx=1)
-        
+        leftcolumn.pack( side="left", fill='both', padx=1)
         
-        rightcolumn = Frame(self.pop,bg = "white")
-        firstrow = Frame(rightcolumn,bg = "white")
-        Frame(rightcolumn,bg = "white").pack(side="top", pady=10 , padx=0)
+        ############################     Right  column       ##################
+        rightcolumn = Frame( self.pop, bg = "white")
+        # space above GRO2LAM letters -> pady 
+        Frame( rightcolumn, bg = "white").pack( side="top", pady = 5)
+        # GRO2LAM letters 
+        Label( rightcolumn, text = "GRO2LAM", fg = "Gray13", bg = "white" ,
+                 font = "Verdana 13 bold").pack( side = "top", anchor = 'w')
         
-        Label(firstrow, text="GROTOLAM",
-                 fg = "Gray13", bg = "white" ,
-                 font = "Verdana 13 bold").pack(side="left", pady=20)
-        Label(firstrow, text='',bg = "white").pack(side="left", padx=40)
-        firstrow.pack(side="top", padx=0)
+        # space after GRO2LAM letters -> pady 
+        Label( rightcolumn, bg = "white").pack( side = "top", pady=15)
         
-        secrow = Frame(rightcolumn, bg = "white")
-        Label(secrow, text="v "+__version__.split()[2],
-                 fg = "Gray13",bg = "white" ,
-                 font = "Verdana 10").pack(side="left")
-        Label(secrow, text='',bg = "white").pack(side="left", padx=75)
-        secrow.pack(side="top", padx=0)
+        # version_row
+        Label( rightcolumn, text = "v " + __version__.split()[2],
+                 fg = "Gray13", bg = "white", font = "Verdana 10"
+             ).pack( side = "top" , anchor = 'w', padx = 10)
         
         # lets create space to do some stuff ...
-        Frame(rightcolumn,bg = "white").pack(side="top", pady=20 , padx=0)
-        
-        thirdrow = Frame(rightcolumn,bg = "white")
+        Frame( rightcolumn,bg = "white").pack(side="top", pady=18)
         
-        Label(thirdrow, text="2018 Python version by",
-                 fg = "Gray13", bg = "white",
-                 font = "Verdana 10").pack(side="left")
-        Label(thirdrow, text='',bg = "white").pack(side="left", padx=16)
-        thirdrow.pack(side="top", padx=0)
-        
-        fourthrow = Frame(rightcolumn,bg = "white")
-        Label(fourthrow, text="Hernan Chavez Thielemann",
-                 fg = "Gray13", bg = "white",
-                 font = "Verdana 10").pack(side="left")
-        Label(fourthrow, text='',bg = "white").pack(side="left", padx=1)
-        fourthrow.pack(side="top", padx=0)
-        
-        fifthrow = Frame(rightcolumn,bg = "white")
-        Label(fifthrow, bg = "white",
-              image= self.im_smlogo).pack(side="left", fill="both", padx=10)
-        fifthrow.pack(side="top", padx=0)
-        
-        sixthrow = Frame(rightcolumn,bg = "white")
-        href = Label(sixthrow, bg = "white", font = "Verdana 10",
-                     text="Small web page", fg="blue",
-                     cursor="hand2")
-        f = Font(href, href.cget("font"))
-        f.configure(underline = True)
-        href.configure(font=f)
-        href.pack(side="left")
-        href.bind("<Button-1>", self.callback)
-        
-        Label(sixthrow, text='',bg = "white").pack(side="left", padx=40)
-        sixthrow.pack(side="top", padx=0)
-        
-        lastrow = Frame(rightcolumn,bg = "white")
-        self.bottom_button_row(lastrow)
-        rightcolumn.pack(side="right", fill='both', padx=5)
-
-    def callback(self, event):
-        open_new(r"http://www.polito.it/small")
+        ####   copyright
+        Label( rightcolumn, text = __copyright__ , fg = "Gray13", bg = "white",
+               font = "Verdana 10").pack( side = "top" , anchor = 'w')
         
+        mail_ = link_row( rightcolumn , "E-mail me!", email_me)
         
-    def bottom_button_row(self, _row_):
+        Label( rightcolumn, bg = "white", image= self.im_smlogo
+             ).pack( side = "top" , anchor = 'w')
         
-        b2 = Button(_row_, text='Close',bg = "white", command= self.exit_pop)
-        b2.pack(side="right", padx=10, pady=4)
-        b1 = Button(_row_, text='Licence',bg="white",command= self._licence_)
-        b1.pack(side="right", padx=10, pady=20)
-        _row_.pack(side="bottom")
+        href = link_row( rightcolumn , "Small-lab web page", polito_callback)
         
-    def openlicence(self):
+        _row_ = Frame(rightcolumn,bg = "white")
+        b2 = Button(_row_, text='Close', bg = "white", command= self.exit_pop)
+        b2.pack( side = "right", padx=10, pady=4)
+        b1 = Button(_row_, text='Licence', bg="white", command= self._licence_)
+        b1.pack( side = "right", padx=10, pady=20)
+        _row_.pack( side = "bottom")
         
-        print 'opening licence file'
+        rightcolumn.pack( side = "left", padx=10)
 
-    def exit_pop(self):
+    def exit_pop( self):
         self.pop.grab_release() # to return to normal
         self.pop.destroy()
+        
+
+def openlicence_hook():
+        show( 'opening licence file')
+    
 
 class WarningPopUp2(PromptPopUp):
 
@@ -802,5 +774,36 @@ def bottom_button_row(self):
         b2.pack(side="right", padx=10, pady=4)
         _br_.pack(side="bottom", expand = True)
         
+        
+def test_launch_about():
+    
+    MasterWin = Tk()
+    w = 460
+    h = 570
+    ws = MasterWin.winfo_screenwidth() # width of the screen
+    hs = MasterWin.winfo_screenheight() # height of the screen
+    x = int( (ws/6) - (w/2))
+    if x <100:
+        x = 100
+    y = int( (hs/3) - (h/2))
+    if y< 40:
+        y = 40
+        
+    MasterWin.MAINVERTEX = [ws, hs, w, h, x, y]
+    MasterWin.geometry('{:d}x{:d}+{:d}+{:d}'.format( *MasterWin.MAINVERTEX[2:]))
+    
+    show( 'Launching about')
+    title_txt = ' '*10+'ABOUT GROTOLAM'
+    pop = AboutPopUp(master = MasterWin,
+                     title = title_txt,
+                     licence = openlicence_hook
+                    )
+    
+    MasterWin.mainloop()
+    MasterWin.destroy()
+    
+if __name__ == "__main__":
+    test_launch_about()
+    
 # vim:tw=80
 
diff --git a/lib/gui/run_gui.py b/lib/gui/run_gui.py
index 3b852cd..f332529 100644
--- a/lib/gui/run_gui.py
+++ b/lib/gui/run_gui.py
@@ -7,12 +7,13 @@
 #///    Packages and globals definitions are here   ///
 #------------------------------------------------------
 
-from Tkinter import Frame, Button, Label
+from lib.gui.tk_lib import Frame, Button, Label
+from lib.gui.tk_lib import create_entry, get_entriesvalue, bottom_hline_deco
+from lib.gui.custom_row import createfileentry
 
-from tk_lib import create_entry, get_entriesvalue, bottom_hline_deco
-from custom_row import createfileentry
 from lib.misc.data import check_vars
 from lib.misc.file import check_file, run_command
+from lib.misc.display import show
 
 #------------------------------------------------------
 '''///////////////        Class       /////////////'''
@@ -83,7 +84,7 @@ def createWidgets(self):
         
         if self.master.test:
             self.after(2000, self.master.master.destroy )
-            print '\n'+'-'*20 + '  Test ended  ' + '-'*20
+            show( '\n'+'-'*20 + '  Test ended  ' + '-'*20)
     
     def build_finalbuttons(self):
         '''    Final Buttons    '''
@@ -108,7 +109,7 @@ def run_script( _file_, machine='lammps', _cores_='1'):
     
     core_flag = min(check_vars([_cores_],['int'],'Run aborted!'))
     if check_file(_file_) and core_flag:
-        print _file_
+        show( _file_)
         command = '{} -echo both -in {}'.format( machine, _file_)
         if int(_cores_)>1:
             command = 'mpirun -np {} '.format( _cores_) + command
diff --git a/lib/gui/script_gui.py b/lib/gui/script_gui.py
index a5747d1..526f778 100644
--- a/lib/gui/script_gui.py
+++ b/lib/gui/script_gui.py
@@ -7,19 +7,17 @@
 #///    Packages and globals definitions are here   ///
 #------------------------------------------------------
 
-from Tkinter import Frame, Button, Label
-
-from popup import PromptPopUp, PromptPopUp_wck
-
-from tk_lib import bottom_hline_deco, format_dec, Drop_Down_List
-from tk_lib import createmenubar, create_entry, get_entriesvalue
-
-from custom_row import File_Entry, createfileentry
+from lib.gui.tk_lib import Frame, Button, Label
+from lib.gui.tk_lib import bottom_hline_deco, format_dec, Drop_Down_List
+from lib.gui.tk_lib import createmenubar, create_entry, get_entriesvalue
+from lib.gui.popup import PromptPopUp, PromptPopUp_wck
+from lib.gui.custom_row import File_Entry, createfileentry
 
 from lib.handling.lammps import write_lammps_input
 
 from lib.misc.data import check_vars
-from lib.misc.warn import print_dec_g, pop_wrg_1
+from lib.misc.display import show, show_in_green
+from lib.misc.warn import pop_wrg_1
 
 #------------------------------------------------------
 '''///////////////        Class       /////////////'''
@@ -78,9 +76,9 @@ def create_main_section(self):
                      '1000','10:1', '1000', '299:300', '100']
         _def_dataname_ = './data.gro2lam'
         
-        if self._convertdata_ <> None:
+        if self._convertdata_ != None:
             _def_dataname_ = self._convertdata_['filename']
-        elif _mainpage_ <> []:
+        elif _mainpage_ != []:
             _def_dataname_ = _mainpage_[0]
             _defvals_ = _mainpage_[1:] 
         
@@ -129,7 +127,7 @@ def build_finalbuttons(self):
             _comb_rl_ = 'No'
             lj_12_13_14_ = '0.0:0.0:0.0'
             co_12_13_14_ = '0.0:0.0:0.0'
-            if self._convertdata_ <> None:
+            if self._convertdata_ != None:
                 _aux_here_ = self._convertdata_['defaults']
                 buckorlj, comb_rule, _flg_, f14_LJ, f14_Co = _aux_here_
                 if int( buckorlj) == 2:
@@ -186,7 +184,7 @@ def build_finalbuttons(self):
         
     def write_script(self):
     
-        print 'Writing the lammps script'
+        show( 'Writing the lammps script')
         _mainpage_ = get_entriesvalue( self.s_entry_c)
         
         _flag_ = True # data ok flag
@@ -205,7 +203,7 @@ def write_script(self):
             _flag_ = write_lammps_input( _script_setup_, self._convertdata_)
             
             if _flag_:
-                print_dec_g( 'Lammps script done!')
+                show_in_green( 'Lammps script done!')
                 self.master.swapbody(3)
     
     def further_config_script( self ):
@@ -251,7 +249,7 @@ def further_config_script( self ):
                         #, 'none',
                         'zero']
         
-        if self._convertdata_ <> None:
+        if self._convertdata_ != None:
             buckorlj = int( self._convertdata_['defaults'][0])
             if buckorlj == 1:
                 _pair_style_ = _pair_style_[ :3] + [_pair_style_[ -1]]
@@ -276,7 +274,7 @@ def further_config_script( self ):
         _defvals_ = []
         
         for _ad_ in range(len(self._container_['advanced'][0])):
-            if self._container_['advanced'][1][_ad_] <> '':
+            if self._container_['advanced'][1][_ad_] != '':
                 _def_ = self._container_['advanced'][0][_ad_]
                 _dfli_ = self._container_['advanced'][1][_ad_]
                 _defvals_.append([ _def_, _dfli_])
@@ -297,8 +295,8 @@ def further_config_script( self ):
         pop.wait_window()
         
         _, bnty_len, anty_len = self.check_datafile()
-        print bnty_len, anty_len
-        if self.master._aux_ <> []:
+        show( bnty_len, anty_len)
+        if self.master._aux_ != []:
             _advanced_ = self.master._aux_
             
             _entries_ = [ 'int', '', '', 'float', 'float', 'float', '',
@@ -314,7 +312,7 @@ def further_config_script( self ):
                                      'Advanced settings not saved!'))
             if _flag_:
                 self._container_['advanced'][0] = _advanced_
-                print_dec_g('Advanced settings saved')
+                show_in_green('Advanced settings saved')
         self.fcb.config(bg = 'gray86')#, width = 145)
     
     def config_restrain( self ):
@@ -337,7 +335,7 @@ def config_restrain( self ):
                 ck_init = [0]
                 second_c = None
                 
-                if self._convertdata_<> None:
+                if self._convertdata_!= None:
                     _mol_niifi_ = self._convertdata_['atomsdata'][1]
                     for mt in range(len( _mol_niifi_)):
                         g_names.append(_mol_niifi_[mt][0])
@@ -348,8 +346,8 @@ def config_restrain( self ):
                 
             else: 
                 g_names, d_ids, kxyz_c, rest_ens, ck_init = _defvals_[0]
-                #print g_names,'\n', d_ids,'\n', kxyz_c,'\n', ck_init 
-                if _defvals_[1]<> None:
+                #show( g_names,'\n', d_ids,'\n', kxyz_c,'\n', ck_init 
+                if _defvals_[1]!= None:
                     second_c = _defvals_[1]
                 else:
                     second_c = None
@@ -372,16 +370,16 @@ def config_restrain( self ):
             
             pop.wait_window()
             
-            if self.master._aux_ <> []:
+            if self.master._aux_ != []:
                 sim_len = 0
                 for x in  self._container_['advanced'][0][-4].split('-'):
-                    if x.strip(' ') <> 'R':
+                    if x.strip(' ') != 'R':
                         sim_len +=1
                 _res_flag_= [[],[]]
                 _restrain_ = self.master._aux_[:]
                 _, _d_ids, _kxyz_c, _runs_c, _res_flag_[0] = _restrain_[0]
                 
-                if _restrain_[1] <> None:
+                if _restrain_[1] != None:
                     _, au2, au3, au4, au5 = _restrain_[1][:] 
                     _restrain_aux = _d_ids+ au2+ _kxyz_c+ au3+ _runs_c+ au4
                     _res_flag_[1] = au5
@@ -395,15 +393,15 @@ def config_restrain( self ):
                 
                 _aux_ = check_vars( _restrain_aux, _entries_,
                                          'Restrain groups not saved!')
-                #print _aux_
+                #show( _aux_
                 _flag_ = min( _aux_)
                 if _flag_:
                     
                     if max(max(_res_flag_)):
                         self._container_['restrain'] = _restrain_
-                        print_dec_g('Restrain data saved')
+                        show_in_green('Restrain data saved')
                     else:
-                        print ('Creating 0 groups, Restraining 0 atoms')
+                        show( 'Creating 0 groups, Restraining 0 atoms')
             self.resb.config(bg = 'gray86')#, width = 45)
 
     def check_datafile(self, _bflag_=None):
@@ -413,7 +411,7 @@ def check_datafile(self, _bflag_=None):
         '''
         max_at_index, bond_types, angle_types = 0, 0, 0
         _flag_ = True
-        if self._convertdata_<> None:
+        if self._convertdata_!= None:
             _numbers_ = self._convertdata_['numbers']
             max_at_index = _numbers_['total'][0]
             _, bond_types, angle_types, _, _ = _numbers_['type']
@@ -434,10 +432,10 @@ def check_datafile(self, _bflag_=None):
                                     break
             except IOError:
                 pop_wrg_1('Data file not found!')
-                print ('Try performing a conversion first!')
+                show( 'Try performing a conversion first!')
                 _flag_ = False
                 
-        if _bflag_<>None:
+        if _bflag_ != None:
             return _flag_
         return max_at_index, bond_types, angle_types
         
diff --git a/lib/gui/tk_lib.py b/lib/gui/tk_lib.py
index 9107eb9..d33422c 100644
--- a/lib/gui/tk_lib.py
+++ b/lib/gui/tk_lib.py
@@ -2,12 +2,35 @@
 #    By Hernan Chavez Thielemann
 __author__ = 'Hernan Chavez Thielemann <hchavezthiele at gmail dot com>'
 
-from Tkinter import Entry, Button, Frame, Label, Scrollbar, StringVar, Menu
-from Tkinter import Listbox, IntVar, Checkbutton, Widget, Spinbox
-from Tkinter import X, Y, SUNKEN, VERTICAL, END
+from lib.misc.version import __python_version__
+from lib.misc.display import show
+
+from webbrowser import open_new
+
+if __python_version__ > 3:
+    '''between trying and knowing, I rather prefer the last one'''
+    from tkinter import Tk, Frame, Label, TclError, PhotoImage, Toplevel
+    from tkinter import Entry, Button, Scrollbar, StringVar, Menu
+    from tkinter import Listbox, IntVar, Checkbutton, Widget, Spinbox
+    from tkinter.font import Font
+else:
+    from Tkinter import Tk, Frame, Label, TclError, PhotoImage, Toplevel
+    from Tkinter import Entry, Button, Scrollbar, StringVar, Menu
+    from Tkinter import Listbox, IntVar, Checkbutton, Widget, Spinbox
+    from tkFont import Font
 
 '''     Home made ttk library     '''
 
+#================================
+#        Variables definition
+#================================
+X = 'x'
+Y = 'y'
+SUNKEN = 'sunken'
+VERTICAL = 'vertical'
+END = 'end'
+'''     Home made ttk library     '''
+
 #================================
 #        Classes definition
 #================================
@@ -40,7 +63,7 @@ def create_entry( _main_row_, _desc_txt, _def_txt_, t_len, ckbc=None, fn=None):
     ''' creates a tkinter entry '''
     _row_ = Frame(_main_row_)
     
-    if ckbc<>None and fn<>None:
+    if ckbc != None and fn != None:
         label0 = Label(_row_, width=3, text="  ", anchor='w')
         cvar = IntVar()
         label1 = Checkbutton(_row_, width=0
@@ -116,7 +139,7 @@ def get_entriesvalue(entries_container):
 
 def bottom_hline_deco(row, func=None):
     ''' adds a sunken line in the bottom of a tkinter function'''
-    if func <> None:
+    if func != None:
         func()
     line = Frame( row ,height=2, bd=1, relief= SUNKEN)
     line.pack(fill=X, padx=1, pady=5)
@@ -261,9 +284,39 @@ def createmenubar( _root_window_, _listofentriesdicts_):
                 else:
                     sub_menu.add_separator(  )
         else:
-             print ' - incomplete data - '
+             show( ' - incomplete data - ')
+    
+
+def link_row( in_frame_, w_txt, link_fo, px=0, py=0):
+    
+    label_obj = Label( in_frame_, bg = "white", font = "Verdana 10",
+                 text = w_txt, fg="blue", cursor = "hand2")
+    font_here = Font( label_obj, label_obj.cget( "font"))
+    font_here.configure( underline = True)
+    label_obj.configure( font = font_here)
+    label_obj.pack( side = "top" , anchor = 'w', padx = px, pady = py)
+    label_obj.bind( "<Button-1>", link_fo)
+    
+    return label_obj
+
+def polito_callback( push_event):
+    open_new(r"http://www.polito.it/small")
+    
+
+def email_me( push_event):
+    ''' You can find this link in the About pop-up inn the help button'''
+    show( 'Writing an email to Hernan')
+    email_address = "herch{0:}vezt{1:}gm{0:}il{2:}com".format('a', '@', '.')
+    subject = "GRO2LAM email"
+    open2send_email( email_address, subject)
+    
+
+def open2send_email( email, subject):
+    open_new("mailto:?to=" + email  + "&subject=" + subject )
+        
+    
 
 def testprint():
-    print '>command33<'
+    show( '>command33<')
 
 # vim:tw=80
diff --git a/lib/handling/gromacs.py b/lib/handling/gromacs.py
index d854ee1..373da6f 100644
--- a/lib/handling/gromacs.py
+++ b/lib/handling/gromacs.py
@@ -4,6 +4,7 @@
 
 
 from lib.misc.warn import wrg_1, wrg_3, pop_err_1, pop_wrg_1
+from lib.misc.display import show
 from lib.misc.file import check_file, debugger_file, fileseeker
 from lib.misc.geometry import rotate, arcos, raiz
 from lib.misc.data import isnot_num
@@ -24,10 +25,10 @@ def extract_gromacs_data( _data_files_, _autoload_):
     data_container  =   {}
     data_container['define'] = {}
     #_solvated_f_, _autoload_= _ck_buttons_
-    print 'Autoload: {}\n'.format( _autoload_)
+    show( 'Autoload: {}\n'.format( _autoload_))
     
     if not _autoload_:
-        print filename_ff # or not
+        show( filename_ff) # or not
         
     _sidemol_f_ = False
     
@@ -69,7 +70,7 @@ def extract_gromacs_data( _data_files_, _autoload_):
                     
         cero = 1e-12
         if Ar[1][0] < cero or Ar[2][0] < cero or Ar[2][1] < cero:
-            print('Your triclinic cell will be rotated to make it!')
+            show( 'Your triclinic cell will be rotated to make it!')
             # y rotation
             a_tor_y = -arcos( (Ar[0][0])/(raiz(Ar[0][0]*Ar[0][0]+Ar[2][0]*Ar[2][0])) )
             Ar = rotate( Ar, a_tor_y, 'y')
@@ -129,7 +130,7 @@ def extract_gromacs_data( _data_files_, _autoload_):
     
     if pure_side_mol_flag:
         startstrings = startstrings[-3:]
-        print wrg_3( 'Using pure side molecule scheme')
+        show( wrg_3( 'Using pure side molecule scheme'))
         #n_atoms     =   0
         #n_bonds     =   0
         #n_angles    =   0
@@ -149,8 +150,8 @@ def extract_gromacs_data( _data_files_, _autoload_):
         if not ok_flag:
             
             if startstrings[ti] not in exclusions_:
-                print wrg_3( 'Not ok flag in <extract_gromacs_data> top file' +
-                         'section, in ' + s_str_)
+                show( wrg_3( 'Not ok flag in <extract_gromacs_data> top file' +
+                         'section, in ' + s_str_) )
                 return {}, [ ok_flag, _sidemol_f_]
             else:
                 ok_flag = True
@@ -244,17 +245,17 @@ def extract_gromacs_data( _data_files_, _autoload_):
     gromosff_flag = False
     data_container[ 'define'][ 'improper'] = {}
     aux_here = {}
-    print( section.split('.')[1])
+    show( section.split('.')[1])
     if filename_nb != filename_ff and filename_nb != filename_bon:
-        print(" Is it GROMOS there ?? ")
+        show( " Is it GROMOS there ?? ")
         aux_here = get_gromos_define( filename_bon)
         
     else:
-        print('no gromos check')
+        show( 'no gromos check')
         
     for key_ in aux_here.keys():
         if aux_here[ key_] != {}:
-            print ( 'GROMOS ' + key_ + ' kind detected!')
+            show( 'GROMOS ' + key_ + ' kind detected!')
             data_container[ 'define'][ key_].update( aux_here[ key_])
             gromosff_flag = True
             
@@ -262,14 +263,14 @@ def extract_gromacs_data( _data_files_, _autoload_):
             if 'dihedraltypes' == key_+'types' and dihe_g_data != []:
                 rewrite_flag = False
                 for gd_ in range( len( dihe_g_data)):
-                    #print dihe_g_data[gd_][2]
+                    #show( dihe_g_data[gd_][2]
                     if dihe_g_data[gd_][2].isdigit():
                         if not rewrite_flag:
-                            print('Dihedral with 2 atoms re-formating to 4: ')
+                            show( 'Dihedral with 2 atoms re-formating to 4: ')
                         rewrite_flag = True
                         dihe_g_data[gd_] = ( [ 'X',] + dihe_g_data[ gd_][:2]
                                             + [ 'X',] + dihe_g_data[ gd_][2:])
-                        print (dihe_g_data[ gd_]) 
+                        show( dihe_g_data[ gd_]) 
                 if rewrite_flag:
                     data_container[ 'dihedraltypes'] = dihe_g_data
                 
@@ -285,12 +286,12 @@ def extract_gromacs_data( _data_files_, _autoload_):
                 if len( data_aux[i][-1].split('.')) < 2:
                     if  not data_aux[i][-1].isdigit():
                         aux = data_aux[i][:-1] + cddd[ data_aux[i][-1]]
-                        #print aux
+                        #show( aux
                         data_container[ s_str_][i] = aux
                     
     
     
-        #print data_container['define']['bond']['gb_33']
+        #show( data_container['define']['bond']['gb_33']
     # Search for impropers in TOP and BON, using crossreference if needed
     data_container = split_define_dihe_impr( data_container)
     
@@ -319,7 +320,7 @@ def extract_gromacs_data( _data_files_, _autoload_):
             improp_x_mol = len( sidemol['data'][sb]['impropers'])
             
             sm_quantity = sidemol['num'][sb]
-            #print(sm_quantity, bonds_x_mol,sm_quantity * bonds_x_mol)
+            #show((sm_quantity, bonds_x_mol,sm_quantity * bonds_x_mol)
             side_bonds_n    += sm_quantity * bonds_x_mol
             side_angles_n   += sm_quantity * angles_x_mol
             side_dihed_n    += sm_quantity * dihedr_x_mol
@@ -329,8 +330,8 @@ def extract_gromacs_data( _data_files_, _autoload_):
         n_anglesnew =   n_angles + side_angles_n
         n_dihednew  =   n_dihedrals + side_dihed_n
         n_impropnew =   n_impropers + side_improp_n
-        #print n_bonds, side_bonds_n, n_bonds + side_bonds_n
-        #print n_angles, side_angles_n, n_angles + side_angles_n
+        #show( n_bonds, side_bonds_n, n_bonds + side_bonds_n
+        #show( n_angles, side_angles_n, n_angles + side_angles_n
         
         ###   A_03
         # tester in case is an asigment for define ore something like that
@@ -344,12 +345,12 @@ def extract_gromacs_data( _data_files_, _autoload_):
                     datacont = sidemol['data'][sb][cont_k+'s']# in its cont-key
                     for dc in range( len( datacont)):# lets look their content
                         if isnot_num( datacont[dc][-1]):#
-                            #print( '{} {} {}'.format( cont_k+'s', dc,
+                            #show(( '{} {} {}'.format( cont_k+'s', dc,
                             #                         datacont[dc][-1])) 
                             aux = datacont[dc][:-1] + cddd[ datacont[dc][-1]]
                             sidemol['data'][sb][cont_k+'s'][dc] = aux
                         #else:
-                            #print datacont[dc]
+                            #show( datacont[dc]
         
         #######################################################################
         ###   A_04
@@ -369,9 +370,9 @@ def extract_gromacs_data( _data_files_, _autoload_):
                 _charge_[a_opls_tag] = a_charge
                 _conv_dict_[ a_elem_tag] = a_opls_tag
                 
-        print '='*45+'\n'+'='*5+'  Charges found: '
-        print _charge_
-        print _conv_dict_
+        show( '='*45+'\n'+'='*5+'  Charges found: ')
+        show( repr(_charge_) )
+        show( repr(_conv_dict_) )
         
         data_container['S_charge'] = _charge_
         data_container['S_translation'] = _conv_dict_
@@ -437,12 +438,12 @@ def extract_gromacs_data( _data_files_, _autoload_):
                     smol_extra_impropertypes.append( aux_here + _smim_[4:])
                 
             if len( _smd_.keys()) > 5:
-                print ('Uuupa!! This thing is not implemented yet' +
-                       ' as side mol part')
+                show( ('Uuupa!! This thing is not implemented yet' +
+                       ' as side mol part'))
                 a_key = [ 'atoms', 'bonds', 'angles', 'dihedrals', 'impropers']
                 for ky in _smd_.keys():
                     if ky not in a_key:
-                        print ('-- > this key : ' + ky)
+                        show( '-- > this key : ' + ky)
                 
                 
             # ---------   !!!!    Update the info    !!!!
@@ -454,7 +455,7 @@ def extract_gromacs_data( _data_files_, _autoload_):
                                         data_container['dihedraltypes'])
         data_container['impropertypes'] = ( smol_extra_impropertypes + 
                                         data_container['impropertypes'])
-        #print(data_container['bondtypes'])
+        #show((data_container['bondtypes'])
     else:
         n_bondsnew  = n_bonds
         n_anglesnew = n_angles
@@ -473,8 +474,8 @@ def extract_gromacs_data( _data_files_, _autoload_):
             poss = 4
         for i in range( len ( data_container[ nice_list[it] ])):
             _aux_set_here.add( data_container[ nice_list[it] ][i][ poss ])
-        #print( nice_list[it][:4])
-        #print(_aux_set_here)
+        #show(( nice_list[it][:4])
+        #show((_aux_set_here)
         data_container[ nice_list[it][:4]+'_kinds'] = _aux_set_here
     
     
@@ -487,12 +488,12 @@ def extract_gromacs_data( _data_files_, _autoload_):
     data_container['numbers']['total'] = [n_atomsnew, n_bondsnew,
                                           n_anglesnew, n_dihednew, n_impropnew
                                          ]
-    #print( 'here', data_container['numbers']['total'])
+    #show(( 'here', data_container['numbers']['total'])
     #exit('')
     data_container['numbers']['type'] = [n_atomtypes, n_bondstypes,
                                          n_anglestypes, n_dihedraltypes,
                                          n_impropertypes]
-    print 'Ending gromacs data parsing\n'
+    show( 'Ending gromacs data parsing\n')
     return data_container, [ ok_flag, _sidemol_f_]
 
 def sidemol_data( _file_top_, data_container):
@@ -529,7 +530,7 @@ def sidemol_data( _file_top_, data_container):
     
     if sm_flag:
         #////////======= Side molecule >>> file <<< search   ==========////////
-        print ('\nLoading side molecule files: ' )
+        show( '\nLoading side molecule files: ' )
         _sm_files_ = []
         root_dir = '/'.join( _file_top_.split('/')[:-1]+[''])
         ok_flag = False
@@ -549,13 +550,13 @@ def sidemol_data( _file_top_, data_container):
                             new_filename = k_line.split('"')[1].lstrip('.')
                         except IndexError:
                             auxt = wrg_1( 'Format error with {}')
-                            print( auxt.format( k_line.split()[-1] ) ) 
+                            show( auxt.format( k_line.split()[-1] ) ) 
                         
                         new_filename = new_filename.lstrip('/').split('/')[-1]
                         po_file = fileseeker( root_dir, new_filename)
                         if po_file != []:
                             _sm_files_.append( po_file[0])
-                            print( 'SM_file : {}'.format(_sm_files_[-1]))
+                            show( 'SM_file : {}'.format(_sm_files_[-1]))
                             ok_flag *= check_file( po_file[0],
                                                   content='[ atoms ]')
                     else:
@@ -566,7 +567,7 @@ def sidemol_data( _file_top_, data_container):
         if ok_flag:
             for sm in sidemol['tag']:
                 aux_data, aux_flag = sidemol_data_gatherer( _sm_files_, sm)
-                #print aux_data
+                #show( aux_data
                 ok_flag *= aux_flag
                 sidemol['data'].append( aux_data)
                 
@@ -574,7 +575,7 @@ def sidemol_data( _file_top_, data_container):
         data_container['sidemol'] = sidemol
         
     else:
-        print ('No side molecule files detected!' )
+        show( 'No side molecule files detected!' )
         ok_flag = True
         
     return data_container, ok_flag, sm_flag
@@ -583,7 +584,7 @@ def sidemol_data_gatherer( _sm_files_, _sm_):
     ''' collects all the data related with one kind of side molecule
         the data types are specified in startstrings
     '''
-    print( '\nSearching for: {}'.format( _sm_ ))#, ' in: ' ,_sm_files_
+    show( '\nSearching for: {}'.format( _sm_ ))#, ' in: ' ,_sm_files_
     _flag_ = True
     _file_ = ''
     _sm_data_c_ = {}
@@ -596,7 +597,7 @@ def sidemol_data_gatherer( _sm_files_, _sm_):
             for j_line in sm_data:
                 j_line = j_line.split(';')[0].strip()
                 
-                #print j_line, read_flag, j_line.startswith(sm)
+                #show( j_line, read_flag, j_line.startswith(sm)
                 
                 if j_line.startswith('['):
                     if j_line.startswith('[ moleculetype ]'):
@@ -618,7 +619,7 @@ def sidemol_data_gatherer( _sm_files_, _sm_):
         pop_err_1('Error!! side molecule {} not found in itp -- '.format( _sm_))
         _flag_ = False
     else:
-        print( 'Success!, found in : {}\n'.format( _file_))
+        show( 'Success!, found in : {}\n'.format( _file_))
         tag_str = [ 'atoms', 'bonds', 'angles', 'dihedrals','fin']
         _sm_data_c_ = { x:[] for x in tag_str if x != 'fin'}
         read_flag = False
@@ -640,20 +641,20 @@ def sidemol_data_gatherer( _sm_files_, _sm_):
                                 i = tag_str.index( j_line[1])
                                 cd_tag = tag_str[i]
                                 iner_flag = True
-                                print( '** Gathering {} data'.format( cd_tag))
+                                show( '** Gathering {} data'.format( cd_tag))
                             elif j_line[1] == 'moleculetype':
                                 break
                             else:
                                 txt_s = '> {} not considered in {}'
-                                print txt_s.format( j_line[1], _sm_)
+                                show( txt_s.format( j_line[1], _sm_) )
                                 iner_flag = False
                         else :
                             cd_tag = tag_str[i]
-                            print( '* Gathering {} data'.format( cd_tag))
+                            show( '* Gathering {} data'.format( cd_tag))
                             iner_flag = True
                             
                     elif iner_flag:
-                        #print j_line
+                        #show( j_line
                         _sm_data_c_[ cd_tag].append( j_line)
                             
                             
@@ -661,13 +662,13 @@ def sidemol_data_gatherer( _sm_files_, _sm_):
                     read_flag = True
                     
         # todo add a new check in case of empty container
-        #print _sm_data_c_
+        #show( _sm_data_c_
         
         ######### Split impropers and dihedrals
         _sm_data_c_, _ = split_dihedral_improper( _sm_data_c_)
         
-        #if _sm_data_c_['impropers'] <>[]:
-        #    print _sm_data_c_['impropers']
+        #if _sm_data_c_['impropers'] != []:
+        #    show( _sm_data_c_['impropers']
         
     return _sm_data_c_, _flag_
 
@@ -713,7 +714,7 @@ def split_dihedral_improper( _data_container_):
                 def_impr_extra.append( _dihedrals_data_[i])
             
         else:
-            print 'Problem #008 here #split_dihedral_improper'
+            show( 'Problem #008 here #split_dihedral_improper')
             dihe_err = "#008_" + dh_dict_kind[ dihe_funct]
             
             if dihe_err != dh_bf_err:
@@ -749,7 +750,7 @@ def split_define_dihe_impr( _data_container_, smt_flag = False):
     _dihe_type_data_ = _data_container_['dihedraltypes']
     im_type_ = []
     dh_type_ = []
-    #print _dihe_type_data_
+    #show( _dihe_type_data_
     for i in range( len ( _dihe_type_data_)):
         #  Dihedral line format:
         #  ai  aj  ak  al  funct   c0  c1  c2  c3  c4  c5
@@ -759,7 +760,7 @@ def split_define_dihe_impr( _data_container_, smt_flag = False):
             else:
                 dh_type_.append( _dihe_type_data_[i])
         except IndexError as _err:
-            print( _err)
+            show( _err)
             exit( _dihe_type_data_[i] )
     #====================================================
     ''' ========  Dihedral "#define" data  ========== ''' 
@@ -798,7 +799,7 @@ def split_define_dihe_impr( _data_container_, smt_flag = False):
             #### TODO Flag
             ## First case with coefs in the top file... c0  c1  c2  c3  c4  c5
             if len( _dhi_) > 6:
-                print ('Coefficients in the top file are not supported yet' +
+                show( 'Coefficients in the top file are not supported yet' +
                         '... or maybe they are '+ u'\u00AC'*2)
                 new_dihedraltypes['-'.join(a_tag)] = (a_tag + _dhi_[4:])
             
@@ -863,11 +864,11 @@ def get_gro_fixed_line( _filename_):
                 if at == at_num:
                     read_flag = False
                     g_names.append(_buffer + [at])
-                    
+                
             elif j_line.startswith(';'):
                 pass
             elif at_num == 0:
-                j_line = indata.next()
+                j_line = next( indata)
                 at_num = int( j_line)
                 read_flag = True
             elif at == at_num:
@@ -904,15 +905,14 @@ def get_topitp_line( _filename_, _ss_):
         read_flag = False
         ok_flag = True
         tag_not_found = True
-        if _verbose_:
-            print  _ss_
+        show(  _ss_, v = 3)
         for j_line in indata:
             # I just whant to read once the flag is on
             j_line_s0 = j_line.split(';')[0].split()
             if read_flag and j_line_s0:
-                #if _verbose_: ### is beter to store and print outside the
+                #if _verbose_: ### is beter to store and show( outside the
                 # cycle with just one if 
-                #print j_line_s0
+                #show( j_line_s0
                 _line_ = j_line_s0
                 # getting out comments and empty lines
                 if len( _line_) <0: 
@@ -920,25 +920,25 @@ def get_topitp_line( _filename_, _ss_):
                     
                 elif _line_[0][0] == '#':
                     if _line_[0] == '#include':
-                        print( wrg_3( _line_[1] + ' skipped this time'))
+                        show( wrg_3( _line_[1] + ' skipped this time'))
                     elif _line_[0] == '#define':
                         _define_[_line_[1]] = _line_[2:]
                     else:
-                        print wrg_3( str(_line_) + '  ??')
+                        show( wrg_3( str(_line_) + '  ??') )
                         
                 elif _line_[0][0] == '[':
-                    print( ' '.join(_line_) + 'Checked!')
+                    show( ' '.join(_line_) + 'Checked!')
                     if  ' '.join(_line_) != _ss_ :
                         read_flag = False
-                    #print 'exit here 424'
+                    #show( 'exit here 424'
                     
                 elif len( _line_) > 0:
                     content_line.append( _line_)
                 else:
-                    print('Ups... please raise an issue at GitHub ;)')
+                    show( 'Ups... please raise an issue at GitHub ;)')
             elif j_line.lstrip().startswith( _ss_):
                 if _verbose_:
-                    print( _ss_+' found!')
+                    show( _ss_ + ' found!')
                 read_flag = True
                 tag_not_found = False
     
@@ -974,12 +974,12 @@ def get_gromos_define( _bondedfile_):
                             aux_dic = { _line_[1] : _line_[2:]}
                             _define_[ _dedic_[ _line_[1][1]]].update( aux_dic)
                         else:
-                            print('Whojojoooh...')
+                            show( 'Whojojoooh...')
                     
                 elif _line_[0] == '#include':
-                    print(wrg_3( _line_[1] + ' skipped!'))
+                    show( wrg_3( _line_[1] + ' skipped!'))
                 else:
-                    print(wrg_3( str(_line_) + '  ??'))
+                    show( wrg_3( str(_line_) + '  ??'))
     
     return _define_
 
@@ -1006,10 +1006,10 @@ def get_ffldfiles( _topfile_):
         # if there is at least one itp, lets parse it
         # first seek for further includes itp
         aux_file_cont = [_topfile_, '']
-        print '----- Loading :'
+        show( '----- Loading :')
         i = 1
         aux_file_cont[i] =  root_folder + ff_file
-        print aux_file_cont[i]
+        show( aux_file_cont[i])
         
         root_folder = '/'.join( aux_file_cont[i].split('/')[:-1]+[''])
         try:
@@ -1018,7 +1018,7 @@ def get_ffldfiles( _topfile_):
                     if k_line.startswith('#include'):
                         i+=1
                         aux_file_cont.append( root_folder+k_line.split('"')[1])
-                        print aux_file_cont[i]
+                        show( aux_file_cont[i] )
                         if i==3:
                             break
                 
@@ -1039,7 +1039,7 @@ def get_ffldfiles( _topfile_):
                     file_cont[-1] = file_cont[0] 
                 
             else:
-                print ('Using :' +file_cont[-1]+' for ' + _di_) 
+                show( 'Using :' + file_cont[-1] + ' for ' + _di_) 
     else:
         pop_err_1('Force field files #include not found!')
         nonerr_flag *= False
@@ -1060,26 +1060,26 @@ def ck_forcefield( _the_file_, _secondoption_ = None):
     podria pedirse solo este archivo y 
     de aqui sacar la iformacion de los otros dos....
     '''
-    _flag_ = False
-    comb_rule = -1
+    found_flag_ = False
+    comb_rule = [ 0, -1]
     with open( _the_file_, 'r')  as indata:
         for j_line in indata:
             line_c = j_line.split()
             if j_line.startswith('[ defaults ]'):
-                _flag_ = True
-            if len(line_c)>1 and 'fudgeQQ'==line_c[-1]:
-                j_line = indata.next()
-                comb_rule= j_line.split()
-                print('---> comb_rule {}'.format(comb_rule[1]))
+                found_flag_ = True
+            if len( line_c) > 1 and 'fudgeQQ'==line_c[-1]:
+                j_line = next( indata)
+                comb_rule = j_line.split()
+                show( '---> comb_rule {}'.format( comb_rule[1]))
                 
-    if not _flag_ and _secondoption_ != None:
-        comb_rule, _flag_, _the_file_ = ck_forcefield( _secondoption_)
+    if not found_flag_ and _secondoption_ != None:
+        comb_rule, found_flag_, _the_file_ = ck_forcefield( _secondoption_)
         
-    if comb_rule < 0 or not _flag_:
-        pop_err_1('forcefield.itp file is missing or incomplete')
-        comb_rule, _flag_, _the_file_ = [ 0, 0, '']
+    if int( comb_rule[1]) < 0 or not found_flag_:
+        pop_err_1( 'forcefield.itp file is missing or incomplete')
+        comb_rule, found_flag_, _the_file_ = [ [0,0], 0, '']
     
-    return comb_rule, _flag_, _the_file_
+    return comb_rule, found_flag_, _the_file_
 
 def seek_for_directive( _list_of_files_, _directive_):
     ''' search for a certain directive in a bunch of files
@@ -1104,9 +1104,7 @@ def seek_for_directive( _list_of_files_, _directive_):
     
     return content_file
     
-    
-    
-    
+
 def get_top_groups( _mtype_container_, _group_):
     
     _mtype_ = _mtype_container_
@@ -1125,8 +1123,9 @@ def get_top_groups( _mtype_container_, _group_):
             
             _buffer += [mt+1]
                 
-    print''
-    print 'Group characterized as: {} with ids {} to {}'.format(*_buffer)
+    show( '')
+    show( 'Group characterized as: {} with ids {} to {}'.format(*_buffer) )
+    
     return _buffer
 
 
diff --git a/lib/handling/lammps.py b/lib/handling/lammps.py
index aba141c..958b39f 100644
--- a/lib/handling/lammps.py
+++ b/lib/handling/lammps.py
@@ -3,7 +3,8 @@
 __author__ = 'Hernan Chavez Thielemann <hchavezthiele at gmail dot com>'
 
 from lib.misc.file import write_file, make_dir
-from lib.misc.warn import print_dec_g, pop_wrg_1, pop_err_1, wrg_1
+from lib.misc.display import show_in_green, show
+from lib.misc.warn import pop_wrg_1, pop_err_1, wrg_1 
 
 from sys import exit
 
@@ -16,7 +17,7 @@ def extract_lammps_data(_data_file_,_ck_buttons_, _forcefield_):
 
 def write_lammps_data( _topodata_, df_name, _config_):
     ''' Write a lammps data file'''
-    print_dec_g ('Writing Lammps data file...')
+    show_in_green('Writing Lammps data file...')
     ####---------------  Unpacking data  ----------------####
     atomstyle, _, _autoload_, root_folder = _config_
     atsty = [ 'atomic', 'angle', 'full', 'charge', 'bond', 'molecular']
@@ -25,9 +26,9 @@ def write_lammps_data( _topodata_, df_name, _config_):
     _content_=''
     if atomstyle in atsty:
         nam = ''.join([ chr( ord(l) - 32) for l in atomstyle])
-        print_dec_g(style_str.format(nam))
+        show_in_green(style_str.format(nam))
         
-        print '\n'+'='*14+'  Still in BETA here  '+'='*14+'\n'
+        show( '\n' + '='*14 + '  Still in BETA here  ' + '='*14 + '\n' ) #   !!
         _content_, _flag_ = write_lammps_data_auto( _topodata_,
                                                     df_name,
                                                     _config_[:-1]
@@ -37,13 +38,13 @@ def write_lammps_data( _topodata_, df_name, _config_):
         if _flag_:
             path_dir = make_dir( root_folder, 'g2l_dir')
             write_file( df_name, _content_, path_dir)
-            print_dec_g ('Successful writing!!')
+            show_in_green ('Successful writing!!')
         elif _autoload_:
             pop_err_1(errmsg + '\nAutoload failed!')
         else:
             pop_err_1(errmsg)
     else: #if atomstyle == 'Angle':
-        print '\n\nExit'
+        show( '\n\nExit')
         exit(('Error 037!!  -  Atom style {} '
               +'not implemented yet').format(atomstyle))
     
@@ -104,7 +105,7 @@ def write_lammps_data_auto( _topodata_, data_name, _config_):
     _text_ +=(' {:.4f} {:.4f} xlo xhi\n {:.4f} {:.4f} ylo yhi\n'
               +' {:.4f} {:.4f} zlo zhi\n').format(*_box_)
     
-    if _box_apend <> [0,0,0]:
+    if _box_apend != [0,0,0]:
         _text_ +=(' {:.4f} {:.4f} {:.4f} xy xz yz\n').format( *_box_apend)
     
     
@@ -124,16 +125,16 @@ def write_lammps_data_auto( _topodata_, data_name, _config_):
         if not float( _atom_mass_) and _sidemol_f_: #meaning that is 0
             
             smols = sidemol['tag']
-            #print 'here1', smols
+            #show( 'here1', smols
             for ad in range( len(smols)):
                 _smda_ = sidemol['data'][ad]['atoms']
-                #print _smda_
+                #show( _smda_
                 at_ = 0
                 for at_ in range( len( _smda_)):
-                    #print _smda_[at_], at_
+                    #show( _smda_[at_], at_
                     if _smda_[at_][1] == _atom_type_ and len(_smda_[at_]) == 8:
                         _atom_mass_ = _smda_[at_][7]
-                        print ( ('Mass for {} not found in atomtypes, taking' 
+                        show( ('Mass for {} not found in atomtypes, taking' 
                                  + ' {} as substitute').format( _atom_type_, 
                                                                 _atom_mass_)
                               )
@@ -181,7 +182,7 @@ def write_lammps_data_auto( _topodata_, data_name, _config_):
     elif atomstyle =='atomic':
         atom_shape = ' {0} {2} {4:7.4f} {5:7.4f} {6:7.4f} # {7}\n'
         
-    #print( dicts[ 0])
+    #show(( dicts[ 0])
     base_atoms_n = len( known_atoms)
     for i in range( base_atoms_n):
         aty = known_atoms[i][1]
@@ -207,7 +208,7 @@ def write_lammps_data_auto( _topodata_, data_name, _config_):
     _residue_buffer_    =   '00'
     _mol_buffer_        =   'nk_type'
     side_at_v           =   []
-    #print dicts[0]
+    #show( dicts[0]
     if _sidemol_f_ == 1:
         
         #charge = _topodata_['S_charge']
@@ -225,31 +226,31 @@ def write_lammps_data_auto( _topodata_, data_name, _config_):
         sm_m            =   1 # side mol multiplier
         for i in range( len ( sidemol['tag'])): 
             sm_at_num = len(sidemol['data'][ i]['atoms'])
-            #print (sidemol['tag'][i],sidemol['num'][i],sm_at_num)
+            #show( (sidemol['tag'][i],sidemol['num'][i],sm_at_num)
             sms_tags +=  [sidemol['tag'][i]]*sidemol['num'][i]*sm_at_num
             
         
-        print( 'Sidemols: ' + str( sum(sidemol['num'][:] ) ) )
+        show( 'Sidemols: ' + str( sum(sidemol['num'][:] ) ) )
         
         
-    #print( len(sms_tags), len(side_at_v), side_at_v[-1])
-        #print n_atoms, base_atoms_n, len( _mtype_), side_at_v, _atype_
+    #show(( len(sms_tags), len(side_at_v), side_at_v[-1])
+        #show( n_atoms, base_atoms_n, len( _mtype_), side_at_v, _atype_
     for i in side_at_v:
         _res_n_ = _mol_[i]
         #aty = conv_dict[_atype_[i]] # _atype_ is the atag in TOP
-        #print _atype_[i], aty
+        #show( _atype_[i], aty
         # meaning new residue or new molecule
-        clause1 = (_res_n_ <> _residue_buffer_)
+        clause1 = (_res_n_ != _residue_buffer_)
         #if clause1:
-            #print(clause1,_mol_[i], _residue_buffer_, multi_residue,prev_mult_res,sm_m)
+            #show((clause1,_mol_[i], _residue_buffer_, multi_residue,prev_mult_res,sm_m)
         if multi_residue and clause1: # mol number
             _residue_buffer_ = _res_n_
             multi_residue -= 1
             if not multi_residue and not sm_m:
                 _mol_buffer_ = 'nk_type'
                 
-            #print ji_, sm_aty[ ji_], _residue_buffer_
-            #print _mol_[i], bool(multi_residue)
+            #show( ji_, sm_aty[ ji_], _residue_buffer_
+            #show( _mol_[i], bool(multi_residue)
                 
         elif clause1: # mol number
             
@@ -258,9 +259,9 @@ def write_lammps_data_auto( _topodata_, data_name, _config_):
             
             # New molecule kind test
             # entering here first time too
-            #print(sm_m)
-            if sms_tags[ i] <> _mol_buffer_:
-                #print (i, sms_tags[ i], _mol_buffer_ )
+            #show((sm_m)
+            if sms_tags[ i] != _mol_buffer_:
+                #show( (i, sms_tags[ i], _mol_buffer_ )
                 # I can change mol-buffer for sms_tags[i-1], but case = 0
                 _mol_buffer_ = sms_tags[i]
                 _sm_i += 1
@@ -272,7 +273,7 @@ def write_lammps_data_auto( _topodata_, data_name, _config_):
                     sm_m -= 1
                     multi_residue = prev_mult_res
                     if sm_m == 0:
-                        print ('here 0')
+                        show( 'here 0')
                         
             else:
                  exit('**** Unhandled case 001')
@@ -281,28 +282,28 @@ def write_lammps_data_auto( _topodata_, data_name, _config_):
             _smol_tag_ = sidemol['tag'][ _sm_i]
             #sm_qty = sidemol['num'][ _sm_i]
             
-            #print i+1, sms_tags[i], _mtype_[i], _smol_tag_
-            #print '\n'+'here ', _sm_i
+            #show( i+1, sms_tags[i], _mtype_[i], _smol_tag_
+            #show( '\n'+'here ', _sm_i
                 
             # meaning also a new type of molecule ???
             if _smol_tag_ not in sm_cont.keys():
                 
                 new_smol_str = '** New side molecule : {} 1st atom : {} **'
-                print( '\n' + new_smol_str.format( _smol_tag_, _atype_[i]))
+                show( '\n' + new_smol_str.format( _smol_tag_, _atype_[i]))
                 ###########     New side mol data     ########
                 sm_data = sidemol['data'][  _sm_i]
                 ##############################################
                 first_atom_check = _atype_[i] == sm_data['atoms'][0][4]
-                #print( sm_data['atoms'][0][4])
+                #show(( sm_data['atoms'][0][4])
                 # In the weird case of a missing residue tag in the GRO file
                 # happens often
                 if _mtype_[i] == '':
-                    print ('\n')
+                    show( '\n')
                     pop_wrg_1('Missing residue tag in .gro file.\n'
                               + 'atom #'+ str( i+1)+ '. Asuming: ' +_smol_tag_)
                     
                     atom_x_ress = len( sm_data['atoms']) * sidemol['num'][ _sm_i]
-                    #print atom_x_ress, len( sm_data['atoms']), sidemol['num'][ _sm_i]
+                    #show( atom_x_ress, len( sm_data['atoms']), sidemol['num'][ _sm_i]
                     _mol_buffer_ = _smol_tag_
                     
                     if first_atom_check:
@@ -319,8 +320,8 @@ def write_lammps_data_auto( _topodata_, data_name, _config_):
                 
                 aux_buffer = ''
                 if not first_atom_check:
-                    print('xx/0   Atoms order mismatch '
-                         +'between residue definition and gro file!')
+                    show( ('xx/0   Atoms order mismatch '
+                         +'between residue definition and gro file!') )
                     
                     exit('xx/0   Error /lammps.py 001!')
                     
@@ -331,38 +332,38 @@ def write_lammps_data_auto( _topodata_, data_name, _config_):
                         
                         # multiple residues per side mol???
                         #lets see how many
-                        if ath[2] <> aux_buffer:
-                            #print ath[2]
-                            if aux_buffer <> '':
+                        if ath[2] != aux_buffer:
+                            #show( ath[2]
+                            if aux_buffer != '':
                                 multi_residue += 1
                             aux_buffer = ath[2]
                         
                     if multi_residue:
-                        print ( '** It is a multi residual structure.\n** '
-                               +'With: '+str( multi_residue +1)+' residues.')
+                        show( ( '** It is a multi residual structure.\n** '
+                               +'With: '+str( multi_residue +1)+' residues.') )
                         #multi_residue -= 1
                         prev_mult_res = multi_residue
                     
                 else:
-                    print('xx/0   Residue: ' + _mtype_[i])
+                    show( 'xx/0   Residue: ' + _mtype_[i] )
                     fst_line = (' ').join( sm_data['atoms'][0]) 
-                    print('xx/0   First line: ' + fst_line)
-                    print('xx/0   Molecules order mismatch '
-                         +'between [ molecules ] definition and gro file!')
-                    exit('xx/0   Error /lammps.py 002!')
+                    show( 'xx/0   First line: ' + fst_line )
+                    show( ('xx/0   Molecules order mismatch '
+                         +'between [ molecules ] definition and gro file!') )
+                    exit( 'xx/0   Error /lammps.py 002!')
                 try:
-                    print 'hop!',sm_data['atoms'][0][3]
+                    show( 'hop!', sm_data['atoms'][0][3] )
                     sm_cont[ _smol_tag_] = {}
                     sm_cont[ _smol_tag_]['data']    = sm_data
                     sm_cont[ _smol_tag_]['charge']  = sm_charge
                     sm_cont[ _smol_tag_]['aty']     = sm_aty
                     
                 except UnboundLocalError as Err_here:
-                    print('Upa mala cosa!')
+                    show( 'Upa mala cosa!')
                     exit( Err_here[0])
                 ######################################################
             else:
-                #print('-> {} +1'.format( _smol_tag_))
+                #show(('-> {} +1'.format( _smol_tag_))
                 sm_data     = sm_cont[ _smol_tag_]['data']
                 sm_charge   = sm_cont[ _smol_tag_]['charge']
                 sm_aty      = sm_cont[ _smol_tag_]['aty']
@@ -371,22 +372,22 @@ def write_lammps_data_auto( _topodata_, data_name, _config_):
             ''' \\\\\\\\\\\-------------------------------------////////// '''
             # This should be done in gromacs.py?? , in a near future
             # changes acording to the molecule kind topology
-            #print sm_data
+            #show( sm_data
             dirtv_data = sm_data[ 'bonds']
-            #print('here')
-            #print dirtv_data
-            #print( ji_, len(sm_aty))
+            #show(('here')
+            #show( dirtv_data
+            #show(( ji_, len(sm_aty))
             for xx in range( len( dirtv_data)):
                 _row = dirtv_data[ xx]
                 
-                #print _row
+                #show( _row
                 i1 =    int( _row[0])
                 i2 =    int( _row[1])
                 aty1 = sm_aty[ i1 - 1 ] 
                 aty2 = sm_aty[ i2 - 1 ]
                 
                 sm_bonds.append([ aty1 + '-' + aty2, i + i1, i + i2])
-                #print sm_bonds[-1]
+                #show( sm_bonds[-1]
                 #if sm_bonds[-1][0]=='HC-HC':
                 #    exit('opls_116-opls_116')
             dirtv_data = sm_data[ 'angles']
@@ -430,8 +431,8 @@ def write_lammps_data_auto( _topodata_, data_name, _config_):
                 sm_impropers.append([ imp_str, i + i1, i + i2, i + i3, i + i4])
         
         
-        #print ji_, sm_aty[ ji_]#, sm_aty
-        #print sm_aty[ ji_]
+        #show( ji_, sm_aty[ ji_]#, sm_aty
+        #show( sm_aty[ ji_]
         _text_ += atom_shape.format(i+1, _res_n_,
                                     dicts[0][ sm_aty[ ji_]],
                                     sm_charge[ ji_],
@@ -443,7 +444,7 @@ def write_lammps_data_auto( _topodata_, data_name, _config_):
     ###########################################################################
     '''==========-----------   4th - Chemical topology   ---------=========='''
     #=========================================================================#
-    print('\n> Checking chemical topology-coefficients')
+    show( '\n> Checking chemical topology-coefficients')
     ''' Building auxiliar atom tag1_tag2 data dictionary -/- OPLS Case '''
     xf = 1
     aat_ddic = {}
@@ -452,7 +453,7 @@ def write_lammps_data_auto( _topodata_, data_name, _config_):
             aat_ddic[ atom_info[i][0]] = atom_info[i][1]
     
     ####################        ------BONDS------          ####################
-    print('>> bonds')
+    show( '>> bonds')
     if _asty_d_[ atomstyle] >= 2:
         known_bonds = _topodata_['bonds']
         base_bonds_n = len (known_bonds)
@@ -461,10 +462,10 @@ def write_lammps_data_auto( _topodata_, data_name, _config_):
         
         a_g_d = {} #aux_goofy_dic
         for i in range(base_bonds_n):
-            # print known_bonds[i][0], known_bonds[i][1]
+            # show( known_bonds[i][0], known_bonds[i][1]
             at1 = int(known_bonds[i][0])
             at2 = int(known_bonds[i][1])
-            #print at1, at2
+            #show( at1, at2
             to_print = '{} {}'.format(at1, at2)
             try:
                 at1_tag = known_atoms[at1-1][xf]
@@ -479,34 +480,34 @@ def write_lammps_data_auto( _topodata_, data_name, _config_):
                     _bond_ty_ = dicts[1][ at1_tag + '-'+ at2_tag]
                     
                 except KeyError:
-                    print at1_tag + '-'+ at2_tag
-                    print('Bond type not found for: ' + to_print)
-                    print 'Exception {}-{}\n'.format( known_atoms[at1-1][4], 
-                                               known_atoms[at2-1][4])
+                    show( at1_tag + '-'+ at2_tag)
+                    show( 'Bond type not found for: ' + to_print)
+                    show( 'Exception {}-{}\n'.format( known_atoms[at1-1][4], 
+                                               known_atoms[at2-1][4]) )
                     _bond_ty_ = 0
             
             _text_ += bond_shape.format( i+1, _bond_ty_, at1, at2)
                         
         if _sidemol_f_ == 1:
-            #print dicts[1]
+            #show( dicts[1]
             # better way to do this is trough corrds ---------  <WFS>
             for i in range( n_bonds - base_bonds_n):
-                # print sm_bonds[i]
+                # show( sm_bonds[i]
                 try:
                     _bond_ty_ = dicts[1][ sm_bonds[i][0]]
                 except KeyError:
                     # OPLS ??
                     try:
                         at1_tag, at2_tag = sm_bonds[i][0].split('-')
-                        #print(at1_tag, at2_tag)
+                        #show((at1_tag, at2_tag)
                         at1_tag = aat_ddic[ at1_tag]
                         at2_tag = aat_ddic[ at2_tag]
                         _bond_ty_ = dicts[1][ at1_tag + '-'+ at2_tag]
                         
                     except KeyError:
                         
-                        print('Bond type not found for:' + sm_bonds[i][0])
-                        print('Or ' + at1_tag + '-'+ at2_tag)
+                        show( 'Bond type not found for:' + sm_bonds[i][0] )
+                        show( 'Or ' + at1_tag + '-'+ at2_tag )
                         exit( dicts[1] )
                         
                 _text_ += bond_shape.format(i+1 + base_bonds_n,
@@ -516,8 +517,8 @@ def write_lammps_data_auto( _topodata_, data_name, _config_):
     
     
     ####################        ------ANGLES------      #######################
-    print('>> angles')
-    if _asty_d_[ atomstyle] >= 3:
+    show( '>> angles')
+    if _asty_d_[ atomstyle] >= 3 and n_angles > 0:
         known_angles = _topodata_['angles']
         base_angles_n = len(known_angles)
         _text_ +='\n Angles\n\n'
@@ -540,9 +541,9 @@ def write_lammps_data_auto( _topodata_, data_name, _config_):
                     at3_tag = aat_ddic[ at3_tag]
                     _angle_ty_ = dicts[2][at1_tag + '-'+ at2_tag+ '-'+ at3_tag]
                 except KeyError as Er_here:
-                    print 'Error ----- '+Er_here.args[0]
+                    show( 'Error ----- '+Er_here.args[0] )
                     
-            #print angle_t, _angle_ty_
+            #show( angle_t, _angle_ty_
             _text_ += angle_shape.format( i+1, _angle_ty_, at1, at2, at3)
             
         if _sidemol_f_ == 1:
@@ -568,8 +569,8 @@ def write_lammps_data_auto( _topodata_, data_name, _config_):
     
     
     ####################        ------DIHEDRAL------       ####################
-    print('>> dihedrals')
-    if _asty_d_[ atomstyle] >= 4:
+    show( '>> dihedrals')
+    if _asty_d_[ atomstyle] >= 4 and n_dihedrals > 0:
         known_dihedrals = _topodata_['dihedrals']
         base_dihedrals_n = len(known_dihedrals)
         if base_dihedrals_n or ( n_dihedrals > base_dihedrals_n):
@@ -615,10 +616,10 @@ def write_lammps_data_auto( _topodata_, data_name, _config_):
                                                          'X-X']
                                 except KeyError as Er_here:
                                     err_str = Er_here.args[0]
-            if _dihe_ty_ == '0' and err_str <> '':
-                print 'Atoms {}-{}-{}-{} : '.format( at1, at2, at3, at4),
-                print (at1_tag +'-'+ at2_tag+'-'+ at3_tag+'-'+ at4_tag)
-                print ( 'Error dihedral ----- '+ err_str +' not found!')
+            if _dihe_ty_ == '0' and err_str != '':
+                show( 'Atoms {}-{}-{}-{} : '.format( at1, at2, at3, at4),
+                     ( at1_tag + '-' + at2_tag + '-' + at3_tag + '-'+ at4_tag))
+                show( 'Error dihedral ----- ' + err_str +' not found!')
             
             _text_+= dihedral_shape.format( i+1, _dihe_ty_, at1, at2, at3, at4)
             
@@ -650,8 +651,8 @@ def write_lammps_data_auto( _topodata_, data_name, _config_):
                             except KeyError as Er_here:
                                 err_str += ' // ' + Er_here.args[0]
                         
-                if _dihe_ty_ == '0' and err_str <> '':
-                    print( 'Atoms {}-{}-{}-{} '.format( *aux_here) + ' '
+                if _dihe_ty_ == '0' and err_str != '':
+                    show( 'Atoms {}-{}-{}-{} '.format( *aux_here) + ' '
                            + err_str + '  //  ' + 
                            wrg_1( 'Dihedral > '+ sm_dihedrals[i][0] +
                                  ' < not found!')
@@ -666,9 +667,9 @@ def write_lammps_data_auto( _topodata_, data_name, _config_):
                                                )
     
     ###################        ------IMPROPERS------       ####################
-    print('>> impropers')
+    show( '>> impropers')
     # TODO SECTION
-    if _asty_d_[ atomstyle] >= 4:
+    if _asty_d_[ atomstyle] >= 4 and n_impropertypes > 0:
         known_impropers = _topodata_['impropers']
         base_impropers_n = len( known_impropers)
         if base_impropers_n or ( n_impropers > base_impropers_n):
@@ -711,11 +712,11 @@ def write_lammps_data_auto( _topodata_, data_name, _config_):
                             if _opt_ == options[-1]:
                                 err_str = Er_here.args[0]
 
-            if _impr_ty_ == '0' and err_str <> '':
-                print 'Atoms {}-{}-{}-{} : '.format( at1, at2, at3, at4),
-                print (at1_tag +'-'+ at2_tag+'-'+ at3_tag+'-'+ at4_tag)
-                print ( 'Error improper dihedral ---- '+ err_str +' not found!')
-                
+            if _impr_ty_ == '0' and err_str != '':
+                show( 'Atoms {}-{}-{}-{} : '.format( at1, at2, at3, at4),
+                     (at1_tag +'-'+ at2_tag+'-'+ at3_tag+'-'+ at4_tag) )
+                show( ( 'Error improper dihedral ---- '+ err_str +' not found!')
+                    )
             
             _text_+= improper_shape.format( i+1, _impr_ty_, at1, at2, at3, at4)
             
@@ -742,10 +743,10 @@ def write_lammps_data_auto( _topodata_, data_name, _config_):
                         except KeyError as Er_here:
                             err_str += Er_here.args[0]
                 
-                if _impr_ty_ == '0' and err_str <> '':
-                    print( wrg_1('Improper dihedral > ' + err_str
+                if _impr_ty_ == '0' and err_str != '':
+                    show( wrg_1('Improper dihedral > ' + err_str
                                              + ' < not found!') 
-                             )
+                         )
                 
                 _text_ += improper_shape.format( i+1 + base_impropers_n,
                                                  _impr_ty_,
@@ -763,20 +764,20 @@ def side_mol_topology_builder( data_in):
     # This should be done in gromacs.py?? , in a near future
     # changes acording to the molecule kind topology
     tag_str = [ 'bonds', 'angles', 'dihedrals', 'impropers']
-    print conv_dict
+    show( conv_dict, v = 2 )
     for _st_ in tag_str:# range( len( tag_str)):
-        #print _st_, i
+        #show( _st_, i
         
         for xx in range( len( sm_data[ _st_])):
             index_c = sm_data[ _st_][ xx]
-            #print index_c
+            #show( index_c
             i1 = int(index_c[0])
             i2 = int(index_c[1])
-            print _atype_[ i + i1-1],  _atype_[ i + i2-1]
+            show( _atype_[ i + i1-1],  _atype_[ i + i2-1] )
             
             aty1 = conv_dict[ _atype_[ i + i1-1]]
             aty2 = conv_dict[ _atype_[ i + i2-1]]
-            print aty1, aty2
+            show( aty1, aty2 )
             
             if _st_ == 'bonds':
                 sm_bonds.append([ aty1+'-'+aty2, i1, i2])
@@ -809,7 +810,7 @@ def side_mol_topology_builder( data_in):
                 sm_impropers.append( [ impr_str, i1, i2, i3, i4])
                 
             else:
-                print aty, _mol_[i], _residue_buffer_, len( _atype_)
+                show( aty, _mol_[i], _residue_buffer_, len( _atype_) )
                 exit('')
         
     
@@ -845,7 +846,7 @@ def write_lammps_potentials( _topodata_, atomstyle = 'full'):
         for x in range( n_atomtypes):
             _A_ = float( atom_info[x][ -1 + minr])
             _B_ = float( atom_info[x][ -2 + minr])
-            if _A_ <> 0 and _B_ == 0 or _A_ == 0 and _B_ <> 0:
+            if _A_ != 0 and _B_ == 0 or _A_ == 0 and _B_ != 0:
                 regular_se = False
                 pairtypename = 'hybrid/overlay'
                 wr_str = ('Using pair style lennard/mdf!!\nThis style can '
@@ -855,7 +856,7 @@ def write_lammps_potentials( _topodata_, atomstyle = 'full'):
                 break
     
     for x in range( n_atomtypes):
-        #print atom_info[x]
+        #show( atom_info[x]
         atom_type_d[atom_info[x][0]] = x+1
         
         _A_ = float( atom_info[x][ -1 + minr])
@@ -918,7 +919,7 @@ def write_lammps_potentials( _topodata_, atomstyle = 'full'):
         bondtypename.append( BondDataBase[ int( bo)- 1])
     
     txt_p_b = ''
-    if bondtypename <> []:
+    if bondtypename != []:
         txt_p_b ='\n Bond Coeffs #{}\n\n'.format( bondtypename[0]) # bond_style hybrid
     
     bondtypes_d = {}
@@ -935,7 +936,7 @@ def write_lammps_potentials( _topodata_, atomstyle = 'full'):
                           float( bty[i][3])*10 ]
             
         elif _bi_ == 2: # bond_style  G96 in gromacs
-            print bty[i]
+            show( bty[i] )
             info_cont = [ i+1, _bdb_[ _bi_ - 1],
                           float( bty[i][4])/ 100/ 4.186/4,
                           float( bty[i][3])*10 ]
@@ -975,7 +976,7 @@ def write_lammps_potentials( _topodata_, atomstyle = 'full'):
         angletypename.append( AngleDataBase[ int( an)- 1])
     
     txt_p_a = ''
-    if angletypename <> []:
+    if angletypename != []:
         txt_p_a ='\n Angle Coeffs #{}\n\n'.format( angletypename[0]) 
     
     angletypes_d = {}
@@ -1036,7 +1037,7 @@ def write_lammps_potentials( _topodata_, atomstyle = 'full'):
             type9_flag = True
     
     txt_p_d = ''
-    if dihedtypename <> []:
+    if dihedtypename != []:
         txt_p_d ='\n Dihedral Coeffs #{}\n\n'.format( dihedtypename[0])
     
     rb_warning = (' Ryckaert-Bellemans angle style conversion in Fourier form' +
@@ -1048,11 +1049,11 @@ def write_lammps_potentials( _topodata_, atomstyle = 'full'):
         pop_wrg_1(' Type 9 Proper dihedral (multiple) angle style conversion'
                   + ' in Fourier form can only be used if LAMMPS was built'
                   + ' with the USER-MISC package!!')
-        print( '='*20 + '  MuXInG  ' + '='*20 )
+        show( '='*20 + '  MuXInG  ' + '='*20 )
         aux_type9 = []
         aux_type9_types = []
         new_dty = []
-        print('Dihedrals before muxing: {}'.format( n_dihedraltypes))
+        show( 'Dihedrals before muxing: {}'.format( n_dihedraltypes))
         for i in range( n_dihedraltypes):
             if int( dty[i][4]) == 9:
                 aux_tag = '{}-{}-{}-{}'.format( *dty[i][:4])
@@ -1069,7 +1070,7 @@ def write_lammps_potentials( _topodata_, atomstyle = 'full'):
         dty = new_dty + aux_type9
         n_dihedraltypes = len( dty) #########Check cases when this is not valid?
         
-        print('Dihedrals after muxing: {}'.format( n_dihedraltypes))
+        show( 'Dihedrals after muxing: {}'.format( n_dihedraltypes))
     ###== end Type 9
     
     dihedraltypes_d = {} # types dictionary
@@ -1157,7 +1158,7 @@ def write_lammps_potentials( _topodata_, atomstyle = 'full'):
         imprtypename.append( ImproperDataBase[ int( im)- 1])
     
     txt_p_i = ''
-    if imprtypename <> []:
+    if imprtypename != []:
         txt_p_i ='\n Improper Coeffs #{}\n\n'.format( imprtypename[0])
             
     impropertypes_d = {} # types dictionary
@@ -1224,7 +1225,7 @@ def write_lammps_potentials( _topodata_, atomstyle = 'full'):
         dicts = [atom_type_d]
         txt_p_ = txt_p_p
     else:
-        print ('\nWeird thing, it is supposed impossible to reach this place\n')
+        show( '\nWeird thing, it is supposed impossible to reach this place\n')
         _flag_ = False
         
         
@@ -1252,7 +1253,7 @@ def write_lammps_input(  _simconfig_, _topodata_= None, in_name= 'in.gro2lam'):
     npt_pstart, npt_pend = npt_pss.split(':')
     npt_tstart, npt_tend = npt_tss.split(':')
     
-    #print (data_file, timestep, nve_steps, nvt_steps, nvt_tstart, nvt_tend,
+    #show( (data_file, timestep, nve_steps, nvt_steps, nvt_tstart, nvt_tend,
     #nvt_tdamp, npt_steps, npt_pstart, npt_pend, npt_tdamp, npt_tdamp,
     #npt_ystart, npt_yend)
     
@@ -1271,9 +1272,9 @@ def write_lammps_input(  _simconfig_, _topodata_= None, in_name= 'in.gro2lam'):
     group_lines = ''
     torestrain = []
     ens_af = []
-    if _simconfig_[2] <> []:
+    if _simconfig_[2] != []:
         g_names, g_aids, k_xyz_c, runs_c, ch_init = _simconfig_[2][0][:]
-        if _simconfig_[2][1] <> None:
+        if _simconfig_[2][1] != None:
             ####### ------------------------------------ Super interesante!!
             ##              este es uno de esos casos donde no es posible 
             ##              utilizar += a pesar de desligar con [:] ... 
@@ -1284,10 +1285,10 @@ def write_lammps_input(  _simconfig_, _topodata_= None, in_name= 'in.gro2lam'):
             k_xyz_c = k_xyz_c + aux3
             runs_c  = runs_c  + aux4
             ch_init = ch_init + aux5
-        print'\n'
+        show( '\n')
         for re in range(len(g_names)):
             if ch_init[re]==1:
-                print 'Restraining group '+g_names[re]+' in '+runs_c[re]
+                show( 'Restraining group ' + g_names[re] + ' in ' + runs_c[re])
                 groupinfo = [g_names[re], g_aids[re]]
                 group_lines += 'group {} id {}\n'.format( *groupinfo)
                 
@@ -1298,7 +1299,7 @@ def write_lammps_input(  _simconfig_, _topodata_= None, in_name= 'in.gro2lam'):
                     for e in ens:
                         if e not in ens_af:
                             ens_af.append(e)
-        if group_lines <> '':
+        if group_lines != '':
             group_lines +='\n'
     
     mix_value = {'1':'geometric', '2':'arithmetic',
@@ -1313,13 +1314,13 @@ def write_lammps_input(  _simconfig_, _topodata_= None, in_name= 'in.gro2lam'):
     #===================================================
     ####------------      TOPO DATA         --------####
     
-    print '\n'+data_file + '\n'
-    if _topodata_ <> None:
+    show( '\n'+data_file + '\n')
+    if _topodata_ != None:
         atomstyle_o, _solvated_, _autoload_, root_folder = _topodata_['config']
         _, comb_rule, _, _, _ = _topodata_['defaults']
         
     else:
-        print '**** Without _topodata_ !!'
+        show( '**** Without _topodata_ !!')
         root_folder = './'
         atomstyle_o = ''
         comb_rule = ''
@@ -1328,17 +1329,17 @@ def write_lammps_input(  _simconfig_, _topodata_= None, in_name= 'in.gro2lam'):
     _aux_her = get_style_info( data_file)
     
     atomstyle, pairstyle, bondstyle, anglstyle, dihestyle, imprstyle = _aux_her
-    if atomstyle_o <> '' and atomstyle_o <> atomstyle[0]:
+    if atomstyle_o != '' and atomstyle_o != atomstyle[0]:
         pop_wrg_1( 'Incongruence between atom styles!')
-        print atomstyle_o, atomstyle[0]
+        show( atomstyle_o, atomstyle[0])
     
     ## ---------------- MIXING RULE
     if f_comb_rule in mix_value.values():
-        mix_value_s=' mix '+f_comb_rule
-    elif f_comb_rule=='from_gromacs' and _topodata_<>None:
-        mix_value_s=' mix '+mix_value[comb_rule]
+        mix_value_s = ' mix ' + f_comb_rule
+    elif f_comb_rule == 'from_gromacs' and _topodata_ != None:
+        mix_value_s = ' mix ' + mix_value[comb_rule]
     else:
-        print 'Using default mixing rule'
+        show( 'Using default mixing rule')
         mix_value_s = ''
         
 
@@ -1348,7 +1349,7 @@ def write_lammps_input(  _simconfig_, _topodata_= None, in_name= 'in.gro2lam'):
     _dtxt_= '# Generated with Gro2lam\n\n'+'units real\nboundary p p p\n'
     # as I understand lammps default is 3
     #_dtxt_+= '%s %d\n'.format('dimension',dimension)
-    _dtxt_+= 'atom_style '+atomstyle[0]+'\n'
+    _dtxt_ += 'atom_style ' + atomstyle[0] + '\n'
     if pairstyle[0] == '':
         if atomstyle[0] not in ['full', 'charge]']: # no charges
             if 'coul' in pairwiseint:
@@ -1375,28 +1376,28 @@ def write_lammps_input(  _simconfig_, _topodata_= None, in_name= 'in.gro2lam'):
     # options like full and non bonded interactions could be asked by the user
     _dsc_txt=['pair_style {} {}'.format( pairwiseint, lj_rcutoff)]
     _dsc_txt.append(' {}\n'.format( c_rcutoff))
-    if bondstyle[0] <> '':
+    if bondstyle[0] != '':
         _dsc_txt.append( 'bond_style '+' '.join( bondstyle)+'\n')
-    if anglstyle[0] <> '':
+    if anglstyle[0] != '':
         _dsc_txt.append( 'angle_style '+' '.join( anglstyle)+'\n')
-    if dihestyle[0] <> '':
+    if dihestyle[0] != '':
         _dsc_txt.append( 'dihedral_style '+' '.join( dihestyle)+'\n')
-    if imprstyle[0] <> '':
+    if imprstyle[0] != '':
         _dsc_txt.append( 'improper_style '+' '.join( imprstyle)+'\n')
     _dtxt_+= ''.join(_dsc_txt[:_asty_d_[atomstyle[0]]])+'\n'
     
     
-    if 'data' in data_file:
-        _dtxt_+= 'read_data {}\n'.format( data_file)
+    if 'restart' in data_file:
+        _dtxt_+= 'read_restart {}\n'.format( data_file.split('/')[-1])
     else:
-        _dtxt_+= 'read_restart {}\n'.format( data_file)
+        _dtxt_+= 'read_data {}\n'.format( data_file.split('/')[-1])
     
     #===================================================
     ####--------------   NEIGHBOR LIST   -----------####
     
     _dtxt_+= '\nneighbor {} bin\n'.format( neighbordistance)
     
-    if lrsolver <> '' and atomstyle[0] in ['full','charge']:
+    if lrsolver != '' and atomstyle[0] in ['full','charge']:
         if 'long' in pairwiseint:
             _dtxt_+= 'kspace_style {} {}\n'.format( lrsolver, lrerror)
         
@@ -1409,11 +1410,11 @@ def write_lammps_input(  _simconfig_, _topodata_= None, in_name= 'in.gro2lam'):
             aux_txt = ' {} {} {} coul {}'.format( *aux_here1)
             sp_bon_3 = [ aux_txt] + aux_here2[1:]
                             
-    elif lrsolver <> '':
+    elif lrsolver != '':
         sp_bon_3 = lj12_13_14.split(':')
         
     
-    if lrsolver <> '':
+    if lrsolver != '':
         _dtxt_+= 'special_bonds lj{} {} {}\n'.format( *sp_bon_3)
         
     _dtxt_+= 'pair_modify shift no tail yes'+mix_value_s+'\n'
@@ -1442,16 +1443,16 @@ def write_lammps_input(  _simconfig_, _topodata_= None, in_name= 'in.gro2lam'):
     ####---------------     SHAKE       ------------####
     shake_bn, shake_an
     shake_txt = 'fix shake_name1 all shake {} 20 0'.format( shake_tol)
-    print shake_bn, shake_an
-    if shake_bn <> '0' or shake_an <> '0':
+    show( shake_bn, shake_an)
+    if shake_bn != '0' or shake_an != '0':
         
-        if shake_bn <> '0':
+        if shake_bn != '0':
             shake_txt += ' b'
             shake_bn = shake_bn.split('-')
             for bn in range(len(shake_bn)):
                 shake_txt += ' '+shake_bn[bn]
 
-        if shake_an <> '0':
+        if shake_an != '0':
             shake_txt += ' a'
             shake_an = shake_an.split('-')
             for an in range(len(shake_an)):
@@ -1469,7 +1470,7 @@ def write_lammps_input(  _simconfig_, _topodata_= None, in_name= 'in.gro2lam'):
     for en in range(len(ensembles)):
         
         
-        if ens_af<>[] and en in ens_af: #       RESTRAIN 
+        if ens_af != [] and en in ens_af: #       RESTRAIN 
             
             for re in range(len(torestrain)):
                 if en in torestrain[re][2]:
@@ -1488,13 +1489,13 @@ def write_lammps_input(  _simconfig_, _topodata_= None, in_name= 'in.gro2lam'):
                         tounfix= [ tounfix[0]+ [unr], tounfix[1]+ [name_2uf]]
             _dtxt_ += '\n'
             
-        if ensembles[en]=='NVE' and nve_steps <> '' and nve_steps.isdigit():
+        if ensembles[en]=='NVE' and nve_steps != '' and nve_steps.isdigit():
             steps = int(nve_steps)
             nve_frm = 'fix nve_name1 all nve\nrun {}\nunfix nve_name1\n\n'
             _dtxt_ += nve_frm.format(steps)
             curr_time += steps*timestep
         
-        elif ensembles[en]=='NPT' and npt_steps <> '' and npt_steps.isdigit():
+        elif ensembles[en]=='NPT' and npt_steps != '' and npt_steps.isdigit():
             steps = int(npt_steps)
             npt_frm = 'fix npt_name1 all npt temp {} {} {} {} {} {} {}\n'
             _dtxt_ += npt_frm.format( npt_tstart, npt_tend, npt_tdamp,
@@ -1502,7 +1503,7 @@ def write_lammps_input(  _simconfig_, _topodata_= None, in_name= 'in.gro2lam'):
             _dtxt_+= 'run {}\nunfix npt_name1\n\n'.format( steps)
             curr_time += steps*timestep
         
-        elif ensembles[en]=='NVT' and nvt_steps <> '' and nvt_steps.isdigit():
+        elif ensembles[en]=='NVT' and nvt_steps != '' and nvt_steps.isdigit():
             steps = int(nvt_steps)
             nvt_frm = 'fix nvt_name1 all nvt temp {} {} {}\n'
             _dtxt_ += nvt_frm.format( nvt_tstart, nvt_tend, nvt_tdamp )
@@ -1524,18 +1525,18 @@ def write_lammps_input(  _simconfig_, _topodata_= None, in_name= 'in.gro2lam'):
             _dtxt_+= emin_frm.format( e_min_tol, f_min_tol )
             
             
-        if tounfix <> [ [], []] and en in tounfix[0]: #       UNFIX RESTRAIN
+        if tounfix != [ [], []] and en in tounfix[0]: #       UNFIX RESTRAIN
             for unf in range( len( tounfix[0])):
                 if tounfix[0][unf] == en:
                     _dtxt_ += 'unfix ' + tounfix[1][unf] + '\n'
     
-    print ('Writing Lammps input script...')
-    if root_folder <> './':
+    show( 'Writing Lammps input script...')
+    if root_folder != './':
         _folder_ = '/'.join( data_file.split('/')[:-1]+[''])
     else:
         _folder_ = root_folder
     write_file( in_name , _dtxt_, _folder_)
-    print_dec_g ('Successful writing!!')
+    show_in_green( 'Successful writing!!')
     
     
      #-------------------- here would be optimum some further check
@@ -1566,16 +1567,16 @@ def get_style_info( lammps_datafile):
                 elif len (line_c) > 3 and  line_c[2] == 'xlo':
                     break
             sty_qty[ 0] = sty_qty[ 5]
-            print 'Quantities : ',
+            show( 'Quantities : ', end = ' ')
             for st in range( len( styles)):
-                print styles[st] + ' : ' + str( sty_qty[st]),
-            print '\n'
+                show( styles[st] + ' : ' + str( sty_qty[st]), end = ' ')
+            show( '\n')
             read_flag = False
             reading_flag = False
             index = 0
             for k_line in indata:
                 line_c =  k_line.split()
-                #print line_c
+                #show( line_c
                 # in cases with hybrid type, it should read from the second 
                 # position in the data line of bond, angle, dihedral and impr
                 if read_flag:
@@ -1587,14 +1588,15 @@ def get_style_info( lammps_datafile):
                         while not sty_qty[ index]:
                             sty_container[ index].append( '')
                             index += 1
-                        print '... done!'
+                        show( '... done!')
                     elif len( line_c) < 1 or line_c[0][0] == '#':
                         pass
                     else:
                         reading_flag = True
                         new_sty = line_c[1]
                         if new_sty not in sty_container[ index]:
-                            print styles[ index] + 'Style : ' , new_sty,
+                            show( styles[ index] + 'Style : ' , new_sty,
+                                 end = ' ')
                             sty_container[ index].append( new_sty)
                     
                 # cutting out the crap, normal empty and commented lines
@@ -1604,7 +1606,7 @@ def get_style_info( lammps_datafile):
                 elif  ( styles[ index] == line_c[0] and 
                       ( index == 5 or 'Coeffs' == line_c[1]) ):
                     
-                    print index, styles[ index], k_line.rstrip()
+                    show( index, styles[ index], k_line.rstrip() )
                     aux_cont = k_line.split('#')
                     
                     if len( aux_cont) > 1:
@@ -1624,12 +1626,12 @@ def get_style_info( lammps_datafile):
                     
     except IOError:
         pop_wrg_1( 'Data file not found!')
-        print ( 'Maybe try performing a conversion first! ;)')
+        show( 'Maybe try performing a conversion first! ;)')
         _flag_ = False
     
-    print('\n')
+    show( '\n')
     sty_container = [sty_container[5],] + sty_container[:5]
-    print sty_container
+    show( sty_container, v = 2)
     return sty_container
 
 if __name__ == '__main__':
diff --git a/lib/misc/data.py b/lib/misc/data.py
index d532f73..a418c8c 100644
--- a/lib/misc/data.py
+++ b/lib/misc/data.py
@@ -2,13 +2,14 @@
 #    By Hernan Chavez Thielemann
 __author__ = 'Hernan Chavez Thielemann <hchavezthiele at gmail dot com>' 
 
-from warn import pop_err_1
+from lib.misc.warn import pop_err_1
+#from display import show
 
 def check_vars( _vars_, _varstype_, _endtext_ = 'Input script not created!'):
     
     _flag_c_ = []
     for v in range(len(_vars_)):
-        #print _vars_[v], _varstype_[v]
+        #show( _vars_[v], _varstype_[v]
         err_txt = ''
         if _varstype_[v] =='int':
             try:
@@ -22,7 +23,7 @@ def check_vars( _vars_, _varstype_, _endtext_ = 'Input script not created!'):
                 err_txt = ' {} is not a valid type of integer'.format( aux_here)
                 
         elif _varstype_[v] =='float':
-            #print _vars_[v]
+            #show( _vars_[v]
             try:
                 float(_vars_[v])
                 _flag_c_.append(1)
@@ -45,7 +46,7 @@ def check_vars( _vars_, _varstype_, _endtext_ = 'Input script not created!'):
             
         elif _varstype_[v] == 'int:':
             _aux_ = _vars_[v].split(':')
-            if len( _aux_) <> 2:
+            if len( _aux_) != 2:
                 err_txt = ( 'Please separate two integers with one colon (:)'
                            +' instead of {}'.format(_vars_[v]))
                 _flag_c_ += [0,0]
@@ -54,7 +55,7 @@ def check_vars( _vars_, _varstype_, _endtext_ = 'Input script not created!'):
             
         elif _varstype_[v] =='float:':
             _aux_ = _vars_[v].split(':')
-            if len( _aux_) <> 2:
+            if len( _aux_) != 2:
                 err_txt = ( 'Please separate two floats with one colon (:)'
                            +' instead of {}'.format(_vars_[v]))
                 _flag_c_ += [0,0]
@@ -64,7 +65,7 @@ def check_vars( _vars_, _varstype_, _endtext_ = 'Input script not created!'):
                         
         elif _varstype_[v] == 'float:xyz':
             _aux_ = _vars_[v].split(':')
-            if len( _aux_) <> 2:
+            if len( _aux_) != 2:
                 err_txt = ( 'Please separate the integer and "xyz" with one'
                            +'colon (:) instead of {}'.format(_vars_[v]))
                 _flag_c_ += [0,0]
@@ -86,7 +87,7 @@ def check_vars( _vars_, _varstype_, _endtext_ = 'Input script not created!'):
             
             _aux_v_ = [x.strip(' ') for x in _vars_[v].split(':')]
             _min, _max = _varstype_[v][1:]
-            if len( _aux_v_) <> 3:
+            if len( _aux_v_) != 3:
                 err_txt = ( 'Please separate three floats with two colons (:)'
                            +' instead of {}'.format(_vars_[v]))
                 _flag_c_ += [0,0,0]
@@ -100,7 +101,7 @@ def check_vars( _vars_, _varstype_, _endtext_ = 'Input script not created!'):
         elif type(_varstype_[v]) == list and _varstype_[v][0] == '<int<':
             _min, _max = _varstype_[v][1:]
             _flag_ = check_vars( [_vars_[v]], ['int'], _endtext_)
-            #print _flag_c_, _flag_
+            #show( _flag_c_, _flag_
             if _flag_[0]:
                 if _min <= int(_vars_[v]) <= _max:
                     _flag_c_ += _flag_
@@ -111,12 +112,12 @@ def check_vars( _vars_, _varstype_, _endtext_ = 'Input script not created!'):
                     _flag_c_ += [0]
             else:
                 _flag_c_ += _flag_
-            #print _flag_c_
+            #show( _flag_c_
             
         elif type(_varstype_[v]) == list and _varstype_[v][0] == '<float<':
             _min, _max = _varstype_[v][1:]
             _flag_ = check_vars( [_vars_[v]], ['float'], _endtext_)
-            #print _flag_c_, _flag_
+            #show( _flag_c_, _flag_
             if _flag_[0]:
                 if _min <= float(_vars_[v]) <= _max:
                     _flag_c_ += _flag_
@@ -129,7 +130,7 @@ def check_vars( _vars_, _varstype_, _endtext_ = 'Input script not created!'):
                 _flag_c_ += _flag_
         
         elif type(_varstype_[v]) == list and _varstype_[v][0]==list:
-            #print _varstype_[v]
+            #show( _varstype_[v]
             _sep_ = _varstype_[v][1]
             _aux_v_ = [x.strip(' ') for x in _vars_[v].split(_sep_)]
             
@@ -139,7 +140,7 @@ def check_vars( _vars_, _varstype_, _endtext_ = 'Input script not created!'):
                 else:
                     err_txt += ( 'Unknown value <{}>.\n'.format( _v_))
                     _flag_c_.append(0)
-            if err_txt <> '':
+            if err_txt != '':
                 err_txt += '\nAllowed values are:\n'
                 for av in _varstype_[v][2:-1]:
                     err_txt += '"{}" or '.format(av)
@@ -203,13 +204,13 @@ def check_vars( _vars_, _varstype_, _endtext_ = 'Input script not created!'):
                 err_txt = 'Not valid number, check {}!'.format(_vars_[v])
                 _flag_c_ += [0]
         else:
-            #print _varstype_[v]
-            print 'Unhandled variable check for: ', _vars_[v]
+            #show( _varstype_[v]
+            show( 'Unhandled variable check for: ' , _vars_[v])
             _flag_c_.append(1)
         
         '''**** Message ****'''
-        if err_txt<>'':
-            pop_err_1( err_txt+'\n'+_endtext_)
+        if err_txt != '':
+            pop_err_1( err_txt + '\n' + _endtext_)
     
         
     return _flag_c_
diff --git a/lib/misc/display.py b/lib/misc/display.py
new file mode 100644
index 0000000..52fd462
--- /dev/null
+++ b/lib/misc/display.py
@@ -0,0 +1,127 @@
+#!/usr/bin/python
+#    By Hernan Chavez Thielemann
+#    this came out of a blend of warn and verbose 
+
+from sys import stdout
+
+__verbose__ = [True,0,0,0]
+
+def in_red_(_text_):
+    return '\033[91m' + _text_ + '\033[0m'
+
+def in_green_(_text_):
+    return '\033[92m' + _text_ + '\033[0m'
+
+def in_yellow_(_text_):
+    return '\033[93m' + _text_ + '\033[0m'
+
+def yellow_deco(func):
+    def func_wrapper(name):
+        return "\033[93m{0}\033[0m".format(func(name))
+    return func_wrapper
+#@yellow_deco
+
+def show_in_green(*text):
+    '''  this should be replaced with a colour_show with sep = "\n"
+        now this looks like nonsense ;)
+    '''
+    n_text = ''
+    if text != None:
+        for line in text[:-1]:
+            n_text += line +'\n'
+        n_text += text[-1]
+    show(  in_green_( n_text ))
+    
+    
+def show( *multi_print, **kargs):
+    ''' show in terminal function
+        to handle basically two things:
+            verbose levels (under development)
+            print under py2 and py3 with the same behaviour
+                e.g. to avoid : print( var1, var2) >< print var1, var2
+    '''
+    # verbose level set in each show command, else is printed
+    try:
+        v = kargs['v'] - 1
+    except KeyError:
+        v = 0
+    
+    if __verbose__[ v]:
+        
+        try:
+            end = kargs['end']
+        except KeyError:
+            end = '\n'
+        try:
+            sep = kargs['sep']
+        except KeyError:
+            sep = ' '
+        try:
+            deep = kargs['deep']
+        except KeyError:
+            deep = 0
+        #ln = len( multi_print)
+        try:
+            string = sep.join( multi_print)
+        except TypeError:
+            #if not deep:
+            #    for item in multi_print:
+            #        show( item, str(type(item)), deep = 1 )
+            #    show( 'other', type(multi_print), deep = 1 )
+            multi_print = [ str( item ) for item in multi_print]
+            string = sep.join( multi_print)
+        #string = '{}'.format( elem)
+        
+        stdout.write( string + end)
+        stdout.flush()
+        
+        
+def set_verbose( degrees ):
+    ''' four degrees of verbose deep should be more than enough'''
+    global __verbose__
+    
+    def_max = 4
+    if degrees < 1:
+        degrees = 1
+    elif degrees > def_max:
+        def_max = degrees
+    
+    __verbose__ = [ True]*degrees + [ False]*( def_max - degrees)
+    
+def test_show():
+    set_verbose( 2)
+    a = 'misc txt here\n oops new line in the same var'
+    b = "another var\n with\n multiple\n lines"
+    c = 'last var'
+    show( )
+    show( a)
+    show( a , b, c)
+    show( a , b, c, sep = '')
+    show( a , b, c, sep = '\n')
+    show( a , b, c, end = ' - ')
+    show( a , b, c, sep = ' - ', end = '\n')
+    show('-------------------')
+    
+    _in_file_ = './lib/./gui/main_gui.py'
+    act_size = 7074
+    show( ["{ '", _in_file_, "':", act_size, '}'])
+    show('-------------------')
+    show( "{ '", _in_file_, "':", act_size, '}', v = 2)
+    show('-------------------')
+    show( "{ '", _in_file_, "':", act_size, '}',sep='', v = 2)
+    
+    aux_dict = {'./lib/./gui/main_gui.py':7074,'./lib/i.py':7,'./lib/mobi.dic':1010}
+    show( aux_dict)
+    show( a , b, c, sep = ' & ', end = '\n', v = 3)
+    show('-------------------')
+
+## test
+if __name__ == '__main__':
+    #print wrg_1('')
+    #print wrg_2('')
+    #print wrg_3('')
+    #print dec_g('deco')
+    test_show()
+    pass
+
+# vim:tw=80
\ No newline at end of file
diff --git a/lib/misc/file.py b/lib/misc/file.py
index 6bb66ed..034e958 100644
--- a/lib/misc/file.py
+++ b/lib/misc/file.py
@@ -3,63 +3,83 @@
 __author__ = 'Hernan Chavez Thielemann <hchavezthiele at gmail dot com>' 
 
 
-from lib.misc.warn import wrg_1, wrg_3, pop_wrg_1 #
-from os import getcwd, walk, system 
-from os.path import join
+from lib.misc.warn import wrg_1, wrg_3, pop_wrg_1
+from lib.misc.display import show #
+from os import getcwd, walk, system, stat
+from os.path import join, sep
 from subprocess import Popen, PIPE
 
-def check_file(_in_file_, content = True, string = ''):
+
+p_sep = sep
+
+def check_file(_in_file_, content = True, string2chk = '', size = 0):
     '''Check the correctness of the given input file address'''
     
-    flag = False
+    ok_flag = False
+    show( _in_file_, v = 4)
     try:
-        _auxf = open(_in_file_,"r")
-        su = 0
-        strcheck = False
-        if content:
-            if string == '':
-                for line in _auxf:
-                    su += len(line.rstrip('\n'))
-            else:
-                for line in _auxf:
-                    line = line.rstrip('\n')
-                    if string in line:
-                        strcheck = True
-                    
-        _auxf.close()
+        with open( _in_file_, "r") as _auxf:
+            su = 0
+            strcheck = False
+            if content:
+                if string2chk == '':
+                    try:
+                        prev_line = ''
+                        for line in _auxf:
+                            su += len( line.rstrip('\n'))
+                            prev_line = line
+                    except UnicodeDecodeError:
+                        show( prev_line, line, sep = '\n', v = 3)
+                else:
+                    for line in _auxf:
+                        line = line.rstrip('\n')
+                        if string2chk in line:
+                            strcheck = True
         
-        if content and su == 0 and string == '':
+        if content and su < 2 and string2chk == '':
             pop_wrg_1(' File {} is empty -- '.format(_in_file_))
-        elif string <> '' and content and not strcheck:
-            pop_wrg_1(' Section {} not found in file {} -- '.format( string, _in_file_))
+        elif content and string2chk != '' and not strcheck:
+            aux_txt = ' Section {} not found in file {} -- '
+            pop_wrg_1( aux_txt.format( string2chk, _in_file_))
         else:
-            flag = True
+            act_size = stat( _in_file_).st_size
+            if size:
+                if act_size == size:
+                    ok_flag = True
+                else:
+                    show( '<{}:{}-{}>'.format(_in_file_, act_size, size),
+                         sep='', v =3)
+            else:
+                show( "Size of '", _in_file_, "': ", act_size, sep='', v = 3)
+                ok_flag = True
             
     except IOError:
-        pop_wrg_1(' File {} not found -- '.format(_in_file_))
-    return flag
+        pop_wrg_1(' File {} not found -- '.format( _in_file_))
+        
+    return ok_flag
 
 def check_file_list(files_list, extensions=['*']):
     '''Check for input files integrity and extension'''
-    finam=''
+    finam = ''
+    _flag = True
     try:
         for x in range(len(files_list)):
-            finam= files_list[x].split('/')[-1]
-            fi= open (files_list[x],'r')
+            finam = files_list[x].split('/')[-1]
+            fi = open (files_list[x],'r')
             fi.close()
-            ext= finam.split('.')[-1]
-            if extensions<>['*'] and not ext in extensions:
-                print extensions
-                print wrg_3(' Invalid format: < '+ext+' >')
-                return False
-        return True
+            ext = finam.split('.')[-1]
+            if extensions != ['*'] and not ext in extensions:
+                show( extensions)
+                show( wrg_3( ' Invalid format: < ' + ext + ' >') )
+                _flag = False
+                break
     except IOError:
-        if finam== '':
-            print wrg_3(" Select a file --- ")
+        if finam == '':
+            show( wrg_3( " Select a file --- "))
         else: 
-            print wrg_3(' No such file or directory: '+ finam)
-        return False
-
+            show( wrg_3(' No such file or directory: ' + finam))
+        _flag = False
+    return _flag
 def check_in_file( _file_, *args, **kwargs):
     ''' Checks if some args are in a file or not '''
     
@@ -71,9 +91,9 @@ def check_in_file( _file_, *args, **kwargs):
             _min, _max = [ int(x) for x in kwargs['slce'].split(':')]
             for k_line in indata:
                 for a in range( len( args)):
-                    #print args[a]
+                    #show( args[a]
                     line_c = k_line[ _min: _max].strip(' ')
-                    #print line_c
+                    #show( line_c
                     if args[a] == line_c:
                         _flags_[a] = 1
                         if min(_flags_):
@@ -93,7 +113,7 @@ def check_in_file( _file_, *args, **kwargs):
         else:
             for k_line in indata:
                 for a in range(len(args)):
-                    #print args[a]
+                    #show( args[a]
                     line_c = k_line.split(args[a])
                     if len(line_c)>1:
                         _flags_[a] = 1
@@ -107,7 +127,7 @@ def make_dir( _path_, _name_):
         path_dir = _path_ + _name_
     else:
         path_dir = _path_ + '/' + _name_
-    print( 'Creating folder {}'.format( path_dir))
+    show( 'Creating folder {}'.format( path_dir))
     # TODO: hadle the case in which the folder already exists
     #try:
     #    write_file( 'test.txt', content=' ', path_dir):
@@ -115,9 +135,10 @@ def make_dir( _path_, _name_):
     #    out_file.close()
     #    
     #except IOError:
+    extrad  = ''
     if len(path_dir.split(' ')) > 1:
-        path_dir = '"' + path_dir + '"'
-    system( 'mkdir ' + path_dir)
+        extrad = '"'
+    system( 'mkdir ' + extrad + path_dir + extrad)
     return path_dir + '/'
 
 def write_xfile(filename='test.txt',  content=''):
@@ -126,7 +147,7 @@ def write_xfile(filename='test.txt',  content=''):
     system('chmod +x '+filename)
 
 def move_file( _file_, _folder_):
-    system('mv ' + _file_ + ' ' + _folder_)
+    system('mv ' + _file_ + ' ' + '"' + _folder_ + '"')
 
 def write_file( filename='test.txt', content='', _folder_ = None):
     '''classic file maker'''
@@ -134,7 +155,7 @@ def write_file( filename='test.txt', content='', _folder_ = None):
     out_file.write( content)
     out_file.close()
     
-    if _folder_ <> None:
+    if _folder_ != None:
         move_file( filename, _folder_)
     
     
@@ -148,7 +169,7 @@ def write_list2file( filename, listofstrings):
                 plaintext+= strings+"\n"
             write_file( filename, plaintext)
         else:
-            print wrg_3(" Wrong format list --- "), listofstrings
+            show( wrg_3(" Wrong format list --- "), listofstrings)
 
 def write_listoflist2file( filename, listoflistsofstrings):
         '''datafile maker'''
@@ -163,11 +184,11 @@ def write_listoflist2file( filename, listoflistsofstrings):
                 listofstrings.append(row_text)
             write_list2file( filename, listofstrings)
         else:
-            print wrg_3(" Wrong format list --- ")
+            show( wrg_3(" Wrong format list --- ") )
 
 def debugger_file( char_str,  container ):
     ''' debugger'''
-    print char_str
+    show( char_str)
     write_listoflist2file( char_str+'.txt', container )
 
 def run_command(__command__):
@@ -193,12 +214,12 @@ def fileseeker( path = getcwd(), word = 'data', notw = None):
         # with remove(list_of_files[fs])
         _file_ = list_of_files[fs]
         if '/' in _file_ and word in _file_.split('/')[-1]:
-            if notw <> None and notw in _file_.split('/')[-1]:
+            if notw != None and notw in _file_.split('/')[-1]:
                 pass
             else:
                 files.append( _file_)
     if files == []:
-        print wrg_3( " No file(s) found with " + word + " criterion---")
+        show( wrg_3( " No file(s) found with " + word + " criterion---"))
         
     return files
  
diff --git a/lib/misc/geometry.py b/lib/misc/geometry.py
index 030ab00..de72736 100644
--- a/lib/misc/geometry.py
+++ b/lib/misc/geometry.py
@@ -3,6 +3,7 @@
 
 from math import pi, cos, sin, sqrt, acos
 from sys import exit
+from lib.misc.display import show
 
 def tensprod(T1, T2):
     
@@ -24,8 +25,8 @@ def tensprod(T1, T2):
     #for i in range( c):
     #    print T1[i]
     #print T1[0][0]*T2[0][0] + T1[0][2]*T2[2][0]
-    if len( T2) <> v:
-        print('Tensors length')
+    if len( T2) != v:
+        show( 'Tensors length')
         pT = [ [], [], []]
     else:
         for j in range( len( T2[0])):
@@ -69,7 +70,7 @@ def rotate( t2r, angle, axis = 'x'):
              [ sin(angle),  cos(angle), 0         ],
              [ 0         ,  0         , 1         ]]
     else:
-        print 'Error! invalid axis : {}'.format( axis)
+        show( 'Error! invalid axis : {}'.format( axis))
     
     return tensprod( R, t2r)
     
@@ -84,13 +85,13 @@ def rotate( t2r, angle, axis = 'x'):
     index = range(3)
     for i in range(3):
         Ar[i][i] = tensorvec[i]
-        print Ar
+        show( Ar)
         for j in range(3):#.remove(i):
-            if i <> j:
+            if i != j:
                 #
                 Ar[i][j] = tensorvec[k+3]
                 k += 1
-    print Ar
+    show( Ar)
     #Ar = [ [4.0465, 0, 0],[-0.73454, 2.63083, 0], [-1.45599, -0.54283, 68.36656 ]]
     #A = [ 4.0465, -0.73454, -1.45599, 0, 2.63083, -0.54283, 0, 0, 68.36656 ]
     
@@ -114,7 +115,7 @@ def rotate( t2r, angle, axis = 'x'):
     Ar = rotate( Ar, a_tor_x)
     
     for i in range( len( Ar)):
-        print Ar[i]
+        show( Ar[i])
     
     
 # vim:tw=80
diff --git a/lib/misc/version.py b/lib/misc/version.py
index 1a6bd52..4658388 100644
--- a/lib/misc/version.py
+++ b/lib/misc/version.py
@@ -1,12 +1,22 @@
 #!/usr/bin/python
 
-# pitiful file, but neats
+# pitiful file, but neats and logs
+from sys import version as python_version
 
-__curr_ver__ = '1.25 10 Aug 2020'.split()
+# 2 modify
+__curr_ver__ = '1.30 06 Jan 2021'.split()
 
-__old_vers__ = ['1.24 20 Feb 2020', '1.23 02 Oct 2019', '1.21 20 Mar 2019', '1.2 27 Jan 2019',
-                '1.12 20 Dic 2018', '1.11 20 Jul 2018', '1.1 07 May 2018','1.0 13 Apr 2018']
+# 2 update
+__old_vers__ = [ '1.25 10 Aug 2020', '1.24 20 Feb 2020', '1.23 02 Oct 2019',
+                 '1.21 20 Mar 2019', '1.2 27 Jan 2019', '1.12 20 Dic 2018',
+                 '1.11 20 Jul 2018', '1.1 07 May 2018','1.0 13 Apr 2018']
+
+# 2 rarely modify
+__copyright__ = ( 'Copyright ' + u'\u00a9' + ' 2018-' + __curr_ver__[-1] + ','
+                + '\nHernan Chavez Thielemann')
 
 __url__ = 'https://github.com/hernanchavezthielemann/GRO2LAM'
 
-__version__ = ' GRO2LAM version {} ({} {} {})'.format(*__curr_ver__)
+__version__ = ' GRO2LAM version {} ({} {} {})'.format( *__curr_ver__)
+
+__python_version__ = float( python_version[ 0:3])
\ No newline at end of file
diff --git a/lib/misc/warn.py b/lib/misc/warn.py
index 2d8bbd4..5169502 100644
--- a/lib/misc/warn.py
+++ b/lib/misc/warn.py
@@ -1,75 +1,54 @@
 #!/usr/bin/python
 #    By Hernan Chavez Thielemann
+#    this came from a blend of warn and verbose 
 
 from lib.gui.popup import message_box
+from lib.misc.display import show
 
-def ye_deco(func):
-    def func_wrapper(name):
-        return "\033[93m{0}\033[0m".format(func(name))
-    return func_wrapper
-
-#@ye_deco
 def wrg_3(*text):
-    '''  yellow warning scheme'''
+    '''  Yellow warning scheme'''
     n_text = ''
-    if text <>None:
+    if text != None:
         for line in text[:-1]:
             n_text += line +'\n'
         n_text +=text[-1]
     return '\033[93mWarning!!\033[0m -- '+ n_text
     
 def wrg_2(*text):
-    '''  warning scheme'''
+    '''  Green warning scheme'''
     n_text = ''
-    if text <>None:
+    if text != None:
         for line in text[:-1]:
             n_text += line +'\n'
         n_text +=text[-1]
     return '\033[92m'+'Warning!!'+'\033[0m -- ' + n_text
     
 def wrg_1(*text):
-    '''  warning scheme'''
+    '''  Red!! warning scheme'''
     n_text = ''
-    if text <>None:
+    if text != None:
         for line in text[:-1]:
             n_text += line +'\n'
         n_text +=text[-1]
     return '\033[91m'+'Warning!!'+'\033[0m -- ' + n_text
 
-def print_dec_g(*text):
-    '''  warning scheme'''
-    n_text = ''
-    if text <>None:
-        for line in text[:-1]:
-            n_text += line +'\n'
-        n_text +=text[-1]
-    print '\033[92m'+n_text+'\033[0m '
-    
-def green(_text_):
-    return '\033[92m'+_text_+'\033[0m '
-
 def pop_wrg_1(*_text_, **kwargs):
-    
+    ''' Warning pop-up'''
     if '_i_' in kwargs.keys():
         _i_ = kwargs['_i_']
     else:
         _i_ = -1
         
-    print wrg_1(*_text_)
+    show( wrg_1(*_text_) )
     message_box(_text_[_i_], 'Warning', icon='warning')
     
 def pop_err_1(*_text_, **kwargs):
-    
+    ''' error pop-up'''
     if '_i_' in kwargs.keys():
         _i_ = kwargs['_i_']
     else:
         _i_ = -1
         
-    print wrg_1(*_text_)
+    show( wrg_1(*_text_) )
     message_box(_text_[_i_], 'Error', icon='error')
-## test
-
-#print wrg_1('')
-#print wrg_2('')
-#print wrg_3('')
-#print dec_g('deco')
+    
\ No newline at end of file
diff --git a/setup b/setup
index 8ca9817..2d2c51b 100644
--- a/setup
+++ b/setup
@@ -1,49 +1,66 @@
 #!/usr/bin/python
+#    By Hernan Chavez Thielemann
+__author__ = 'Hernan Chavez Thielemann <hchavezthiele at gmail dot com>' 
 
+### My homemade setup
+# If you think/know that there is a better option for this please contact me.
+# I am positively sure that there is something that I am skipping here ;)
 
 from os import system, getcwd
 from sys import exit, argv
-
+from lib.misc.display import in_green_, show, set_verbose, in_red_
+from lib.misc.file import write_xfile, fileseeker, check_file, p_sep
+from lib.misc.version import __version__, __python_version__
 
 _okflag_ = True
 _exit_txt_ = ''
 _test_ = False
+_size_check_ = True
 
-if len(argv)>1 and argv[1] == '-t':
-    _test_ = True
+if len( argv) > 1:
+    if argv[1] == '-t':
+        _test_ = True
+        set_verbose( 4)
+    if argv[1] == '-f':
+        _size_check_ = False
+        set_verbose( 4)
     
 try:
-    import Tkinter
+    if __python_version__ > 3:
+        tk_install_txt = 3
+        from tkinter import Tk
+    else:
+        tk_install_txt = ''
+        from Tkinter import Tk
+        
 except ImportError:
-    print 'Package python-tk not found!'
+    show( 'Package python-tk not found!' )
     sel = ''
-    print ('You need to install the python-tk package,'
+    show( 'You need to install the python-tk package,'
            +' would you like to do it now?')
     while sel not in ['yes','no','y','n']:
         sel = raw_input('Please write "yes" or "no": ')
     if sel in ['y','yes']:
-        su =''
+        su = ''
         while su not in ['yes','no','y','n']:
             su = raw_input('Do you possess sudo rights?,'
                            +' Please answer "yes" or "no": ')
         if su in ['y','yes']:
-            system('sudo apt-get install python-tk')
+            system('sudo apt-get install python{}-tk'.format( tk_install_txt))
+            
             exit('\nPlease run again the setup\n' +
-                  '"./GRO2LAM-27ene19/python2 setup"')
+                  '"./GRO2LAM-ddmmmYY/python setup"')
         else:
             _exit_txt_ += ('Without sudo rights you cannot install tkinter.'
                            +'\n - Please contact your system administrator\n')
             _okflag_ = False
     else:
-        _exit_txt_ += ( 'In order to use GRO2LAM GUI you need to install'
-                       +' the python-tk package first\n')
+        _exit_txt_ += ( 'In order to use GRO2LAM GUI you need to install the '
+                       + 'python{}-tk package first\n'.format( tk_install_txt))
         _okflag_ = False
 
 if _okflag_:
-    from lib.misc.file import write_xfile, fileseeker, check_file
-    from lib.misc.warn import green
-    from lib.grotolam import __version__
-    print 'Checking instalation files for {}'.format(__version__)
+    show( 'Checking instalation files for {}'.format( __version__) )
     
     commands_container = ['mkdir ./lib/docs','mv README.md ./lib/docs',
                           'mv setup ./lib', 'mv LICENSE ./lib/docs',
@@ -51,78 +68,112 @@ if _okflag_:
     
     for com_txt in commands_container:
         system( com_txt)
+        
+# The voodoo call *~* >o
+_run_s_ = ( ')(rehcnual_malotorg    \n:"__niam__" == __eman__ fi\n'
+           + 'rehcnual_malotorg tropmi malotorg.bil morf\n'
+           + 'nohtyp/nib/rsu/!#', 'just to play with the reader ;)')
     
-_run_s_ = (')(rehcnual_malotorg    \n:"__niam__" == __eman__ fi\n'
-           +'rehcnual_malotorg tropmi malotorg.bil morf\n'+
-           'nohtyp/nib/rsu/!#', 'just to play with the reader ;)')
-           
-    
-_run_t_ = (')(tset_rehcnual_malotorg\ntset_rehcnual_malotorg tropmi'
+_run_t_ = ( ')(tset_rehcnual_malotorg\ntset_rehcnual_malotorg tropmi'
            +' malotorg.bil morf\nnohtyp/nib/rsu/!#', 'same here ;)')
 
 
-files = {'./grotolam.py', './__init__.py', './setup',
-
-         './gui/popup.py', './gui/__init__.py', './gui/main_gui.py',
-         './gui/script_gui.py', './gui/run_gui.py','./gui/custom_row.py',
-         './gui/conversion_gui.py', './gui/tk_lib.py', '/gui/img/README.md',
+files = { 'grotolam.py', '__init__.py', 'setup',
+          
+          'gui/img/README_img.md','gui/popup.py', 'gui/__init__.py',
+          'gui/script_gui.py', 'gui/run_gui.py','gui/custom_row.py',
+          'gui/conversion_gui.py', 'gui/tk_lib.py', 'gui/main_gui.py',
          
-          './gui/img/logo.ppm', './gui/img/gear.ppm',
-         './gui/img/logo2.ppm', './gui/img/small_logo2.ppm', 
-         './gui/img/file.ppm', './gui/img/help.ppm',
+          'gui/img/logo.ppm', 'gui/img/gear.ppm',
+          'gui/img/logo2.ppm', 'gui/img/small_logo2.ppm', 
+          'gui/img/file.ppm', 'gui/img/help.ppm',
          
-         './docs/README.md', './docs/LICENSE',
+          'docs/README.md', 'docs/LICENSE',
          
-         './misc/file.py', './misc/warn.py', './misc/version.py',
-         './misc/data.py', './misc/geometry.py', './misc/__init__.py',
+          'misc/file.py', 'misc/warn.py', 'misc/version.py', 
+          'misc/data.py', 'misc/geometry.py', 'misc/__init__.py',
+          'misc/display.py',
          
-         './handling/lammps.py',
-         './handling/__init__.py', './handling/gromacs.py'}
+          'handling/lammps.py',
+          'handling/__init__.py', 'handling/gromacs.py'}
 
+file_sizes = { '__init__.py': 45, 'file.ppm': 1861, 'setup': 6565,  
+              'README_img.md': 739, 'grotolam.py': 1407, 'logo2.ppm': 85389,
+              'logo.ppm': 153519, 'conversion_gui.py': 12417, 'file.py': 7597,
+              'data.py': 8895, 'LICENSE': 1086, 'README.md': 20765,
+              'custom_row.py': 4441, 'warn.py': 1363, 'gromacs.py': 45807,
+              'main_gui.py': 7177,'geometry.py': 3686, 'lammps.py': 69726,
+              'script_gui.py': 18621, 'tk_lib.py': 11589, 'gear.ppm': 1957,
+              'popup.py': 27086, 'version.py': 719, 'display.py': 3432,
+              'help.ppm': 1039, 'run_gui.py': 4573, 'small_logo2.ppm': 51134
+             }
+if not _size_check_:
+    file_sizes = {}
+    
 if _okflag_:
-    # This gathers all the files downstream /lib
-    act_files = fileseeker(getcwd()+'/lib', '', '.pyc')
+    # This gathers all the files downstream /lib, avoiding pyc extension
+    act_files = fileseeker( getcwd() + '/lib', '', '.pyc')
+    show( act_files, v = 4)
     
-    # so the we can compare
-    if len(files) == len(act_files):
-        print 'File check sum ok!'
-        fi_o_flag = 1
+    # so then we can compare
+    if len( files) == len( act_files):
+        show( 'File check sum ok!')
+        fi_o_flag = []
+        file_not_found = []
+        file_corrupt = []
         _eq_files_ = []
-        for _fi_ in act_files:
-                _eq_files_.append(_fi_.split('/lib/')[-1])
+        for _fi_a in act_files:
+                _eq_files_.append( _fi_a.split('/lib/')[-1])
                 #print _eq_files_[-1]
-        for _fi_ in files:
-            if _fi_.lstrip('./') in _eq_files_:
-                if check_file('./lib/'+_fi_):
-                    print ('File {} '+green('ok')).format(_fi_)
+        
+        for _fi_l in files:
+            if _fi_l in _eq_files_:
+                c_flag = True
+                file_name = _fi_l.split( p_sep )[-1]
+                if file_name.split('.')[-1] == 'ppm':
+                    c_flag = False
+                try:
+                    fsz_ = file_sizes[ file_name]
+                except KeyError:
+                    fsz_ = 0
+                
+                f_f_name = './lib/'+ _fi_l #full_file_name
+                if check_file( f_f_name , content = c_flag, size = fsz_):
+                    show( ( 'File {} ' + in_green_('ok') ).format( f_f_name) )
+                    fi_o_flag.append( 1)
                 else:
-                    fi_o_flag *= 0
+                    fi_o_flag.append( 0)
+                    file_corrupt.append( { _fi_l : fsz_} )
+                    show( ( 'File {} ' + in_red_('corrupt') ).format( f_f_name) )
             else:
-                print 'File {} not found'.format(_fi_)
-                fi_o_flag *= 0
+                show( 'File {} not found'.format( f_f_name) )
+                fi_o_flag.append( 0)
+                file_not_found.append( _fi_l)
                 
-        if fi_o_flag:
-            print 'File check done!'
+        if min( fi_o_flag):
+            show( 'File check done!')
         else:
-            _exit_txt_ += 'Warning!! some file not found!\n'
             _okflag_ = False
+            _exit_txt_ += 'Warning!! some file not found or corrupt!\n'
+            _exit_txt_ += str( file_not_found) + '/corrupt:' + str( file_corrupt)
+            
     else:
         _okflag_ = False
-        print len(files), len(act_files)
+        show(  len( files), len(act_files) )
         for fi in act_files:
-            print fi
-        _exit_txt_ += 'Warning!! there is something rare in these files\n'
+            show( fi)
+        _exit_txt_ += 'Warning!! there is a mismatch in the files\n'
         
         
         
 if not _okflag_:
     for fi in act_files:
-        print fi
-    exit(_exit_txt_)
+        show( fi, v = 4)
+    exit( _exit_txt_)
 
 if _test_:
-    write_xfile('run', _run_t_[0][::-1])
+    write_xfile( 'run', _run_t_[0][::-1])
 else:
-    write_xfile('run', _run_s_[0][::-1])
+    write_xfile( 'run', _run_s_[0][::-1])
 
-print 'Exit status OK!'
+show( 'Exit status OK!')