fix typo causing unit length errors

This commit is contained in:
erkn 2022-10-11 00:48:44 +02:00
parent 0706a0eca3
commit aeadaf0b2f

View file

@ -802,9 +802,11 @@ void write_mcpl_source_bank(mcpl_outfile_t file_id, bool surf_source_bank)
p.position[1]=_site->r.y;
p.position[2]=_site->r.z;
//mcpl requires that the direction vector is unit length
//which is also the case in openmc
p.direction[0]=_site->u.x;
p.direction[1]=_site->u.y;
p.direction[2]=_site->u.x;
p.direction[2]=_site->u.z;
//mcpl stores kinetic energy in MeV
p.ekin=_site->E*1e-6;