Visualization to Various 3-Dimensional
MHD Models of Earth's Magnetosphere


We will show how to use the graphics programs to make VRML files of various 3-Dimensional MHD models of Earthfs Magnetosphere. The output binary files produced by various MHD simulation codes by using modified leap-frog scheme are used as input.
execution of VRML graphics program
1. f77 -c -O zvrsubb.f 2. f77 -O zvrmagb.f zvrsubb.o
where file must be defined in open statement like c open(10,file='earthb10.data', c 1 access='sequential',form='unformatted') open(11,file='earthb11.data', 1 access='sequential',form='unformatted') c
3. a.out & : output is written in fort.10 4. mv fort.10 fort.102
1. f77 -c -O zvrsubb.f 5. f77 -O zvrcrob.f zvrsubb.o 6. a.out & : output is written in fort.10 7. mv fort.10 fort.101 8. cat fort.101 fort.102 > zvrml01.wrl
Summary of parameters in MHD Simulation Code
(nx,ny,nz) = (180,60,60)
grid number without boundary
(nx,ny,nz) = (180,60,60)
grid number without boundary
nxp=30
parameter to determine earth position
(xxl,yyl,zzl) = (90.5,30.5,30.5)Re
length in each direction
hx=xxl/float(nx+1)=0.5Re
grid interval in x direction
hy=yyl/float(ny+1)=0.5Re
grid interval in y direction
hz=zzl/float(nz+1)=0.5Re
grid interval in z direction

1-dimensional array variable f(i1)=f(i,j,k,m)
i1=i+n1*(j-1)+n2*(k-1)+n3*(m-1)
m=1 : rho, plasma density
m=2 : Vx, x component of velocity
m=3 : Vy, y component of velocity
m=4 : Vz, z component of velocity
m=5 : P, plasma pressure
m=6 : Bx, x component of magnetic field
m=7 : By, y component of magnetic field
m=8 : Bz, z component of magnetic field

5-a. Half volume model of earth's magnetosphere with IMF By and Bz components
5-b. Quarter volume model of earth's magnetosphere
5-c. Half volume model of earth's magnetosphere with dipole tilt
5-d. Whole volume model of earth's magnetosphere

return to file list (English)
return to file list (Japanese)