GEANT321/cgpack/cgsaan.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.1  1995/10/24 10:19:44  cernlib
   6 * Geant
   7 *
   8 *
   9 #include "geant321/pilot.h"
  10 *CMZ :  3.21/02 29/03/94  15.41.32  by  S.Giani
  11 *-- Author :
  12       SUBROUTINE CGSAAN(A1,A2,NA,NB,IREP)
  13 ************************************************************************
  14 *                                                                      *
  15 *     Name: CGSAAN                                                     *
  16 *     Author: E. Chernyaev                       Date:    01.02.89     *
  17 *                                                Revised:              *
  18 *                                                                      *
  19 *     Function: Set parameters of rotation                             *
  20 *                                                                      *
  21 *     References: none                                                 *
  22 *                                                                      *
  23 *     Input: A1 - intitial angle                                       *
  24 *            A2 - end angle                                            *
  25 *            NA - 1-st number of steps                                 *
  26 *            NB - 2-nd number of steps                                 *
  27 *                                                                      *
  28 *     Output: IREP - reply (-1 if error in parameters)                 *
  29 *                                                                      *
  30 *     Errors: none                                                     *
  31 *                                                                      *
  32 ************************************************************************
  33 #include "geant321/cgcaan.inc"
  34       SAVE IFIRST
  35       DATA      IFIRST/1/
  36 *-
  37       IF (IFIRST .NE. 0)                RAD = ATAN(1.0D+00)*4./180.
  38       IFIRST = 0
  39       IREP   =-1
  40       IF (NA .LE. 0)                    GOTO 999
  41       IF (NB .LE. 0)                    GOTO 999
  42       IF (A1 .GE. A2)                   GOTO 999
  43       DA     = A2 - A1
  44       IF (DA .GT. 360.)                 GOTO 999
  45       IF (DA .EQ. 360.)                 IFULL = 1
  46       IF (DA .LT. 360.)                 IFULL = 0
  47       NASTP(1) = NA
  48       NASTP(2) = NB
  49       COSI     = COS(A1*RAD)
  50       SINI     = SIN(A1*RAD)
  51       COSDA(1) = COS(DA*RAD/NA)
  52       SINDA(1) = SIN(DA*RAD/NA)
  53       COSDA(2) = COS(DA*RAD/NB)
  54       SINDA(2) = SIN(DA*RAD/NB)
  55       IREP   = 0
  56   999 RETURN
  57       END