Changeset 811 for codes/icosagcm/trunk/src

Timestamp:

02/22/19 18:09:12 (5 years ago)

Author:

adurocher

Message:

trunk+devel : fixed cross product

File:

• codes/icosagcm/trunk/src/sphere/vector.f90 (modified) (2 diffs)

codes/icosagcm/trunk/src/sphere/vector.f90

@@ -23 +23 @@
    END FUNCTION dot_product
    
-   FUNCTION cross_product(V1,V2)
+   FUNCTION cross_product(Va,Vb)
    IMPLICIT NONE
      REAL(rstd)     ::cross_product(3)
-     REAL(rstd),INTENT(IN) :: V1(3)
-     REAL(rstd),INTENT(IN) :: V2(3)
+     REAL(rstd),INTENT(IN) :: Va(3)
+     REAL(rstd),INTENT(IN) :: Vb(3)
+     REAL(rstd) :: V1(3),V2(3)
+     V1=.5*(Va+Vb)
+     V2=Vb-Va
      cross_product(1)=V1(2)*V2(3)-V1(3)*V2(2)
      cross_product(2)=V1(3)*V2(1)-V1(1)*V2(3)
@@ -33 +36 @@
    END FUNCTION cross_product
 
-   SUBROUTINE cross_product2(Va,Vb,cross_product)
+   SUBROUTINE cross_product2(V1,V2,cross_product_res)
    IMPLICIT NONE
-     REAL(rstd),INTENT(OUT):: cross_product(3)
-     REAL(rstd),INTENT(IN) :: Va(3), Vb(3)
-     REAL(rstd) :: V1(3),V2(3)
-     V1=.5*(Va+Vb)
-     V2=Va-Vb
-     cross_product(1)=V1(2)*V2(3)-V1(3)*V2(2)
-     cross_product(2)=V1(3)*V2(1)-V1(1)*V2(3)
-     cross_product(3)=V1(1)*V2(2)-V1(2)*V2(1)
+     REAL(rstd),INTENT(OUT):: cross_product_res(3)
+     REAL(rstd),INTENT(IN) :: V1(3), V2(3)
+     cross_product_res=cross_product(V1,V2)
    END SUBROUTINE cross_product2