Changeset 745

Timestamp:

10/03/18 23:12:28 (6 years ago)

Author:

dubos

Message:

devel : fix accuracy issue with cross_product2 => -xHost now conserves mass on Irene

File:

• codes/icosagcm/devel/src/sphere/vector.f90 (modified) (1 diff)

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

@@ -28 +28 @@
      REAL(rstd),INTENT(IN) :: V1(3)
      REAL(rstd),INTENT(IN) :: V2(3)
-    
+     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)     
+   END FUNCTION cross_product
+
+   SUBROUTINE cross_product2(Va,Vb,cross_product)
+   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)
-     
-    END FUNCTION cross_product
-
-
-   SUBROUTINE cross_product2(V1,V2,cross_product)
-   IMPLICIT NONE
-     REAL(rstd),INTENT(OUT):: cross_product(3)
-     REAL(rstd),INTENT(IN) :: V1(3)
-     REAL(rstd),INTENT(IN) :: V2(3)
-    
-     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)
-     
-    END SUBROUTINE cross_product2
+   END SUBROUTINE cross_product2
 
     FUNCTION arc(lon,lat, lonc,latc)