Initial commit from Polytechnique Montreal

1 files changed, 112 insertions, 0 deletions

diff --git a/Donjon/src/SIMLIB.f b/Donjon/src/SIMLIB.f
new file mode 100644
index 0000000..1cda07e
--- /dev/null
+++ b/Donjon/src/SIMLIB.f
@@ -0,0 +1,112 @@
+*DECK SIMLIB
+      SUBROUTINE SIMLIB(IMPX,MODE,KPMAP,IPLIB,NTOT,NIS,IFMIX,HFOLLO,
+     > RFOLLO)
+*
+*-----------------------------------------------------------------------
+*
+*Purpose:
+* Put/get number densities of particularized isotopes in the microlib
+*
+*Copyright:
+* Copyright (C) 2017 Ecole Polytechnique de Montreal
+*
+*Author(s): 
+* A. Hebert
+*
+*Parameters: input
+* IMPX    print parameter.
+* MODE    transfert mode (=1: get from KPMAP; =2: put to KPMAP).
+* KPMAP   HCYCLE subdirectory in the fuelmap.
+* IPLIB   pointer to the microlib.
+* NTOT    number of fuel bundles.
+* NIS     number of particularized isotopes.
+* IFMIX   fuel mixture assigned to each fuel bundle.
+* HFOLLO  character*8 names of the particularized isotopes.
+*
+*Parameters: input/output
+* RFOLLO  number densities of the particularized isotopes.
+*
+*-----------------------------------------------------------------------
+*
+      USE GANLIB
+*----
+*  SUBROUTINE ARGUMENTS
+*----
+      TYPE(C_PTR) KPMAP,IPLIB
+      INTEGER IMPX,MODE,NIS,NTOT,IFMIX(NTOT)
+      REAL RFOLLO(NTOT,NIS)
+      CHARACTER*8 HFOLLO(NIS)
+*----
+*  LOCAL VARIABLES
+*----
+      PARAMETER(NSTATE=40)
+      INTEGER ISTATE(NSTATE)
+      CHARACTER*12 HCYCL
+*----
+*  ALLOCATABLE ARRAYS
+*----
+      INTEGER, ALLOCATABLE, DIMENSION(:) :: IMIX,IVB
+      REAL, ALLOCATABLE, DIMENSION(:) :: DENS
+      CHARACTER(LEN=12), ALLOCATABLE, DIMENSION(:) :: HUSE
+*
+      IF(.NOT.C_ASSOCIATED(IPLIB)) THEN
+        CALL XABORT('SIMLIB: MICROLIB LCM OBJECT MISSING AT RHS.')
+      ENDIF
+      CALL LCMGET(IPLIB,'STATE-VECTOR',ISTATE)
+      NBMIX=ISTATE(1)
+      NBISO=ISTATE(2)
+      IF(NTOT.GT.NBMIX) CALL XABORT('SIMLIB: NBMIX OVERFLOW.')
+      ALLOCATE(HUSE(NBISO),DENS(NBISO),IMIX(NBISO),IVB(NBMIX))
+      CALL LCMGTC(IPLIB,'ISOTOPESUSED',12,NBISO,HUSE)
+      CALL LCMGET(IPLIB,'ISOTOPESMIX',IMIX)
+      CALL LCMGET(IPLIB,'ISOTOPESDENS',DENS)
+      IVB(:NBMIX)=0
+      IBM=0
+      DO ITOT=1,NTOT
+        IF(IFMIX(ITOT).EQ.0) CYCLE
+        IBM=IBM+1
+        IVB(IBM)=ITOT
+      ENDDO
+      CALL LCMGTC(KPMAP,'ALIAS',12,HCYCL)
+      IF(MODE.EQ.1) THEN
+*       recover number densities from KCYCLE directory
+        IF(IMPX.GE.0) WRITE(6,'(/34H SIMLIB: recover number densities ,
+     >  5Hfrom ,A,11H directory.)') HCYCL
+        CALL LCMGET(KPMAP,'FOLLOW',RFOLLO)
+        DO ISO=1,NBISO
+          IBM=IMIX(ISO)
+          ITOT=IVB(IBM)
+          IF(ITOT.EQ.0) CALL XABORT('SIMLIB: MISSING FUEL BUNDLE(1).')
+          DO JSO=1,NIS
+            IF(HUSE(ISO)(:8).EQ.HFOLLO(JSO)) THEN
+              DENS(ISO)=RFOLLO(ITOT,JSO)
+              GO TO 10
+            ENDIF
+          ENDDO
+   10     CONTINUE
+        ENDDO
+        CALL LCMPUT(IPLIB,'ISOTOPESDENS',NBISO,2,DENS)
+      ELSE IF(MODE.EQ.2) THEN
+*       put number densities in KCYCLE directory
+        IF(IMPX.GE.0) WRITE(6,'(/33H SIMLIB: put number densities in ,
+     >  A,11H directory.)') HCYCL
+        RFOLLO(:NTOT,:NIS)=0.0
+        DO ISO=1,NBISO
+          IBM=IMIX(ISO)
+          ITOT=IVB(IBM)
+          IF(ITOT.EQ.0) CALL XABORT('SIMLIB: MISSING FUEL BUNDLE(2).')
+          DO JSO=1,NIS
+            IF(HUSE(ISO)(:8).EQ.HFOLLO(JSO)) THEN
+              RFOLLO(ITOT,JSO)=DENS(ISO)
+              GO TO 20
+            ENDIF
+          ENDDO
+   20     CONTINUE
+        ENDDO
+        CALL LCMPUT(KPMAP,'FOLLOW',NTOT*NIS,2,RFOLLO)
+      ELSE
+        CALL XABORT('SIMLIB: INVALID VALUE OF MODE.')
+      ENDIF
+      DEALLOCATE(IVB,IMIX,DENS,HUSE)
+      RETURN
+      END