Commit 0686e162 authored by Alain O' Miniussi's avatar Alain O' Miniussi
Browse files

Merge branch 'master' of gitlab.oca.eu:DISC/fargOCA

Conflicts:
	src/diskPhysic.cpp
parents e6d94516 c70465b5
......@@ -790,11 +790,12 @@ namespace fargOCA {
real o1an = pow(o1, 2);
real o2an = pow(o2, 2);
// smoothed and continuous opacity law for regions 1, 2, and 3.
opacity=pow(pow(o1an*o2an/(o1an+o2an),2)+pow(o3/(1+1.e22/pow(temp,10)),4),0.25);
opacity = pow((pow(o1an*o2an/(o1an+o2an),2)
+pow(o3/(1+1.e22/pow(temp,10)),4)),
0.25);
} else if ( temp < t456*pow(rho,2.381e-2)){
// to avoid overflow
real ts4 = 1.e-4*temp;
// disjoint opacity laws for 3, 4, and 5.
real o4 = bk4*pow(cbrt(rho),2)/pow(ts4,9);
real o5 = bk5*pow(cbrt(rho),2)*pow(ts4,3);
......@@ -802,8 +803,8 @@ namespace fargOCA {
real o4an = pow(o4,4);
real o3an = pow(bk3*ts4,4);
// smoothed and continuous opacity law for regions 3, 4, and 5.
opacity = pow((o4an*o3an/(o4an+o3an))
+pow(o5/(1+6.561e-5/pow(ts4,8)),4),
opacity = pow(((o4an*o3an/(o4an+o3an))
+pow(o5/(1+6.561e-5/pow(ts4,8)),4)),
0.25);
} else if (temp < t678*pow(rho,2.267e-1) || rho <= 1.0e-10){
// to avoid overflow
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment