GEANT321/peanut/frhonc.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.1  1995/10/24 10:22:01  cernlib
   6 * Geant
   7 *
   8 *
   9 #include "geant321/pilot.h"
  10 *CMZ :  3.21/02 29/03/94  15.41.46  by  S.Giani
  11 *-- Author :
  12 *$ CREATE FRHONC.FOR
  13 *COPY FRHONC
  14 *
  15 *=== frhonc ===========================================================*
  16 *
  17       FUNCTION FRHONC (R)
  18
  19 #include "geant321/dblprc.inc"
  20 #include "geant321/dimpar.inc"
  21 #include "geant321/iounit.inc"
  22 *
  23 #include "geant321/nucgeo.inc"
  24 *
  25 #include "geant321/nucstf.inc"
  26 *
  27       IF ( R .LE. RADIU0 ) THEN
  28          FRHONC = RHCORE (R)
  29       ELSE IF ( R .LE. RADIU1 ) THEN
  30          FRHONC = RHSKIN (R)
  31       ELSE IF ( R .LE. RADTOT ) THEN
  32          FRHONC = RHHALO (R)
  33       ELSE
  34          FRHONC = 0.D+00
  35       END IF
  36       RETURN
  37       END