c -- context.h --
c 
      real ddd0(nvar,1-nbdyf:npx+nbdyf, 1-nbdyf:npy+nbdyf, 1-nbdyf:npz+nbdyf+1)
      real ddd1(nvar,1-nbdyf:npx+nbdyf, 1-nbdyf:npy+nbdyf, 1-nbdyf:npz+nbdyf+1)
c
      common /context/ ddd0, ddd1