Files to duplicate mass/age underflow with DSED models

indivWD.csh

#!/bin/csh
#
# Paper has been submitted, but found a bug in wdEvol.c (had to change "log" to "log10"), 
# so rerunning everything
#
# rerunning Renedo models to see if the slight mass adjustments helped

set nProc    = 2		# number of processors to run on (must be >= 2)
set thin     = 050 		# starting with this 
set seed     = 001
set IFMR     = 1	# all IFMRs
set WDModel  = (2 3) 		# Renedo WD model

set dir  = "./photfiles"
set dir2 = "./bin"
set code = "mpiMcmc"	# BigStep and log->log10 fix, then improved Renedo mass values

cd $dir
set photfiles = J*new
cd ../

foreach model ($WDModel)
   foreach photfile ($photfiles)

      set outfile = "$dir/mcmc.$photfile:r.m3.w$model.i$IFMR.s$seed.t$thin.zeus"

      echo "$dir2/$code --photFile $dir/$photfile --outputFileBase $outfile --wdModel $model --ifmr $IFMR --seed $seed --thin $thin"
      $dir2/$code --photFile $dir/$photfile --outputFileBase $outfile --wdModel $model --ifmr $IFMR --seed $seed --thin $thin

   end
end

indivWD.yaml

general:
    seed:               73    # initialize the random number generator
    verbose:            0     # 0 is quiet, 1 is verbose, 2 is very verbose

    files:
        photFile:       "photfiles/test"
        outputFileBase: "photfiles/mcmc.test"
        scatterFile:    ""    # Empty string.  This is primarily a CLI option
        modelDirectory: "models/"

    main_sequence:
        filterSet:      2     # 0 = UBVRIJHK, 1 = ACS, 2 = SDSS + JHK
        msRgbModel:     3     # 0 = Girardi, 1 = Chaboyer-Dotter w/He sampling, 2 = Yale-Yonsei, 3 = DSED

    white_dwarfs:
        wdModel:        1     # 0 = Wood, 1 = Montgomery, 2 = Althaus, 3 = Renedo
        M_wd_up:        8.0   # The maximum mass for a WD-producing star
        ifmr:           1     # 0 = Weidemann, 1 = Williams, 2 = Salaris lin, 3 = Salaris pw lin, 4+ = tunable
        carbonicity:    0.38

    brown_dwarfs:
        bdModel:        1     # 0 = None, 1 = Baraffe

    cluster:
        priors:
            Fe_H:       0.0
            distMod:    4.0
            Av:         0.0
            Y:          0.29

        sigmas:
            Fe_H:       0.3
            distMod:    2.5
            Av:         0.0
            Y:          0.0

        logClusAge:     9.62 # not a prior, just a starting value

        minMag:        -10.0  # these variables can be used to read a subset of the photometry
        maxMag:         30.0
        index:          1     # with 0 being the first filter in the dataset

makeIsochrone:
    # No makeIsochrone-specific settings yet

scatterCluster:
     brightLimit:       0     # apparant mags, can remove bright stars, e.g. RGB
     faintLimit:        30    # apparant mags, can remove faint stars, e.g. faint MS and WDs
     relevantFilt:      1     # 0=bluest band available
     limitS2N:          1     # use to remove objects with overly low signal-to-noise
     exposures:         [0, 1, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0]   # exposure times in hours, e.g., (1 1 1 0 0 ...) for UBV 

makeCMD:
    # No makeCMD-specific settings yet

mpiMcmc:
    burnIter:           10000
    runIter:            10000
    thin:               1

sampleWDMass:
    # No sampleWDMass-specific settings yet

simCluster:
    nStars:             100
    percentBinary:      0     # percent binaries (drawn randomly)
    percentDB:          0     # percent of WDs that have He atmospheres (drawn randomly)
    nFieldStars:        0
    nBrownDwarfs:       0

J0003-0111.DA.grizJHK.new

id g r i z J H K sigg sigr sigi sigz sigJ sigH sigK mass1 massRatio stage1 CMprior useDuringBurnIn
1 19.32 18.79 18.56 18.52 17.63 17.31 17.35 0.02 0.02 0.02 0.04 0.04 0.04 0.05 4.0 0.0 3 1 1