// We create a integer array SizeX=7, SizeY=5
  TArray<int_4> ia(7,5);
  // We fill it with a sequence of numbers starting at 10., with step = 2.
  ia = RegularSequence(10., 2.);
  cout << " ----- Array IA = \n " << ia << endl;
  TArray<int_4> ib = ia(Range(0,3), Range(3,4), Range(0));
  cout << " ----- Array IB IA(Range(0,3), Range(3,2)) = \n" << ib << endl;
  // sub array extraction, Range(2,4) : starting position=2 , End=4
  TArray<int_4> ic = ia(Range(2,3),Range(1,3),Range(0));
  cout << " ----- Array IC IA(Range(2,3),Range(1,3)) = \n " << ic << endl;
  // we set the sub-array to zero, this should reflect in the original array
  // sub-arrays share their data with parent array
  ic = 0;
  cout << " ----- Array IC Apres (=0) = \n " << ic << endl;
  cout << " ----- Array IB Apres IC=0 = \n " << ib << endl;
  cout << " ----- Array IA Apres IC=0 = \n " << ia << endl;
  cout << " >>>>>> Writing in  arrt.ppf <<<<<<< " << endl;
  POutPersist pos("arrt.ppf");
  // We write  the three arrays in the stream
  pos << ia << ib << ic;