Télécharger pacalc.eso

Retour à la liste

Numérotation des lignes :

pacalc
  1. C PACALC SOURCE PV090527 23/04/28 07:45:29 11659
  2. c Sous-Programme du module MISTRAL1
  3. C --------------------------------------------------------------------------
  4. SUBROUTINE PACALC (DXMI,DXMA,XMAX, II, DYMA,DYPRMA,
  5. & YP, DXA,YPA,NP, DX)
  6. C --------------------------------------------------------------------------
  7. C Determination du pas de calcul DX en fonction :
  8. C - du minimum et du maximum imposes : DXMI et DXMA,
  9. C - de l'intervalle maximal de X considere XMAX,
  10. C - des increments maximaux autorises DYMA(I) des Y(I),
  11. C - des increments relatifs maximaux autorises DYPRMA(I)
  12. C de leurs derivees YP(I),
  13. C - des valeurs de ces derivees YP(I) en X,
  14. C et eventuellement, si NP est superieur ou egal a 1 :
  15. C - des derivees YPA(I) au X precedent,
  16. C - du pas de calcul precedent DXA,
  17. C II etant le nombre de variables Y considerees pour determiner DX.
  18. C --------------------------------------------------------------------------
  19. IMPLICIT INTEGER(I-N)
  20. IMPLICIT REAL*8 (A-H, O-Z)
  21. -INC CCREEL
  22. DIMENSION DYMA(1:*),DYPRMA(1:*),YP(1:*),YPA(1:*)
  23. DX = DXMA
  24. IF (DXMI.GE.DXMA) RETURN
  25. IF (NP.EQ.0) THEN
  26. DO I = 1,II
  27. IF (ABS(YP(I))*DX.GT.DYMA(I)) THEN
  28. DX = DYMA(I)/ABS(YP(I))
  29. END IF
  30. END DO
  31. ELSE
  32. DO I = 1,II
  33. AYPI = ABS(YP(I))
  34. IF (AYPI*DX.GT.DYMA(I)) THEN
  35. AYSI = ABS((YP(I)-YPA(I))/DXA)
  36. DYPIMA = MAX(DYMA(I)/XMAX,AYPI*DYPRMA(I))
  37. IF (AYSI*DX.GT.DYPIMA) THEN
  38. DXYI = DYMA(I)/AYPI
  39. DXYPI = DYPIMA/(AYSI+XPETIT)
  40. DX = MAX(DXYI,DXYPI)
  41. END IF
  42. END IF
  43. END DO
  44. END IF
  45. DX = MAX(DX,DXMI)
  46. RETURN
  47. END
  48.  
  49.  
  50.  
  51.  

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