Changes between Version 9 and Version 10 of u/bliu/polarMap

Timestamp:

04/02/18 15:29:22 (8 years ago)

Author:

Baowei Liu

Comment:

—

u/bliu/polarMap

@@ -35 +35 @@
 || python code || [http://www.pas.rochester.edu/~bliu/pnStudy/PolarizationMap/polMap_const0.py python code] ||
 || data || [http://www.pas.rochester.edu/~bliu/pnStudy/PolarizationMap/polMap_const0.okc data from Visit] ||
+
+== MHD TEST ==
+
+|| const Bz everywhere || [[Image(http://www.pas.rochester.edu/~bliu/pnStudy/PolarizationMap/test/Bz_0.6_m1.png, width=250)]] || ||
+|| const aux_z and Bz everywhere || [[Image(http://www.pas.rochester.edu/~bliu/pnStudy/PolarizationMap/test/Bz_0.6.png, width=250)]] || same mag as 1 ||
+|| Bz in Outflow only ||  [[Image(http://www.pas.rochester.edu/~bliu/pnStudy/PolarizationMap/test/Bz_Jet0.6.png, width=250)]] || 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. 
+
+{{{
+
+      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
+}}}