sigcoh.h what is cw

[ameliajo]

assuming that wal2 = 1/M int p(w)/w coth(w/kbT) dw, then we want that

What is in the GA Report = What is in the sigcoh.h code right now
-h^2 tau^2 / 2 * wal2 = - tau^2 * t2 * wint
h^2 / 2 * wal2 = t2 * wint
h^2 / 2 * wal2 = hbar * cw * amsc * wal2 / ( 2 * amu * amsc )
h^2 = hbar * cw / amu
h^2 = h * cw / (amu * 2pi)
cw = 2 pi amu h

So in order for the code to match the report, the above statement has to be true. But it is not. It's very not true. (cw = 0.658e-15, by the way). Interestingly enough, this almost equals h / 2pi

So maybe the equations are wrong? or my math, that's possible too.