I need to get the area vector at each face on my boundary domain for which i am using F_AREA(a,face,t) where a:- a[ND_ND] this is where the value gets stored Face:- is face_t type variable. t:- face thread.

Now ND_ND would be 3 for 3D , I am printing out a[0] , a[1], a[2], a[3] would the area vector be a[1] a[2] a[3] ? or would it be a[0] a[1] a[2] ? or something else !

Thanks, Ajay

a[0] a[1] and a[2] give you the face area in the x, y and z directions, respectivly. Be careful about your signs - read the manual for more info.