Sources Esope | Cast3M

jacmrs
C JACMRS    SOURCE    CHAT      05/01/13    00:48:11     5004
CCC
C **********************************************************************
CCC
      SUBROUTINE JAC_MRSMAC (KMAT,NDIMK,SIG,NDIMS,VINTCON,VINTCAP,LCON,
     .                  lcap,xmat,nescri,ues,nnumer,deltax,kerre)
      IMPLICIT INTEGER(I-N)
 
      integer   ndims,kerre,nescri,nnumer,ues,ndimk,
     .          j,i,ndimx,kmod
      real*8    sig(ndims),vintcon,vintcap,Kmat(*),xmat(*),deltax,x(7),
     .          lcon,lcap,lam
 
      kerre=0
      do i=1,7
          x(i)=0.D0
        enddo
      call CaracmateMAC(xmat)
 
      if ((lcon.le.0.D0).and.(lcap.le.0.D0)) then
         do i=1,ndims
           do j=1,ndims
             Kmat(ndimk*(i-1)+j)=0.D0
           enddo
           Kmat(ndimk*(i-1)+i)=1.D0
         enddo
        return
      else if ((lcon.gt.0.D0).and.(lcap.gt.0.D0)) then
         write(*,*)' Modelo MRSLMAC: no tiene corner'
         return
      else if ((lcon.gt.0.D0).and.(lcap.lt.0.D0)) then
        kmod=21
        lam=lcon
      else
        kmod=22
        lam=lcap
      endif
      do i=1,ndims
          x(i)=sig(i)
        enddo
      ndimx=ndims+2
      x(ndims+1)=vintcon
      x(ndims+2)=vintcap
      call jac21(Kmat,ndimk,x,ndimx,lam,ndims,
     .           kmod,nescri,ues,nnumer,deltax,kerre)
      return
      end