Group 1. run title and other preliminaries
TEXT(Gala Example 1D                   
TITLE
 
    ** This example appears to work satisfactorily for
       GALA=F but the pressure blows up when GALA=T
 
       Correction in Ground activated by USOURC=T
 
       1-D, X-direction, FIXFLU inlet, FIXP outlet
 
    Group 2. transience; time-step specification
STEADY=F;LSTEP=1; TLAST=0.1; GRDPWR(T,LSTEP,TLAST,1.0)
    Group 3. x-direction grid specification
CARTES=T; NX=10;XULAST=0.05; GRDPWR(X,NX,XULAST,1.0)
    Group 3. y-direction grid specification
NY=1;YVLAST=1.0; GRDPWR(Y,NY,YVLAST,1.0)
    Group 4. z-direction grid specification
NZ=1;ZWLAST=1.0; GRDPWR(Z,NZ,ZWLAST,1.0)
    Group 7. variables stored, solved & named
#solvel
    Group 8. terms (in differential equations) & devices
GALA=F;USOURC=F
    Group 9. properties of the medium (or media)
RHO1=1000.0;ENUL=1.E-4
    Group 13. boundary conditions and special sources
  ** Inlet
PATCH(GATE2,WEST,1,1,1,NY,1,NZ,1,LSTEP)
COVAL(GATE2,P1,FIXFLU,1000.0*0.001); COVAL(GATE2,U1,ONLYMS,0.001)
   ** Exit
PATCH(EXIT,EAST,NX,NX,1,NY,1,NZ,1,LSTEP)
COVAL(EXIT,P1,FIXP,0.0)
    Group 15. termination of sweeps
LSWEEP=10
RESREF(P1)=0.001*1.0*1.E-12; RESREF(U1)=1000.*0.001*0.001*1.0*1.E-12
    Group 16. termination of iterations
LITER(P1)=15;LITER(U1)= 5
    Group 22. spot-value print-out
IXMON=1; IYMON=1; IZMON=1
SPEDAT(SET,GXMONI,TRANSIENT,L,F)
    Group 23. field print-out and plot control
NTPRIN=LSTEP; TSTSWP=1; OUTPUT(P1,Y,Y,Y,Y,Y,Y)
OUTPUT(U1,Y,Y,Y,Y,Y,Y)