Ticket #41 (new defect - FPE)

Opened 6 months ago

Last modified 6 months ago

very deep dusty model that does not conserve energy in the spectrum

Reported by: peter Assigned to: nobody
Priority: major Milestone:
Component: radiative transfer Version:
Keywords: FIR energy leak Cc:

Description (Last modified by peter)

This is a bug report from 14 Sep 2004. It describes a problem with RT where energy is not conserved. The original description is as follows (the "other problems" referred to below are probably solved by now, I think they were convergence problems in the molecular solver):


I have dug up the test case that shows the RT problem. You need to change GrnVryDepth to run it. The model runs for a long time, and it runs into lots of other problems which muddles the issue. The problem is that if you integrate over the transmitted spectrum, it has far less energy than the incident spectrum. Since this is a closed geometry, the energy in the spectrum should be conserved.


I have not modified anything from the original input file / code. I also did not test whether it still produces the same behavior with the current code.

Attachments

bug.in (1.7 kB) - added by peter on 12/01/07 21:09:14.
input script
GrnVryDpth.c (1.8 kB) - added by peter on 12/01/07 21:09:54.
modified version of GrnVryDpth?()

Change History

12/01/07 21:09:14 changed by peter

  • attachment bug.in added.

input script

12/01/07 21:09:54 changed by peter

  • attachment GrnVryDpth.c added.

modified version of GrnVryDpth?()

12/01/07 21:11:07 changed by peter

  • description changed.