Télécharger la plaquette Cast3M


fr en

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