Changes between Version 13 and Version 14 of u/bliu/polarMap

Timestamp:

04/02/18 16:15:57 (8 years ago)

Author:

Baowei Liu

Comment:

—

u/bliu/polarMap

@@ -104 +104 @@
 }}}
 
-3. 
+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
 
 {{{