Sources Esope | Cast3M

vecyim
C VECYIM    SOURCE    CHAT      05/01/13    04:05:29     5004
CCC
C **********************************************************************
CCC
      SUBROUTINE VECYIMAC (X,NDIMX,VECN,NMODEL)
      IMPLICIT INTEGER(I-N)
      integer ndimx,nmodel,ndims,ndimv,npcon,npcap,kdummy
      real*8 x(ndimx),vecn(ndimx)
      ndims=ndimx-2
      ndimv=2
      npcon=ndimx-1
      npcap=ndimx
      kdummy=nmodel
      if (nmodel.eq.22) kdummy=23
      call vyisigMAC(x,ndims,x(npcon),ndimv,vecn,kdummy)
      if (nmodel.eq.21) then
        call vyicoMACvar(x,ndims,x(npcon),vecn(npcon))
        vecn(npcap)=0.D0
      else if (nmodel.eq.22) then
        call vyicaMACvar(x,ndims,x(npcap),x(npcon),vecn(npcap))
        call vyicaMACvarco(x,ndims,x(npcap),x(npcon),vecn(npcon))
      endif
      return
      end