Set magnetic field variables
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| integer, | intent(in) | :: | ndir |
function var_set_bfield(ndir) result(iw) integer, intent(in) :: ndir integer :: iw(ndir), idir if (allocated(iw_mag)) call errormsg("Error: set_mag was already called") allocate(iw_mag(ndir)) do idir = 1, ndir nwflux = nwflux + 1 nwfluxbc = nwfluxbc + 1 nw = nw + 1 iw_mag(idir) = nwflux iw(idir) = nwflux write(cons_wnames(nwflux),"(A1,I1)") "b", idir write(prim_wnames(nwflux),"(A1,I1)") "b", idir end do end function var_set_bfield