Improve calculation of geometry info for Nek BPs

[jeremylt]

I found a minor bug in the geodatq function when converting BP3/4 for libCEED. Specifically, the variable tmp had incorrect values between applying intp_rstd and loc_grad3, leading to minor variations in the correctness of the geometric factors. Creating a variable local to the subroutine than using the one in the common variable space corrected the issue.

The easiest way to verify this error is to insert a write statement for each value of jacmq in geodatq.

Also, I included the determinant of the Jacobian in the geometric factors on BP1/2.

I discovered these small changes because in libCEED we compute the RHS directly rather than calculating the RHS by applying the operator to the true solution. (see the MFEM and PETSc examples) If it is desirable, I add add something like that here.

I do not believe these changes have any performance impact, just a corectness impact.

[jeremylt]

Please let me know if any changes are required for this PR.

[v-dobrev]

@jeremylt , for this repository, feel free to push directly to master.