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.

_______________________________________________
Mailing list: https://launchpad.net/~yade-users
Post to     : [email protected]
Unsubscribe : https://launchpad.net/~yade-users
More help   : https://help.launchpad.net/ListHelp

Reply via email to