diff --git a/schainpy/scripts/schain_paralel.py b/schainpy/scripts/schain_paralel.py
new file mode 100644
index 0000000..d1e2e98
--- /dev/null
+++ b/schainpy/scripts/schain_paralel.py
@@ -0,0 +1,93 @@
+from mpi4py import MPI
+import datetime
+import os, sys
+#import timeit
+
+path = os.path.split(os.getcwd())[0]
+sys.path.append(path)
+
+from controller import *
+
+def conversion(x1,x2):
+    a=[x1,x2]
+    for x in a:
+        m,s = divmod(x,60)
+        h,m = divmod(m,60)
+        if x==x1:
+            startime= str("%02d:%02d:%02d" % (h, m, s))
+        if x==x2:
+            endtime =str("%02d:%02d:%02d" % (h, m, s))
+    return startime,endtime
+
+
+
+def loop(startime,endtime,rank):
+    desc = "HF_EXAMPLE"+str(rank)
+    path= "/home/alex/Documents/hysell_data/pdata/sp1_f0"
+    figpath= "/home/alex/Pictures/pdata_plot"+str(rank)
+
+    filename = "hf_test"+str(rank)+".xml"
+    
+    controllerObj = Project()
+    
+    controllerObj.setup(id = '191', name='test01'+str(rank), description=desc)   
+
+    readUnitConfObj = controllerObj.addReadUnit(datatype='SpectraReader',
+                                                 path=path,
+                                                 startDate='2015/09/26',
+                                                 endDate='2015/09/26',
+                                                 startTime=startime,
+                                                 endTime=endtime,
+                                                 online=0,
+                                                 #set=1426485881,
+                                                 delay=10,
+                                                 walk=1
+                                                 #timezone=-5*3600
+                                                 )
+
+    #opObj11 = readUnitConfObj.addOperation(name='printNumberOfBlock')
+     
+    procUnitConfObj1 = controllerObj.addProcUnit(datatype='SpectraProc', inputId=readUnitConfObj.getId())
+    
+    opObj11 = procUnitConfObj1.addOperation(name='SpectraPlot', optype='other')
+    opObj11.addParameter(name='id', value='1000', format='int')
+    opObj11.addParameter(name='wintitle', value='HF_Jicamarca_Spc', format='str')
+    #opObj11.addParameter(name='channelList', value='0', format='intlist') 
+    opObj11.addParameter(name='zmin', value='-120', format='float')
+    opObj11.addParameter(name='zmax', value='-70', format='float')
+    opObj11.addParameter(name='save', value='1', format='int')
+    opObj11.addParameter(name='figpath', value=figpath, format='str')
+    
+
+    print "Escribiendo el archivo XML"
+    controllerObj.writeXml(filename)
+    print "Leyendo el archivo XML"
+    controllerObj.readXml(filename)
+    
+    controllerObj.createObjects()
+    controllerObj.connectObjects()
+    
+    #timeit.timeit('controllerObj.run()', number=2)
+    
+    controllerObj.run()
+
+    
+
+def parallel():
+    
+    comm = MPI.COMM_WORLD
+    rank = comm.Get_rank()
+    size = comm.Get_size()  
+    totalStartTime = time.time()
+    print "Hello world from process %d/%d"%(rank,size)
+    # First just for one day :D!
+    num_hours = 4/size
+    time1,time2 = rank*num_hours*3600,(rank+1)*num_hours*3600-60
+    #print time1,time2
+    startime,endtime =conversion(time1,time2)
+    print startime,endtime
+    loop(startime,endtime,rank)
+    print "Total time %f seconds" %(time.time() -totalStartTime)
+
+if __name__=='__main__':
+    parallel()
\ No newline at end of file