Working f90 massively parallel source codes for the CRAY T3E.
The codes heavily employs the CRAY-specific SHMEM routines.
All codes use Poisson asynchrony for the discrete-events (spin-flips) and a conservative synchronization scheme [(Lubachevsky,  J. Comput. Phys. 75, 103 (1988)] for the parallel discrete-event simulations.
No warranties, no regrets. Use them at your own risk.
Appropriate acknowledgement of the source and the story behind it is appreciated:
"Parallelization of a Dynamic Monte Carlo Algorithm: A Partially Rejection-Free Conservative Approach",
G. Korniss, M.A. Novotny, and P.A. Rikvold,  J. Comput. Phys.153, 488 (1999), (paper on the implementation).
"From Massively Parallel Algorithms and Fluctuating Time Horizons to Non-equilibrium Surface Growth",
G. Korniss, Z. Toroczkai, M.A. Novotny,  and P.A. Rikvold, Phys. Rev. Lett. 84, 1351 (2000), (paper on scalability).


Continuous-time "Metropolis" algorithm for the two-dimensional kinetic Ising model after an instateneous field reversal. The  code can be used to study hysteresis and metastability in large spatially-extended bistable systems.
par_ctmp.f90


Partially rejection-free  (n-fold-way-like) version of the same algorithm as above.
par_ctnf.f90


Continuous-time "Metropolis" algorithm for the two-dimensional kinetic Ising model in a periodic square-wave external field. The  code can be used to study the dynamic phase transition in large spatially-extended bistable systems.
par_dpt_ctmp_J.f90


The "KISS" random number generator by G. Marsaglia and A. Zaman (1993).
kiss.f90