if(:ask:.eq.f) then
 goto proceed
endif
LSTEP
MESG( WHAT LSTEP DO YOU WANT?
READVDU(LSTEP,INT,LSTEP)
LSTEP
 
    GROUP 3. X-direction grid specification
NX=12
NX
MESG( WHAT NX DO YOU WANT? X IS CIRCUMFERENTIAL FOR SPHERE
READVDU(NX,INT,NX)
NX
 
    GROUP 4. Y-direction grid specification
NY=12
NY
MESG( WHAT NY  DO YOU WANT? Y IS AZIMUTHAL FOR SPHERE
READVDU(NY,INT,NY)
NY
 
    GROUP 5. Z-direction grid specification
NZ=3
NZ
MESG( WHAT NZ  DO YOU WANT? Z IS RADIAL FOR SPHERE
READVDU(NZ,INT,nz)
NZ
label proceed
STEADY=F
GRDPWR(T,LSTEP,10.0,1.0)
#unigrid
 
    GROUP 6. Body-fitted coordinates or grid distortion
    The grid will be specified by way of coding in GROUND,
    and will vary with time.
char(movebfc);movebfc=$B005
#movebfc
 
    GROUP 8. Terms (in differential equations) & devices
TERMS(U1,Y,Y,N,Y,Y,Y);TERMS(V1,Y,Y,N,Y,Y,Y);TERMS(W1,Y,Y,N,Y,Y,Y)
 
    GROUP 9. Properties of the medium (or media)
REAL(GAMA); GAMA=1.4; RHO1=COMPRESS; RHO1B=1./GAMA; PRESS0=1.E5
RHO1A=1./PRESS0**rho1b
  DRH1DP=COMPRESS
    GROUP 11. Initialization of variable or porosity fields
FIINIT(P1)=0.0;FIINIT(U1)=0.0;FIINIT(V1)=0.0;FIINIT(W1)=0.0
FIINIT(CONI)=0.0;FIINIT(CONJ)=0.0;FIINIT(CONK)=0.0
    GROUP 15. Termination of sweeps
lsweep=50;selref=t;resfac=0.1

patch(fixu,east,1,nx,1,ny,1,nz,1,lstep)
coval(fixu,u1,1.e15,0.0)
patch(fixv,north,1,nx,1,ny,1,nz,1,lstep)
coval(fixv,u1,1.e15,0.0)
    GROUP 16. Termination of iterations
LITER(P1)=50
 
    GROUP 17. Under-relaxation devices
  RELAX(U1,FALSDT,0.0001);RELAX(V1,FALSDT,0.0001);RELAX(W1,FALSDT,0.1)
RELAX(P1,LINRLX,0.5)
 
    GROUP 19. Data communicated by satellite to GROUND
 
   
    Set the initial sphere size:
SPEDAT(GRIDS, CASE,  C, SPHERE)
SPEDAT(GRIDS, RADINN, R, 0.95)
SPEDAT(GRIDS, RADOUT, R, 1.00)
    Set the time factors
SPEDAT(GRIDS, TFCRIN, R, 1.00)
SPEDAT(GRIDS, TFCROU, R, 1.01)
SPEDAT(SET,GXMONI,TRANSIENT,L,F) 
    GROUP 22. Spot-value print-out
TSTSWP=-1;IXMON=NX/2;IYMON=NY/2;IZMON=NZ/2
 
    GROUP 23. Field print-out and plot control
NPLT=1;NXPRIN=1;NTPRIN=1
PATCH(LONGPLOT,PROFIL,1,10,1,1,1,1,1,LSTEP)
PLOT(LONGPLOT,P1,0.0,0.0);PLOT(LONGPLOT,U1,0.0,0.0)
 
CSG1=PHI;CSG2=XYZ;IDISPA=1
 
IPLTF=3;XZPR=T