!elte(iout,*) "rms_nacc before contact"
call contact(.false.,ncont,icont,co)
frac=contact_fract(ncont,ncont_ref,icont,icont_ref)
!elte(iout,*) "rms_nacc before contact"
call contact(.false.,ncont,icont,co)
frac=contact_fract(ncont,ncont_ref,icont,icont_ref)