Numerical Simulation of Space Plasmas (I)
Example of Frotran 77 Programs

by Ling-Hsiao Lyu / lyu@jupiter.ss.ncu.edu.tw


 

Lectures Programs

Program Descriptions

1 GETASCII.f This program shows the binary structures of the ASCII code.
1 ascii_test.f Example of Fortran Program to determine the ASCII code of a given character and vice versa.
1 integer_test.f

Example of Fortran Program to determine the binary structure of a given integer number.

1 real_test.f Example of Fortran Program to determine the binary structure of a given real number.
1 maximum_test.f Example of Fortran Program to determine the extrema of the real number & integer number that can be resolved by the current computing system.
(Note that the extrema are not unique but are machine dependent.)
1 GETU.f This single-precision subroutine determines machine-dependent relative error relative to 1.
1 DGETU.f This double-precision subroutine determines machine-dependent relative error relative to 1.
2 Fast Fourier Transform This is the Fast Fourier transform code originally published at IBM/SSP. This code has been modified by Prof. John Dawson's group in UCLA to fit in the prticle simulation codes developed at UCLA. Ling-Hsiao Lyu made additional modification in response to the default setting in most of the modern FORTRAN compliers.
4 ran.f This function obtains a machine-dependent uniform random number.
This function is good for 32 bits computer.
This function is modified from IBM/SSP subroutine RANDU
4 ran_func.f This subroutine obtains a random number of a given function FUNC(x).
4 GAUSS.f This subroutine obtains a random number of a normal distribution function with a given mean and standard deviation based on the method described on page 77 in IBM/SSP.
6 1-D electrostatic particle code This code originally come from UCLA. This is a second-order finite-size particle code. Spatial deposiition is constructed based on the dipole approximation. Time integration is constructed based on the leap-frog method.
It has been modified by Ling-Hsiao Lyu in response to the default setting in the modern FORTRAN compliers. The input parameters are good for study electron-electron, or electron-ion, or ion-ion two-stream instabilities.

 


Contact information:

886-3-4227151 ext. 65771 lyu@jupiter.ss.ncu.edu.tw