0425
实际需求:
给定一个 temp, top , example.mdp 自动流程化调参
- gmx_mpi make_ndx -f npt2.gro -o index.ndx
- for aa in $(seq 1 $num_replica)
- do
- b=
`sed -n "${aa}p" ${file_dir}/temp`
- mkdir -p ${input_dir}/md${aa}
- sed "s/xxxx/${b}/g" ${file_dir}/example.mdp > ${input_dir}/md${aa}/md.mdp
- gmx_mpi grompp -f ${input_dir}/md${aa}/md.mdp -c ${file_dir}/npt2.gro -p ${file_dir}/topol.top -n ${file_dir}/index.ndx -o ${input_dir}/md${aa}/md.tpr -maxwarn 1
- cd ${input_dir}
- mpirun -np 100 gmx_mpi mdrun -v -deffnm md -multidir md1 md2 md3 md4 md5 md6 md7 md8 md9 md10 md11 md12 md13 md14 md15 md16 md17 md18 md19 md20 -replex 1000
- mpirun -np 100 gmx_mpi mdrun -v -deffnm md -multidir md1 md2 md3 md4 md5 md6 md7 md8 md9 md10 md11 md12 md13 md14 md15 md16 md17 md18 md19 md20 -replex 100
- /home/databank/lywu/apo-proteins/hiv-1/tune_temp/grep_exchange/test
- lywu 76932 0.0 0.0 113288 1452 ? S 01:17 0:00 sh tune_temp.sh
0426
- (base) [lywu@localhost 0425]$ ps 97715
- PID TTY STAT TIME COMMAND
- 97715 pts/125 S 0:00 sh tune_temp.sh