Télécharger kaam.eso

Retour à la liste

Numérotation des lignes :

  1. C KAAM SOURCE CHAT 05/01/13 00:50:19 5004
  2. SUBROUTINE KAAM(FN,HR,PG,PGSQ,AM,K,NEL,NP,NPG,IES,IAXI,RPG)
  3. C************************************************************************
  4. C CALCUL DE LA MATRICE M ( P*DIV(U) ) --> AM(NEL,NP,IES)
  5. C
  6. C************************************************************************
  7. IMPLICIT INTEGER(I-N)
  8. IMPLICIT REAL*8 (A-H,O-Z)
  9.  
  10. C***
  11. DIMENSION FN(NP,NPG),HR(IES,NP,NPG),PG(NPG),PGSQ(NPG),RPG(NPG)
  12. DIMENSION AM(NEL,NP,IES)
  13. C? DIMENSION PORO(NPR,*)
  14. C? DIMENSION LE(NP),IPADL(*)
  15. C***
  16.  
  17. IF(IAXI.EQ.0)CALL INITD(RPG,NPG,1.D30)
  18.  
  19. DO 20 N=1,IES
  20. IAX=0
  21. IF((3-IAXI).EQ.N)IAX=1
  22. IF(IES.EQ.3)IAX=0
  23.  
  24. IF(IAX.EQ.0)THEN
  25.  
  26. DO 2 J=1,NP
  27. SJ=0.D0
  28. C? NF=IPADL(LE(J))
  29. C? NFP=IPOR*(NF-1)+1
  30. DO 1 L=1,NPG
  31. 1 SJ=SJ+HR(N,J,L)*PGSQ(L)
  32. AM(K,J,N)=SJ
  33. C? AM(K,J,N)=SJ*PORO(N,NFP)
  34. 2 CONTINUE
  35.  
  36. ELSE
  37.  
  38. DO 12 J=1,NP
  39. SJ=0.D0
  40. C? NF=IPADL(LE(J))
  41. C? NFP=IPOR*(NF-1)+1
  42. DO 11 L=1,NPG
  43. 11 SJ=SJ+(HR(N,J,L)+FN(J,L)/RPG(L))*PGSQ(L)
  44. AM(K,J,N)=SJ
  45. C? AM(K,J,N)=SJ*PORO(N,NFP)
  46. 12 CONTINUE
  47.  
  48. ENDIF
  49. 20 CONTINUE
  50.  
  51. C
  52. C
  53. RETURN
  54. 1001 FORMAT(20(1X,I5))
  55. 1002 FORMAT(10(1X,1PE11.4))
  56. END
  57.  
  58.  

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