Projecting distance for Unit Vector

[andykol]

Hello,
I am calculating Unit vector for non orthogonal Cartesian grid. To calculate Unit Vector I need to project eg. y and x distance of east face of cell on i and j axis resp. I am using following formula-
Unit vector= Ue
Surface area= Se=sqrt(dx2+dy2)

Ue=(dy/se)i-(dx/se)j

Can anybody please tell me how I can calculate projection in Fortran?

 

[mikrom]

Can anybody please tell me how I can calculate projection in Fortran?

IMHO, the same way as you calculate it on the paper.
But how to calculate projection is topic from Linear Algebra - google:

http://en.wikipedia.org/wiki/Vector_projection

http://en.wikipedia.org/wiki/Gram-Schmidt_process

etc...

 

[gummibaer]

andy,

sorry, I do not understand your problem.

You have x and y to be the coordinates in your non-orthogonal system with i,j axis.

Then x and y are the projections of your vector on these axis. And you have your unit vector as given in your formula.

Or could you post just a small sketch of what you are after ?

Norbert