Update 4 files

- /15_Sampler/solutions/lwr_pin_r2.inp
- /15_Sampler/solutions/pincell.inp
- /15_Sampler/starters/lwr_pin_r2_starter.inp
- /15_Sampler/starters/pincellsampler.inp

15_Sampler/solutions/lwr_pin_r2.inp

+=sampler
+
+' ------ SAMPLER parameters -------
+
+read parameters
+  n_samples = 500
+  library = "xn56v7.1"
+  perturb_xs = yes
+  force_run = no
+  run_cases = no
+end parameters
+
+' ------ SAMPLER responses -------
+
+read response[lambda]
+  type = grep
+  regexp = ":xsdrn.lambda:"
+end response
+
+' ------ SAMPLER analysis -------
+
+read analysis[analysis_id]
+  type    = r2
+  targets = c1:lambda       end
+  sources = h-1 h-2 o-16 o-17 cr-50 cr-52 cr-53 cr-54 fe-54 fe-56 fe-57 fe-58 zr-90 zr-91 zr-92 zr-94 zr-96 u-234 sn-112 sn-114 sn-115 sn-116 sn-117 sn-118 sn-119 sn-120 sn-122 sn-124 hf-174 hf-176 hf-177 hf-178 hf-179 hf-180 u-235 u-236 u-238 end
+end analysis
+
+' ------ T-XSDRN input -------
+
+read case[c1]
+sequence=t-xsdrn parm=(bonami)
+lwr pin
+xn56v7.1
+
+read comp
+ uo2     1 den=10.257 1.0 800 
+              92234 0.0263 
+              92235 3.1 
+              92236 0.0143 
+              92238 96.8594 end
+ zirc4   2 den=6.56      1.0    585 end
+ h2o     3 den=0.70045   1.0    585 end
+end comp
+
+read model
+lwr pin
+
+read parameters
+  sn=16
+  epsglobal=1e-5
+end parameters
+
+read materials
+  mix=1 com='fuel' end
+  mix=2 com='clad' end
+  mix=3 com='moderator' end
+end materials
+
+read geom
+  geom=cylinder
+  rightBC=white
+  zoneIDs 1 2 3 end zoneids
+  zoneDimensions 0.4100 0.4750 0.7406 end zoneDimensions
+  zoneIntervals 3r10 end zoneIntervals
+end geom
+
+end model
+end sequence
+end case
+
+end
+

15_Sampler/solutions/pincell.inp

+=%sampler
+read parameters
+  n_samples=10
+  library="v7-56"
+  perturb_xs=yes
+  perturb_decay=yes
+  perturb_yields=yes
+  run_cases=yes
+end parameters
+
+read case [c1]
+sequence=polaris
+title "W17x17 pincell 80 GwD/MTU 8 wt%"
+lib "v7-56"
+system PWR
+% increase pitch to account for GT and gap
+opt ESSM NumAzim=1 NumPolar=1
+opt KEFF PnOrder=0 TransportCorrection=no
+         EigSolver='power' PartitionUpscatter=no
+         UpscatterSolver='gauss_seidel' UpscatterInnerSolver='source'
+
+geometry W17 : ASSM 1 1.32077
+% compositions and materials
+comp c_f50 : UOX enr=8.0
+mat FUEL.1 : c_f50 dens=10.26
+pin F : 0.4096 0.418 0.475
+      : FUEL.1 GAP.1 CLAD.1
+mesh FUEL : nr=1
+% state variables
+state ALL : temp=590
+COOL : dens=0.68 boron=800
+FUEL : temp=900
+CLAD : temp=700 TUBE : temp=580
+pow 40
+bu 1 5 10 20 30 40 50 60 70 80
+end sequence
+
+sequence=shell
+mv ft71f001 discharge_iso
+end sequence
+
+sequence=origen
+case
+{
+  lib{ file=end7dec }
+  mat{ load{ file=discharge_iso pos=150 } }
+  time=[10L 1e-5 5000]
+  save=yes
+}
+end sequence
+
+sequence=opus
+data="ft71f001"
+library="end7dec"
+nrank=40
+units=watts
+end sequence
+
+sequence=shell
+rm ft71f001
+end sequence
+
+sequence=origen
+case
+{
+  lib{ file=end7dec }
+  mat{ load{ file=discharge_iso pos=154 } }
+  time=[10L 1e-5 5000]
+  save=yes
+}
+end sequence
+
+sequence=opus
+data="ft71f001"
+library="end7dec"
+nrank=40
+units=watts
+end sequence
+end case
+
+'response extracted from *.samplerfiles/c1_pert_*0.plt
+read response[dh]
+type=opus_plt
+ndataset=0
+nuclides=total end
+end response
+
+'response extracted from *.samplerfiles/c1_pert_*1.plt
+read response[hbdh]
+type=opus_plt
+ndataset=1
+nuclides=total end
+end response
+
+read analysis[pearson_steps2]
+  type    = pearson_corr
+  targets= c1:dh.total(last)  end
+  sources = c1:hbdh.total(last)  end
+end analysis
+read analysis[pearson_set]
+  type = correlation_matrix
+'  targets= c1:dh.total(all) end
+  sources = c1:hbdh.total(all)  end
+end analysis
+
+end
\ No newline at end of file

15_Sampler/starters/lwr_pin_r2_starter.inp

+=sampler
+
+' ------ SAMPLER parameters -------
+
+read parameters
+  n_samples = ???
+  library = "xn56v7.1"
+  perturb_xs = yes
+  force_run = no
+  run_cases = no
+end parameters
+
+' ------ SAMPLER responses -------
+
+
+' ------ SAMPLER analysis -------
+
+
+' ------ T-XSDRN input -------
+
+read case[c1]
+sequence=t-xsdrn parm=(bonami)
+lwr pin
+xn56v7.1
+
+read comp
+ uo2     1 den=10.257 1.0 800 
+              92234 0.0263 
+              92235 3.1 
+              92236 0.0143 
+              92238 96.8594 end
+ zirc4   2 den=6.56      1.0    585 end
+ h2o     3 den=0.70045   1.0    585 end
+end comp
+
+read model
+lwr pin
+
+read parameters
+  sn=16
+  epsglobal=1e-5
+end parameters
+
+read materials
+  mix=1 com='fuel' end
+  mix=2 com='clad' end
+  mix=3 com='moderator' end
+end materials
+
+read geom
+  geom=cylinder
+  rightBC=white
+  zoneIDs 1 2 3 end zoneids
+  zoneDimensions 0.4100 0.4750 0.7406 end zoneDimensions
+  zoneIntervals 3r10 end zoneIntervals
+end geom
+
+end model
+end sequence
+end case
+
+end
+

15_Sampler/starters/pincellsampler.inp

+=%sampler
+read parameters
+  n_samples=
+  library="v7-56"
+  perturb_xs=
+  perturb_decay=
+  perturb_yields=
+  run_cases=
+end parameters
+
+read case [c1]
+sequence=polaris
+title "W17x17 pincell 80 GwD/MTU 8 wt%"
+lib
+system PWR
+% increase pitch to account for GT and gap
+opt ESSM NumAzim=1 NumPolar=1
+opt KEFF PnOrder=0 TransportCorrection=no
+         EigSolver='power' PartitionUpscatter=no
+         UpscatterSolver='gauss_seidel' UpscatterInnerSolver='source'
+
+geometry W17 : ASSM 1 1.32077
+% compositions and materials
+comp c_f50 : UOX enr=8.0
+mat FUEL.1 : c_f50 dens=10.26
+pin F : 0.4096 0.418 0.475
+      : FUEL.1 GAP.1 CLAD.1
+mesh FUEL : nr=1
+% state variables
+state ALL : temp=590
+COOL : dens=0.68 boron=800
+FUEL : temp=900
+CLAD : temp=700 TUBE : temp=580
+pow 40
+bu 1 5 10 20 30 40 50 60 70 80
+end sequence
+
+'decay 5000 days
+sequence=shell
+mv ft71f001 discharge_iso
+end sequence
+
+sequence=origen
+case
+{
+  lib{ file=end7dec }
+  mat{ load{ file=discharge_iso pos=150 } }
+  time=[10L 1e-5 5000]
+  save=yes
+}
+end sequence
+
+sequence=opus
+data="ft71f001"
+library="end7dec"
+nrank=40
+units=watts
+end sequence
+
+'decay 5000 days
+sequence=shell
+rm ft71f001
+end sequence
+
+sequence=origen
+case
+{
+  lib{ file=end7dec }
+  mat{ load{ file=discharge_iso pos=154 } }
+  time=[10L 1e-5 5000]
+  save=yes
+}
+end sequence
+
+sequence=opus
+data="ft71f001"
+library="end7dec"
+nrank=40
+units=watts
+end sequence
+end case
+
+'response extracted from *.samplerfiles/c1_pert_*0.plt
+read response[dh]
+type=opus_plt
+ndataset=0
+nuclides=total end
+end response
+
+'response extracted from *.samplerfiles/c1_pert_*1.plt
+read response[hbdh]
+
+end response
+
+'investigate correlation between dh and hbdh at 5000 days values
+read analysis[pearson_steps2]
+  type    = pearson_corr
+  targets= c1:dh.total(last)  end
+  sources =
+end analysis
+
+'is decay heat at different decay related?
+read analysis[pearson_set]
+  type = correlation_matrix
+  sources =
+end analysis
+
+end