The maps were created using the projections object in Astrobear, using a "field" for inverse Beta. We are choosing inverse beta because the denominator will not be allowed to get arbitrarily small due to pressure protections, whereas the magnetic field pressure could approach zero. Once inverse beta was defined for each cell using the fields object, we summed these through 2 dimensions, down the barrel and down the y-axis.

Modifications to the code

In fields.f90, added a new integer parameter, "Inv_Beta", and a case in the function "GetField":

CASE(Inv_Beta)

GetField = (half*sum(q(iBx:iBz)**2))/(Press(q))

which is used in problem.f90, as:

! For integrated inverse beta maps along the flow axis:

CALL CreateProjection(projection)
Projection%dim=1
Projection%Field(1)%iD=Inv_Beta
Projection%Field(1)%component=BOTHCOMP
Projection%field(1)%name='Inv_Beta_downx'