Merge pull request #1 from tsuiyong/master

Added the %VTK card for parameter visualization with ParaView
imronuke · Dec 6, 2023 · d4c80d5 · d4c80d5
2 parents 1a6322d + 5c9f578
commit d4c80d5
Show file tree

Hide file tree

Showing 28 changed files with 8,711 additions and 17 deletions.
diff --git a/docs/images/MOX.png b/docs/images/MOX.png
diff --git a/docs/images/adpres1.png b/docs/images/adpres1.png
diff --git a/docs/images/adpres2.png b/docs/images/adpres2.png
diff --git a/docs/images/fortran.png b/docs/images/fortran.png
diff --git a/docs/images/geom_1.png b/docs/images/geom_1.png
diff --git a/docs/images/geom_2.png b/docs/images/geom_2.png
diff --git a/docs/images/neacrp.png b/docs/images/neacrp.png
diff --git a/docs/images/vertice.png b/docs/images/vertice.png
diff --git a/docs/vtk.md b/docs/vtk.md
@@ -0,0 +1,94 @@
+---
+title: %VTK
+theme: _config.yml
+filename: vtk
+---
+
+# %VTK Card
+
+This card is considered to display the three dimensional physical output parameters by a visual method. Once %VTK card is defined, after running the case, a *.vtk file will be generated, and it can be dealt with the open source code ParaView.
+
+Different from the definition of nodal geometry, %VTK card focus on each point of every node.
+
+For example, the following figure is the discretized mesh of case NEACRP displayed by ParaView.
+
+<img src="./images/neacrp.png" alt="alt text" title="Discretized mesh of case NEACRP" style="zoom: 50%;" />
+
+## %GEOM card of case NEACRP
+```
+%GEOM
+9 9 18               !nx, ny, nz
+10.803   8*21.606    !x-direction assembly size in cm
+  1         8*2      !x-direction assembly divided 2x2
+8*21.606 10.803      !y-direction assembly size in cm
+  8*2       1        !y-direction assembly divided 2x2
+30.  7.7  11.0  15.0  10*30.0  2*12.8  8.0  30. !z-direction size in cm (from bottom to top)
+18*1                 !z-direction assembly divided into 20 cm each
+3                    !np number of planar type
+1  2  15*3  1        !planar assignment (from bottom to top)
+! Planar_Reg 1
+ 1  1  1  1  1  1  1  1  2
+ 1  1  1  1  1  1  1  1  2
+ 1  1  1  1  1  1  1  3  2
+ 1  1  1  1  1  1  1  2  0
+ 1  1  1  1  1  1  3  2  0
+ 1  1  1  1  1  3  2  0  0
+ 1  1  1  1  3  2  0  0  0
+ 1  1  3  2  2  0  0  0  0
+ 2  2  2  0  0  0  0  0  0
+! Planar_Reg 2
+ 4  5  4  5  4  5  4  6  2
+ 5  4  5  4  5  4  6  6  2
+ 4  5  4  5  4  5  6  3  2
+ 5  4  5  4  5  6  6  2  0
+ 4  5  4  5  4  6  3  2  0
+ 5  4  5  6  6  3  2  0  0
+ 4  6  6  6  3  2  0  0  0
+ 6  6  3  2  2  0  0  0  0
+ 2  2  2  0  0  0  0  0  0
+! Planar_Reg 3
+ 4  9  4  9  4  7  4  6  2
+ 9  4  8  4  8  4 11  6  2
+ 4  8  4  8  4  8  6  3  2
+ 9  4  8  4  8 10  6  2  0
+ 4  8  4  8  4  6  3  2  0
+ 7  4  8 10  6  3  2  0  0
+ 4 11  6  6  3  2  0  0  0
+ 6  6  3  2  2  0  0  0  0
+ 2  2  2  0  0  0  0  0  0
+! Boundary conditions (east), (west), (north), (south), (bottom), (top)
+0 2 2 0 0 0
+```
+
+As described in %GEOM card, each node is divided into 2 parts in x and y directions except for the first node. Therefore, the actual node distribution of each planar for calculation is depicted as follows:
+
+<img src="./images/vertice.png" style="zoom:50%;" />
+
+## How to define the %VTK card?
+
+The %VTK card describes the points distribution on one planar which constitute the nodes and they have nothing to do with the material.
+
+The basic rule to define %VTK card is if one point constitute a node, then the value is 1, otherwise, similar to %GEOM card, it is 0, as shown in above figure. For example, the %VTK card of case NEACRP is as follows:
+```
+1   1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	1
+1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	1
+1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	1
+1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	1
+1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	1
+1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	1
+1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	0	0
+1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	0	0
+1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	0	0
+1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	0	0
+1	1	1	1	1	1	1	1	1	1	1	1	1	1	0	0	0	0
+1	1	1	1	1	1	1	1	1	1	1	1	1	1	0	0	0	0
+1	1	1	1	1	1	1	1	1	1	1	1	0	0	0	0	0	0
+1	1	1	1	1	1	1	1	1	1	1	1	0	0	0	0	0	0
+1	1	1	1	1	1	1	1	1	1	0	0	0	0	0	0	0	0
+1	1	1	1	1	1	1	1	1	1	0	0	0	0	0	0	0	0
+1	1	1	1	1	1	0	0	0	0	0	0	0	0	0	0	0	0
+1	1	1	1	1	1	0	0	0	0	0	0	0	0	0	0	0	0
+```
+## Example: Power distribution of Part 2 of MOX benchmark at full power
+
+<img src="/home/cuiyong/Program_files/KOMODO/docs/images/MOX.png" style="zoom: 50%;" />
diff --git a/smpl/VTK/static/NEACRP/A2 b/smpl/VTK/static/NEACRP/A2
@@ -0,0 +1,103 @@
+! CRITICAL BORON CONCENTRATION SEARCH CALCULATION
+! FOR PWR NEACRP BENCHMARK
+%MODE
+BCSEARCH
+
+! CASE CARD
+%CASE
+NEACRP Critical Boron Concentration Search - Case A2
+10.803 CM NODE SIZE
+
+! XSEC CARD
+%XSEC
+FILE ./neacrp_xsec
+
+! GEOMETRY CARD
+%GEOM
+FILE ./neacrp_geom
+
+! CRITICAL BORON SEARCH CARD
+%CBCS
+FILE ./neacrp_cbcs
+
+! FUEL TEMPERATURE CARD
+%FTEM
+FILE ./neacrp_ftem
+
+! MODERATOR TEMPERATURE CARD
+%MTEM
+FILE ./neacrp_mtem
+
+! COOLANT DENSITY CARD
+%CDEN
+FILE ./neacrp_cden
+
+! VTK OUTPUT
+%VTK
+1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	1
+1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	1
+1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	1
+1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	1
+1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	1
+1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	1
+1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	0	0
+1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	0	0
+1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	0	0
+1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	0	0
+1	1	1	1	1	1	1	1	1	1	1	1	1	1	0	0	0	0
+1	1	1	1	1	1	1	1	1	1	1	1	1	1	0	0	0	0
+1	1	1	1	1	1	1	1	1	1	1	1	0	0	0	0	0	0
+1	1	1	1	1	1	1	1	1	1	1	1	0	0	0	0	0	0
+1	1	1	1	1	1	1	1	1	1	0	0	0	0	0	0	0	0
+1	1	1	1	1	1	1	1	1	1	0	0	0	0	0	0	0	0
+1	1	1	1	1	1	0	0	0	0	0	0	0	0	0	0	0	0
+1	1	1	1	1	1	0	0	0	0	0	0	0	0	0	0	0	0
+
+! CONTROL CARD
+%CROD
+7    228                      ! Number of CR banks and max number of steps
+37.7 1.5942237                ! Zero step pos. (cm) and cm/step (total 228 steps)
+100.00  200.00  100.00  200.00  200.00  200.00  200.00           ! CR Bank pos.
+ 1  0  2  0  0  0  3  0  0
+ 0  4  0  0  0  6  0  0  0
+ 2  0  5  0  6  0  6  0  0
+ 0  0  0  4  0  0  0  0  0
+ 0  0  6  0  7  0  0  0  0
+ 0  6  0  0  0  0  0  0  0
+ 3  0  6  0  0  0  0  0  0
+ 0  0  0  0  0  0  0  0  0
+ 0  0  0  0  0  0  0  0  0
+! CX changes
+!  sigtr        siga        nu*sigf        kappa*sigf      sigs_g1   sigs_g2
+ 3.73220E-03  2.47770E-03  -1.02786E-04  -1.21448E-15  0.0  -3.19253E-03
+-2.19926E-02  2.55875E-02  -2.82319E-03  -3.70238E-14  0.0   0.00000E+00    !COMP 1
+ 3.73220E-03  2.47770E-03  -1.02786E-04  -1.21448E-15  0.0  -3.19253E-03
+-2.19926E-02  2.55875E-02  -2.82319E-03  -3.70238E-14  0.0   0.00000E+00    !COMP 2
+ 3.73220E-03  2.47770E-03  -1.02786E-04  -1.21448E-15  0.0  -3.19253E-03
+-2.19926E-02  2.55875E-02  -2.82319E-03  -3.70238E-14  0.0   0.00000E+00    !COMP 3
+ 3.73220E-03  2.47770E-03  -1.02786E-04  -1.21448E-15  0.0  -3.19253E-03
+-2.19926E-02  2.55875E-02  -2.82319E-03  -3.70238E-14  0.0   0.00000E+00    !COMP 4
+ 3.73220E-03  2.47770E-03  -1.02786E-04  -1.21448E-15  0.0  -3.19253E-03
+-2.19926E-02  2.55875E-02  -2.82319E-03  -3.70238E-14  0.0   0.00000E+00    !COMP 5
+ 3.74092E-03  2.42926E-03  -1.22634E-04  -1.47557E-15  0.0  -3.14239E-03
+-1.67503E-02  2.56478E-02  -3.28086E-03  -4.30444E-14  0.0   0.00000E+00    !COMP 6
+ 3.73220E-03  2.47770E-03  -1.02786E-04  -1.21448E-15  0.0  -3.19253E-03
+-2.19926E-02  2.55875E-02  -2.82319E-03  -3.70238E-14  0.0   0.00000E+00    !COMP 7
+ 3.73220E-03  2.47770E-03  -1.02786E-04  -1.21448E-15  0.0  -3.19253E-03
+-2.19926E-02  2.55875E-02  -2.82319E-03  -3.70238E-14  0.0   0.00000E+00    !COMP 8
+ 3.73220E-03  2.47770E-03  -1.02786E-04  -1.21448E-15  0.0  -3.19253E-03
+-2.19926E-02  2.55875E-02  -2.82319E-03  -3.70238E-14  0.0   0.00000E+00    !COMP 9
+ 3.73220E-03  2.47770E-03  -1.02786E-04  -1.21448E-15  0.0  -3.19253E-03
+-2.19926E-02  2.55875E-02  -2.82319E-03  -3.70238E-14  0.0   0.00000E+00    !COMP 10
+ 3.73220E-03  2.47770E-03  -1.02786E-04  -1.21448E-15  0.0  -3.19253E-03
+-2.19926E-02  2.55875E-02  -2.82319E-03  -3.70238E-14  0.0   0.00000E+00    !COMP 11
+1
+
+! THERMAL-HYDRAULIC CARD
+%THER
+100.                                   ! Percent power in %
+693.75e6                               ! Reactor thermal power for quarter geometry in Watt
+559.15  82.12102                       ! Inlet coolant temp. (Kelvin) and Fuel Assembly Mass flow rate (kg/s)
+4.11950E-03  6.8E-05  5.71E-04  1.2655E-2 ! Fuel meat rad., gap thickness, cladding thickness and pin pitch (m)
+264  25                                ! Number of fuel pin and guide tubes
+0.019                                  ! FRACTION OF HEAT DEPOSITED IN COOLANT
diff --git a/smpl/VTK/static/NEACRP/B2 b/smpl/VTK/static/NEACRP/B2
@@ -0,0 +1,107 @@
+! CRITICAL BORON CONCENTRATION SEARCH CALCULATION
+! FOR PWR NEACRP BENCHMARK
+%MODE
+BCSEARCH
+
+! CASE CARD
+%CASE
+NEACRP Critical Boron Concentration Search - Case B2
+10.803 CM NODE SIZE
+
+%ITER
+500 2 1.e-5 1.e-5 5 30 20 60
+
+! XSEC CARD
+%XSEC
+FILE ./neacrp_xsec
+
+! GEOMETRY CARD
+%GEOM
+FILE ./neacrp_geom
+
+! CRITICAL BORON SEARCH CARD
+%CBCS
+FILE ./neacrp_cbcs
+
+! FUEL TEMPERATURE CARD
+%FTEM
+FILE ./neacrp_ftem
+
+! MODERATOR TEMPERATURE CARD
+%MTEM
+FILE ./neacrp_mtem
+
+! COOLANT DENSITY CARD
+%CDEN
+FILE ./neacrp_cden
+
+! VTK OUTPUT
+%VTK
+1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	1
+1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	1
+1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	1
+1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	1
+1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	1
+1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	1
+1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	0	0
+1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	0	0
+1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	0	0
+1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	1	0	0
+1	1	1	1	1	1	1	1	1	1	1	1	1	1	0	0	0	0
+1	1	1	1	1	1	1	1	1	1	1	1	1	1	0	0	0	0
+1	1	1	1	1	1	1	1	1	1	1	1	0	0	0	0	0	0
+1	1	1	1	1	1	1	1	1	1	1	1	0	0	0	0	0	0
+1	1	1	1	1	1	1	1	1	1	0	0	0	0	0	0	0	0
+1	1	1	1	1	1	1	1	1	1	0	0	0	0	0	0	0	0
+1	1	1	1	1	1	0	0	0	0	0	0	0	0	0	0	0	0
+1	1	1	1	1	1	0	0	0	0	0	0	0	0	0	0	0	0
+
+! CONTROL CARD
+%CROD
+7   228                       ! Number of CR banks and max number of steps
+37.7 1.5942237                ! Zero step pos. (cm) and cm/step (total 228 steps)
+150.00  200.00  150.00  200.00  200.00  200.00  200.00           ! CR Bank pos.
+ 1  0  2  0  0  0  3  0  0
+ 0  4  0  0  0  6  0  0  0
+ 2  0  5  0  6  0  6  0  0
+ 0  0  0  4  0  0  0  0  0
+ 0  0  6  0  7  0  0  0  0
+ 0  6  0  0  0  0  0  0  0
+ 3  0  6  0  0  0  0  0  0
+ 0  0  0  0  0  0  0  0  0
+ 0  0  0  0  0  0  0  0  0
+! CX changes
+!  sigtr        siga        nu*sigf        kappa*sigf      sigs_g1   sigs_g2
+ 3.73220E-03  2.47770E-03  -1.02786E-04  -1.21448E-15  0.0  -3.19253E-03
+-2.19926E-02  2.55875E-02  -2.82319E-03  -3.70238E-14  0.0   0.00000E+00    !COMP 1
+ 3.73220E-03  2.47770E-03  -1.02786E-04  -1.21448E-15  0.0  -3.19253E-03
+-2.19926E-02  2.55875E-02  -2.82319E-03  -3.70238E-14  0.0   0.00000E+00    !COMP 2
+ 3.73220E-03  2.47770E-03  -1.02786E-04  -1.21448E-15  0.0  -3.19253E-03
+-2.19926E-02  2.55875E-02  -2.82319E-03  -3.70238E-14  0.0   0.00000E+00    !COMP 3
+ 3.73220E-03  2.47770E-03  -1.02786E-04  -1.21448E-15  0.0  -3.19253E-03
+-2.19926E-02  2.55875E-02  -2.82319E-03  -3.70238E-14  0.0   0.00000E+00    !COMP 4
+ 3.73220E-03  2.47770E-03  -1.02786E-04  -1.21448E-15  0.0  -3.19253E-03
+-2.19926E-02  2.55875E-02  -2.82319E-03  -3.70238E-14  0.0   0.00000E+00    !COMP 5
+ 3.74092E-03  2.42926E-03  -1.22634E-04  -1.47557E-15  0.0  -3.14239E-03
+-1.67503E-02  2.56478E-02  -3.28086E-03  -4.30444E-14  0.0   0.00000E+00    !COMP 6
+ 3.73220E-03  2.47770E-03  -1.02786E-04  -1.21448E-15  0.0  -3.19253E-03
+-2.19926E-02  2.55875E-02  -2.82319E-03  -3.70238E-14  0.0   0.00000E+00    !COMP 7
+ 3.73220E-03  2.47770E-03  -1.02786E-04  -1.21448E-15  0.0  -3.19253E-03
+-2.19926E-02  2.55875E-02  -2.82319E-03  -3.70238E-14  0.0   0.00000E+00    !COMP 8
+ 3.73220E-03  2.47770E-03  -1.02786E-04  -1.21448E-15  0.0  -3.19253E-03
+-2.19926E-02  2.55875E-02  -2.82319E-03  -3.70238E-14  0.0   0.00000E+00    !COMP 9
+ 3.73220E-03  2.47770E-03  -1.02786E-04  -1.21448E-15  0.0  -3.19253E-03
+-2.19926E-02  2.55875E-02  -2.82319E-03  -3.70238E-14  0.0   0.00000E+00    !COMP 10
+ 3.73220E-03  2.47770E-03  -1.02786E-04  -1.21448E-15  0.0  -3.19253E-03
+-2.19926E-02  2.55875E-02  -2.82319E-03  -3.70238E-14  0.0   0.00000E+00    !COMP 11
+1
+
+! THERMAL-HYDRAULIC CARD
+%THER
+100.                                   ! Percent power in %
+693.75e6                               ! Reactor thermal power for quarter geometry in Watt
+559.15  82.12102                       ! Inlet coolant temp. (Kelvin) and Fuel Assembly Mass flow rate (kg/s)
+4.11950E-03  6.8E-05  5.71E-04  1.2655E-2 ! Fuel meat rad., gap thickness, cladding thickness and pin pitch (m)
+264  25                                ! Number of fuel pin and guide tubes
+0.019                                  ! FRACTION OF HEAT DEPOSITED IN COOLANT
+1
diff --git a/smpl/VTK/static/NEACRP/neacrp_cbcs b/smpl/VTK/static/NEACRP/neacrp_cbcs
@@ -0,0 +1,25 @@
+1200.2  ! Boron concentration ref. in ppm
+! CX change per unit ppm change of Boron concentration
+!  sigtr          siga         nu*sigf    kappa*sigf  sigs_g1   sigs_g2
+ 6.11833E-08  1.87731E-07   0.00000E+00   0.00000E+00  0.0   7.91457E-10
+ 5.17535E-06  1.02635E-05   0.00000E+00   0.00000E+00  0.0   0.00000E+00    !COMP 1
+ 0.00000E+00  0.00000E+00   0.00000E+00   0.00000E+00  0.0   0.00000E+00
+ 7.76184E-04  8.44695E-05   0.00000E+00   0.00000E+00  0.0   0.00000E+00    !COMP 2
+ 0.00000E+00  0.00000E+00   0.00000E+00   0.00000E+00  0.0   0.00000E+00
+ 7.76184E-04  8.44695E-05   0.00000E+00   0.00000E+00  0.0   0.00000E+00    !COMP 3
+ 3.47809E-08  1.28505E-07  -1.12099E-09  -1.76188E-20  0.0  -1.08590E-07
+-9.76510E-06  7.08807E-06  -2.43045E-06  -3.19085E-17  0.0   0.00000E+00    !COMP 4
+ 3.53826E-08  1.26709E-07  -1.67880E-09  -2.49965E-20  0.0  -1.06951E-07
+-8.50169E-06  6.82311E-06  -2.72445E-06  -3.57680E-17  0.0   0.00000E+00    !COMP 5
+ 3.59838E-08  1.24986E-07  -2.21038E-09  -3.20225E-20  0.0  -1.05374E-07
+-7.46251E-06  6.59798E-06  -2.95883E-06  -3.88451E-17  0.0   0.00000E+00    !COMP 6
+ 3.37806E-08  1.19869E-07  -1.71323E-09  -2.49965E-20  0.0  -1.00873E-07
+-6.73744E-06  6.29310E-06  -2.55359E-06  -3.35223E-17  0.0   0.00000E+00    !COMP 7
+ 3.32495E-08  1.17585E-07  -1.72421E-09  -2.54896E-20  0.0  -9.88578E-08
+-6.19725E-06  6.11904E-06  -2.48880E-06  -3.26704E-17  0.0   0.00000E+00    !COMP 8
+ 3.27201E-08  1.15319E-07  -1.73502E-09  -2.56049E-20  0.0  -9.68489E-08
+-5.68220E-06  5.94711E-06  -2.42240E-06  -3.17976E-17  0.0   0.00000E+00    !COMP 9
+ 3.43859E-08  1.18186E-07  -2.24335E-09  -3.20225E-20  0.0  -9.93312E-08
+-5.86898E-06  6.08443E-06  -2.77657E-06  -3.64509E-17  0.0   0.00000E+00    !COMP 10
+ 3.38559E-08  1.15917E-07  -2.25369E-09  -3.24873E-20  0.0  -9.73291E-08
+-5.38345E-06  5.91697E-06  -2.70780E-06  -3.55476E-17  0.0   0.00000E+00    !COMP 11
diff --git a/smpl/VTK/static/NEACRP/neacrp_cden b/smpl/VTK/static/NEACRP/neacrp_cden
@@ -0,0 +1,26 @@
+! CDEN CARD FILE
+0.7464  0.7125  ! Average Coolant Density and Coolant Density Ref. in g/cm3
+! CX change per Coolant Density Changes
+!  sigtr          siga         nu*sigf  kappa*sigf sigs_g1   sigs_g2
+7.45756E-02  2.07688E-04  0.00000E+00  0.00000E+00  0.0  3.71310E-02
+5.33634E-01  7.58421E-03  0.00000E+00  0.00000E+00  0.0  0.00000E+00    !COMP 1
+0.00000E+00  0.00000E+00  0.00000E+00  0.00000E+00  0.0  0.00000E+00
+0.00000E+00  0.00000E+00  0.00000E+00  0.00000E+00  0.0  0.00000E+00    !COMP 2
+0.00000E+00  0.00000E+00  0.00000E+00  0.00000E+00  0.0  0.00000E+00
+0.00000E+00  0.00000E+00  0.00000E+00  0.00000E+00  0.0  0.00000E+00    !COMP 3
+1.35665E-01  1.55185E-03  9.20694E-04  1.02392E-14  0.0  2.93195E-02
+9.92628E-01  2.52662E-02  2.47746E-02  3.25255E-13  0.0  0.00000E+00    !COMP 4
+1.35748E-01  1.61491E-03  9.64160E-04  1.08141E-14  0.0  2.92696E-02
+9.81985E-01  2.86667E-02  3.14993E-02  4.13542E-13  0.0  0.00000E+00    !COMP 5
+1.35827E-01  1.68015E-03  1.01410E-03  1.14771E-14  0.0  2.92154E-02
+9.72267E-01  3.19571E-02  3.81097E-02  5.00328E-13  0.0  0.00000E+00    !COMP 6
+1.31033E-01  1.68397E-03  9.81951E-04  1.08141E-14  0.0  2.82489E-02
+9.34697E-01  3.14240E-02  3.51588E-02  4.61715E-13  0.0  0.00000E+00    !COMP 7
+1.29379E-01  1.71972E-03  9.88437E-04  1.11322E-14  0.0  2.78895E-02
+9.18171E-01  3.24715E-02  3.63251E-02  4.77078E-13  0.0  0.00000E+00    !COMP 8
+1.27682E-01  1.74989E-03  9.95175E-04  1.12209E-14  0.0  2.75202E-02
+9.01293E-01  3.35945E-02  3.74499E-02  4.91900E-13  0.0  0.00000E+00    !COMP 9
+1.31116E-01  1.75528E-03  1.03522E-03  1.14771E-14  0.0  2.81877E-02
+9.24925E-01  3.49853E-02  4.20693E-02  5.52387E-13  0.0  0.00000E+00    !COMP 10
+1.29463E-01  1.79499E-03  1.04291E-03  1.18534E-14  0.0  2.78259E-02
+9.08456E-01  3.61032E-02  4.33215E-02  5.68857E-13  0.0  0.00000E+00    !COMP 11
diff --git a/smpl/VTK/static/NEACRP/neacrp_ftem b/smpl/VTK/static/NEACRP/neacrp_ftem
@@ -0,0 +1,26 @@
+! FTEM CARD FILE
+891.19  891.45  ! Average Fuel Temperature and fuel temperature Ref. in Kelvin
+! CX change per Fuel Temperature Changes
+!  sigtr          siga         nu*sigf     kappa*sigf sigs_g1   sigs_g2
+ 0.00000E+00   0.00000E+00   0.00000E+00   0.00000E+00  0.0   0.00000E+00
+ 0.00000E+00   0.00000E+00   0.00000E+00   0.00000E+00  0.0   0.00000E+00    !COMP 1
+ 0.00000E+00   0.00000E+00   0.00000E+00   0.00000E+00  0.0   0.00000E+00
+ 0.00000E+00   0.00000E+00   0.00000E+00   0.00000E+00  0.0   0.00000E+00    !COMP 2
+ 0.00000E+00   0.00000E+00   0.00000E+00   0.00000E+00  0.0   0.00000E+00
+ 0.00000E+00   0.00000E+00   0.00000E+00   0.00000E+00  0.0   0.00000E+00    !COMP 3
+-3.09197E-05   3.49709E-05   6.40134E-07   7.15412E-18  0.0  -2.75536E-05
+-0.000137292  -3.71806E-05  -5.63037E-05  -7.39188E-16  0.0   0.00000E+00    !COMP 4
+-3.08607E-05   3.51798E-05   9.97431E-07   1.18685E-17  0.0  -2.76766E-05
+-0.000117481  -3.77039E-05  -6.04155E-05  -7.93170E-16  0.0   0.00000E+00    !COMP 5
+-3.09165E-05   3.53841E-05   1.41847E-06   1.74269E-17  0.0  -2.78390E-05
+-0.000101337  -3.77558E-05  -0.000063096  -8.28363E-16  0.0   0.00000E+00    !COMP 6
+-3.13746E-05   3.48699E-05   9.45431E-07   1.18685E-17  0.0  -2.73550E-05
+-0.000108271  -3.72748E-05  -5.79662E-05  -7.60849E-16  0.0   0.00000E+00    !COMP 7
+-3.15503E-05   3.47274E-05   9.26078E-07   1.08935E-17  0.0  -2.72381E-05
+-0.000105521  -3.71808E-05  -5.71108E-05  -7.49575E-16  0.0   0.00000E+00    !COMP 8
+-3.17281E-05   3.46026E-05   9.05802E-07   1.06166E-17  0.0  -2.71169E-05
+-0.000102525  -3.70201E-05  -5.61543E-05  -7.36969E-16  0.0   0.00000E+00    !COMP 9
+-3.14192E-05   3.50637E-05   1.35642E-06   1.74269E-17  0.0  -2.75049E-05
+-9.38886E-05  -3.71403E-05  -6.05052E-05  -7.94252E-16  0.0   0.00000E+00    !COMP 10
+-3.15908E-05   3.49119E-05   1.33336E-06   1.62769E-17  0.0  -2.73835E-05
+-9.17126E-05  -3.69909E-05  -5.96284E-05  -7.82716E-16  0.0   0.00000E+00    !COMP 11