anchor.diff:  
     integer(kind=st), parameter :: levmax=32, submax=32
   
  -  real(kind = lg)  , dimension(levmax,submax) :: ntic
  -  real(kind = lg)  , dimension(levmax,submax) :: ntoc
  +  integer(kind = lg)  , dimension(levmax,submax) :: ntic
  +  integer(kind = lg)  , dimension(levmax,submax) :: ntoc
     real(kind = lg)  , dimension(levmax,submax) :: time
     integer(kind=st) , dimension(levmax,submax) :: calls
     character(len=32),dimension(submax)         :: subname
  
positive.code:        if (flag) then
            nbsub = nbsub + 1
            subname(nbsub)=TRIM(string)
            call cpu_time(ntic(lev,nbsub))
            time(lev,nbsub)  = 0._8
            calls(lev,nbsub) = 0
         endif
  
negative.code:        !- if yes -> cpu_time(tic)
         do ns=1, nbsub
            if (TRIM(string) == subname(ns)) then
               !call cpu_time(ntic(lev,ns))
               call system_clock(ntic(lev,ns))
               flag = .false.
               exit
            endif
  