From 93fd8246fbe3e660e4a27fdfa3ead18006d98203 Mon Sep 17 00:00:00 2001
From: Joe Taylor <joe@princeton.edu>
Date: Thu, 15 Feb 2024 14:12:45 -0500
Subject: [PATCH] Ann new file getpp3.f90.

---
 lib/superfox/getpp3.f90 | 22 ++++++++++++++++++++++
 1 file changed, 22 insertions(+)
 create mode 100644 lib/superfox/getpp3.f90

diff --git a/lib/superfox/getpp3.f90 b/lib/superfox/getpp3.f90
new file mode 100644
index 000000000..38d96357d
--- /dev/null
+++ b/lib/superfox/getpp3.f90
@@ -0,0 +1,22 @@
+subroutine getpp3(s3,workdat,p)
+
+  use sfox_mod
+  real s3(NQ,NN)
+  integer workdat(NN)
+  integer a(NN)
+
+!  a(1:NN)=workdat(NN:1:-1)
+  a=workdat
+
+  psum=0.
+  do j=1,NN
+     i=a(j)+1
+     x=s3(i,j)
+     s3(i,j)=0.
+     psum=psum + x
+     s3(i,j)=x
+  enddo
+  p=psum/NN
+
+  return
+end subroutine getpp3