TEXT(Solve For Velocities And Pressure
73 start of ....solvel
mesg(solvel macro has been loaded
GROUP 7. Variables stored, solved & named
SOLUTN(1,Y,Y,Y,P,P,P); OUTPUT(1,Y,Y,Y,Y,Y,Y)
BOOLEAN(NXNYNZ1);NXNYNZ1=T
IF(.NOT.ONEPHS) THEN
SOLUTN(9,Y,Y,Y,P,P,P); OUTPUT(9,Y,Y,Y,Y,Y,Y)
SOLUTN(10,Y,Y,Y,P,P,P); OUTPUT(10,Y,Y,Y,Y,Y,Y)
ENDIF
IF(NX.GT.1) THEN
NXNYNZ1=F
SOLUTN(3,Y,Y,Y,P,P,P); OUTPUT(3,Y,Y,Y,Y,Y,Y)
IF(.NOT.ONEPHS) THEN
+ SOLUTN(4,Y,Y,Y,P,P,P); OUTPUT(4,Y,Y,Y,Y,Y,Y)
ENDIF
ENDIF
IF(NY.GT.1) THEN
NXNYNZ1=F
SOLUTN(5,Y,Y,Y,P,P,P); OUTPUT(5,Y,Y,Y,Y,Y,Y)
IF(.NOT.ONEPHS) THEN
+ SOLUTN(6,Y,Y,Y,P,P,P); OUTPUT(6,Y,Y,Y,Y,Y,Y)
ENDIF
ENDIF
IF(NZ.GT.1) THEN
NXNYNZ1=F
+ SOLUTN(7,Y,Y,Y,P,P,P); OUTPUT(7,Y,Y,Y,Y,Y,Y)
IF(.NOT.ONEPHS) THEN
+ SOLUTN(8,Y,Y,Y,P,P,P); OUTPUT(8,Y,Y,Y,Y,Y,Y)
ENDIF
ENDIF
IF(NXNYNZ1) THEN
mesg(nx=ny=nz=1, so no velocities stored.
ENDIF
73 end of ....solvel