read(12,*) status11 read(12,*) PC_dispersion read(12,*) polar_flag read(12,*) analytic_flag do i=1, ncp read(12,*) v00_tmp(i) read(12,*) m_tmp(i) read(12,*) u0_tmp(i,i) read(12,*) mw_tmp(i) read(12,*) e_ii_tmp(i) read(12,*) N_sites_tmp(i) !these are the polar parameters. Set them to zero if no polarity read(12,*) q_tmp(i) read(12,*) xq_tmp(i) read(12,*) mu_tmp(i) read(12,*) xd_tmp(i) end do write(6,*) ' ' write(6,*) 'Enter the association epsilon (Well depth) interaction parameters below' do j = 1, ncp do i = 1, N_sites_tmp(j) do k = 1, ncp do l = 1, N_sites_tmp(k) !write(6,*) 'interaction between site',i,' on molecule', j,' and' !write(6,*) 'site',l, ' on molecule', k read(12,*) fa_tmp(j,k,i,l,2) end do end do end do end do write(6,*) ' ' write(6,*) 'Enter the association kappa (Well width) interaction parameters below' do j = 1, ncp do i = 1, N_sites_tmp(j) do k = 1, ncp do l = 1, N_sites_tmp(k) !write(6,*) 'interaction between site',i,' on molecule', j, ' and' !write(6,*) 'site',l, ' on molecule', k read(12,*) kappa_tmp(j,k,i,l) end do end do end do end do !write(6,*) ' ' !write(6,*) 'Enter the interaction parameters Kij below' do i = 1, ncp do j = i+1, ncp !write(6,*) 'For i =',i, 'j=',j read(12,*) kij_tmp(i,j) kij_tmp(j,i) = kij_tmp(i,j) end do end do