Numérotation des lignes :

bb302
C BB302     SOURCE    MAGN      10/05/31    21:15:00     6679      SUBROUTINE BB302(XREF,XCOPG,XPOPG,     $FFPGV,DFFPGV,$     FFPGP,DFFPGP,     \$     X,Y,PG,FN,GR,FM,GM,ND,NP,MP,NPG)      IMPLICIT INTEGER(I-N)      IMPLICIT REAL*8 (A-H,O-Z)C************************************************************************CC     CALCULE LES FONCTIONS DE FORME D'UN : TRI4   P1 + bulle P1C                                              ou  P1 + bulle P0C     ^ etaC     |C   a |n3C     |\C     | \     a=sqrt(2)C     |n4\C     |___\ _____>ksiC     0   aC     n1  n2CCC************************************************************************* Entrées      INTEGER ND,NP,MP,NPG*      REAL*8 XREF(ND,NP)      REAL*8 XCOPG(ND,NPG)      REAL*8 XPOPG(NPG)      REAL*8 FFPGV(NP,NPG)      REAL*8 DFFPGV(NP,ND,NPG)      REAL*8 FFPGP(MP,NPG)      REAL*8 DFFPGP(MP,ND,NPG)* Sorties      REAL*8 X(NPG),Y(NPG)      REAL*8 PG(NPG)      REAL*8 FN(NP,NPG),GR(ND,NP,NPG)      REAL*8 FM(MP,NPG),GM(ND,MP,NPG)* Travail      INTEGER IPG,IP,ID      REAL*8 UNSVOL,USVUSDC***      A=SQRT(2.D0)      XREF(1,1)=0.D0      XREF(2,1)=0.D0      XREF(1,2)=A      XREF(2,2)=0.D0      XREF(1,3)=0.D0      XREF(2,3)=A      XREF(1,4)=A/3.D0      XREF(2,4)=A/3.D0 * magn veut que les volumes de ses éléments fasse 1      UNSVOL=2.D0      USVUSD=UNSVOL**(1.D0/DBLE(ND))* Recopie des points de Gauss      DO 1 IPG=1,NPG         X(IPG)=XCOPG(1,IPG)*USVUSD         Y(IPG)=XCOPG(2,IPG)*USVUSD         PG(IPG)=XPOPG(IPG)*UNSVOL 1    CONTINUE* Recopie des fns. de forme vitesse      DO 3 IPG=1,NPG         DO 32 IP=1,NP            FN(IP,IPG)=FFPGV(IP,IPG)            DO 322 ID=1,ND               GR(ID,IP,IPG)=DFFPGV(IP,ID,IPG) 322        CONTINUE 32      CONTINUE 3    CONTINUE* Recopie des fns. de forme pression      DO 5 IPG=1,NPG         DO 52 IP=1,MP            FM(IP,IPG)=FFPGP(IP,IPG)            DO 522 ID=1,ND               GM(ID,IP,IPG)=DFFPGP(IP,ID,IPG) 522        CONTINUE 52      CONTINUE 5    CONTINUEC     WRITE(6,101)C     WRITE(6,1002)FMC     WRITE(6,1002)GMC     WRITE(6,1002)FNC     WRITE(6,1002)GRC     WRITE(6,101)      RETURN 1002 FORMAT(10(1X,1PD11.4)) 1001 FORMAT(20(1X,I5)) 101  FORMAT(1X,'... SUBPB302 ... FM,GM,FN,GR ',9(10H..........)/)C      END        

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