Télécharger jac_j2.eso

Retour à la liste

Numérotation des lignes :

  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
  16. call carac_mate_j2(xmat)
  17. if (lam.le.0.D0) then
  18. do i=1,ndims
  19. do j=1,ndims
  20. Kmat(ndimk*(i-1)+j)=0.D0
  21. enddo
  22. Kmat(ndimk*(i-1)+i)=1.D0
  23. enddo
  24. return
  25. endif
  26. if ((kiso.eq.0.D0).and.(velo.eq.0.D0)) then
  27. ndimv=0
  28. call jac0J2(Kmat,ndimk,sig,ndims,ndimv,lam,
  29. . nescri,ues,kerre)
  30. return
  31. else
  32. ndimx=ndims+1
  33. do i=1,ndims
  34. x(i)=sig(i)
  35. enddo
  36. x(ndimx)=vint
  37. void=0.D0
  38. call jac1J2(Kmat,ndimk,x,ndimx,lam,ndims,
  39. . nescri,ues,kerre)
  40. endif
  41. return
  42. end
  43.  
  44.  
  45.  

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