Télécharger valmel.eso

Retour à la liste

Numérotation des lignes :

valmel
  1. C VALMEL SOURCE OF166741 24/10/07 21:15:51 12016
  2. SUBROUTINE VALMEL (IPMELV,IPMINT,IPRES)
  3. *______________________________________________________________________
  4. *
  5. * CALCUL DES MELVALS AU SUPPORT DE IPMINT
  6. *
  7. * IPMELV POINTEUR SUR UN SEGMENT MELVAL A RECALCULER
  8. * (SUPPOSE ACTIF)
  9. *
  10. * IPMINT POINTEUR SUR LE SEGMENT D'INTEGRATION DU SUPPORT OU L'ON
  11. * SOUHAITE PROJETER LE MELVAL
  12. * (SUPPOSE ACTIF EN ENTREE ET EN SORTIE SANS CHANGMENT)
  13. *
  14. * IPRES POINTEUR SUR LE MELVAL RESULTAT
  15. * CREE PUIS ACTIF EN SORTIE
  16. *
  17. *_____________________________________________________________________
  18.  
  19. IMPLICIT INTEGER(I-N)
  20. IMPLICIT REAL*8(A-H,O-Z)
  21.  
  22. -INC PPARAM
  23. -INC CCOPTIO
  24.  
  25. -INC SMCHAML
  26. -INC SMINTE
  27.  
  28. MELVAL=IPMELV
  29. N1PTEL=VELCHE(/1)
  30. N2PTEL=IELCHE(/1)
  31. IF(N2PTEL.NE.0)THEN
  32. IF(N2PTEL.EQ.1)THEN
  33. SEGINI,MELVA1=MELVAL
  34. IPRES=MELVA1
  35. RETURN
  36. ELSE
  37. MOTERR(1:8)=' '
  38. CALL ERREUR(124)
  39. IPRES=IPMELV
  40. RETURN
  41. ENDIF
  42. ELSE
  43.  
  44. IF (N1PTEL.EQ.1) THEN
  45. SEGINI,MELVA1=MELVAL
  46. IPRES=MELVA1
  47. RETURN
  48. ELSE
  49. *
  50. * ON RECUPERE LE NOMBRE DE POINTS SUPPORT
  51. *
  52. MINTE=IPMINT
  53. C* SEGACT MINTE <- ACTIF EN ENTREE
  54. NBPGAU=SHPTOT(/3)
  55. NBNO =SHPTOT(/2)
  56. if (velche(/1).gt.1.and.velche(/1).ne.nbno) then
  57. interr(1)=nbno
  58. interr(2)=velche(/1)
  59. call erreur(1139)
  60. endif
  61.  
  62. N1PTEL=NBPGAU
  63. N1EL =VELCHE(/2)
  64. N2PTEL=0
  65. N2EL =0
  66. SEGINI MELVA1
  67. IPRES=MELVA1
  68. DO IB=1,N1EL
  69. DO IGAU=1,NBPGAU
  70. r_z = 0.D0
  71. DO INBNO=1,NBNO
  72. IGMN=MIN(VELCHE(/1),INBNO)
  73. r_z = r_z + VELCHE(IGMN,IB)*SHPTOT(1,INBNO,IGAU)
  74. ENDDO
  75. MELVA1.VELCHE(IGAU,IB) = r_z
  76. ENDDO
  77. ENDDO
  78. ENDIF
  79. ENDIF
  80.  
  81. C* SEGDES MELVAL <- ACTIF EN SORTIE
  82. C* SEGDES MINTE <- ACTIF EN SORTIE (NON MODIFIE)
  83.  
  84. RETURN
  85. END
  86.  
  87.  
  88.  

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