anchor.diff:  #ifdef FLUKA
     use mod_fluka
   #endif
  -
  +#ifdef HDF5
  +  use hdf5_output
  +#endif
   #ifdef ROOT
     use iso_c_binding
     use root_output
  
positive.code:           if(plost(jj).ne.0) then
  #ifdef FLUKA
              if( fluka_uid(j).eq.plost(jj).or. fluka_gen(j).eq.plost(jj) ) then
  #else
              if ( (     do_coll .and. (  ipart(j) .eq. plost(jj) )) .or. &
                   (.not.do_coll .and. ( nlostp(j) .eq. plost(jj) ))       ) then
  #endif
                lparID=.true.
              end if
  
              jjx=jj+1 !points to the last zero
            end if
  
negative.code:         this_name = trim(adjustl(bez(ix))) // C_NULL_CHAR
  #if defined(FLUKA)
          call ApertureCheckWriteLossParticleF(turn, i, ix, this_name, len_trim(this_name), slos, &
            fluka_uid(j), fluka_gen(j), fluka_weight(j), &
            xlos(1)*c1m3, ylos(1)*c1m3, xlos(2)*c1m3, ylos(2)*c1m3, ejfvlos*c1m3, (ejvlos-e0)*c1e6, &
            -c1m3 * (sigmvlos/clight) * (e0/e0f), naalos, nzzlos)
  #else
          call ApertureCheckWriteLossParticle(turn, i, ix, this_name, len_trim(this_name), slos, plost(j),&
            xlos(1)*c1m3, ylos(1)*c1m3, xlos(2)*c1m3, ylos(2)*c1m3, ejfvlos*c1m3, (ejvlos-e0)*c1e6, &
            -c1m3 * (sigmvlos/clight) * (e0/e0f), naalos, nzzlos)
  #endif
        end if
  #endif
  