Version 6 (modified by 3 years ago) ( diff ) | ,
---|
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.
! Physical bounds of boxes box_0=64000*reshape((/-1,-1,-1,1,1,1/),(/3,2/)) box_4=4000*reshape((/-1,-1,-1,1,1,1/),(/3,2/)) box_7=750*reshape((/-1,-1,-1,1,1,1/),(/3,2/)) ! index of boxes in each level's index space ibox_0=nint(box_0-spread(GxBounds(:,1),2,2))/levels(0)%dx + spread((/1,0/),1,3) if (maxlevel >= 4) then ibox_4=nint(box_4-spread(GxBounds(:,1),2,2))/levels(4)%dx + spread((/1,0/),1,3) if (maxlevel >= 7) then ibox_7=nint(box_7-spread(GxBounds(:,1),2,2))/levels(7)%dx + spread((/1,0/),1,3) end if end if CASE(0) open(unit=195,file='NewLargeGrid_density_g_cm3_grid.dat',status='old',form='formatted',action='read') open(unit=196,file='NewLargeGrid_v_x_cm_s_grid.dat',status='old',form='formatted',action='read') open(unit=197,file='NewLargeGrid_v_y_cm_s_grid.dat',status='old',form='formatted',action='read') open(unit=198,file='NewLargeGrid_v_z_cm_s_grid.dat',status='old',form='formatted',action='read') open(unit=199,file='NewLargeGrid_u_erg_g_grid.dat',status='old',form='formatted',action='read') ibox=ibox_0 CASE(4) open(unit=195,file='NewMedGrid_density_g_cm3_grid.dat',status='old',form='formatted',action='read') open(unit=196,file='NewMedGrid_v_x_cm_s_grid.dat',status='old',form='formatted',action='read') open(unit=197,file='NewMedGrid_v_y_cm_s_grid.dat',status='old',form='formatted',action='read') open(unit=198,file='NewMedGrid_v_z_cm_s_grid.dat',status='old',form='formatted',action='read') open(unit=199,file='NewMedGrid_u_erg_g_grid.dat',status='old',form='formatted',action='read') ibox=ibox_4 CASE(7) open(unit=195,file='NewSmallGrid_density_g_cm3_grid.dat',status='old',form='formatted',action='read') open(unit=196,file='NewSmallGrid_v_x_cm_s_grid.dat',status='old',form='formatted',action='read') open(unit=197,file='NewSmallGrid_v_y_cm_s_grid.dat',status='old',form='formatted',action='read') open(unit=198,file='NewSmallGrid_v_z_cm_s_grid.dat',status='old',form='formatted',action='read') open(unit=199,file='NewSmallGrid_u_erg_g_grid.dat',status='old',form='formatted',action='read') ibox=ibox_7 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
(z,x,y) | ![]() |
Attachments (2)
- oldData.png (2.9 MB ) - added by 3 years ago.
-
new_AMR7_500.png
(1.2 MB
) - added by 3 years ago.
old data with new AMR7 as 500
Note:
See TracWiki
for help on using the wiki.