|
| |
|
Главная Журналы CONSTRUCT THE X-DIRECTION GRID LI=NCVLX+2 XU(2)=0. XU(LI)=XL L2=LI-1 FCVLX=FLOAT(NCVLX) DO 20 1=3,L2 DD=FLOAT(1-2)/FCVLX IFIPOWERX.GT.0.) THEN XU(I)=XL*DD**POWERX ELSE XU (I) =XL* ( I . - ( 1 . -DD) ** (-POWERX) ) ENDIF 2 0 CONTINUE CONSTRUCT THE Y-DIRECTION GRID MI=NCVLY+2 YV(2)=0. YV(MI)=YL M2=MI-I FCVLY=FLOAT(NCVLY) DO 30 J=3,M2 DD=FLOAT(J-2)/FCVLY IF(POWERY.GT.O.) THEN YV(J)=YL*DD**POWERY ELSE YV(J)=YL*(I.-(I.-DD)**(-POWERY)) ENDIF 30 CONTINUE RETURN ENTRY ZGRID CONSTRUCT THE GRID ZONE-BY-ZONE CONSIDER THE X DIRECTION XU(2)=0. 12 = 2 DO IOO NZ=I,NZX FCVLX=FLOAT(NCVX(NZ)) ILAST=I2 I1=ILAST+I I2=ILAST+NCVX(NZ) DO IOO 1=11,12 DD=FLOAT(I-ILAST)/FCVLX IF(POWRX(NZ).GT.O.) THEN XU(I)=XU(ILAST)+XZONE(NZ)*DD**POWRX(NZ) ELSE XU(I)=XU(ILAST)+XZONE(NZ)*(I.-(I.-DD)*(-POWRX(NZ))) ENDIF 100 CONTINUE LI = I2 CONSIDER THE Y DIRECTION YV(2)=0. J2 = 2 DO 110 NZ=1,NZY FCVLY=FLOAT(NCVY(NZ)) JLAST=J2 J1=JLAST+1 J2=JLAST+NCVY(NZ) DO 110 J=J1,J2 DD=FLOAT(J-JLAST)/FCVLY IF(POWRY(NZ).GT.O.) THEN YV(J)=YV(JLAST)+YZONE(NZ)*DD**POWRY(NZ) ELSE YV(J)=YV(JLAST)+YZONE(NZ)*(1.-(1.-DD)**(-POWRY(NZ))) ENDIF 110 CONTINUE M1 = J2 RETURN ENTRY PRINT DO 999 IUNIT=IU1, IU2 COME HERE TO ARRANGE THE PRINTOUT OF TWO-DIMENSIONAL FIELDS IF(KPGR.NE.O) THEN CREATE PRINTOUT FOR GRID С WRITE(lUNIT,1) 1 FORMAT( ) IBEG=1 IEND=L1 IREP=(IEND-IBEG+7)/7 DO 200 K=1,IREP INCR=MIN(6,lEND-IBEG) ISTOP=IBEG+INCR WRITE(lUNIT,2) (I,I=IBEG,ISTOP) 2 FORMAT(/2X, I =,2X,7(14,5X) ) IF(M0DE.EQ.3) THEN WRITE(lUNIT,3) (X(I),I=IBEG,ISTOP) 3 FORMATdX, TH =,1P7E9.2) ELSE WRITE(lUNIT,4) (X(I),I=IBEG,ISTOP) 4 FORMAT(2X,X =,1P7E9.2) ENDIF IBEG=IST0P+1 200 CONTINUE С WRITE(lUNIT,1) JBEG=1 JEND=MI JREP=(JEND-JBEG+7)/7 DO 210 K=I,JREP INCR=MIN(6,JEND-JBEG) JSTOP=JBEG+INCR WRITE(lONIT,5) {J,J=JBEG,JSTOP) 5 FORMAT(/2X,J =,2X,7(14,5X)) WRITE(lUNIT,6) (У(J),J=JBEG,JSTOP) 6 FORMAT(2X,Y =,1P7E9.2) JBEG=JST0P+1 210 CONTINUE ENDIF CREATE PRINTOUT FOR THE VALUES OF DEPENDENT VARIABLES DO 220 N=1,NFMAX IF(KPRINT(N).NE.O) THEN WRITE(lUNIT,7) TITLE(N) 7 FORMAT(/IX,6(IH*),3X,A18,3X,6(IH*)/9Х,20(1H-)) IBEG=1 JBEG=1 IEND=L1 JEND=MI IREP=(IEND-IBEG+7)/7 DO 230 K=l,IREP INCR=MIN(6,lEND-IBEG) ISTOP=IBEG+INCR WRITE(lUNIT,8) (I,I=IBEG,ISTOP) 8 FORMAT(/ I =,I6,6I9) WRITE(lUNIT, 9) 9 FORMAT( J) DO 240 J=JEND,JBEG,-1 WRITE(I UNIT,10) J, (F(I,J,N) ,I = IBEG,ISTOP) 10 FORMAT(IX,12,3X,1P7E9.2) 24 0 CONTINUE IBEG=IST0P+1 230 CONTINUE ENDIF 220 CONTINUE 999 CONTINUE RETURN ENTRY PLOT 0PEN(UNIT=8,FILE=PLOTF) 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 [ 95 ] 96 97 98 99 |