Télécharger resdpr.eso

Retour à la liste

Numérotation des lignes :

  1. C RESDPR SOURCE CB215821 16/04/21 21:18:11 8920
  2.  
  3. SUBROUTINE RESIDUO_DPRAL(X,R,N,NDIMX,NMODEL,XTRI,SIG,VECM)
  4. IMPLICIT INTEGER(I-N)
  5. integer n,i,j,ndims,ndimv,nmodel,ndimx
  6. real*8 x(n),r(n),void,sig(4),vecm(4),xtri(ndimx)
  7. integer augla
  8. real*8 c
  9. common /auglagrang1/ augla
  10. common /auglagrang2/ c
  11. call zzero(sig,4)
  12. call zzero(vecm,4)
  13. void=0.D0
  14. ndimv=1
  15. ndims=3
  16. call der_enerelas_dpral(x,sig,nmodel)
  17. if (ndimx.eq.3) then
  18. call yieldd(sig,ndims,void,ndimv,r(n),nmodel)
  19. call vflsig(sig,ndims,void,ndimv,vecm,nmodel)
  20. else if (ndimx.eq.4) then
  21. call yieldd(sig,ndims,x(ndimx),ndimv,r(n),nmodel)
  22. call vflsig(sig,ndims,x(ndimx),ndimv,vecm,nmodel)
  23. call vflvar(sig,ndims,x(ndimx),ndimv,vecm(ndimx),nmodel)
  24. endif
  25. do i=1,ndimx
  26. r(i)=x(i)-xtri(i)+ABS(x(n)+c*r(n))*vecm(i)
  27. enddo
  28. C
  29. c VISCOSIDAD FICTICIA AÑADIDA
  30. c r(n)=r(n)-1.D-8
  31. c
  32. return
  33. end
  34.  
  35.  
  36.  
  37.  
  38.  
  39.  

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