Esope source files

   1 : C CAMARH    SOURCE    CHAT      05/01/12    21:48:13     5004           
   2 : c MODELO RHMC
   3 : CCC
   4 : C **********************************************************************
   5 : CCC
   6 :       SUBROUTINE CARAC_MATE_RHMC (XMAT)
   7 :       IMPLICIT INTEGER(I-N)
   8 :       real*8    xmat(*)
   9 :       real*8    you,xnu,coh,phi,ar,tmax,aap,aan,bbp,bbn
  10 :       common /elasdata/ you,xnu
  11 :       common /rhmcdata/ coh,phi,ar,tmax,aap,aan,bbp,bbn
  12 :       real*8                  kap,xmu,xmu1,alp1,xmu2,alp2,xmu3,alp3
  13 :       logical                 ogden
  14 :       common /elasdata_ogden/ kap,xmu,xmu1,alp1,xmu2,alp2,xmu3,alp3,
  15 :      .                        ogden
  16 :       ogden=.false.      
  17 :       you    = xmat( 1)
  18 :       xnu    = xmat( 2)
  19 : c en anteriores versiones el dato era realmente coh*sqrt(2/3)
  20 : c pq se multiplicaba luego por ese valor para usarlo en todos lados !!
  21 :       coh    = xmat( 5) 
  22 : c *dsqrt(2.D0/3.D0)
  23 :       phi    = xmat( 6)
  24 :       ar     = 0.05D0*cos(phi)
  25 :       tmax   = 25.D0*atan(1.D0)/45.D0
  26 :       aap    = cos(tmax)/3.D0*(3.D0+tan(tmax)*tan(3.D0*tmax)+
  27 :      .         sin(phi)*(tan(3.D0*tmax)-3.D0*tan(tmax))/sqrt(3.D0))
  28 :       aan    = cos(tmax)/3.D0*(3.D0+tan(tmax)*tan(3.D0*tmax)-
  29 :      .         sin(phi)*(tan(3.D0*tmax)-3.D0*tan(tmax))/sqrt(3.D0))
  30 :       bbp    = (sin(phi)*cos(tmax)/sqrt(3.D0)+sin(tmax))
  31 :      .         /3.D0/cos(3.D0*tmax)
  32 :       bbn    = (sin(phi)*cos(tmax)/sqrt(3.D0)-sin(tmax))
  33 :      .         /3.D0/cos(3.D0*tmax)
  34 :       return
  35 :       end
  36 :  
  37 :  
  38 :