diff options
Diffstat (limited to 'Donjon/src/FLPRNT.f')
| -rw-r--r-- | Donjon/src/FLPRNT.f | 272 |
1 files changed, 272 insertions, 0 deletions
diff --git a/Donjon/src/FLPRNT.f b/Donjon/src/FLPRNT.f new file mode 100644 index 0000000..bdbdc22 --- /dev/null +++ b/Donjon/src/FLPRNT.f @@ -0,0 +1,272 @@ +*DECK FLPRNT
+ SUBROUTINE FLPRNT(IPMAP,NCH,NB,NX,NY,NZ,POWB,PBNM,ICHM,IBNM,POWC,
+ 1 PCHM,IPCH,BAVG,BFACT,CAVG,CFACT,IMPX)
+*
+*-----------------------------------------------------------------------
+*
+*Purpose:
+* Print the bundle and channel powers over the fuel lattice.
+*
+*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.
+* NY number of elements along y-axis.
+* NZ number of elements along z-axis.
+* POWB bundle powers in kW.
+* PBNM maximum bundle power.
+* ICHM maximum-power channel number.
+* IBNM maximum-power bundle number.
+* POWC channel powers in kW.
+* PCHM maximum channel power.
+* IPCH maximum-power channel number.
+* BAVG average bundle power.
+* BFACT bundle power-form factor.
+* CAVG average channel power.
+* CFACT channel power-form factor.
+* IMPX printing index: 0 = no print
+* 1 = minimal printing
+* 2 = channel power only
+* 3 = bundle power by plane only
+* 10 = bundle power by channel
+* any added values of 2, 3 and 10: 5,12,13,15
+*
+*-----------------------------------------------------------------------
+*
+ USE GANLIB
+*----
+* SUBROUTINE ARGUMENTS
+*----
+ TYPE(C_PTR) IPMAP
+ INTEGER NCH,NB,NX,NY,NZ,ICHM,IBNM,IPCH,IMPX
+ REAL POWB(NCH,NB),POWC(NCH),PBNM,PCHM
+ DOUBLE PRECISION BAVG,CAVG,BFACT
+*----
+* LOCAL VARIABLES
+*----
+ PARAMETER(IOUT=6)
+ REAL RADB(NX,NY,NB),RADC(NX,NY)
+ INTEGER MIX(NX*NY,NZ),NAMX(NX),NAMY(NY)
+ CHARACTER TEXT*12,CHANX*2,CHANY*2,TEXT1A*17,TEXT2A*17,TEXT3A*17,
+ 1 TEXT1B*17,TEXT2B*17,TEXT3B*17
+*
+ MIX(:NX*NY,NZ)=0
+ NAMX(:NX)=0
+ NAMY(:NY)=0
+ CALL LCMGET(IPMAP,'BMIX',MIX)
+* CHANNEL NAMES
+ CALL LCMGET(IPMAP,'XNAME',NAMX)
+ CALL LCMGET(IPMAP,'YNAME',NAMY)
+*----
+* BUNDLE POWERS OVER EACH CHANNEL
+*----
+ IF(IMPX.GE.10) WRITE(IOUT,1009)
+ IEL=0
+ ICH=0
+ JCM=0
+ ICM=0
+ JBM=0
+ IBM=0
+ DO 11 J=1,NY
+ DO 10 I=1,NX
+ IEL=IEL+1
+ DO 5 K=1,NZ
+ IF(MIX(IEL,K).NE.0)GOTO 6
+ 5 CONTINUE
+ GO TO 10
+ 6 ICH=ICH+1
+ IF(ICH.EQ.IPCH)THEN
+ JCM=J
+ ICM=I
+ ENDIF
+ IF(ICH.EQ.ICHM)THEN
+ JBM=J
+ IBM=I
+ ENDIF
+ IF(IMPX.GE.10) THEN
+ WRITE(TEXT,'(A9,I3.3)')'CHANNEL #',ICH
+ WRITE(CHANX,'(A2)') (NAMX(I))
+ WRITE(CHANY,'(A2)') (NAMY(J))
+ WRITE(IOUT,1000)TEXT,CHANY,CHANX,POWC(ICH)
+ IF(PBNM.LT.1.)THEN
+ WRITE(IOUT,'(1X,1P,12E11.4)')(POWB(ICH,IB),IB=1,NB)
+ ELSE IF(PBNM.LT.1000.)THEN
+ WRITE(IOUT,'(1X,12F11.3)')(POWB(ICH,IB),IB=1,NB)
+ ELSE
+ WRITE(IOUT,'(1X,12F11.1)')(POWB(ICH,IB),IB=1,NB)
+ ENDIF
+ ENDIF
+ 10 CONTINUE
+ 11 CONTINUE
+*
+ WRITE(TEXT1A,'(A6,I2,A8)') '(A4,',(NX/2),'(A9,2X))'
+ WRITE(TEXT2A,'(A4,I2,A6)') '(A3,',(NX/2),'F11.3)'
+ WRITE(TEXT3A,'(A8,I2,A6)') '(A3,1P,',(NX/2),'E11.3)'
+ WRITE(TEXT1B,'(A4,I2,A8)') '(A4,',NX-(NX/2),'(A9,2X))'
+ WRITE(TEXT2B,'(A4,I2,A6)') '(A3,',NX-(NX/2),'F11.3)'
+ WRITE(TEXT3B,'(A8,I2,A6)') '(A3,1P,',NX-(NX/2),'E11.3)'
+ IF((IMPX.LT.3).OR.((IMPX.GE.10).AND.(IMPX.LT.13)))GOTO 50
+*----
+* BUNDLE POWERS PER RADIAL PLANE
+*----
+ RADB(:NX,:NY,:NB)=0.0
+ WRITE(IOUT,1010)
+ DO IB=1,NB
+ IEL=0
+ ICH=0
+ DO 25 J=1,NY
+ DO 20 I=1,NX
+ IEL=IEL+1
+ DO 15 K=1,NZ
+ IF(MIX(IEL,K).NE.0)GOTO 16
+ 15 CONTINUE
+ GO TO 20
+ 16 ICH=ICH+1
+ RADB(I,J,IB)=POWB(ICH,IB)
+ 20 CONTINUE
+ 25 CONTINUE
+ ENDDO
+ DO IB=1,NB
+ WRITE(IOUT,1011)IB
+ WRITE(IOUT,TEXT1A)' ',(NAMX(I),I=1,(NX/2))
+ WRITE(IOUT,*)' '
+ DO 30 J=1,NY
+ WRITE(CHANY,'(A2)') (NAMY(J))
+ IF(INDEX(CHANY,'-').EQ.1)GOTO 30
+ IF(PBNM.GE.1.)THEN
+ WRITE(IOUT,TEXT2A)CHANY,(RADB(I,J,IB),I=1,(NX/2))
+ ELSE
+ WRITE(IOUT,TEXT3A)CHANY,(RADB(I,J,IB),I=1,(NX/2))
+ ENDIF
+ 30 CONTINUE
+ WRITE(IOUT,*)' '
+ WRITE(IOUT,TEXT1B)' ',(NAMX(I),I=(NX/2+1),NX)
+ WRITE(IOUT,*)' '
+ DO 40 J=1,NY
+ WRITE(CHANY,'(A2)') (NAMY(J))
+ IF(INDEX(CHANY,'-').EQ.1)GOTO 40
+ IF(PBNM.GE.1.)THEN
+ WRITE(IOUT,TEXT2B)CHANY,(RADB(I,J,IB),I=(NX/2+1),NX)
+ ELSE
+ WRITE(IOUT,TEXT3B)CHANY,(RADB(I,J,IB),I=(NX/2+1),NX)
+ ENDIF
+ 40 CONTINUE
+ ENDDO
+ 50 IF((IMPX.EQ.0).OR.(IMPX.EQ.1).OR.(IMPX.EQ.3).OR.(IMPX.EQ.4)
+ 1 .OR.(IMPX.EQ.10).OR.(IMPX.EQ.11).OR.(IMPX.EQ.13).OR.(IMPX.EQ.14))
+ 2 GOTO 90
+*----
+* CHANNEL POWERS IN RADIAL PLANE
+*----
+ RADC(:NX,:NY)=0.0
+ WRITE(IOUT,1013)
+ IEL=0
+ ICH=0
+ DO 65 J=1,NY
+ DO 60 I=1,NX
+ IEL=IEL+1
+ DO 55 K=1,NZ
+ IF(MIX(IEL,K).NE.0)GOTO 56
+ 55 CONTINUE
+ GO TO 60
+ 56 ICH=ICH+1
+ RADC(I,J)=POWC(ICH)
+ 60 CONTINUE
+ 65 CONTINUE
+ WRITE(IOUT,TEXT1A)' ',(NAMX(I),I=1,(NX/2))
+ WRITE(IOUT,*)' '
+ DO 70 J=1,NY
+ WRITE(CHANY,'(A2)') (NAMY(J))
+ IF(INDEX(CHANY,'-').EQ.1)GOTO 70
+ IF(PBNM.GE.1.)THEN
+ WRITE(IOUT,TEXT2A)CHANY,(RADC(I,J),I=1,(NX/2))
+ ELSE
+ WRITE(IOUT,TEXT3A)CHANY,(RADC(I,J),I=1,(NX/2))
+ ENDIF
+ 70 CONTINUE
+ WRITE(IOUT,*)' '
+ WRITE(IOUT,TEXT1B)' ',(NAMX(I),I=(NX/2+1),NX)
+ WRITE(IOUT,*)' '
+ DO 80 J=1,NY
+ WRITE(CHANY,'(A2)') (NAMY(J))
+ IF(INDEX(CHANY,'-').EQ.1)GOTO 80
+ IF(PBNM.GE.1.)THEN
+ WRITE(IOUT,TEXT2B)CHANY,(RADC(I,J),I=(NX/2+1),NX)
+ ELSE
+ WRITE(IOUT,TEXT3B)CHANY,(RADC(I,J),I=(NX/2+1),NX)
+ ENDIF
+ 80 CONTINUE
+*----
+* FINAL INFORMATION
+*----
+ 90 IF((IBM.EQ.0).OR.(JBM.EQ.0)) CALL XABORT('FLPRNT: INVALID POWERS')
+ WRITE(IOUT,1002)
+ WRITE(CHANX,'(A2)') (NAMX(IBM))
+ WRITE(CHANY,'(A2)') (NAMY(JBM))
+ IF(PBNM.LT.1000.)THEN
+ WRITE(IOUT,1003)PBNM,CHANY,CHANX,IBNM
+ ELSE
+ WRITE(IOUT,1016)PBNM,CHANY,CHANX,IBNM
+ ENDIF
+ IF(BAVG.LT.1000.)THEN
+ WRITE(IOUT,1005)BAVG
+ ELSE
+ WRITE(IOUT,1017)BAVG
+ ENDIF
+ FACT=1./REAL(BFACT)
+ IF((ICM.EQ.0).OR.(JCM.EQ.0)) CALL XABORT('FLPRNT: INVALID POWERS')
+ WRITE(IOUT,1006)BFACT,FACT
+ WRITE(CHANX,'(A2)') (NAMX(ICM))
+ WRITE(CHANY,'(A2)') (NAMY(JCM))
+ IF(PCHM.LT.10000.)THEN
+ WRITE(IOUT,1004)PCHM,CHANY,CHANX
+ ELSE
+ WRITE(IOUT,1018)PCHM,CHANY,CHANX
+ ENDIF
+ IF(CAVG.LT.10000.)THEN
+ WRITE(IOUT,1007)CAVG
+ ELSE
+ WRITE(IOUT,1019)CAVG
+ ENDIF
+ FACT=1./CFACT
+ WRITE(IOUT,1008)CFACT,FACT
+ RETURN
+*
+ 1000 FORMAT(/5X,A12,5X,'NAME:',1X,A2,A2,5X,'CHANNEL POWER =',1X,1P,
+ 1 E11.4,'kW')
+ 1002 FORMAT(/5X,5('--o--',6X)/)
+ 1003 FORMAT(/1X,'MAXIMUM BUNDLE POWER =',1X,F9.1,1X,'kW',
+ 1 3X,'=>',3X,'CHANNEL:',1X,A2,A2,3X,'BUNDLE #',I2.2)
+ 1004 FORMAT(/1X,'MAXIMUM CHANNEL POWER =',1X,F9.1,1X,'kW',
+ 1 3X,'=>',3X,'CHANNEL:',1X,A2,A2)
+ 1005 FORMAT(1X,'AVERAGE POWER OVER ALL BUNDLES',
+ 1 1X,'=',1X,F9.1,1X,'kW')
+ 1006 FORMAT(1X,'BUNDLE-POWER FORM FACTOR',2X,'=>',2X,
+ 1 'AVG/MAX =',1X,F8.4,3X,'(MAX/AVG = ',F8.4,')')
+ 1007 FORMAT(1X,'AVERAGE POWER OVER ALL CHANNELS',
+ 1 1X,'=',1X,F9.1,1X,'kW')
+ 1008 FORMAT(1X,'CHANNEL-POWER FORM FACTOR',2X,'=>',2X,
+ 1 'AVG/MAX =',1X,F8.4,2X,'(MAX/AVG = ',F8.4,')'/)
+ 1009 FORMAT(/20X,'** BUNDLE POWERS OVER EACH',
+ 1 1X,'CHANNEL (kW) **'/)
+ 1010 FORMAT(//20X,'** BUNDLE POWERS PER RADIAL',
+ 1 1X,'PLANE **'/)
+ 1011 FORMAT(//1X,'BUNDLE POWERS',1X,'(kW)',1X,
+ 1 '=>',1X,'RADIAL PLANE',1X,'#',I2.2/)
+ 1013 FORMAT(//20X,'** CHANNEL POWERS IN RADIAL',1X,'PLANE (kW) **'/)
+ 1016 FORMAT(/1X,'MAXIMUM BUNDLE POWER =',1X,F9.1,1X,'kW',
+ 1 3X,'=>',3X,'CHANNEL:',1X,A2,A2,3X,'BUNDLE #',I2.2)
+ 1017 FORMAT(1X,'AVERAGE POWER OVER ALL BUNDLES',
+ 1 1X,'=',1X,F9.1,1X,'kW')
+ 1018 FORMAT(/1X,'MAXIMUM CHANNEL POWER =',1X,F9.1,1X,'kW',
+ 1 3X,'=>',3X,'CHANNEL:',1X,A2,A2)
+ 1019 FORMAT(1X,'AVERAGE POWER OVER ALL CHANNELS',
+ 1 1X,'=',1X,F9.1,1X,'kW')
+ END
|