Télécharger jacmrs.eso

Retour à la liste

Numérotation des lignes :

  1. C JACMRS SOURCE CHAT 05/01/13 00:48:11 5004
  2. CCC
  3. C **********************************************************************
  4. CCC
  5. SUBROUTINE JAC_MRSMAC (KMAT,NDIMK,SIG,NDIMS,VINTCON,VINTCAP,LCON,
  6. . lcap,xmat,nescri,ues,nnumer,deltax,kerre)
  7. IMPLICIT INTEGER(I-N)
  8.  
  9. integer ndims,kerre,nescri,nnumer,ues,ndimk,
  10. . j,i,ndimx,kmod
  11. real*8 sig(ndims),vintcon,vintcap,Kmat(*),xmat(*),deltax,x(7),
  12. . lcon,lcap,lam
  13.  
  14. kerre=0
  15. do i=1,7
  16. x(i)=0.D0
  17. enddo
  18. call CaracmateMAC(xmat)
  19.  
  20. if ((lcon.le.0.D0).and.(lcap.le.0.D0)) then
  21. do i=1,ndims
  22. do j=1,ndims
  23. Kmat(ndimk*(i-1)+j)=0.D0
  24. enddo
  25. Kmat(ndimk*(i-1)+i)=1.D0
  26. enddo
  27. return
  28. else if ((lcon.gt.0.D0).and.(lcap.gt.0.D0)) then
  29. write(*,*)' Modelo MRSLMAC: no tiene corner'
  30. return
  31. else if ((lcon.gt.0.D0).and.(lcap.lt.0.D0)) then
  32. kmod=21
  33. lam=lcon
  34. else
  35. kmod=22
  36. lam=lcap
  37. endif
  38. do i=1,ndims
  39. x(i)=sig(i)
  40. enddo
  41. ndimx=ndims+2
  42. x(ndims+1)=vintcon
  43. x(ndims+2)=vintcap
  44. call jac21(Kmat,ndimk,x,ndimx,lam,ndims,
  45. . kmod,nescri,ues,nnumer,deltax,kerre)
  46. return
  47. end
  48.  
  49.  
  50.  
  51.  

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