gkowal/amun-code
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

Make the satellite orbit really eliptical.

Browse Source
This commit is contained in:
Grzegorz Kowal 2011-05-09 18:49:28 -03:00
parent 09b6955cd4
commit 22a0e1b2ed
1 changed files with 3 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

9
src/problem.F90
View File

@ -1409,10 +1409,10 @@ module problem
! !
integer :: i, j, k integer :: i, j, k
real :: dx, dy, dz, dr, ang = 0.0d0, sn = 0.0d0, cs = 1.0d0, xsh, ysh real :: dx, dy, dz, dr, ang = 0.0d0, sn = 0.0d0, cs = 1.0d0, xsh, ysh
real :: dist, rad, amp, asat, bsat, xsat
real :: dnamb, pramb real :: dnamb, pramb
real :: dnstar, prstar, dnvstar, rc real :: dnstar, prstar, dnvstar, rc
real :: dnsat , prsat , dnvsat , rs, xs, ys, zs real :: dnsat , prsat , dnvsat , rs, xs, ys, zs
real :: asat, bsat, xsat
#ifdef ADI #ifdef ADI
real :: ekin, ekstar, eksat real :: ekin, ekstar, eksat
#ifdef MHD #ifdef MHD
@ -1480,11 +1480,8 @@ module problem
asat = dsat / (1.0d0 - esat) asat = dsat / (1.0d0 - esat)
xsat = asat * esat xsat = asat * esat
bsat = sqrt(asat * asat - xsat * xsat) bsat = sqrt(asat * asat - xsat * xsat)
rad = 0.5d0 * (asat + bsat) xsh = asat * cs - xsat
amp = 0.5d0 * (asat - bsat) ysh = bsat * sn
dist = asat - amp * sn * sn
xsh = dist * cs - xsat
ysh = dist * sn
! reset update ! reset update
! !