anchor.diff:  
       write(*,*) "\x1B[32mrun_particle:\x1B[0m end particle calculation."
   
  +#if MPI
  +    call MPI_FINALIZE ( ierr )
  +#endif
  +
   end program
  
positive.code:         ! load parameters from namelist
          call surface_physics_par_load(surface%par,trim(nml_file))
          ! open file for CRMSD output
          write(output,"(A,I3.3,A)") trim(prefix), i,  ".out"
          open(2,file=trim(output),form='formatted')
  
  
          do k=1,nloop ! re-iterate 'nloop' times
  
negative.code:     allocate(state%melt(ntime))
      allocate(state%acc(ntime))
  
  #if MPI
      do i = i0+my_id,i1,num_procs
  #else
      do i = i0, i1
  #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"
  