first example: NettlesOHT

Fixed many errors:
- coh sub elem node numbering order in xbrick
- coh3d elem Nmatrix was wrong previously
- in xbrick, nodecnc was taken as dof cnc, caused huge error; a separate
dofcnc array is constructed locally
- wedge elem 2 point integration has hourglass mode, cause crazy disp.
of nodes; changed to 9 pnt intg
- include residual stiffness Kres in parameter module and implemented in
interface and lamina modules
- all modules' local dynamic arrays are now deallocated upon exit.
- preprocessing module preproc_libnee now writes max 90 char. per line
for 2 lines when writing elem nodal cnc (previously it writes on 1 line
only, which exceeds the F90 line length limit of 132 sometimes)
- in uexternaldb: glb clock updated before each increment
- in output module: integer length was too small (i5) for the size of
mesh, changed all to i10

model must constrain all directions, in 3D, remember to constrain z dir,
otherwise the simulation wont converge

output module has problems, it causes the program to hault for no
reason. it needs to be changed

xbrick edge status partiton algorithm can be improved:
when a tip, ref, or trans elem has one/two edge that is broken (eg val
>= cohcrack), then it should be repartitioned