: dens0,aln1,aln2      ! dens[1e19 m-3] = dens0*(1 - psin**aln1)**aln2
: te0,dte0,alt1,alt2   ! temp[keV] = (te0-dte0)*(1 - psin**alt1)**alt2 + dte0
: zeffan