.5 4 4 : dens0,aln1,aln2 ! dens[1e19 m-3] = dens0*(1 - psin**aln1)**aln2! dens[1e19 m-3] = dens0*(1 - psin**aln1)**aln2 .1 0 8 8 : te0,dte0,alt1,alt2 ! temp[keV] = (te0-dte0)*(1 - psin**alt1)**alt2 + dte0temp[keV] = (te0-dte0)*(1 - psin**alt1)**alt2 + dte0 2 : zeffan