commit ebb87790813a741ac8c65f92ce47a8c27ace30c4
author Adrià Vilanova Martínez <me@avm99963.com> Sat Jun 04 20:07:20 2022 +0200
committer Adrià Vilanova Martínez <me@avm99963.com> Sat Jun 04 20:07:20 2022 +0200
tree c428936c1b22ffb40d927545721f7b7a8ef7b2b8
parent 61897898c5b660cdec6e1d3a06804f88fede4a89

Minor changes to fenomens p4

Change-Id: I2eafd58ea0a4dbb32f5f2e94ca05d9a472949a0d

quad10/fenomens/lab/p4/MC-2.f90

diff --git a/quad10/fenomens/lab/p4/MC-2.f90 b/quad10/fenomens/lab/p4/MC-2.f90
index 3a2260d..48c6a11 100644
--- a/quad10/fenomens/lab/p4/MC-2.f90
+++ b/quad10/fenomens/lab/p4/MC-2.f90
@@ -3,8 +3,8 @@
   use iso_fortran_env
   implicit none
   integer*4 :: L
-  real*8 :: MAGNE, ENERG, TEMP, E, DIFE, DELTA, M, SUMA, V1, V2, V3, V4, W(-8:8)
-  integer*4 :: I, J, IMC, MCTOT, MCINI, MCD, IPAS, N, SEED, NSEED, SEED0, PBC
+  real*8 :: magne, energ, TEMP, E, DIFE, DELTA, M, SUMA, W(-8:8)
+  integer*4 :: I, J, IMC, MCTOT, MCINI, MCD, IPAS, N, SEED, NSEED, SEED0, pbc
   integer :: genrand_int31
   character(100) :: NOM
 
@@ -43,23 +43,22 @@
     print *, "Seed: ", SEED
 
     ! Inicialitzem la matriu d'spins aleatòriament
-    call WRITECONFIGSEED(S, L, SEED)
+    call generateSpinMatrix(S, L, SEED)
 
-    E = ENERG(S, L)
-    M = MAGNE(S, L)
+    E = energ(S, L)
+    M = magne(S, L)
     !write (12, *) "0", E, M
 
     ! Iterem amb el mètode de Montecarlo
     do IMC = 1, MCTOT
       do IPAS = 1, N
         ! LOS LOOPS HACEN COSAS
-        I = MOD(genrand_int31(), L) + 1
-        J = MOD(genrand_int31(), L) + 1
-        V1 = S(PBC(I - 1, L), J)
-        V2 = S(PBC(I + 1, L), J)
-        V3 = S(I, PBC(J - 1, L))
-        V4 = S(I, PBC(J + 1, L))
-        SUMA = V1 + V2 + V3 + V4
+        I = mod(genrand_int31(), L) + 1
+        J = mod(genrand_int31(), L) + 1
+        SUMA = S(PBC(I - 1, L), J) + &
+               S(PBC(I + 1, L), J) + &
+               S(I, PBC(J - 1, L)) + &
+               S(I, PBC(J + 1, L))
         DIFE = 2*SUMA*S(I, J)
 
         if (DIFE > 0) then