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