||January 20, 2003 12:30
It sounds like you want to use a structured rectangular mesh. I suggest you partition your problem space into blocks. The region either side of the discontinuity would lie in different blocks. You could then partition judiciously-chosen blocks by logically-rectangular (structured) meshes compatible at the block interfaces. If using FEM, you could assemble each rectangular block separately, and couple the blocks with constraints before solving. If your graphics program for plotting results is limited to rectangular meshes, then plot each block separately on the same graph.