MHD with stress form

[guo-yichen]

In an MHD setup where stress form is enabled (i.e., ifmhd = .true. and ifstrs = .true.), during the boundary condition setting stage, the stsmask subroutine is invoked.

Nek5000/core/subs2.f

Lines 700 to 701 in 4eb31b9

	CALL SETHMSK (HVMASK,HFMASK,IFLD,NEL)
	CALL SETMLOG (HVMASK,HFMASK,IFLD,NEL)

It seems that SETHMSK subroutine only accounts for ifield=1(velocity). Should it also include for magnetic field?
https://github.com/Nek5000/Nek5000/blob/4eb31b918ed6c2b137bcd6b2d63f069f790c3c6f/core/subs2.f#L746C18-L746C18

When SETMLOG is executed,

Nek5000/core/subs2.f

Line 939 in 4eb31b9

CALL LFALSE (IFMSCR(1,1,IFLD),NTOTC)

it appears the sizes of IFMSFC may not be taking the magnetic field into consideration, as shown by the sizes in the GEOM.

Nek5000/core/GEOM

Lines 81 to 84 in 4eb31b9

	$ ,iflmsf(0:1),iflmse(0:1),iflmsc(0:1)
	$ ,ifmsfc(2*ldim,lelt,0:1)
	$ ,ifmseg(12,lelt,0:1)
	$ ,ifmscr(8,lelt,0:1)

Discussed with Yu-Hsiang and Paul