9.21.2021 New Phantom data from Orsola

1. Dimensions

A small correction.
The resolution of the three boxes is:

Grid 1 (r+g+b): 512x512x299

Grid 2 (g+b): 512x512x512

Grid 3 (b): 512x512x512

Grid 1:

# SPLASH: A visualisation tool for SPH data (c)2004-2021 Daniel Price and contributors
# binary_08000_density_grid.dat produced using "splash to grid" on file binary_08000
#
# time:
#   2.0201106E+01
#
# xmin,xmax,ymin,ymax,zmin,zmax:
#  -1.2000000E+04   1.2000000E+04  -1.2000000E+04   1.2000000E+04  -7.0000000E+03   7.0000000E+03
#
# file contains:
# density interpolated to 3D cartesian   x,y,z grid 
#
# written in the form: 
#   do k=1,nz
#      do j=1,ny
#         write(*,*) (dat(i,j,k),i=1,nx)
#      enddo
#   enddo
#
# grid dimensions:
# nx   ny   nz
         512         512         299

Grid 2:

# SPLASH: A visualisation tool for SPH data (c)2004-2021 Daniel Price and contributors
# binary_08000_density_grid.dat produced using "splash to grid" on file binary_08000
#
# time:
#   2.0201106E+01
#
# xmin,xmax,ymin,ymax,zmin,zmax:
#  -5.0000000E+03   5.0000000E+03  -5.0000000E+03   5.0000000E+03  -5.0000000E+03   5.0000000E+03
#
# file contains:
# density interpolated to 3D cartesian   x,y,z grid 
#
# written in the form: 
#   do k=1,nz
#      do j=1,ny
#         write(*,*) (dat(i,j,k),i=1,nx)
#      enddo
#   enddo
#
# grid dimensions:
# nx   ny   nz
         512         512         512

Grid 3:

# SPLASH: A visualisation tool for SPH data (c)2004-2021 Daniel Price and contributors
# binary_08000_density_grid.dat produced using "splash to grid" on file binary_08000
#
# time:
#   2.0201106E+01
#
# xmin,xmax,ymin,ymax,zmin,zmax:
#  -2.5000000E+03   2.5000000E+03  -2.5000000E+03   2.5000000E+03  -2.5000000E+03   2.5000000E+03
#
# file contains:
# density interpolated to 3D cartesian   x,y,z grid 
#
# written in the form: 
#   do k=1,nz
#      do j=1,ny
#         write(*,*) (dat(i,j,k),i=1,nx)
#      enddo
#   enddo
#
# grid dimensions:
# nx   ny   nz
         512         512         512

2. scales

Phantom units	Conversion factor
Length	6.960E+10 cm
Mass	1.989E+33 g
Density	5.901E+00g/cm3
Velocity	4.367E+07 cm/s

3. Figures

9.13.2021 Old data mapping and Simulations

1. Phantom data

Large grid:

# SPLASH: A visualisation tool for SPH data (c)2004-2014 Daniel Price
# time:
#   1.5868375E+01
#
# file contains:
# density [g/cm\u3\d] interpolated to 3D grid 
#
# written in the form: 
#   do k=1,nz
#      do j=1,ny
#         write(*,*) (dat(i,j,k),i=1,nx)
#      enddo
#   enddo
#
# grid dimensions:
# nx    ny    nz
         128         128         128

Medium grid:

#
# grid dimensions:
# nx    ny    nz
         128         128         128

Small grid

# grid dimensions:
# nx    ny    nz
         192         192         192

---

Details

​density large grid

​density medium grid

​density small grid

2. Parameters and Mapping code in Astrobear

!=============================================================================================
! Parameters Related to Space
!=============================================================================================
nDim     = 3                            ! number of dimensions for this problem (1-3)
GmX      = 128,128,128                  ! Base grid resolution [x,y,z]
MaxLevel = 7                            ! Maximum level for this simulation (0 is fixed grid)
LastStaticLevel = 0                     ! Use static AMR for levels through LastStaticLevel [-1]
GxBounds = -64000d0,-64000d0,-64000d0,64000d0,64000d0,64000d0   ! Problem boundaries in computational units,format: 
                                        ! (xlower, ylower, zlower, xupper, yupper, zupper)
                                        ! For 2D problems, set zlower and zupper to 0.d0.

    select case(info%level)
    case(0)
        do k=1,128
           do j=1,128
              read(195,*),a(1,:)
              read(196,*),a(2,:)
              read(197,*),a(3,:)
              read(198,*),a(4,:)
              read(199,*),a(5,:)
              if (j.ge.57.and.j.le.72) then
                  do i=1,16
                      Info%q(i,j-56,k,1)=a(1,56+i)/rscale
                      Info%q(i,j-56,k,2)=a(1,56+i)*a(2,56+i)/rscale/velscale
                      Info%q(i,j-56,k,3)=a(1,56+i)*a(3,56+i)/rscale/velscale
                      Info%q(i,j-56,k,4)=a(1,56+i)*a(4,56+i)/rscale/velscale
                      Info%q(i,j-56,k,5)=a(1,56+i)*(a(5,56+i)+0.5*(a(2,56+i)&
                      *a(2,56+i)+a(3,56+i)*a(3,56+i)+a(4,56+i)*a(4,56+i)))/pscale
                  end do
              end if
           end do
        end do
    case(4)
        do k=1,128
           do j=1,128
              read(195,*),a(1,:)
              read(196,*),a(2,:)
              read(197,*),a(3,:)
              read(198,*),a(4,:)
              read(199,*),a(5,:)
              Info%q(1:128,j,k,1)=a(1,:)/rscale
              Info%q(1:128,j,k,2)=a(1,1:128)*a(2,1:128)/rscale/velscale
              Info%q(1:128,j,k,3)=a(1,1:128)*a(3,1:128)/rscale/velscale
              Info%q(1:128,j,k,4)=a(1,1:128)*a(4,1:128)/rscale/velscale
              Info%q(1:128,j,k,5)=a(1,1:128)*(a(5,1:128)+0.5*(a(2,1:128)*a(2,1:128)+a(3,1:128)*a(3,1:128)+a(4,1:128)*a(4,1:128)))/pscale
           end do
        end do
    case(7)
        do k=1,192
           do j=1,192
              read(195,*),b(1,:)
              read(196,*),b(2,:)
              read(197,*),b(3,:)
              read(198,*),b(4,:)
              read(199,*),b(5,:)
              Info%q(1:192,j,k,1)=b(1,:)/rscale
              Info%q(1:192,j,k,2)=b(1,1:192)*b(2,1:192)/rscale/velscale
              Info%q(1:192,j,k,3)=b(1,1:192)*b(3,1:192)/rscale/velscale
              Info%q(1:192,j,k,4)=b(1,1:192)*b(4,1:192)/rscale/velscale
              Info%q(1:192,j,k,5)=b(1,1:192)*(b(5,1:192)+0.5*(b(2,1:192)*b(2,1:192)+b(3,1:192)*b(3,1:192)+b(4,1:192)*b(4,1:192)))/pscale
           end do
        end do

---

Details

​problem.f90

​global.data

​physics.data

​problem.data

3. output from AstroBear

Imag(oldData.png, width=800)