Bug Report: Incorrect variable usage in find_bounds subroutine for EXTENDED_BOUNDS == 1

[Jason-JunjieWU]

I found a bug in the find_bounds subroutine when EXTENDED_BOUNDS == 1. In the turbulence equation block, the code currently uses U_C, but it should use U_CT to correctly handle the primitive/turbulent variables.

! Change from:
UTEMP(k,...) = U_C(ILOCAL_RECON3(I)%IHEXL(1,IQ))%VAL(1,...)

! To:
UTEMP(k,...) = U_CT(ILOCAL_RECON3(I)%IHEXL(1,IQ))%VAL(1,...)

And I have some small questions about the extended boundary. When is it best to enable EXTENDED_BOUNDS? Is it specifically recommended for MUSCL schemes higher than 2nd-order, or perhaps when using the MOGE limiter?