• Code version for 3D JetClump Module with MHD enabled: 6aeed49a12538e0661c85eb44ddd0ed102aba781

• Data files for JetClump Module

​global.data

​physics.data

​problem.data

---

• 3D Results @ z/2 for frame 8

logRho
Bx
By
Bz

• Extract polarization data using visit: ​Python code running on Visit from Martin

• ​Xmdv Data from Python code

• ​Matlab drawing polarization map

• Result:
  

---

• Constant Magnetic Field test

Bz
Polarization Map
python code	​python code
data	​data from Visit

---

---

MHD TEST

const Bz everywhere
const aux_z and Bz everywhere		same mag as 1
const aux_y and By everywhere	;	unzero Bz at frame 1
Bz in Outflow only		similar as 1&2

MHD TEST CODE

1. const Bz everywhere

         DO i=1-rmbc, mx+rmbc  ; x = (xl+(REAL(i,xPrec)-half)*dx)
         DO j=1-rmbc, my+rmbc  ; y = (yl+(REAL(j,xPrec)-half)*dy)
         DO k=1-zrmbc,mz+zrmbc ; z = (zl+(REAL(k,xPrec)-half)*dz)
   
             q(i,j,k,iBx)=0d0
             q(i,j,k,iBy)=0d0
             q(i,j,k,iBz)=sigmaB
             q(i,j,k,iE)=q(i,j,k,iE)+half*sigmaB**2
   
             CALL ConvertInternaltoTotalEnergy(q(i,j,k,:))
         END DO
         END DO
         END DO
   

2. const aux along z direction and Bz everywhere

      if ( nDim .eq. 3 .AND. lMHD .AND. lTestConstBz) then
         aux(1:mx+1, 1:my, 1:mz, 1) = 0d0
         aux(1:mx, 1:my+1, 1:mz, 2) = 0d0
         DO i=1, mx
           DO j=1, my
             DO k=1, mz+1
               pos=CellPos(Info, i, j, k)
               aux(i,j,k,3)=sigmaB !constant along z
             END DO
           END DO
         END DO
       else
         aux(:,:,:,:)=0d0
      end if
      
     ....

      DO i=1-rmbc, mx+rmbc  ; x = (xl+(REAL(i,xPrec)-half)*dx)
      DO j=1-rmbc, my+rmbc  ; y = (yl+(REAL(j,xPrec)-half)*dy)
      DO k=1-zrmbc,mz+zrmbc ; z = (zl+(REAL(k,xPrec)-half)*dz)

          q(i,j,k,iBx)=0d0
          q(i,j,k,iBy)=0d0
          q(i,j,k,iBz)=sigmaB
          q(i,j,k,iE)=q(i,j,k,iE)+half*sigmaB**2

          CALL ConvertInternaltoTotalEnergy(q(i,j,k,:))
      END DO
      END DO
      END DO

3. Const aux along y direction and By everywhere

         if ( nDim .eq. 3 .AND. lMHD .AND. lTestConstBy) then
            aux(1:mx+1, 1:my, 1:mz, 1) = 0d0
            aux(1:mx, 1:my, 1:mz+1, 3) = 0d0
            DO i=1, mx
              DO j=1, my+1
                DO k=1, mz
                  aux(i,j,k,2)=sigmaB !constant along z
                END DO
              END DO
            END DO
         else
            aux(:,:,:,:)=0d0
         end if
        ....
   
       
         DO i=1-rmbc, mx+rmbc  ; x = (xl+(REAL(i,xPrec)-half)*dx)
         DO j=1-rmbc, my+rmbc  ; y = (yl+(REAL(j,xPrec)-half)*dy)
         DO k=1-zrmbc,mz+zrmbc ; z = (zl+(REAL(k,xPrec)-half)*dz)
   
             q(i,j,k,iBx)=0d0
             q(i,j,k,iBy)=sigmaB
             q(i,j,k,iBz)=0d0
             q(i,j,k,iE)=q(i,j,k,iE)+half*sigmaB**2
   
             CALL ConvertInternaltoTotalEnergy(q(i,j,k,:))
         END DO
         END DO
         END DO
   
   

4. const Bz in outflow only

      DO i=1-rmbc, mx+rmbc  ; x = (xl+(REAL(i,xPrec)-half)*dx)
      DO j=1-rmbc, my+rmbc  ; y = (yl+(REAL(j,xPrec)-half)*dy)
      DO k=1-zrmbc,mz+zrmbc ; z = (zl+(REAL(k,xPrec)-half)*dz)

                 qjet(iBx) = 0d0 
                qjet(iBy) = 0d0 
                qjet(iBz) = sigmaB !sigmaFACTOR
                qjet(iE)=qjet(iE)+half*qjet(iBz)**2 
            ...
           qnew=qwind+qjet
           ...
            q(i,j,k,1:NrHydroVars)=q(i,j,k,1:NrHydroVars)+(qnew(1:NrHydroVars)-q(i,j,k,1:NrHydroVars))
            CALL ConvertInternaltoTotalEnergy(q(i,j,k,:))

      END DO
      END DO
      END DO