accel1.F File Reference

#include "implicit_f.inc"
#include "com08_c.inc"
#include "param_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	accel1 (a, ff, a2, a1, a0, as, vs, skew)

Function/Subroutine Documentation

accel1()

subroutine accel1	(		a,
			ff,
			a2,
			a1,
			a0,
			as,
			vs,
			skew )

Definition at line 28 of file accel1.F.

C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "com08_c.inc"
#include      "param_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      my_real a(3),a0(3,2),a1(3,2),a2(3,2),as(3),vs(3),skew(lskew),ff
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER J
      my_real pi1,pi8,pi38,spi8,spi38,c0,c1,c2,c3,c4,c5,c6,c7,c8,c9,x1,x2,x3,y1,y2,y3,z1,z2,z3,d,dd,d2,dp,e,g,f
 
      as(1) = a(1)*skew(1) + a(2)*skew(2) + a(3)*skew(3)
      as(2) = a(1)*skew(4) + a(2)*skew(5) + a(3)*skew(6)
      as(3) = a(1)*skew(7) + a(2)*skew(8) + a(3)*skew(9)
      vs(1) = vs(1)+as(1)*dt12
      vs(2) = vs(2)+as(2)*dt12
      vs(3) = vs(3)+as(3)*dt12
      IF(ff==zero)RETURN
C
      f = min(ff,zep4/dt2)
C-----------------------------------------------
C INITIALIALISATION DES COEFFICIENTS DU FILTRE
C-----------------------------------------------
      pi1   = two*atan2(one,zero)
      pi8  = pi1*one_over_8
      pi38 = three*pi8
      spi8  = sin(pi8)
      spi38 = sin(pi38)
 
      d  = tan(pi1*f*dt2)
 
      dd = d*d
      d2 = two*d
      dp = one + dd
      e  = d2*spi8
      g  = e + dp
      g  = one/g
 
      c0 = dd * g
      c1 = two* c0
      c2 = c0
      c3 = two * g - c1
      c4 = (e - dp) * g
 
      e  = d2*spi38
      g  = e + dp
      g  = one/g
 
      c5 = dd * g
      c6 = two * c5
      c7 = c5
      c8 = two * g - c6
      c9 = (e - dp) * g
C-----------------------------------------------
C FILTERING
C-----------------------------------------------
      DO j=1,3
        x1 = a0(j,2)
        x2 = a0(j,1)
        x3 = a(j)
        y1 = a1(j,2)
        y2 = a1(j,1)
        y3 = c0 * x3 + c1 * x2 + c2 * x1 + c3 * y2 + c4 * y1
        z1 = a2(j,2)
        z2 = a2(j,1)
        z3 = c5 * y3 + c6 * y2 + c7 * y1 + c8 * z2 + c9 * z1
 
        a0(j,2) = x2
        a0(j,1) = x3
        a1(j,2) = y2
        a1(j,1) = y3
        a2(j,2) = z2
        a2(j,1) = z3
      ENDDO
 
      as(1) = a2(1,1)*skew(1) + a2(2,1)*skew(2) + a2(3,1)*skew(3)
      as(2) = a2(1,1)*skew(4) + a2(2,1)*skew(5) + a2(3,1)*skew(6)
      as(3) = a2(1,1)*skew(7) + a2(2,1)*skew(8) + a2(3,1)*skew(9)
 
      RETURN