root / src / sinangle.f90 @ 3
Historique | Voir | Annoter | Télécharger (469 octet)
| 1 |
FUNCTION SinAngle(v1x,v1y,v1z,norm1,v2x,v2y,v2z,norm2) |
|---|---|
| 2 |
|
| 3 |
use Path_module, only : Pi,KINT, KREAL |
| 4 |
|
| 5 |
IMPLICIT NONE |
| 6 |
|
| 7 |
real(KREAL) :: v1x,v1y,v1z,norm1 |
| 8 |
real(KREAL) :: v2x,v2y,v2z,norm2 |
| 9 |
real(KREAL) :: SinAngle |
| 10 |
|
| 11 |
if (norm1<=1d-6) STOP "ERROR in SinAngle, norm1<=1e-6" |
| 12 |
if (norm2<=1d-6) STOP "ERROR in SinAngle, norm2<=1e-6" |
| 13 |
Sinangle=sqrt(1.-((v1x*v2x+v1y*v2y+v1z*V2z)/(norm1*norm2))**2) |
| 14 |
|
| 15 |
|
| 16 |
END FUNCTION SINANGLE |