Esope source files

   1 : C CMP1DE    SOURCE    CHAT      05/01/12    22:11:00     5004           
   2 : CCC
   3 : C **********************************************************************
   4 : CCC
   5 :       SUBROUTINE CARAC_MATE_POWDER1_DENSI(XMAT,XDENSI)
   6 :       IMPLICIT INTEGER(I-N)
   7 :       real*8    xmat(*),xdensi,eta
   8 :       real*8            you,xnu,sigy0,siginf,kiso,velo,cpar,mpar
   9 :       common /elasdata/ you,xnu
  10 :       common /miehdata/         sigy0,kiso,siginf,velo,cpar,mpar
  11 :       real*8                  kap,xmu,xmu1,alp1,xmu2,alp2,xmu3,alp3
  12 :       logical                 ogden
  13 :       common /elasdata_ogden/ kap,xmu,xmu1,alp1,xmu2,alp2,xmu3,alp3,
  14 :      .                        ogden
  15 :       real*8               nn1,nn2,eta0,aa1,aa2
  16 :       common /poder1_data/ nn1,nn2,eta0,aa1,aa2
  17 : 
  18 :       ogden  =.false.
  19 :       you    = xmat( 1)
  20 :       xnu    = xmat( 2)
  21 :       sigy0  = xmat( 5)
  22 :       nn1    = xmat( 6)
  23 :       nn2    = xmat( 7)
  24 :       eta0   = xmat( 8)
  25 :       eta    = xdensi*eta0
  26 :       aa1    = (1.D0 - eta**2)/(2.D0 + eta**2)
  27 :       if (aa1.lt.1.D-30) aa1=1.D-30
  28 :       aa1    = aa1**nn1
  29 :       if (aa1.lt.1.D-5) aa1=1.D-5
  30 :       aa2    = abs(eta - 0.98*eta0)/(1.D0 - 0.98*eta0)
  31 :       if (eta.lt.eta0)  aa2=0.02D0*eta0/(1.D0 - 0.98D0*eta0)
  32 :       aa2    = aa2**nn2
  33 : 
  34 :       siginf = 0.D0
  35 :       kiso   = 0.D0
  36 :       velo   = 0.D0
  37 :       cpar   = 0.D0
  38 :       mpar   = 1.D0
  39 :       return
  40 :       end
  41 :  
  42 :  
  43 :