From 7dfcc480ba1e19bd3232349fc733caef94034292 Mon Sep 17 00:00:00 2001
From: stainer_t <thomas.stainer@oecd-nea.org>
Date: Mon, 8 Sep 2025 13:48:49 +0200
Subject: Initial commit from Polytechnique Montreal

---
 Utilib/src/PSFREG.f | 45 +++++++++++++++++++++++++++++++++++++++++++++
 1 file changed, 45 insertions(+)
 create mode 100644 Utilib/src/PSFREG.f

(limited to 'Utilib/src/PSFREG.f')

diff --git a/Utilib/src/PSFREG.f b/Utilib/src/PSFREG.f
new file mode 100644
index 0000000..822d9a3
--- /dev/null
+++ b/Utilib/src/PSFREG.f
@@ -0,0 +1,45 @@
+*DECK PSFREG
+      SUBROUTINE PSFREG(ISPSP,NPTS,XYPTS)
+C
+C---------------------------  PSFREG  ---------------------------------
+C
+C  1- PROGRAMME STATISTICS:
+C      NAME     : PSFREG
+C      USE      : FILL REGION
+C                 REPLACES PSPLOT ROUTINE FILRGNC
+C
+C  2- ROUTINE PARAMETERS:
+C    INPUT/OUTPUT
+C      ISPSP    : PSP FILE UNIT                          I
+C      NPTS     : NUMBER OF POINTS                       I
+C      YXPTS    : POSITION (X,Y) OF LINE INTERSECTION    R(2,NPTS)
+C
+C---------------------------   PSFREG  --------------------------------
+C
+      IMPLICIT         NONE
+      INTEGER          ISPSP,NPTS
+      REAL             XYPTS(2,NPTS)
+C----
+C  LOCAL VARIABLES
+C----
+      CHARACTER        NAMSBR*6
+      REAL             CONVER
+      PARAMETER       (NAMSBR='PSFREG',CONVER=72.0)
+      INTEGER          IPT
+      CHARACTER        CMDSTR*132
+      CMDSTR='Np'
+      CALL PSCPUT(ISPSP,CMDSTR)
+      CMDSTR=' '
+      WRITE(CMDSTR,'(2(F8.2,1X),A1)')
+     >  XYPTS(1,1)*CONVER,XYPTS(2,1)*CONVER,'M'
+      CALL PSCPUT(ISPSP,CMDSTR)
+      DO 100 IPT=2,NPTS
+        CMDSTR=' '
+        WRITE(CMDSTR,'(2(F8.2,1X),A1)')
+     >    XYPTS(1,IPT)*CONVER,XYPTS(2,IPT)*CONVER,'L'
+        CALL PSCPUT(ISPSP,CMDSTR)
+ 100  CONTINUE
+      CMDSTR='cf'
+      CALL PSCPUT(ISPSP,CMDSTR)
+      RETURN
+      END
-- 
cgit v1.2.3