[Debian-astro-commits] [gyoto] 140/221: ThinDiskIronLine: don't assume coordinate system is spherical

Thibaut Jean-Claude Paumard thibaut at moszumanska.debian.org
Fri May 22 20:52:41 UTC 2015


This is an automated email from the git hooks/post-receive script.

thibaut pushed a commit to branch master
in repository gyoto.

commit 4900995e52f403f006a1b88fb386e233e330ca61
Author: Thibaut Paumard <paumard at users.sourceforge.net>
Date:   Wed Dec 10 16:49:42 2014 +0100

    ThinDiskIronLine: don't assume coordinate system is spherical
---
 lib/ThinDiskIronLine.C | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

diff --git a/lib/ThinDiskIronLine.C b/lib/ThinDiskIronLine.C
index fc493fe..fb70383 100644
--- a/lib/ThinDiskIronLine.C
+++ b/lib/ThinDiskIronLine.C
@@ -88,7 +88,7 @@ Gyoto::Astrobj::ThinDiskIronLine::~ThinDiskIronLine()
 double ThinDiskIronLine::emission(double nu_em, double /* dsem */,
 				  double *,
 				  double coord_obj[8]) const{
-  double rr=coord_obj[1];
+  double rr=projectedRadius(coord_obj);
   if (rr<cutradius_) return 0.;
   double dfreq=linefreq_/100.;
   /*
@@ -104,7 +104,7 @@ double ThinDiskIronLine::emission(double nu_em, double /* dsem */,
 }
 
 void ThinDiskIronLine::getVelocity(double const pos[4], double vel[4]) {
-  if (pos[1]<cutradius_){
+  if (projectedRadius(pos)<cutradius_){
     //any velocity, emission=0 anyway
     for (int ii=1;ii<4;ii++) vel[ii]=0;
     vel[0] = 1.;

-- 
Alioth's /usr/local/bin/git-commit-notice on /srv/git.debian.org/git/debian-astro/packages/gyoto.git



More information about the Debian-astro-commits mailing list