Problem with Adding a new model parameter in CosmoMC
Posted: January 12 2020
Hi,
As mentioned in my post from last month I'm interested in adding a new model parameter in the recent 2019 version.
Following the Readme, I added a 17'th parameter in the SetForH subroutine of CosmologyParameterizations.f90
right after CMB%fdm = Params(16) , and changed
call this%SetTheoryParameterNumbers(16,last_power_index)
to
call this%SetTheoryParameterNumbers(17,last_power_index).
I also defined the new parameter in Calculator_CAMB.f90, CosmologyTypes.f90 , camb.f90, and model f.90,
and of course modified results.f90 and recfast.f90 to reflect the new physics accounted for by the new model parameter.
In addition, I added the new parameter in the paramnames-file (in the 17'th line between fdm and logA)
and defined its range in the ini-file.
The code compiles OK. However, it stops shortly after it starts without calculating anything.
The error message I get is: "SetTheoryParameterNumbers: parameter numbers do not match".
Obviously, I'm doing something wrong but I can't tell what exactly. Does anyone know how to fix this error ?
Thanks,
Meir
As mentioned in my post from last month I'm interested in adding a new model parameter in the recent 2019 version.
Following the Readme, I added a 17'th parameter in the SetForH subroutine of CosmologyParameterizations.f90
right after CMB%fdm = Params(16) , and changed
call this%SetTheoryParameterNumbers(16,last_power_index)
to
call this%SetTheoryParameterNumbers(17,last_power_index).
I also defined the new parameter in Calculator_CAMB.f90, CosmologyTypes.f90 , camb.f90, and model f.90,
and of course modified results.f90 and recfast.f90 to reflect the new physics accounted for by the new model parameter.
In addition, I added the new parameter in the paramnames-file (in the 17'th line between fdm and logA)
and defined its range in the ini-file.
The code compiles OK. However, it stops shortly after it starts without calculating anything.
The error message I get is: "SetTheoryParameterNumbers: parameter numbers do not match".
Obviously, I'm doing something wrong but I can't tell what exactly. Does anyone know how to fix this error ?
Thanks,
Meir