** geom1 **
IF(IDIR.EQ.1) THEN
 NREGX= 3; IREGX= 1; GRDPWR(X,1,WSL1,1.0)
 IREGX= 2; GRDPWR(X,10,WGAP,1.0); IREGX= 3; GRDPWR(X,1,WSL2,1.0)
 PATCH(SOL1,INIVAL, 1, 1,1,NY,1,NZ,1,LSTEP)
 PATCH(SOL2,INIVAL,NX,NX,1,NY,1,NZ,1,LSTEP)
 PATCH( HOT,WWALL, 1, 1,1,NY,1,NZ,1,LSTEP)
 PATCH(COLD,EWALL,NX,NX,1,NY,1,NZ,1,LSTEP)
ENDIF
IF(IDIR.EQ.2) THEN
 NREGY= 3; IREGY= 1; GRDPWR(Y,1,WSL1,1.0)
 IREGY= 2; GRDPWR(Y,10,WGAP,1.0); IREGY= 3; GRDPWR(Y,1,WSL2,1.0)
 PATCH(SOL1,INIVAL,1,NX, 1, 1,1,NZ,1,LSTEP)
 PATCH(SOL2,INIVAL,1,NX,NY,NY,1,NZ,1,LSTEP)
 PATCH( HOT,SWALL,1,NX, 1, 1,1,NZ,1,LSTEP)
 PATCH(COLD,NWALL,1,NX,NY,NY,1,NZ,1,LSTEP)
ENDIF
IF(IDIR.EQ.3) THEN
 NREGZ= 3; IREGZ= 1; GRDPWR(Z,1,WSL1,1.0)
 IREGZ= 2; GRDPWR(Z,10,WGAP,1.0); IREGZ= 3; GRDPWR(Z,1,WSL2,1.0)
 PATCH(SOL1,INIVAL,1,NX,1,NY, 1, 1,1,LSTEP)
 PATCH(SOL2,INIVAL,1,NX,1,NY,NZ,NZ,1,LSTEP)
 PATCH( HOT,LWALL,1,NX,1,NY, 1, 1,1,LSTEP)
 PATCH(COLD,HWALL,1,NX,1,NY,NZ,NZ,1,LSTEP)
ENDIF