Télécharger caljdc.eso

Retour à la liste

Numérotation des lignes :

  1. C CALJDC SOURCE CHAT 05/01/12 21:46:49 5004
  2. SUBROUTINE CALJDC(FN,GR,PG,XYZ,HR,PGSQ,RPG,ND,IES,NP,NPG,
  3. *IAXI,AIRE)
  4. C************************************************************************
  5. C
  6. C CALJDC 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 ND DIMENSION ESPACE DE DEFINITION DE L'ELEMENT (ALIAS NES)
  17. C IES DIMENSION ESPACE DE TRAVAIL >= ND
  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. C
  28. REAL*8 FN(NP,NPG),GR(ND,NP,NPG),HR(IES,NP,NPG),HHR(162)
  29. REAL*8 PG(NPG),XYZ(IES,NP),PGSQ(NPG),RPG(NPG)
  30. REAL*8 AJ(3,3,64)
  31. C
  32. C
  33. C WRITE(6,*)' *** SUB CALJBC *** ND,IES',ND,IES
  34. IF(ND.EQ.IES)GOTO 20
  35.  
  36. CALL CALJCC(FN,GR,PG,XYZ,HR,PGSQ,RPG,IES,ND,NP,NPG,IAXI,AIRE,AJ,
  37. &HHR)
  38. RETURN
  39. 20 CONTINUE
  40.  
  41. CALL CALJ22(FN,GR,PG,XYZ,HR,PGSQ,RPG,ND,NP,NPG,IAXI,AIRE,AJ)
  42.  
  43. RETURN
  44. 1002 FORMAT(10(1X,1PE11.4))
  45. 1001 FORMAT(20(1X,I5))
  46. END
  47.  
  48.  
  49.  
  50.  

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