anchor.diff:        call write_double(6,delta_E,"delta_E")
         delta_E = dabs(delta_E)
         call save_wavefunction
  -      call ezfio_set_mrcepa0_energy(ci_energy_dressed(1))
  +      call ezfio_set_dress_zmq_energy(ci_energy_dressed(1))
         if (iteration >= n_it_dress_max) then
           exit
         endif
  
positive.code:   print *,  'E+PT2    = ', energy+pt2 
    print *,  '-----' 
  
    call ezfio_set_mrcepa0_energy_pt2(energy(1)+pt2(1))
  
  end 
  
  
negative.code:     enddo
      SOFT_TOUCH psi_coef ci_energy_dressed
      call write_double(6,ci_energy_dressed(1),"Final dress energy")
      call ezfio_set_dress_zmq_energy(ci_energy_dressed(1))
      call save_wavefunction
    else
      E_new = 0.d0
  