gwm17/catima
1
0
Fork 0
mirror of https://github.com/gwm17/catima.git synced 2024-11-23 02:38:51 -05:00
Code Issues Projects Releases Wiki Activity

deff zcheck

Browse Source
This commit is contained in:
hrocho 2018-02-08 00:59:31 +01:00
parent a515091478
commit 1bc7fd14ef
1 changed files with 4 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
calculations.cpp
View File

@ -134,12 +134,14 @@ double bethek_barkas(double zp_eff,double eta, double zt){
double bethek_density_effect(double beta, int zt){ double bethek_density_effect(double beta, int zt){
double gamma = 1/sqrt(1-(beta*beta)); double gamma = 1/sqrt(1-(beta*beta));
double x = log(beta * gamma) / 2.3025851; double x = log(beta * gamma) / 2.3025851;
int i = zt-1; int i;
double del = 0; double del = 0;
if(zt>97){ // check if data are available, if not take highest z data if(zt>97){ // check if data are available, if not take highest z data
i=96; zt=97;
} }
i = zt-1;
if (x < density_effect::x0[i] ){ if (x < density_effect::x0[i] ){
if(density_effect::del_0[i] > 0.)del = density_effect::del_0[i] * pow(10.0,(2.*(x-density_effect::x0[i]))); if(density_effect::del_0[i] > 0.)del = density_effect::del_0[i] * pow(10.0,(2.*(x-density_effect::x0[i])));
} }