docs: def plotprofMulti(self, ini, end, delta, what_specie, xlim1, xlim2,
                        ylim1, ylim2, symbol=None):
  
          '''
          create a movie with mass fractions vs mass coordinate between
          xlim1 and xlim2, ylim1 and ylim2. Only works with instances of
          se.
  
          Parameters
          ----------
          ini : integer
              Initial model i.e. cycle.
          end : integer
              Final model i.e. cycle.
          delta : integer
              Sparsity factor of the frames.
          what_specie : list
              Array with species in the plot.
          xlim1, xlim2 : integer or float
              Mass coordinate range.
          ylim1, ylim2 : integer or float
              Mass fraction coordinate range.
          symbol : list, optional
              Array indicating which symbol you want to use. Must be of
              the same len of what_specie array.  The default is None.
  
          '''
          plotType=self._classTest()
          if plotType=='se':
              for i in range(ini,end+1,delta):
                  step = int(i)
                  #print step
                  if symbol==None:
                      symbol_dummy = '-'
                      for j in range(len(what_specie)):
                          self.plot_prof_1(step,what_specie[j],xlim1,xlim2,ylim1,ylim2,symbol_dummy)
                  else:
                      for j in range(len(what_specie)):
                          symbol_dummy = symbol[j]
                          self.plot_prof_1(step,what_specie[j],xlim1,xlim2,ylim1,ylim2,symbol_dummy)
  
                  #
                  filename = str('%03d' % step)+'_test.png'
                  pl.savefig(filename, dpi=400)
                  print('wrote file ', filename)
                  #
                  pl.clf()
  
          else:
              print('This method is not supported for '+str(self.__class__))
              return