src/d3d/lecture.f

315           READ(IFAC) NUMX,NUMY
318           NUMX = IWORK(1)
321         IF (NUMX.EQ.0) THEN
324             READ(IFAC) NUMX,NUMY
327             NUMX = IWORK(1)
334         NUMNP0 = NUMX*NUMY
335         IF (NUMX.GT.NXYMAX.OR.NUMY.GT.NXYMAX) THEN
338      &           '*** Votre fichier comprend',NUMX,' x',NUMY,' points'
341             PRINT*,'*** Your file contains',NUMX,' x',NUMY,' points'
348           READ(IFAC) (XSU(I),I=1,NUMX),(YSU(I),I=1,NUMY)
351           DO I=1,NUMX
355             YSU(I) = WORK(I+NUMX)
387               THETA = XSU(1+MOD(I-1,NUMX))
390               Z(I) = YSU(1+(I-1)/NUMX)
401               THETA = XSU(1+MOD(I-1,NUMX))
422               THETA = XSU(1+MOD(I-1,NUMX))
423               PHI   = YSU(1+(I-1)/NUMX)
437               THETA = XSU(1+MOD(I-1,NUMX))
438               PHI   = YSU(1+(I-1)/NUMX)
461               RRR   = XSU(1+MOD(I-1,NUMX))
462               THETA = YSU(1+(I-1)/NUMX)
476               THETA = YSU(1+(I-1)/NUMX)
519           DO I=1,NUMX
538               X(I) = XSU(1+MOD(I-1,NUMX))
539               Y(I) = YSU(1+(I-1)/NUMX)
551               XI = XSU(1+MOD(I-1,NUMX))
552               YI = YSU(1+(I-1)/NUMX)
595           CALL LIENTIER('Nombre de points en X ?',0,NUMX)
598           CALL LIENTIER('Number of points in the X direction ?',0,NUMX)
601         NUMX = MAX(NUMX,3)
603         IF (NUMX.GT.NXYMAX.OR.NUMY.GT.NXYMAX) THEN
613         NUMNP = NUMX*NUMY
658         PASX = XL/REAL(NUMX-1)
660         DO I=1,NUMX
667           X(I) = XSU(1+MOD(I-1,NUMX))
668           Y(I) = YSU(1+(I-1)/NUMX)
694           CALL LIENTIER('Nombre de points en U ?',0,NUMX)
697           CALL LIENTIER('Number of points along U ?',0,NUMX)
700         NUMX = MAX(NUMX,3)
702         IF (NUMX.GT.NXYMAX.OR.NUMY.GT.NXYMAX) THEN
712         NUMNP = NUMX*NUMY
745         PASU = UL/REAL(NUMX-1)
747         DO I=1,NUMX
788      &     ,UUU(1+MOD(I-1,NUMX)),VVV(1+(I-1)/NUMX),X(I),ICODE)
800      &     ,UUU(1+MOD(I-1,NUMX)),VVV(1+(I-1)/NUMX),Y(I),ICODE)
812      &     ,UUU(1+MOD(I-1,NUMX)),VVV(1+(I-1)/NUMX),Z(I),ICODE)
839           READ(IFAC,*) NUMX,NUMY
842           READ(IFAC,*) NUMX,NUMY,MAXVAL
844         NUMNP = NUMX*NUMY
845         IF (NUMX.GT.NXYMAX.OR.NUMY.GT.NXYMAX) THEN
848      &           '*** Votre fichier comprend',NUMX,' x',NUMY,' points'
852      &           '*** Your file contains',NUMX,' x',NUMY,' points'
928         XMAX = XMIN+DX*REAL(NUMX-1)
934         DO I=1,NUMX
942         NN = (NUMY/2)*NUMX
945             P = 1 + MOD(I-1,NUMX)
946             Q = 1 + (I-1)/NUMX
947             II = (NUMY-Q)*NUMX + P
952           X(I) = XSU(1+MOD(I-1,NUMX))
953           Y(I) = YSU(1+(I-1)/NUMX)