Sources Esope | Cast3M

sigy
C SIGY      SOURCE    CB215821  16/04/21    21:18:27     8920
 
             FUNCTION SIGY (VINT)
      IMPLICIT INTEGER(I-N)
      real*8 vint,sigy0,kiso,siginf,velo,sigy
      common /j2j2data/ sigy0,kiso,siginf,velo
      if (velo.eq.0.D0) then
         sigy=(sigy0+kiso*vint)*sqrt(2.D0/3.D0)
      else
         sigy=(sigy0+kiso*vint+(siginf-sigy0)*(1.D0-EXP(-velo*vint)))
     .        *sqrt(2.D0/3.D0)
      endif
      return
      end