diff options
| author | stainer_t <thomas.stainer@oecd-nea.org> | 2025-09-08 13:48:49 +0200 |
|---|---|---|
| committer | stainer_t <thomas.stainer@oecd-nea.org> | 2025-09-08 13:48:49 +0200 |
| commit | 7dfcc480ba1e19bd3232349fc733caef94034292 (patch) | |
| tree | 03ee104eb8846d5cc1a981d267687a729185d3f3 /Donjon/src/CVR.f | |
Initial commit from Polytechnique Montreal
Diffstat (limited to 'Donjon/src/CVR.f')
| -rw-r--r-- | Donjon/src/CVR.f | 114 |
1 files changed, 114 insertions, 0 deletions
diff --git a/Donjon/src/CVR.f b/Donjon/src/CVR.f new file mode 100644 index 0000000..3ca5280 --- /dev/null +++ b/Donjon/src/CVR.f @@ -0,0 +1,114 @@ +*DECK CVR
+ SUBROUTINE CVR(NENTRY,HENTRY,IENTRY,JENTRY,KENTRY)
+*
+*-----------------------------------------------------------------------
+*
+*Purpose:
+* Perform reordering of fuel-regions properties in the reactor core,
+* according to the specified voiding pattern.
+*
+*Copyright:
+* Copyright (C) 2007 Ecole Polytechnique de Montreal.
+*
+*Author(s):
+* D. Sekki
+*
+*Parameters: input
+* NENTRY number of data structures transfered to this module.
+* HENTRY name of the data structures.
+* IENTRY data structure type where:
+* IENTRY=1 for LCM memory object;
+* IENTRY=2 for XSM file;
+* IENTRY=3 for sequential binary file;
+* IENTRY=4 for sequential ASCII file.
+* JENTRY access permission for the data structure where:
+* JENTRY=0 for a data structure in creation mode;
+* JENTRY=1 for a data structure in modifications mode;
+* JENTRY=2 for a data structure in read-only mode.
+* KENTRY data structure pointer.
+*
+*Comments:
+* The CVR: module specification is:
+* FMAPV := CVR: FMAP :: (descrcvr) ;
+* where
+* FMAP : name of a read-only \emph{fmap} object,
+* created in the RESINI: module. This object must contain the non-perturbed
+* fuel-cell properties.
+* FMAPV : name of a new \emph{fmap} object,
+* that will contain the modified fuel-type indices and reordered coolant
+* densities according to the specified core-voiding pattern.
+* (descrcvr) : structure describing the input data to the CVR: module.
+*
+*-----------------------------------------------------------------------
+*
+ USE GANLIB
+*----
+* SUBROUTINE ARGUMENTS
+*----
+ INTEGER NENTRY,IENTRY(NENTRY),JENTRY(NENTRY)
+ TYPE(C_PTR) KENTRY(NENTRY)
+ CHARACTER HENTRY(NENTRY)*12
+*----
+* LOCAL VARIABLES
+*----
+ PARAMETER(NSTATE=40)
+ CHARACTER TEXT*12,HSIGN*12
+ INTEGER ISTATE(NSTATE),IGST(NSTATE)
+ DOUBLE PRECISION DFLOT
+ TYPE(C_PTR) IPMAP,JPMAP
+*----
+* PARAMETER VALIDATION
+*----
+ IF(NENTRY.NE.2)CALL XABORT('@CVR: TWO PARAMETERS EXPECTED.')
+ TEXT=HENTRY(1)
+ IF((IENTRY(1).NE.1).AND.(IENTRY(1).NE.2))CALL XABORT('@CVR:'
+ 1 //' LCM OBJECT EXPECTED AT LHS ('//TEXT//').')
+ IF(JENTRY(1).NE.0)CALL XABORT('@CVR: FUEL MAP OBJECT IN CRE'
+ 1 //'ATE MODE EXPECTED AT LHS ('//TEXT//').')
+ TEXT=HENTRY(2)
+ IF((IENTRY(2).NE.1).AND.(IENTRY(2).NE.2))CALL XABORT('@CVR:'
+ 1 //' LCM OBJECT EXPECTED AT RHS ('//TEXT//').')
+ IF(JENTRY(2).NE.2)CALL XABORT('@CVR: FUEL MAP OBJECT IN REA'
+ 1 //'D-ONLY MODE EXPECTED AT RHS ('//TEXT//').')
+ CALL LCMGTC(KENTRY(2),'SIGNATURE',12,HSIGN)
+ IF(HSIGN.NE.'L_MAP')THEN
+ TEXT=HENTRY(2)
+ CALL XABORT('@CVR: SIGNATURE OF '//TEXT//' IS '//HSIGN//
+ 1 '. L_MAP EXPECTED.')
+ ENDIF
+ IPMAP=KENTRY(1)
+ CALL LCMEQU(KENTRY(2),IPMAP)
+*----
+* RECOVER INFORMATION
+*----
+ ISTATE(:NSTATE)=0
+ CALL LCMGET(IPMAP,'STATE-VECTOR',ISTATE)
+ NB=ISTATE(1)
+ NCH=ISTATE(2)
+ NFUEL=ISTATE(7)
+ NPARM=ISTATE(8)
+* FUEL-MAP GEOMETRY
+ JPMAP=LCMGID(IPMAP,'GEOMAP')
+ IGST(:NSTATE)=0
+ CALL LCMGET(JPMAP,'STATE-VECTOR',IGST)
+ IF(IGST(1).NE.7)CALL XABORT('@CVR: ONLY 3-D CART'
+ 1 //'ESIAN GEOMETRY ALLOWED.')
+ NX=IGST(3)
+ NY=IGST(4)
+ NZ=IGST(5)
+* PRINTING INDEX
+ CALL REDGET(ITYP,NITMA,FLOT,TEXT,DFLOT)
+ IF(ITYP.NE.3)CALL XABORT('@CVR: CHARACTER DATA EXPECTED.')
+ IF(TEXT.NE.'EDIT')CALL XABORT('@CVR: KEYWORD EDIT EXPECTED.')
+ CALL REDGET(ITYP,NITMA,FLOT,TEXT,DFLOT)
+ IF(ITYP.NE.1)CALL XABORT('@CVR: INTEGER FOR EDIT EXPECTED.')
+ IMPX=MAX(0,NITMA)
+* READ INPUT DATA
+ CALL CVRDRV(IPMAP,NCH,NB,NFUEL,NPARM,NX,NY,NZ,NVD,IVD,IMPX)
+* UPDATE STATE-VECTOR
+ ISTATE(10)=NVD
+ ISTATE(11)=IVD
+ CALL LCMPUT(IPMAP,'STATE-VECTOR',NSTATE,1,ISTATE)
+ IF(IMPX.GT.1) CALL LCMLIB(IPMAP)
+ RETURN
+ END
|