anchor.diff:      allocate(state%melt(ntime))
       allocate(state%acc(ntime))
   
  +#if MPI
  +    do i = i0+my_id,i1,num_procs
  +#else
       do i = i0, i1
  -        write(nml_file,"(A,I3.3,A)") trim(prefix), i, ".nml"
  +#endif
  +        write(nml_file,"(A,I6.6,A)") trim(prefix), i, ".nml"
           inquire(file=nml_file, exist=file_exists)
           if (.not. file_exists) then
               write(*,*) "file '", trim(nml_file), "' not found, exiting loop"
  
positive.code: 
      write(*,*) "\x1B[32mrun_particle:\x1B[0m end particle calculation."
  
  end program
  
negative.code:         ! load parameters from namelist
          call surface_physics_par_load(surface%par,trim(nml_file))
          ! open file for CRMSD output
          write(output,"(A,I6.6,A)") trim(prefix), i,  ".out"
          open(2,file=trim(output),form='formatted')
          if (mod(i,100) == 0) write(*,*) "\x1B[32mrun_particle:\x1B[0m Load file ", trim(nml_file)
  
  
          do k=1,nloop ! re-iterate 'nloop' times
  