Initial commit from Polytechnique Montreal

1 files changed, 124 insertions, 0 deletions

diff --git a/Donjon/src/CVRPRN.f b/Donjon/src/CVRPRN.f
new file mode 100644
index 0000000..dc60cc5
--- /dev/null
+++ b/Donjon/src/CVRPRN.f
@@ -0,0 +1,124 @@
+*DECK CVRPRN
+      SUBROUTINE CVRPRN(IPMAP,NCH,NB,NX,NY,NZ,MIXNEW,PVALUE,LCOOL,IMPX)
+*
+*-----------------------------------------------------------------------
+*
+*Purpose:
+* Print the fuel-type indices per bundle for each reactor channel.
+*
+*Copyright:
+* Copyright (C) 2007 Ecole Polytechnique de Montreal.
+*
+*Author(s): 
+* D. Sekki
+*
+*Parameters: input/output
+* IPMAP   pointer to fuel-map information.
+* NCH     number of reactor channels.
+* NB      number of fuel bundles per channel.
+* NX      number of elements along x-axis in fuel map.
+* NY      number of elements along y-axis in fuel map.
+* NZ      number of elements along z-axis in fuel map.
+* MIXNEW  updated fuel-type index.
+* PVALUE  structure containing the modified coolant density
+*         values throughout the reactor core.
+* LCOOL   flag with respect to the coolant densities:
+*          =.true. coolant densities were modified;
+*          =.false. coolant densities not provided.
+* IMPX    printing index (=0 for no print).
+*
+*-----------------------------------------------------------------------
+*
+      USE GANLIB
+*----
+*  SUBROUTINE ARGUMENTS
+*----
+      TYPE(C_PTR) IPMAP
+      INTEGER NCH,NB,NX,NY,NZ,MIXNEW(NCH,NB),IMPX
+      REAL PVALUE(NCH,NB)
+      LOGICAL LCOOL
+*----
+*  LOCAL VARIABLES
+*----
+      PARAMETER(IOUT=6)
+      INTEGER MIX(NX*NY*NZ),MIXOLD(NCH,NB),
+     1 NSCH(NCH),NAMX(NX),NAMY(NY),FMIX(NX,NY,NB)
+      CHARACTER TEXT*12,CHANX*2,CHANY*2,FORM1*14,FORM2*14
+*----
+*  RECOVER INFORMATION
+*----
+      MIX(:NX*NY*NZ)=0
+      MIXOLD(:NCH,NB)=0
+      NSCH(:NCH)=0
+      NAMX(:NX)=0
+      NAMY(:NY)=0
+      FMIX(:NX,:NY,:NB)=0
+      CALL LCMGET(IPMAP,'BMIX',MIX)
+      CALL LCMGET(IPMAP,'FLMIX',MIXOLD)
+*     CHANNEL NAMES
+      CALL LCMGET(IPMAP,'XNAME',NAMX)
+      CALL LCMGET(IPMAP,'YNAME',NAMY)
+*     REFUELLING SCHEME
+      CALL LCMLEN(IPMAP,'REF-SCHEME',LENGT,ITYP)
+      IF(LENGT.EQ.0)CALL XABORT('@CVRPRN: MISSI'
+     1 //'NG REF-SCHEME DATA IN FUEL-MAP.')
+      CALL LCMGET(IPMAP,'REF-SCHEME',NSCH)
+      IF((IMPX.NE.2).AND.(IMPX.LT.4))GOTO 20
+*----
+*  PRINTING OVER EACH CHANNEL
+*----
+      IEL=0
+      ICH=0
+      DO 15 J=1,NY
+      DO 10 I=1,NX
+      IEL=IEL+1
+      IF(MIX(IEL).EQ.0)GOTO 10
+      ICH=ICH+1
+      WRITE(TEXT,'(A9,I3.3)')'CHANNEL #',ICH
+      WRITE(CHANX,'(A2)') (NAMX(I))
+      WRITE(CHANY,'(A2)') (NAMY(J))
+      WRITE(IOUT,1001)TEXT,CHANY,CHANX,NSCH(ICH)
+      WRITE(IOUT,1002) (MIXOLD(ICH,IB),IB=1,NB)
+      WRITE(IOUT,1003) (MIXNEW(ICH,IB),IB=1,NB)
+      IF(LCOOL)WRITE(IOUT,1004)(PVALUE(ICH,IB),IB=1,NB)
+   10 CONTINUE
+   15 CONTINUE
+   20 IF((IMPX.NE.3).AND.(IMPX.LT.4))GOTO 50
+*----
+*  PRINTING PER RADIAL PLANE
+*----
+      WRITE(FORM1,'(A4,I2,A8)')'(A4,',NX,'(A3,1X))'
+      WRITE(FORM2,'(A4,I2,A8)')'(A2,',NX,'(I3,1X))'
+      WRITE(IOUT,1005)
+      IEL=0
+      DO IB=1,NB
+        ICH=0
+        DO 35 J=1,NY
+        DO 30 I=1,NX
+        IEL=IEL+1
+        IF(MIX(IEL).EQ.0)GOTO 30
+        ICH=ICH+1
+        FMIX(I,J,IB)=MIXNEW(ICH,IB)
+   30   CONTINUE
+   35   CONTINUE
+      ENDDO
+      DO IB=1,NB
+        WRITE(IOUT,1006)IB
+        WRITE(IOUT,FORM1)' ',(NAMX(I),I=1,NX)
+        WRITE(IOUT,*)' '
+        DO 40 J=1,NY
+        WRITE(CHANY,'(A2)') (NAMY(J))
+        IF(INDEX(CHANY,'-').EQ.1)GOTO 40
+        WRITE(IOUT,FORM2)CHANY,(FMIX(I,J,IB),I=1,NX)
+   40   CONTINUE
+      ENDDO
+   50 RETURN
+*
+ 1001 FORMAT(/10X,'* ',A12,' *',5X,'NAME:',1X,A2,A2,
+     1 5X,'REF-SCHEME:',1X,I2)
+ 1002 FORMAT(2X,'OLD FUEL INDICES :',2X,12(I2,2X))
+ 1003 FORMAT(2X,'NEW FUEL INDICES :',2X,12(I2,2X)/)
+ 1004 FORMAT(2X,'COOLANT DENSITIES:',1X,12(F4.2,1X)/)
+ 1005 FORMAT(//20X,'** FUEL-TYPE INDICES PER RADIAL PLANE **')
+ 1006 FORMAT(//2X,'RADIAL PLANE',1X,'#',I2.2/)
+      END