commit 802583e7b41e076c94a8279782af52dbd5df2f55
author avm99963 <me@avm99963.com> Thu May 12 10:56:56 2022 +0200
committer avm99963 <me@avm99963.com> Thu May 12 10:58:11 2022 +0200
tree f36579c43d929795735c14bc6ee04b36842ea7fa
parent 354187350193313e85cfed50a07c0c3971425a66

Afegir sessió 3 del lab de fenòmens

Change-Id: Ic21d679decc844e3c2899d6da4cca48301dbae0c

quad10/fenomens/lab/p3/energ.f90

diff --git a/quad10/fenomens/lab/p3/energ.f90 b/quad10/fenomens/lab/p3/energ.f90
new file mode 100644
index 0000000..4ee8508
--- /dev/null
+++ b/quad10/fenomens/lab/p3/energ.f90
@@ -0,0 +1,14 @@
+real*8 function ENERG(S, L, PBC)
+  integer*2 :: S(1:L, 1:L)
+  integer*4 :: I, J, L
+  integer*4 :: PBC(0:L+1)
+  real*8 :: ENE
+  ENE = 0.0d0
+  do I = 1, L
+    do J = 1, L
+      ENE = ENE - S(I, J)*S(PBC(I + 1), J) - S(I, J)*S(I, PBC(J + 1))
+    enddo
+  enddo
+  ENERG = ENE
+  return
+endfunction