yade-users team mailing list archive

[Yousef Golabchi <question689817@xxxxxxxxxxxxxxxxxxxxx>]

Question #689817 on Yade changed:
https://answers.launchpad.net/yade/+question/689817

    Status: Answered => Solved

Yousef Golabchi confirmed that the question is solved:
Hi Luc,

Thanks a lot. The code works fine and I am really appreciated you help
me with that. It just needed a tiny change in the lines.

#####

 def recordInteractionStatus():
	global intOnFracPlane
	intOnFracPlane=[]
	for i in O.interactions:
		if i.phys.isOnJoint:
			intOnFracPlane.append(i)
	inFile=open('interactions_timePositionsAperturePressure_'+str(O.iter)+'_'+str(round(O.time,3)),'a')
	for i in intOnFracPlane:
		inFile.write( str(i.geom.contactPoint[0]) + '\t' + str(i.geom.contactPoint[1]) + '\t' + str(i.geom.contactPoint[2]) + '\t' + str(i.phys.crackJointAperture) + '\t' + str(flow.getPorePressure((i.geom.contactPoint[0],i.geom.contactPoint[1],i.geom.contactPoint[2]))) + '\n')
	inFile.close()

#####

Best regards,
Yousef

-- 
You received this question notification because your team yade-users is
an answer contact for Yade.