[FLASH-USERS] normalization confusion

For my one-fluid MHD model I am trying to implement a (supersonic) flow
with
rho=2.5e-20 kg/m^3 (or n=1.5e7 m^-3 since only considering protons)
v = 5.0e5 m/s
T = 2.0e5 K
The pressure of this flow should be
p = n*k_B*T = 1.5e7*1.38e-23*2.0e5 = 4.1e-11 Pa

For the normalization I use reference values rho_0 = rho, v_0 = v and
r_0 = Earth radius. I have found a reference saying that pressure should be normalized with
p' = p/(rho_0*v_0^2)
This gives me p' = 6.6e-3

I can not find a reference explaining how temperature should be
normalized, but to acheive p'=6.6e-3 I can set t'=8.0e-11 in FLASH. This
would be in accordance with Eq (10.5) from FLASH UG:
p = N_a*k_B*rho*T
if using k_B = 1.38e-16 cm^2 g K^-1 s^-2, rho=1 and N_a=6.022e23

I have 4 questions for which I would appreciate any information:
a) Is there any documentation explaining the normalizations used in FLASH?
b) Shouldn't it be possible to use the relation p=n*k_B*T to calculate
what temperature I should set for my flow? I must be missing some
trick.
c) How come the EOS does not use k_B in SI-units?
d) (the short story) How should T be normalized for FLASH?

Thanks in advance
/Andreas

-- 
Andreas Ekenbäck                      andreas@irf.se
Ph.D. Student
Swedish Institute of Space Physics    http://www.irf.se
Tel direct: +46 (0)980 79113