Sources Esope | Cast3M

camarh
C CAMARH    SOURCE    CB215821  16/04/21    21:15:31     8920
c MODELO RHMC
CCC
C **********************************************************************
CCC
      SUBROUTINE CARAC_MATE_RHMC (XMAT)
      IMPLICIT INTEGER(I-N)
      real*8    xmat(*)
      real*8    you,xnu,coh,phi,ar,tmax,aap,aan,bbp,bbn
      common /elasdata/ you,xnu
      common /rhmcdata/ coh,phi,ar,tmax,aap,aan,bbp,bbn
      real*8                  kap,xmu,xmu1,alp1,xmu2,alp2,xmu3,alp3
      logical                 ogden
      common /elasdata_ogden/ kap,xmu,xmu1,alp1,xmu2,alp2,xmu3,alp3,
     .                        ogden
      ogden=.false.
      you    = xmat( 1)
      xnu    = xmat( 2)
c en anteriores versiones el dato era realmente coh*sqrt(2/3)
c pq se multiplicaba luego por ese valor para usarlo en todos lados !!
      coh    = xmat( 5)
c *SQRT(2.D0/3.D0)
      phi    = xmat( 6)
      ar     = 0.05D0*cos(phi)
      tmax   = 25.D0*atan(1.D0)/45.D0
      aap    = cos(tmax)/3.D0*(3.D0+tan(tmax)*tan(3.D0*tmax)+
     .         sin(phi)*(tan(3.D0*tmax)-3.D0*tan(tmax))/sqrt(3.D0))
      aan    = cos(tmax)/3.D0*(3.D0+tan(tmax)*tan(3.D0*tmax)-
     .         sin(phi)*(tan(3.D0*tmax)-3.D0*tan(tmax))/sqrt(3.D0))
      bbp    = (sin(phi)*cos(tmax)/sqrt(3.D0)+sin(tmax))
     .         /3.D0/cos(3.D0*tmax)
      bbn    = (sin(phi)*cos(tmax)/sqrt(3.D0)-sin(tmax))
     .         /3.D0/cos(3.D0*tmax)
      return
      end