anchor.diff:  
       write(*,*) "\x1B[32mrun_particle:\x1B[0m end particle calculation."
   
  +#if MPI
  +    call MPI_FINALIZE ( ierr )
  +#endif
  +
   end program
  
positive.code:         cost_smb = calculate_crmsd(state%smb,vali%smb,ntime)
          cost_alb = calculate_crmsd(state%alb,vali%alb,ntime)
          cost = dsqrt(cost_alb**2+cost_stt**2+cost_smb**2)
          write(2,*) cost
          close(2)
      end do
          
  
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"
  