argiopetech/base

Missing CHAB variables

Closed this issue · 3 comments

argiopetech commented

If loading CHAB models, mpiMcmc still requests two variables from the command line.

    /* open model file, choose model set, and load models */
    if (settings.mainSequence.msRgbModel == CHABHELIUM)
    {
        scanf ("%lf %lf", &ctrl.priorMean[YYY], &ctrl.priorVar[YYY]);

        if (ctrl.priorVar[YYY] < 0.0)
        {
            ctrl.priorVar[YYY] = 0.0;
        }
    }
tedvh commented

Right. These are the subset of models that have a range of helium (Y),
so they need a helium prior value and prior sigma.

On 7/4/13 12:57 PM, Elliot Robinson wrote:

If loading CHAB models, mpiMcmc still requests two variables from the command line.

 /* open model file, choose model set, and load models */
 if  (settings.mainSequence.msRgbModel  ==  MsModel::CHABHELIUM)
 {
     scanf  ("%lf %lf",  &ctrl.priorMean[YYY],  &ctrl.priorVar[YYY]);

     if  (ctrl.priorVar[YYY]  <  0.0)
     {
         ctrl.priorVar[YYY]  =  0.0;
     }
 }

—
Reply to this email directly or view it on GitHub
#27.

Ted von Hippel

Department of Physical Sciences
Embry-Riddle Aeronautical University
600 S. Clyde Morris Boulevard
Daytona Beach, FL 32114-3900
386-226-7751

argiopetech commented

This block of code will fix the issue if anyone needs to use the Chab models between now and the 9.3.0 release. A modified version of this code is already in the cpp_conversion branch and will be in 9.3.0.

    if (settings.mainSequence.msRgbModel == CHABHELIUM)
    {
        ctrl.priorMean[YYY] = settings.cluster.Y;
        ctrl.priorVar[YYY] = settings.cluster.sigma.Y;
    }
    else
    {
        ctrl.priorMean[YYY] = 0.0;
        ctrl.priorVar[YYY] = 0.0;
    }
tedvh commented

Thanks Elliot.

This might be relevant to your work David, so I'm cc'ing you on this.

-Ted

On 7/4/13 4:41 PM, Elliot Robinson wrote:

This block of code will fix the issue if anyone needs to use the Chab models
between now and the 9.3.0 release. A modified version of this code is already in
the cpp_conversion branch and will be in 9.3.0.

 if  (settings.mainSequence.msRgbModel  ==  CHABHELIUM)
 {
     ctrl.priorMean[YYY]  =  settings.cluster.Y;
     ctrl.priorVar[YYY]  =  settings.cluster.sigma.Y;
 }
 else
 {
     ctrl.priorMean[YYY]  =  0.0;
     ctrl.priorVar[YYY]  =  0.0;
 }

—
Reply to this email directly or view it on GitHub
#27 (comment).

Ted von Hippel

Department of Physical Sciences
Embry-Riddle Aeronautical University
600 S. Clyde Morris Boulevard
Daytona Beach, FL 32114-3900
386-226-7751