anchor.diff:    double precision t !< current time
     double precision dt !< time step size
     double precision tFinal !< final time
  -  integer iStep !< current time step number
  -  integer nStep !< total number of time steps
  +  integer iWrite !< index of result output snapshot
   end module
   
  -!> foeuler main program
  -program foeuler
  +!> fons main program
  +program fons
     use moduleBasicDataStruct
     use moduleFileIO
     use moduleGrid
  
positive.code: !----------------------------------------------------------------------------- best with 100 columns
  
  !> global variables for foeuler
  module gVarEuler
    use moduleGrid
    public
    
  
negative.code:   use moduleInterpolation
    use moduleFVMGrad
    use moduleFVMConvect
    use moduleCLIO
    use gVarNS
    double precision,allocatable::rhoNode(:) !< density at node
    double precision,allocatable::uBlock(:,:) !< velocity at block
    double precision,allocatable::uIntf(:,:) !< velocity at interface
  