Minor changes to fenomens p4
Change-Id: I2eafd58ea0a4dbb32f5f2e94ca05d9a472949a0d
diff --git a/quad10/fenomens/lab/p4/generateSpinMatrix.f90 b/quad10/fenomens/lab/p4/generateSpinMatrix.f90
new file mode 100644
index 0000000..0128119
--- /dev/null
+++ b/quad10/fenomens/lab/p4/generateSpinMatrix.f90
@@ -0,0 +1,21 @@
+subroutine generateSpinMatrix(S, L, SEED)
+ use iso_fortran_env
+ implicit none
+ integer*4, intent(in) :: L
+ integer*4, intent(in) :: SEED
+ integer*2, intent(out) :: S(1:64,1:64)
+ integer*4 :: i, j
+ integer(int32) :: genrand_int31
+
+ call init_genrand(SEED)
+
+ do i = 1, L
+ do j = 1, L
+ if (genrand_int31() < 2**30) then
+ S(i, j) = 1
+ else
+ S(i, j) = -1
+ endif
+ enddo
+ enddo
+end subroutine generateSpinMatrix