Ticket #57 (closed physics: fixed)

Opened 3 months ago

Last modified 2 months ago

fine opacity doppler shift not propagated to rt_line_one_tauinc

Reported by: gary Assigned to: nobody
Priority: major Milestone: milestone4
Component: radiative transfer Version:
Keywords: dynamics Cc:

Description

the fine opacities resolve all lines. when dynamics is used the local rest frame line center is

long int ipLineCenter = t->Emis->ipFine + rfield.ipFineConVelShift;

where the last term is the velocity shift wrt the illuminated face

in rt_line_one_tauinc.cpp line 43 is not used:

dTau_total = rfield.fine_opac_zone[t->Emis->ipFine] * (realnum)radius.drad_x_fillfac;

note that this is the <= velocity case. the lvg case, in the following if, does not include the fine opacities at all, presumbably due to large velocity gradient

rethink all this

Change History

03/05/08 07:44:05 changed by gary

  • status changed from new to closed.
  • resolution set to fixed.