Commit 31313a7e authored by Alain O' Miniussi's avatar Alain O' Miniussi
Browse files

Merge commit 'e6d94516' into 597-compile-with-cuda

Conflicts:
	src/radiationManager.cpp
parents a3b6f126 e6d94516
......@@ -187,9 +187,9 @@ namespace fargOCA {
real specificHeat = disk.physic().adiabatic->specificHeat();
kk::parallel_for("init_stell_heating", disk.grid().fullRange(),
LBD(int i, int h, int j) {
real denom = 1+16*dt*PlankOpacity(i,h,j)*SIGMARAD*pow(temperature(i,h,j),3.0)/specificHeat;
real denom = 1+16*dt*PlankOpacity(i,h,j)*SIGMARAD*pow(temperature(i,h,j),3)/specificHeat;
eta1(i,h,j) = (temperature(i,h,j)
+ dt*(12*PlankOpacity(i,h,j)/specificHeat*SIGMARAD*pow(temperature(i,h,j),4.0)
+ dt*(12*PlankOpacity(i,h,j)/specificHeat*SIGMARAD*pow(temperature(i,h,j),4)
+ (stellarRadiation(i,h,j)+qplus(i,h,j)+planetRadiation(i,h,j))/(density(i,h,j)*specificHeat))) / denom;
eta2(i,h,j) = dt*C0*PlankOpacity(i,h,j)/(specificHeat*denom);
});
......
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