weiss_hot_component.param

This file contains the parameters defining a radial grid
for simulating the hot gas component of APM08279 according
to the results of Weiss et al 2007
This file is used as input to the program logrid3.f, which
generates the output file:

weiss_hot_component.logrid

which is how you define the radial shells for the radiative transfer
program molnd_7.f

Here's an example run:

orthanc% ../bin/molnd_7
  Start new (N) or continue old (O) calculation ?
N
  Starting NEW calculation...
  Input file name for radial grid ?
weiss_hot_component.logrid
  Molecular transition data file ?
p-h2o.spec
  Collisional excitation rates file ?
p-h2o.collrates
  Dust parameters file ?
dust.param
  Molecular abundance factor (normally 1.0) ?
1.0
  Dust absorption factor (normally 1.0) ?
1.0
  Output file ?
weiss_hot_component.p-h2o.out
  Use accelerated (1) or normal (0) lambda-iteration ?
1
  Relative fractional population change
  threshold for convergence (i.e. 1.E-3) ?
1.e-3
  Maximum number of iterations (i.e. 30) ?
10
  Solving rates for iteration # 0
  Solving transfer for iteration # 1
  Solving rates for iteration # 1
  Maximum fractional population change is  8.53544669
  This occurs for level 23( 7( 2, 6) ) in radial shell 50( / 50)
  Old population:   0.00529396742
  New population:   0.000554293065
  Population-weighted r.m.s pop. change is:  0.00939823663
  Solving transfer for iteration # 2
  Solving rates for iteration # 2
  Maximum fractional population change is  1.17859415
  This occurs for level 19( 7( 1, 7) ) in radial shell 50( / 50)
  Old population:   0.00187483721
  New population:   0.000860031051
  Population-weighted r.m.s pop. change is:  0.00477847345
....
keep going, this converged after around 50 iterations.

Now you are ready to run the program that generates the
predicted spectra. You first need to create a little file
weiss_hot_component.observ, that specifies the FWHM beam size
for the observation. 


orthanc% ../bin/observ_spec
  Name of file created by molnd ?
weiss_hot_component.p-h2o.out
  File name giving observation info ?
weiss_hot_component.observ
  Abundance factor (normally 1.0) ?
1.
  Dust absorption factor (normally 1.0) ?
1.
  Output file ?
weiss_hot_component.p-h2o.spec.txt
orthanc% 
