MINICERN/mathlib/gen/c/legz.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.1  1996/04/01 15:01:59  mclareni
   6 * Mathlib gen
   7 *
   8 *
   9 #include "gen/pilot.h"
  10       SUBROUTINE LEGZ
  11 #include "legbl.inc"
  12      1,F1,F2,ZV
  13      2,CLOGOK,RGAM
  14       GR=CLOGOK(2.0*ZZ,NFRIG,2)
  15       A=1.5
  16       B=1.0
  17       ZV=RGAM(VV,A,B)*EXP(VV*GR)
  18       A=VV/2.0+1.0
  19       B=VV/2.0+0.5
  20       C=VV+1.5
  21       ACCC=ACC/100.
  22       CALL HYPGM(A,B,C,Z2,F1,ACCC,NCVG)
  23       F1=F1/(2.0*ZZ*ZV)
  24       QQ=PISR*F1
  25       IF(ABS(CVV).GT.0.001) GOTO 9
  26 C   TRDZ EXPECTS GR=CLOGOK(2.0*ZZ,NFRIG,2) BUT DESTROYS CONTENTS
  27       CALL TRDZ
  28       GOTO 80
  29     9 A=-VV/2.0
  30       B=(1.0-VV)/2.0
  31       C=0.5-VV
  32       CALL HYPGM(A,B,C,Z2,F2,ACCC,NCV )
  33       NCVG=NCVG+2*NCV
  34       F2=F2*ZV/(VV+0.5)
  35       PP=(F1*SVV/CVV+F2)/PISR
  36    80 RETURN
  37       END