Télécharger delskm.eso

Retour à la liste

Numérotation des lignes :

  1. C DELSKM SOURCE CB215821 16/04/21 21:16:22 8920
  2.  
  3. SUBROUTINE DETERMINA_LS_KMA(SIG,NESCRI,KMA,IZONE)
  4. IMPLICIT INTEGER(I-N)
  5. integer kma,nescri,izone
  6. real*8 sig(*),xi1,xj2,valor
  7. real*8 nn1,nn2,eta0,aa1,aa2,bb1,bb2,
  8. . pia,pib,sigy0,pcc,rrr
  9. common /poder2_data/ nn1,nn2,eta0,aa1,aa2,bb1,bb2,
  10. . pia,pib,sigy0,pcc,rrr
  11. call Invari_J2(sig,3,xj2)
  12. call Invari_I1(sig,3,xi1)
  13. if (xi1.lt.pia) then
  14. kma=2
  15. valor=rrr**2*xj2+(xi1-pia)**2-(pia-pib)**2
  16. valor=xj2+((xi1-pia)/rrr)**2-aa2*sigy0**2/3.D0
  17. if (nescri.eq.1)
  18. . write(199,*)' CAP ',xi1,xj2,
  19. . pia,pib,rrr,((xi1-pia)/rrr)**2,aa2,valor
  20. izone = 0
  21. else
  22. kma=3
  23. valor=SQRT(2.D0*xj2)+bb1*xi1/3.D0-bb2
  24. if (nescri.eq.1)
  25. . write(199,*)' CONE',xi1,pia,valor,kma
  26. izone = 1
  27. endif
  28. return
  29. end
  30.  
  31.  
  32.  
  33.  
  34.  
  35.  

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