Télécharger caljbc.eso

Retour à la liste

Numérotation des lignes :

  1. C CALJBC SOURCE CHAT 05/01/12 21:46:40 5004
  2. SUBROUTINE CALJBC(FN,GR,PG,XYZ,HR,PGSQ,RPG,NES,IES,NP,NPG,
  3. *IAXI,AIRE)
  4. C************************************************************************
  5. C
  6. C CALJBC APPEL CALJ22
  7. C
  8. C LA DIFFERENCE AVEC CALJAC EST LE RANGEMENT DE HR
  9. C HR(IES,NP)
  10. C
  11. C
  12. C CALCUL DE L'INVERSE DU JACOBIEN AJ=1/J
  13. C CALCUL DE L'ELEMENT D'AIRE SQ=DET(J)
  14. C DANS LES CAS 1D 2D ET 3D
  15. C
  16. C NES DIMENSION ESPACE DE DEFINITION DE L'ELEMENT
  17. C IES DIMENSION ESPACE
  18. C NP NOMBRE DE NOEUDS DE L'ELEMENT
  19. C NPG NOMBRE DE POINTS D'INTEGRATION
  20. C
  21. C XYZ COORDONNEES
  22. C GR GRADIENT
  23. C
  24. C************************************************************************
  25. IMPLICIT INTEGER(I-N)
  26. IMPLICIT REAL*8 (A-H,O-Z)
  27. -INC CCREEL
  28. C
  29. REAL*8 FN(NP,NPG),GR(NES,NP,NPG),HR(IES,NP,NPG)
  30. REAL*8 PG(NPG),XYZ(IES,NP),PGSQ(NPG),RPG(NPG)
  31. REAL*8 AJ(3,3,64)
  32. C
  33. CALL CALJBR
  34. &(FN,GR,PG,XYZ,HR,PGSQ,RPG,NES,IES,NP,NPG,IAXI,AIRE,AJ,SGN)
  35. IF(IAXI.NE.0)THEN
  36. DEUPI=2.D0*XPI
  37. DO 1 L=1,NPG
  38. PGSQ(L)=PGSQ(L)*DEUPI*RPG(L)
  39. 1 CONTINUE
  40. ENDIF
  41.  
  42. RETURN
  43. 1002 FORMAT(10(1X,1PE11.4))
  44. 1001 FORMAT(20(1X,I5))
  45. END
  46.  
  47.  
  48.  
  49.  
  50.  
  51.  

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