3-Dimensional MHD Simulation of Earth's Magnetosphere

1995.7.25  by T. Ogino


main program : ears1651.f

graphics program : earg2621.f
subroutine for graphics : sub3ahm.f

(NX,NY,NZ)=(60,30,30)        : grid number without boundary
(XL,YL,ZL)=(61.0,31.0,31.0)  : length in each direction
 HX=XL/FLOAT(NX+1)=1.0       : grid interval in x direction
 HY=YL/FLOAT(NY+1)=1.0       : grid interval in y direction
 HZ=ZL/FLOAT(NZ+1)=1.0       : grid interval in z direction
 T=0.5*HX*THX                : time interval

 X=0.5*HX*FLOAT(2*I-NX2-1+2*NXP) : x position versus grid number
 Y=0.5*HY*FLOAT(2*J-3)           : y position versus grid number
 Z=0.5*HZ*FLOAT(2*K-3)           : z position versus grid number

  where NXP=NX+2 and NXP=0.

 RO01=5.0E-4 (5/cc)          : mass density of solar wind
 PR01=3.56E-8                : pressure of solar wind
 VSW=0.044   (300km/s)       : speed of solar wind
 BIS=CP(11)=1.5E-4 (5nT)     : z-component of IMF

 EATT                        : resistivity
 RMUU                        : viscosity
 EUD0                        : friction or collision term


execution of main program
1. f77 -O ears1651.f
2. a.out &

  where file must be defined in open statement like

       open(12,file='/usr2/ogdata/mhd3a/mhd4001',access='sequential',
     1         form='unformatted')

or

1. f77 -o ears1651 -O ears1651.f
2. eras1651 &

execution of graphics program

1. f77 -c -O sub3ahm.f
2. f77 -c -O earg2621.f
3. f77 -O earg2621.o sub3ahm.o
4. a.out > zout001.ps &
5. ghostscript zout001.ps
6. lpr -s zout001.ps

or
3. f77 -o earg2651 -O earg2621.o sub3ahm.o
4. earg2651 > zout001.ps & 
5. ghostscript zout001.ps
6. lpr -s zout001.ps