modify GAUSS_DFDZ routine to correctly calculate DFUDX,DFPDX