Esope source files

   1 : C JAC_J2    SOURCE    CHAT      05/01/13    00:48:08     5004           
   2 : CCC
   3 : C **********************************************************************
   4 : CCC
   5 :       SUBROUTINE JAC_J2 (KMAT,NDIMK,SIG,NDIMS,VINT,LAM,XMAT,
   6 :      .                   nescri,ues,kerre)
   7 :       IMPLICIT INTEGER(I-N)
   8 :       integer ndims,kerre,nescri,ues,ndimk,ndimx,ndimv,i,j
   9 :       real*8  sig(ndims),vint,Kmat(*),xmat(*),x(7),void,lam
  10 :       real*8            sigy0,kiso,siginf,velo
  11 :       common /j2j2data/ sigy0,kiso,siginf,velo
  12 :       kerre=0
  13 :       do i=1,7
  14 :         x(i)=0.D0
  15 :       enddo
  15 :       call carac_mate_j2(xmat)
  16 :       if (lam.le.0.D0) then
  17 :          do i=1,ndims
  18 :            do j=1,ndims
  19 :              Kmat(ndimk*(i-1)+j)=0.D0
  20 :            enddo
  21 :            Kmat(ndimk*(i-1)+i)=1.D0
  22 :          enddo
  23 :          return      
  24 :       endif
  25 :       if ((kiso.eq.0.D0).and.(velo.eq.0.D0)) then
  26 :          ndimv=0
  26 :          call jac0J2(Kmat,ndimk,sig,ndims,ndimv,lam,
  27 :      .               nescri,ues,kerre)
  28 :          return
  29 :       else
  30 :          ndimx=ndims+1
  31 :          do i=1,ndims
  32 :               x(i)=sig(i)
  33 :          enddo
  34 :          x(ndimx)=vint
  35 :          void=0.D0
  35 :          call jac1J2(Kmat,ndimk,x,ndimx,lam,ndims,
  36 :      .               nescri,ues,kerre)
  37 :       endif
  38 :       return
  39 :       end
  40 :  
  41 :