Statistiques
| Révision :

root / src / sinangle.f90 @ 7

Historique | Voir | Annoter | Télécharger (460 octet)

1
     FUNCTION SinAngle(v1x,v1y,v1z,norm1,v2x,v2y,v2z,norm2)
2

    
3
       use Path_module, only :  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