Initial commit from Polytechnique Montreal

1 files changed, 280 insertions, 0 deletions

diff --git a/Donjon/data/AFMtest.x2m b/Donjon/data/AFMtest.x2m
new file mode 100644
index 0000000..81e6ac0
--- /dev/null
+++ b/Donjon/data/AFMtest.x2m
@@ -0,0 +1,280 @@
+!*****************************************************************
+!* Test for AFM module development in DONJON4                   **
+!*****************************************************************
+!* Calcul DONJON a l'equilibre du rechargement a deux (2) zones **
+!*                                                              **
+!* Determination du burnup moyen de sortie pour un coeur        **
+!* avec des grappes quelconques provenant de FBMDT0             **
+!*                                                              **
+!* COURS ENE6209                                                **
+!*                                                              **
+!* Options:                                                     **
+!*   1- Rapport B1/B2 (B1 zone centrale, B2 zone peripherique)  **
+!*   2- Position des BC                                         **
+!*****************************************************************
+!* B. Dionne (03/04/2001)                                       **
+!* E. Varin (23/02/2004)                                        **
+!* R. Chambon (18/10/2007)                                      **
+!* A. Landry (13/09/2011)                                       **
+!* A. Hebert (11/11/2011) (pour DONJON4)                        **
+!*****************************************************************
+ PROCEDURE  Pburncal Pcalflu ;
+ PROCEDURE  assertS ;
+
+ MODULE       RESINI: MACINI: POWER: AFM: DELETE: NEWMAC: INIKIN:
+              KINSOL: TRIVAA: DMAC: GREP: DSET: FLPOW: ABORT: END: UTL: ;
+
+ SEQ_ASCII FBMDT0 ;
+
+ LINKED_LIST  G23D2 INDEX MACRO DEVICE RESEAU ;
+ LINKED_LIST  MACRO2 MACRES FLUX TRACK MACOLD DBASE SYSOLD ;
+ LINKED_LIST  MACRO3 MACRES2 KINET SYSTEMP MACOLDP POWER ;
+
+ INTEGER TPS ;
+
+!* Base de donnees AFM
+!**********************
+ DBASE := FBMDT0 ;
+
+!* Declaration des variables neutronique
+!***************************************
+ REAL Rapport BCinout ;
+ INTEGER Maxreg  ;
+ REAL B1 B2     ;
+ STRING       Refl ;
+ REAL         Xfacc Xfacp ;
+ REAL         Keff ;
+ INTEGER      iter := 1 ;
+ REAL         Precf := 1.0E-7 ;
+ STRING    NamCpo1 ;
+ EVALUATE NamCpo1 := "FULFBMDATA" ;
+ REAL Powi Powf ;
+ EVALUATE Powi := 2061.4 ;
+ REAL REFVALUE DELTA ;
+
+!* Parametres de controle de la simulation
+!*****************************************
+ EVALUATE Rapport := 0.8 ; ! Rapport B1 / B2
+ EVALUATE BCinout := 1.0 ; ! BC entrees = 1.0, BC sorties = 0.0
+ EVALUATE Maxreg  := 27700 ; ! nombre de regions dans la geometrie
+
+!* Appel de la procedure principale
+!**********************************
+G23D2 INDEX MACRO DEVICE RESEAU
+           := Pburncal DBASE :: <<Rapport>> <<BCinout>> <<Maxreg>> ;
+
+ ECHO 'START TRANSIENT CALCULATION' ;
+
+ EVALUATE  Refl  := "REFL-SIGF" ;
+ EVALUATE  Xfacp := 1. ;
+ EVALUATE  Xfacc := 2. ;
+
+!*Age pattern (TA -> Instantaneous)
+  RESEAU := RESINI: RESEAU ::
+      CELL
+                               358  91  39
+  67 162 234
+                   315 113 294 196 248 366
+ 309 109 287 191 245 363
+                59  33 208 346  19 219  61
+  25 202 340  14 213  55  31
+           271 322 123 261 103 156 273 323
+ 118 254  97 149 265 318 122 260
+       333   7 184  51 377 135 334   8 185
+  44 371 127 327   1 179  50 376 133
+        176 145 283 226 301  79 173 142 279
+ 232 307  86 177 147 285 228 303  81
+    241 360  93  40  69 164 236 353  88  37
+  76 171 243 361  95  42  71 166 238 355
+    296 198 250 367 311 110 289 193 246 364
+ 316 116 298 200 252 369 313 112 291 195
+210 348  21 221  63  27 204 342  16 215  57
+  35 211 350  23 223  65  29 206 344  18 217
+262 105 158 275 324 119 256  99 151 267 319
+ 125 263 107 160 276 325 121 258 101 153 269
+378 137 336  10 187  46 372 129 329   3 181
+  53 379 139 338  12 189  48 374 131 331   5
+292  72 167 132 270 224 299  77 172 140 277
+ 216 290  70 165 130 268 222 297  75 170 138
+154 229 345  82  30  66 161 233 351  87  36
+  58 152 227 343  80  28  64 159 231 349  85
+102 281 183 239 356 308 108 286 190 244 362
+ 302 100 280 182 237 354 306 106 284 188 242
+    332   6 207  49  24 201 339  13 212  54
+  17 194 330   4 205  47  22 199 337  11
+     90 144 259 314 117 253  96 148 264 317
+ 111 247  89 143 257 312 115 251  94 146
+        321 375 175  43 370 126 326 380 178
+  38 365 120 320 373 174  41 368 124
+        168 134 272 214 288  68 163 128 266
+ 220 295  74 169 136 274 218 293  73
+             83  32  56 150 225 341  78  26
+  62 157 230 347  84  34  60 155
+                357 300  98 278 180 235 352
+ 305 104 282 186 240 359 304
+                     15 192 328   2 203  45
+  20 197 335   9 209  52
+                                141 255 310
+ 114 249  92
+;
+
+  ECHO " +++++++++++++++++++++++++++++++++++++++ " ;
+  ECHO " +++++++++++++++++++++++++++++++++++++++ " ;
+  ECHO " +++ CONDITION INITIALE INSTANTANNEE +++ " ;
+  ECHO " +++++++++++++++++++++++++++++++++++++++ " ;
+  ECHO " +++++++++++++++++++++++++++++++++++++++ " ;
+
+!* Calcul des sections efficaces
+!***************************************************
+
+  MACRES := AFM: DBASE RESEAU  ::
+   MAP
+   INFOR C6_FBM_CRP0
+   DNAME 1 <<NamCpo1>>
+   EDIT 1
+   NPWF
+   REFT 1 <<NamCpo1>> ;
+
+   UTL: MACRES :: DIR STEP UP GROUP STEP AT 1 DIR IMPR H-FACTOR * ;
+
+  ECHO "RECOVER FUEL AFM MACROLIB" ;
+
+  MACRO2 INDEX := MACINI: INDEX MACRO MACRES ;
+
+!* Calcul du flux pour la repartition convergee du burnup
+!********************************************************
+  FLUX MACOLD TRACK SYSOLD INDEX := Pcalflu INDEX MACRO2 DEVICE G23D2 ::
+    <<iter>>  <<Maxreg>> <<Refl>> <<Xfacc>> <<Xfacp>> <<Precf>> ;
+  assertS FLUX :: "K-EFFECTIVE" 1 1.000627 ;
+
+!* Calcul des puissances (normalisees) de canal et de grappe
+!***********************************************************
+  POWER RESEAU := FLPOW: RESEAU FLUX TRACK INDEX
+                  :: EDIT 0 PTOT <<Powi>> ;
+
+  ECHO " +++++++++++++++++++++++++++++++++++++++ " ;
+  ECHO " +++++++++++++++++++++++++++++++++++++++ " ;
+  ECHO " ++ CINETIQUE ++++++++++++++++++++++++++ " ;
+  ECHO " +++++++++++++++++++++++++++++++++++++++ " ;
+  ECHO " +++++++++++++++++++++++++++++++++++++++ " ;
+
+KINET := INIKIN: MACOLD TRACK SYSOLD FLUX ::
+  EDIT 1
+  NDEL 6
+  BETA 0.000250 0.001562 0.001505 0.002755 0.000830 0.000317
+  LAMBDA 0.0129 0.0311 0.134 0.331 1.26 3.21
+  CHID 1.0 1.0 1.0 1.0 1.0 1.0
+       0.0 0.0 0.0 0.0 0.0 0.0
+  NORM POWER-INI <<Powi>>
+;
+
+EVALUATE TPS := 1 ;
+
+REAL TEMPS := 0.0 ;
+REAL LEVEL := 0.5 ;
+REAL ADJPOS := 1.0 ;
+REAL ADJPOSVAR ;
+REAL TEXIT := 1.0 ;
+REAL DT := 0.1 ;
+EVALUATE ADJPOSVAR := DT TEXIT / ;
+
+!* Boucle de temps
+!************************************************************
+REPEAT
+
+ EVALUATE TEMPS := TEMPS DT + ;
+!*mettre a jour les devices: POSITION DE BC #11 ET DE BL
+ IF TEMPS TEXIT < THEN
+    EVALUATE ADJPOS := ADJPOS ADJPOSVAR - ;
+    EVALUATE LEVEL := 0.5 ;
+ ELSE
+    EVALUATE ADJPOS := 0.0 ;
+    EVALUATE LEVEL := 0.5 ;
+ ENDIF ;
+
+ ECHO " +++++++++++++++++++++++++++++++++++++++ " ;
+ ECHO " CINETIQUE AU TEMPS: " TPS ;
+ ECHO " TEMPS: " TEMPS ;
+ ECHO " ADJUSTER #11 LEVEL:" ADJPOS ;
+ ECHO " +++++++++++++++++++++++++++++++++++++++ " ;
+
+ MACRES2 := AFM: DBASE RESEAU  ::
+   MAP
+   INFOR C6_FBM_CRP0
+   DNAME 1 <<NamCpo1>>
+   EDIT 1
+   NPWF
+   REFT 1 <<NamCpo1>> ;
+
+*   ADJUST DEVICE ADJUSTER AND LZC POSITION            *
+********************************************************
+DEVICE := DSET: DEVICE ::
+  EDIT 1
+  ROD-GROUP 1 LEVEL 1.0 END
+  ROD 11 LEVEL <<ADJPOS>> END
+
+  LZC 1 LEVEL <<LEVEL>> END
+  LZC 2 LEVEL <<LEVEL>> END
+  LZC 3 LEVEL <<LEVEL>> END
+  LZC 4 LEVEL <<LEVEL>> END
+  LZC 5 LEVEL <<LEVEL>> END
+  LZC 6 LEVEL <<LEVEL>> END
+  LZC 7 LEVEL <<LEVEL>> END
+  LZC 8 LEVEL <<LEVEL>> END
+  LZC 9 LEVEL <<LEVEL>> END
+  LZC 10 LEVEL <<LEVEL>> END
+  LZC 11 LEVEL <<LEVEL>> END
+  LZC 12 LEVEL <<LEVEL>> END
+  LZC 13 LEVEL <<LEVEL>> END
+  LZC 14 LEVEL <<LEVEL>> END
+ ;
+
+ MACRO3 INDEX := MACINI: INDEX MACRO MACRES2 ;
+
+ MACOLDP INDEX := NEWMAC: INDEX MACRO3 DEVICE ::
+    EDIT 1  XFAC  <<Xfacc>> ;
+
+ MACRES2 MACRO3 := DELETE: MACRES2 MACRO3 ;
+
+ SYSTEMP := TRIVAA: MACOLDP TRACK :: UNIT ;
+
+ KINET := KINSOL: KINET MACOLDP TRACK SYSTEMP MACOLD SYSOLD ::
+   EDIT 1
+   DELTA <<DT>>
+   SCHEME FLUX THETA 0.8
+   PREC THETA 0.8
+  ;
+
+ !* Calcul des puissances (normalisees) de canal et de grappe
+ !***********************************************************
+ POWER := DELETE: POWER ;
+ POWER RESEAU := FLPOW: RESEAU KINET TRACK INDEX ;
+
+ EVALUATE TPS := TPS 1 + ;
+
+ MACOLD SYSOLD := DELETE: MACOLD SYSOLD ;
+ MACOLD := MACOLDP ;
+ SYSOLD := SYSTEMP ;
+ MACOLDP SYSTEMP := DELETE: MACOLDP SYSTEMP ;
+
+UNTIL TPS 10 > ;
+
+!* Fichier de sortie
+!*******************
+GREP: KINET :: GETVAL E-POW 1 >>Powf<< ;
+EVALUATE REFVALUE := 2.357928E3 ;
+EVALUATE DELTA := Powf REFVALUE - REFVALUE / ABS ;
+IF DELTA 1.0E-3 < THEN
+  PRINT "TEST SUCCESSFUL; DELTA=" DELTA ;
+ELSE
+  PRINT "------------" ;
+  PRINT "TEST FAILURE" ;
+  PRINT "------------" ;
+  PRINT "REFERENCE=" REFVALUE " CALCULATED=" Powf ;
+  ABORT: ;
+ENDIF ;
+ECHO "Initial power=" Powi " Final power=" Powf ;
+ECHO "Successful end of AFMtest." ;
+END: ;
+
+QUIT .