MINICERN/mathlib/gen/x/strag.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.1  1996/04/01 15:02:58  mclareni
   6 * Mathlib gen
   7 *
   8 *
   9 #include "gen/pilot.h"
  10       FUNCTION STRAG(PINV,EM)
  11       DIMENSION ENT(15),TAB(15)
  12       DATA  ENT/.0,.001,.00189,.00355,.00675,.0128,.0282,.0460,.0875,.16
  13      +   5,.312,.593,1.012,5.,100000./
  14       DATA  TAB/.1,.01573,.01445,.0134,.0124,.0117,.0112,.0107,.0102,.00
  15      +   97,.0092,.0087,.0082,0.,0./
  16       BI     =SQRT(1.+1./(PINV*EM)**2)-1.
  17       DO 1 I=2,15
  18       IF(ENT(I).GT.BI)    GO TO 2
  19     1 CONTINUE
  20     2 A=(BI-ENT(I-1))*(TAB(I)-TAB(I-1))/(ENT(I)-ENT(I-1))+TAB(I-1)
  21       STRAG=SQRT(.938213/EM)*A
  22       RETURN
  23       END