ScaleRotateVector.F90

Go to the documentation of this file.

! *********************************************************************
! * Rocstar Simulation Suite                                          *
! * Copyright@2015, Illinois Rocstar LLC. All rights reserved.        *
! *                                                                   *
! * Illinois Rocstar LLC                                              *
! * Champaign, IL                                                     *
! * www.illinoisrocstar.com                                           *
! * sales@illinoisrocstar.com                                         *
! *                                                                   *
! * License: See LICENSE file in top level of distribution package or *
! * http://opensource.org/licenses/NCSA                               *
! *********************************************************************
! *********************************************************************
! * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND,   *
! * EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES   *
! * OF MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND          *
! * NONINFRINGEMENT.  IN NO EVENT SHALL THE CONTRIBUTORS OR           *
! * COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER       *
! * LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE,   *
! * Arising FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE    *
! * USE OR OTHER DEALINGS WITH THE SOFTWARE.                          *
! *********************************************************************
!******************************************************************************
!
! Purpose: Scale and rotate vector.
!
! Description: None.
!
! Input: vect = vector to scale and rotate.
!
! Output: vect = scaled and rotated.
!
! Notes: None.
!
!******************************************************************************
!
! $Id: ScaleRotateVector.F90,v 1.3 2008/12/06 08:44:10 mtcampbe Exp $
!
! Copyright: (c) 2002 by the University of Illinois
!
!******************************************************************************

SUBROUTINE scalerotatevector( global,vect )

  USE moddatatypes
  USE modglobal, ONLY : t_global
  USE moderror
  USE modparameters

  IMPLICIT NONE
    
! ... parameters
  REAL(RFREAL), POINTER :: vect(:,:)

  TYPE(t_global), POINTER :: global
  
! ... local variables
  CHARACTER(CHRLEN) :: rcsidentstring

  INTEGER :: ic, ncomp

  REAL(RFREAL) :: anglex, angley, anglez, cx, cy, cz, sx, sy, sz, x, y, z

!******************************************************************************

  rcsidentstring = '$RCSfile: ScaleRotateVector.F90,v $ $Revision: 1.3 $'

  IF ( global%verbLevel > verbose_none ) THEN 
    WRITE(stdout,'(A,1X,A)') solver_name,'Scaling and rotating vector...'
    WRITE(stdout,'(A,3X,A,1X,E15.6)') solver_name,'X-scale:',global%scaleX
    WRITE(stdout,'(A,3X,A,1X,E15.6)') solver_name,'Y-scale:',global%scaleY
    WRITE(stdout,'(A,3X,A,1X,E15.6)') solver_name,'Z-scale:',global%scaleZ
    WRITE(stdout,'(A,3X,A,1X,E15.6)') solver_name,'X-angle:',global%angleX
    WRITE(stdout,'(A,3X,A,1X,E15.6)') solver_name,'Y-angle:',global%angleY
    WRITE(stdout,'(A,3X,A,1X,E15.6)') solver_name,'Z-angle:',global%angleZ
  END IF ! global%verbLevel

  ncomp = SIZE(vect,dim=2)

! scale vector

  DO ic = 1,ncomp
    vect(1,ic) = global%scaleX*vect(1,ic)
    vect(2,ic) = global%scaleY*vect(2,ic)
    vect(3,ic) = global%scaleZ*vect(3,ic)
  END DO ! ic

! rotate vector

  anglex = global%angleX*global%rad
  angley = global%angleY*global%rad
  anglez = global%angleZ*global%rad   

  cx = cos(anglex)
  sx = sin(anglex)
  cy = cos(angley)
  sy = sin(angley)
  cz = cos(anglez)
  sz = sin(anglez)
  
  DO ic = 1,ncomp
    x = vect(1,ic)
    y = vect(2,ic)
    z = vect(3,ic)

    vect(1,ic) = cy*cz*x - (sx*sy*cz + cx*sz)*y - (cx*sy*cz - sx*sz)*z
    vect(2,ic) = cy*sz*x - (sx*sy*sz - cx*cz)*y - (cx*sy*sz + sx*cz)*z
    vect(3,ic) =    sy*x + (           sx*cy)*y + (           cx*cy)*z
  END DO ! ic
  
  IF ( global%verbLevel > verbose_none ) THEN
    WRITE(stdout,'(A,1X,A)') solver_name,'Scaling and rotating vector done.'
  END IF ! global%verbLevel
  
END SUBROUTINE scalerotatevector

!******************************************************************************
!
! RCS Revision history:
!
! $Log: ScaleRotateVector.F90,v $
! Revision 1.3  2008/12/06 08:44:10  mtcampbe
! Updated license.
!
! Revision 1.2  2008/11/19 22:17:23  mtcampbe
! Added Illinois Open Source License/Copyright
!
! Revision 1.1  2004/12/01 16:51:21  haselbac
! Initial revision after changing case
!
! Revision 1.5  2004/01/29 22:52:48  haselbac
! Clean up
!
! Revision 1.4  2003/02/01 00:28:20  haselbac
! Some clean up, added diagnostic output
!
! Revision 1.3  2002/09/05 17:40:20  jblazek
! Variable global moved into regions().
!
! Revision 1.2  2002/07/16 21:34:37  jblazek
! Prefixed screen output with SOLVER_NAME.
!
! Revision 1.1  2002/05/07 18:51:27  haselbac
! Initial revision
!
!******************************************************************************