Numérotation des lignes :

C CAMARH    SOURCE    CB215821  16/04/21    21:15:31     8920c MODELO RHMCCCCC **********************************************************************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

© Cast3M 2003 - Tous droits réservés.
Mentions légales