c twod.inc
      COMMON /TWOD/ VOLd(LLOW:LIHIGH,LLOW:LJHIGH),
     $     BXd(LLOW:LIHIGH,LLOW:LJHIGH),BYd(LLOW:LIHIGH,LLOW:LJHIGH),
     $     BZd(LLOW:LIHIGH,LLOW:LJHIGH)
     $	   ,BX1d(LLOW:LIHIGH,LLOW:LJHIGH),BY1d(LLOW:LIHIGH,LLOW:LJHIGH),
     $     BZ1d(LLOW:LIHIGH,LLOW:LJHIGH),RHOd(LLOW:LIHIGH,LLOW:LJHIGH),
     q 	   VXd(LLOW:LIHIGH,LLOW:LJHIGH),VYd(LLOW:LIHIGH,LLOW:LJHIGH),
     $     VZd(LLOW:LIHIGH,LLOW:LJHIGH),Cd(LLOW:LIHIGH,LLOW:LJHIGH),
     q	   RHO1d(LLOW:LIHIGH,LLOW:LJHIGH),VX1d(LLOW:LIHIGH,LLOW:LJHIGH),
     $     VY1d(LLOW:LIHIGH,LLOW:LJHIGH),VZ1d(LLOW:LIHIGH,LLOW:LJHIGH),
     q	   C1d(LLOW:LIHIGH,LLOW:LJHIGH),
     $	   BId(LLOW:LIHIGH,LLOW:LJHIGH),BI1d(LLOW:LIHIGH,LLOW:LJHIGH),
     $     BJd(Llow:Lihigh,Llow:Ljhigh),BJ1d(Llow:lihigh,llow:ljhigh),
     $     xqd(llow:lihigh,llow:ljhigh),yqd(llow:lihigh,llow:ljhigh),
     $     zqd(llow:lihigh,llow:ljhigh),
     $     xnormdi(llow:lihigh,llow:ljhigh),
     $     ynormdi(llow:lihigh,llow:ljhigh),
     $     znormdi(llow:lihigh,llow:ljhigh),
     $     xnormdj(llow:lihigh,llow:ljhigh),
     $     ynormdj(llow:lihigh,llow:ljhigh),
     $     znormdj(llow:lihigh,llow:ljhigh),
     $     facedi(llow:lihigh,llow:ljhigh),
     $     facedj(llow:lihigh,llow:ljhigh)
	COMMON /TWOD/ rho2d(li,lj), VX2d(LI,LJ),VY2d(LI,LJ),
     q	   VZ2d(LI,LJ),C2d(LI,LJ), BDVX2d(LI,LJ),
     q		BDVY2d(LI,LJ),BDVZ2d(LI,LJ)
	COMMON /TWOD/ FACEd(LIP1,LJ),TRANSd(LIP1,LJ,9),
     $     BZqnd(LIP1,LJ),BXqnd(LIP1,LJ),BYqnd(LIP1,LJ),
     $	   BNqd(LIP1,LJ),bsqqd(lip1,lj),BXqd(LIP1,LJ),BYqd(LIP1,LJ),
     $     BZqd(LIP1,LJ),BNORMd(LIP1,LJ),XNORMd(LIP1,LJp1),
     $	   YNORMd(LIP1,LJp1),TRAN2d(LIP1,LJP1,9),BI2d(LIP1,LJ),
     $	   BJ2d(LI,LJP1),EDGEd(LIP1,LJP1),FACejd(LI,LJP1),
     $     bx0d(lip1,ljp1),by0d(lip1,ljp1),bz0d(lip1,ljp1),
     $       e2(lip1,ljp1)