Initial commit from Polytechnique Montreal

1 files changed, 91 insertions, 0 deletions

diff --git a/Dragon/src/LIBND7.f b/Dragon/src/LIBND7.f
new file mode 100644
index 0000000..730641c
--- /dev/null
+++ b/Dragon/src/LIBND7.f
@@ -0,0 +1,91 @@
+*DECK LIBND7
+      SUBROUTINE LIBND7 (MAXDIL,NGRO,NAMFIL,HNISOR,NDIL,DILUT)
+*
+*-----------------------------------------------------------------------
+*
+*Purpose:
+* Find the dilutions corresponding to a resonant isotope within a
+* library in NDAS format.
+*
+*Copyright:
+* Copyright (C) 2006 Ecole Polytechnique de Montreal
+*
+*Author(s): A. Hebert
+*
+*Parameters: input
+* MAXDIL  maximum number of dilutions.
+* NGRO    number of energy groups.
+* NAMFIL  NDAS library name.
+* HNISOR  library name of the isotope.
+*
+*Parameters: output
+* NDIL    number of finite dilutions.
+* DILUT   dilutions.
+*
+*Reference:
+* Copyright (C) from NDAS Atomic Energy of Canada Limited utility (2006)
+*
+*-----------------------------------------------------------------------
+*
+      USE FSDF
+      IMPLICIT NONE
+*----
+*  Subroutine arguments
+*----
+      CHARACTER NAMFIL*(*),HNISOR*12
+      INTEGER MAXDIL,NGRO,NDIL
+      REAL DILUT(MAXDIL)
+*----
+*  Local variables
+*----
+      INTEGER I,IND,IERR,ISOID,HEADER(16),NISOLB,NF,NTEM,IHEAD(200)
+      REAL RHEAD(200)
+      CHARACTER HSMG*131,TEXT8*8
+      EQUIVALENCE(RHEAD(1),IHEAD(1))
+*----
+*  Read NDAS library parameters
+*----
+      CALL XSDOPN(NAMFIL,IERR)
+      IF(IERR.NE.0) CALL XABORT('LIBND7: XSDOPN could not open Library'
+     >  //' files')
+      CALL XSDBLD(6001,HEADER,IERR)
+      IF(IERR.NE.0) CALL XABORT('LIBND7: XSDBLD could not read library'
+     > //' parameters')
+      IF(NGRO.NE.HEADER(2)) CALL XABORT('LIBND7: Invalid number of e'
+     > //'nergy groups')
+      NISOLB=HEADER(1)
+*----
+*  Read through NDAS file
+*----
+      IND=0
+      DO I=1,NISOLB
+        CALL XSDNAM(I,ISOID,TEXT8,IERR)
+        IF(IERR.NE.0) CALL XABORT('LIBND7: XSDNAM index overflow')
+        IF(TEXT8.EQ.HNISOR(:8)) THEN
+          IND=I
+          GO TO 10
+        ENDIF
+      ENDDO
+      WRITE (HSMG,100) HNISOR,NAMFIL
+      CALL XABORT(HSMG)
+*     Load nuclide header
+   10 CALL XSDISO(7000,6001,IND,RHEAD,IERR)
+      NF=IHEAD(5)
+      IF((NF.GE.1).AND.(NF.LE.3)) THEN
+        CALL XSDRES(IND,IHEAD,IERR)
+        NTEM=IHEAD(1)
+        NDIL=IHEAD(2)
+        IF(NDIL.GT.MAXDIL) CALL XABORT('LIBND7: MAXDIL overflow')
+        DO I=1,NDIL
+          DILUT(I)=RHEAD(2+NTEM+I)
+        ENDDO
+        NDIL=NDIL-1
+      ELSE
+        NDIL=0
+      ENDIF
+      CALL XSDCL()
+      RETURN
+*
+  100 FORMAT(26HLIBND7: Material/isotope ',A12,20H' is missing on NDAS,
+     > 12H file named ,A24,1H.)
+      END