** geom2 **
IF(IDIR.EQ.1) THEN
 GRDPWR(X,10,WGAP,1.0)
 PATCH(IMSWL1,WWALL, 1, 1,1,NY,1,NZ,1,LSTEP)
 PATCH(IMSWL2,EWALL,NX,NX,1,NY,1,NZ,1,LSTEP)
ENDIF
IF(IDIR.EQ.2) THEN
 GRDPWR(Y,10,WGAP,1.0)
 PATCH(IMSWL1,SWALL,1,NX, 1, 1,1,NZ,1,LSTEP)
 PATCH(IMSWL2,NWALL,1,NX,NY,NY,1,NZ,1,LSTEP)
ENDIF
IF(IDIR.EQ.3) THEN
 GRDPWR(Z,10,WGAP,1.0)
 PATCH(IMSWL1,LWALL,1,NX,1,NY, 1, 1,1,LSTEP)
 PATCH(IMSWL2,HWALL,1,NX,1,NY,NZ,NZ,1,LSTEP)
ENDIF