*----
*  TEST CASE TCW612
*  WIMS-AECL 89 GROUPS LIBRARY FILE E6MLIB
*  GENERATE A LIBRARY USING INFO AND OTHER OPTIONS
*
*  REF: None
*
*----
*  Define variables and initialize
*    Coolant properties
*    a) Input
*       TempCool    = Coolant temperature (K)
*       Purity      = D2/(D2+H1) Weight % ratio in Coolant
*    b) Output
*       DensCool    = Coolant Density (g/cm**3)
*       WH1C        = H1 Weight % in Coolant
*       WD2C        = D2 Weight % in Coolant
*       WO16C       = O16 Weight % in Coolant
*    Fuel properties
*    a) Input
*       TempFuel    = Fuel temperature  (K)
*       Enrichment  = U235/(U235+U238) Weight % ratio in Fuel
*       DensFuel    = Fuel Density (g/cm**3)
*    b) Output
*       WU235F      = U235 Weight % in Fuel
*       WU238F      = U238 Weight % in Fuel
*       WO16F       = O16 Weight % in Fuel
*----
REAL
  TempCool Purity TempFuel Enrichment DensFuel   :=
  560.66   99.95  941.29   0.72       10.437501 ;
REAL
 WH1C   WD2C   WO16C  DensCool
 WU235F WU238F WO16F ;
*----
*  Define STRUCTURES and MODULES used
*----
LINKED_LIST
  LIBRARY ISOT ;
MODULE
  LIB: INFO: END: ;
PROCEDURE assertS ;
*----
*  Get Coolant properties
*----
ECHO
  "Input - Coolant temperature (K)             " TempCool ;
ECHO
  "Input - D2/(D2+H1) Weight % ratio in Coolant" Purity ;
INFO: ::
  TMP: <<TempCool>> K
  PUR: <<Purity>>   WGT%
  CALC DENS WATER >>DensCool<<
  LIB: WIMSAECL FIL: E6MLIB
  ISO: 3 H1H2O    DD2O     O16
  CALC WGT% D2O >>WH1C<< >>WD2C<< >>WO16C<<
  ;
ECHO
  "Output - Coolant Density (g/cm**3)" DensCool ;
ECHO
  "Output - H1 Weight % in Coolant   " WH1C ;
ECHO
  "Output - D2 Weight % in Coolant   " WD2C ;
ECHO
  "Output - O16 Weight % in Coolant  " WO16C ;
*----
*  Get Fuel properties
*----
ECHO
  "Input - Fuel temperature (K)                   " TempFuel ;
ECHO
  "Input - U235/(U235+U238) Weight % ratio in Fuel" Enrichment ;
ECHO
  "Input - Fuel Density (g/cm**3)                 " DensFuel ;
INFO: ::
  ENR: <<Enrichment>> WGT%
  LIB: WIMSAECL FIL: E6MLIB
  ISO: 3 U235     U238     O16
  CALC WGT% UO2 >>WU235F<< >>WU238F<< >>WO16F<< 
  ;
ECHO
  "Output - U235 Weight % in Fuel" WU235F ;
ECHO
  "Output - U238 Weight % in Fuel" WU238F ;
ECHO
  "Output - O16 Weight % in Fuel " WO16F ;
*----
*  Microscopic cross sections from file E6MLIB format WIMS-AECL
*----
LIBRARY := LIB: ::
  NMIX 5 CTRA WIMS
  MIXS LIB: WIMSAECL FIL: E6MLIB
  MIX  1  <<TempCool>> <<DensCool>>  O16      = O16      <<WO16C>>
    D2D2O    = DD2O     <<WD2C>>     H1H2O    = H1H2O    <<WH1C>>
  MIX  2  <<TempFuel>> <<DensCool>>  O16      = O16      <<WO16C>>
    D2D2O    = DD2O     <<WD2C>>     H1H2O    = H1H2O    <<WH1C>>
  MIX  3  <<TempFuel>> <<DensFuel>>  O16      = O16      <<WO16F>>
    U235     = U235     <<WU235F>> 1
    U238     = U238     <<WU238F>> 1
  MIX  4 COMB 2 0.5 0 0.5
  MIX  5 COMB 2 0.1 3 0.9
  ;
ISOT := LIBRARY :: STEP UP ISOTOPESLIST STEP AT 17 ; (*U238    0005*)
assertS ISOT :: NG 35 4.984570E+00 ;
ECHO "test TCW612 completed" ;
END: ; 
QUIT "LIST" .