The lifting flux at the boundary is not computed correctly for Navier-Stokes

[amrueda]

The lifting flux for Navier-Stokes is computed with the conservative variables and not converted to the other possible "lifting variables". In MHD, this is done correcty:

fluxo/src/equation/mhd/getboundaryflux.f90

Lines 545 to 562 in c7e0cc9

	DO SideID=1,nBCSides
	!for BR1 and BR2, lifting is in strong form: Flux=Flux-U_master
	#if PP_Lifting_Var==1
	P_m(:,:,:)=U_master(:,:,:,SideID)
	Floc(:,:,:)=Flux(:,:,:,SideID)
	#elif PP_Lifting_Var==2
	!convert BC lifting flux to primitive variables
	CALL ConsToPrimVec(nTotal_face,P_m,U_master(:,:,:,SideID)) !prim_var
	CALL ConsToPrimVec(nTotal_face,Floc,Flux(:,:,:,SideID)) !prim_var
	#elif PP_Lifting_Var==3
	!convert BC lifting flux to entropy variables,
	CALL ConsToEntropyVec(nTotal_face,P_m,U_master(:,:,:,SideID)) !entropy_var
	CALL ConsToEntropyVec(nTotal_face,Floc,Flux(:,:,:,SideID)) !entropy_var
	#endif /*PP_Lifting_Var**/
	DO q=0,PP_N; DO p=0,PP_N
	Flux(:,p,q,SideID)=(Floc(:,p,q)-P_m(:,p,q))*SurfElem(p,q,SideID)
	END DO; END DO
	END DO ! iSide