• Author: M.TROISNE, marc.troisne@millennium.fr
• Date: 2016
• Subject: MCNP6 and ADVANTG code launch script
• Keywords: mcnp6, advantG, variance reduction, cadis, monte carlo

• Use 'DEBUT_ADVANTG' and 'FIN_ADVANTG' to write the ADVANTG input
• Skip a line after 'DEBUT_ADVANTG' and after 'FIN_ADVANTG'
• Write the MCNP input after the ADVANTG input

#!/bin/bash
 
export MCNP_HOME=$SEC_HOME/CODES/MCNP6/MCNP_CODE/MCNP6
export MCNP_BIN=$MCNP_HOME/bin/linux/mcnp6_linux_x86_64_omp
export LD_LIBRARY_PATH=$MCNP_HOME
export DATAPATH=$SEC_HOME/lib/ace6
 
export ADVANTG_HOME=$SEC_HOME/CODES/advantg/
export ADVANTG_BIN=$ADVANTG_HOME/bin/advantg
 
# ---------------------------------------------
# check license for user
# ---------------------------------------------
AUTHORIZED=`grep "$USERNAME" $MCNP_HOME/authorized_users.txt | wc -l`
if [[ $AUTHORIZED == 0 ]]; then
  echo "user <"$USERNAME"> is NOT authorized."
  exit -1 
else
  echo "user <"$USERNAME"> is authorized."
fi
 
# ---------------------------------------------
# split file for adv and mcnp calculations
# ---------------------------------------------
 
function getsection()      
{
 
	local arg1=$tmp   # input file
	infilei=${arg1#*=}	
        local arg2=$2     # output file for MCNP
        if [[ $arg2 == 0 ]]; then
           infilen='out_mcnp.'
        else
           infilen=${arg2#*=}
        fi
        # Split accordingly
	ADV=$(grep DEBUT_ADVANTG $infilei | wc -l)  #look for advantg keyword
	if [[ $ADV == 0 ]]; then
		echo nul >> input_adv.adv
		cp $infilei i.mcnp
	else
		csplit -s --prefix output $infilei /'DEBUT_ADVANTG'/+2 /'FIN_ADVANTG'/
		rm output00
		mv output01 input_adv.adv
		csplit -s --prefix out output02 /'FIN_ADVANTG'/+2
		rm out00 output02
		mv out01 i.mcnp
	fi
} 
 
# ---------------------------------------------
# split file for AdvantG
# ---------------------------------------------
 
function run_advantG()
{
	if [[ $ADV == 0 ]]; then
		echo no_advantg
	else	
		csplit -s --prefix output2 input_adv.adv /'mcnp_input'/
		csplit -s --prefix output21 output201 /'mcnp_input'/+1
		echo 'mcnp_input     i.mcnp' >> tmp
		rm input_adv.adv
		cat output200 tmp output2101 >> input_adv.adv  
		rm tmp
		rm output2*
	fi	
}
 
#: <<'END'
# ---------------------------------------------
# process number in "PID" file
# ---------------------------------------------
echo $$ >> PID
 
tmp=$1                 #mcnp input in argument
shift
INP=$*
var="i=i.mcnp"   #new mcnp from advantG
var2="i.mcnp"    #var for dosunix
 
if [[ ( $INP == ${INP##"inp="} ) && ( $INP == ${INP##"n="} ) && ( $INP == ${INP##"i="} ) ]]; then 
        dos2unix $var2
	echo $INP
	getsection $INP
	if [[ $ADV -gt 0 ]]; then
		echo 'ADVANTG...'	
		run_advantG
		$ADVANTG_BIN -v input_adv.adv    #run AdvantG
		mv output/inp ./i.mcnp           #move new input mcnp with modifications from AdvantG
		mv output/wwinp ./		 #move wwinp file for mcnp calculations
		$MCNP_BIN $var $INP &
		PID_MCNP=$!
		echo $PID_MCNP >> PID
		wait $PID_MCNP
	else			 
                $MCNP_BIN $var $INP &
		PID_MCNP=$!
		echo $PID_MCNP >> PID
		wait $PID_MCNP
	fi
else 
        dos2unix $var2
        getsection $INP
        if [[ $ADV -gt 0 ]]; then
		echo 'ADVANTG...'	
		run_advantG
		$ADVANTG_BIN -v input_adv.adv    #run AdvantG
		mv output/inp ./i.mcnp           #move new input mcnp with modifications from AdvantG
		mv output/wwinp ./	         #move wwinp file for mcnp calculations		
                $MCNP_BIN $var $INP &
		PID_MCNP=$!
		echo $PID_MCNP >> PID
		wait $PID_MCNP
	else
                $MCNP_BIN $var $INP &
		PID_MCNP=$!
		echo $PID_MCNP >> PID
		wait $PID_MCNP
	fi	
fi
 
PRINTED_KSEN=`grep "KSEN" *`
if [ "$PRINTED_GRAF" == "" ]; then
  DIR="$( cd "$( dirname "${BASH_SOURCE[0]}" )" && pwd )"
  python $DIR/convert_mcnp_sens_238_.py -d *o
fi
 
# ---------------------------------------------
# Suppress "PID" file
# ---------------------------------------------
if [ -f "PID" ]; then
	rm -f PID
fi
 
#END
 
...