diff --git a/src/pec.f90 b/src/pec.f90
index 05c0482..5151da9 100644
--- a/src/pec.f90
+++ b/src/pec.f90
@@ -138,6 +138,8 @@ contains
     dpdv=facpds*(dpdv/dvol)
     ajphiv=facjs*(ajphiv/darea)
     ajcd=ajphiv*ratjbv
+    pins=facpds*pins
+    currins=facjs*currins
 
 ! now dpdv is  dP/dV [MW/m^3]
 ! now ajphiv is J_phi=dI/dA [MA/m^2]