ModeCode: an inflation add-on for CAMB and CosmoMC

[Michael Mortonson]

We are pleased to announce that ModeCode, an add-on to CAMB and CosmoMC for general single-field inflationary models, is now available to download here. The paper accompanying this code release is Mortonson, Peiris, & Easther, 1007.4205.

For a given inflationary potential specified by the user, ModeCode numerically integrates the background equations and the evolution of scalar and tensor perturbations to compute the primordial power spectra from inflation. The code is interfaced with CAMB and CosmoMC to facilitate parameter estimation from CMB observations and other data sets. The parameters varied in CosmoMC include not only the parameters describing the inflationary potential but also an additional parameter to account for uncertainty in the post-inflationary reheating epoch. We have included several standard models of inflation in the code (for which chains from various combinations of current CMB data are available on the web site), and additional models can be easily added.

Please see the paper and the code web site or contact us for more information.


Michael Mortonson, Hiranya Peiris, and Richard Easther

[Christophe Ringeval]

Hi,
for those interested in the subject, I am using this post to recap this link to our patches. The onefield patch has been public for 4 years now. It is good to see that numerical methods for inflation + reheating and data analysis are becoming of some interest ;)

http://theory.physics.unige.ch/~ringeva ... /onefield/

Cheers,
chris.

PS: the original paper is there:

astro-ph/0605367

the WMAP7 analysis there:

1004.5525

[Michael Mortonson]

For any ModeCode users, there is an updated version that's compatible with the latest version of CAMB (January 2011) available to download here.

[floriberta rodriguez]

hi!!
I am a student of master. I am working modecode program. By following the steps that come in aid, generates me errors when compiling the cosmomc, errors are as follows:

gfortran -O2 -c modpk_modules.f90
gfortran -O2 -c modpk_potential.f90
gfortran -O2 -c modpk_utils.f90
gfortran -O2 -c modpk_odeint.f90
gfortran -O2 -c modpk_backgrnd.f90
gfortran -O2 -c modpk.f90
gfortran -O2 -c constants.f90
gfortran -O2 -c utils.F90
gfortran -O2 -c subroutines.f90
gfortran -O2 -c inifile.f90
gfortran -O2 -c power_tilt_inflation.f90
gfortran -O2 -c recfast.f90
gfortran -O2 -c reionization.f90
gfortran -O2 -c modules_inflation.f90
gfortran -O2 -c bessels.f90
gfortran -O2 -c equations.f90
equations.f90:1083.26:

subroutine Nu_Intvsq(EV,y, a, nu_i, G11,G30)
1
equations.f90:90.21:

use MassiveNu
2
Error: Procedure 'nu_intvsq' at (1) is already defined at (2)
equations.f90:1083.50:

subroutine Nu_Intvsq(EV,y, a, nu_i, G11,G30)
1
Error: Cannot change attributes of USE-associated symbol nu_intvsq at (1)
equations.f90:1084.30:

type(EvolutionVars) EV
1
Error: Unexpected data declaration statement in CONTAINS section at (1)
equations.f90:1085.35:

integer, intent(in) :: nu_i
1
Error: Unexpected data declaration statement in CONTAINS section at (1)
equations.f90:1086.39:

real(dl), intent(in) :: a, y(EV%nvar)
1
Error: Expected another dimension in array declaration at (1)
equations.f90:1087.41:

real(dl), intent(OUT) :: G11,G30
1
Error: Unexpected data declaration statement in CONTAINS section at (1)
equations.f90:1091.27:

real(dl) aq,q,v, am
1
Error: Unexpected data declaration statement in CONTAINS section at (1)
equations.f90:1092.23:

integer iq, ind
1
Error: Unexpected data declaration statement in CONTAINS section at (1)
equations.f90:1095.28:

am=a*nu_masses(nu_i)
1
Error: Unexpected assignment statement in CONTAINS section at (1)
equations.f90:1096.8:

ind=EV%nu_ix(nu_i)
1
Error: Unclassifiable statement at (1)
equations.f90:1097.17:

G11=0._dl
1
Error: Unexpected assignment statement in CONTAINS section at (1)
equations.f90:1098.17:

G30=0._dl
1
Error: Unexpected assignment statement in CONTAINS section at (1)
equations.f90:1099.14:

if (EV%nq(nu_i)/=nqmax) stop 'Nu_Intvsq nq/=nqmax0'
1
Error: Syntax error in IF-expression at (1)
equations.f90:1100.19:

do iq=1, EV%nq(nu_i)
1
Error: Syntax error in DO statement at (1)
equations.f90:1101.22:

q=nu_q(iq)
1
Error: Unexpected assignment statement in CONTAINS section at (1)
equations.f90:1102.19:

aq=am/q
1
Error: Unexpected assignment statement in CONTAINS section at (1)
equations.f90:1103.47:

v=1._dl/sqrt(1._dl+aq*aq)
1
Error: Unexpected assignment statement in CONTAINS section at (1)
equations.f90:1104.51:

G11=G11+nu_int_kernel(iq)*y(ind+1)*v**2
1
Error: Unexpected assignment statement in CONTAINS section at (1)
equations.f90:1105.18:

if (EV%lmaxnu_tau(nu_i)>2) then
1
Error: Syntax error in IF-expression at (1)
equations.f90:1106.51:

G30=G30+nu_int_kernel(iq)*y(ind+3)*v**2
1
Error: Unexpected assignment statement in CONTAINS section at (1)
equations.f90:1107.15:

end if
1
Error: Expecting END MODULE statement at (1)
equations.f90:1108.12:

ind = ind+EV%lmaxnu_tau(nu_i)+1
1
Error: Unclassifiable statement at (1)
equations.f90:1109.11:

end do
1
Error: Expecting END MODULE statement at (1)
equations.f90:1111.11:

end subroutine Nu_Intvsq
1
Error: Expecting END MODULE statement at (1)
equations.f90:1585.35:

if (HighAccuracyDefault) ep = ep*0.65
1
Error: Symbol 'highaccuracydefault' at (1) has no IMPLICIT type
Fatal Error: Error count reached limit of 25.
make: *** [equations.o] Error 1



someone could guide me with this error ...

[floriberta rodriguez]

hi!!
I am a student of master. I am working modecode program. By following the steps that come in aid, generates me errors when compiling the cosmomc, errors are as follows:

gfortran -O2 -c modpk_modules.f90
gfortran -O2 -c modpk_potential.f90
gfortran -O2 -c modpk_utils.f90
gfortran -O2 -c modpk_odeint.f90
gfortran -O2 -c modpk_backgrnd.f90
gfortran -O2 -c modpk.f90
gfortran -O2 -c constants.f90
gfortran -O2 -c utils.F90
gfortran -O2 -c subroutines.f90
gfortran -O2 -c inifile.f90
gfortran -O2 -c power_tilt_inflation.f90
gfortran -O2 -c recfast.f90
gfortran -O2 -c reionization.f90
gfortran -O2 -c modules_inflation.f90
gfortran -O2 -c bessels.f90
gfortran -O2 -c equations.f90
equations.f90:1083.26:

subroutine Nu_Intvsq(EV,y, a, nu_i, G11,G30)
1
equations.f90:90.21:

use MassiveNu
2
Error: Procedure 'nu_intvsq' at (1) is already defined at (2)
equations.f90:1083.50:

subroutine Nu_Intvsq(EV,y, a, nu_i, G11,G30)
1
Error: Cannot change attributes of USE-associated symbol nu_intvsq at (1)
equations.f90:1084.30:

type(EvolutionVars) EV
1
Error: Unexpected data declaration statement in CONTAINS section at (1)
equations.f90:1085.35:

integer, intent(in) :: nu_i
1
Error: Unexpected data declaration statement in CONTAINS section at (1)
equations.f90:1086.39:

real(dl), intent(in) :: a, y(EV%nvar)
1
Error: Expected another dimension in array declaration at (1)
equations.f90:1087.41:

real(dl), intent(OUT) :: G11,G30
1
Error: Unexpected data declaration statement in CONTAINS section at (1)
equations.f90:1091.27:

real(dl) aq,q,v, am
1
Error: Unexpected data declaration statement in CONTAINS section at (1)
equations.f90:1092.23:

integer iq, ind
1
Error: Unexpected data declaration statement in CONTAINS section at (1)
equations.f90:1095.28:

am=a*nu_masses(nu_i)
1
Error: Unexpected assignment statement in CONTAINS section at (1)
equations.f90:1096.8:

ind=EV%nu_ix(nu_i)
1
Error: Unclassifiable statement at (1)
equations.f90:1097.17:

G11=0._dl
1
Error: Unexpected assignment statement in CONTAINS section at (1)
equations.f90:1098.17:

G30=0._dl
1
Error: Unexpected assignment statement in CONTAINS section at (1)
equations.f90:1099.14:

if (EV%nq(nu_i)/=nqmax) stop 'Nu_Intvsq nq/=nqmax0'
1
Error: Syntax error in IF-expression at (1)
equations.f90:1100.19:

do iq=1, EV%nq(nu_i)
1
Error: Syntax error in DO statement at (1)
equations.f90:1101.22:

q=nu_q(iq)
1
Error: Unexpected assignment statement in CONTAINS section at (1)
equations.f90:1102.19:

aq=am/q
1
Error: Unexpected assignment statement in CONTAINS section at (1)
equations.f90:1103.47:

v=1._dl/sqrt(1._dl+aq*aq)
1
Error: Unexpected assignment statement in CONTAINS section at (1)
equations.f90:1104.51:

G11=G11+nu_int_kernel(iq)*y(ind+1)*v**2
1
Error: Unexpected assignment statement in CONTAINS section at (1)
equations.f90:1105.18:

if (EV%lmaxnu_tau(nu_i)>2) then
1
Error: Syntax error in IF-expression at (1)
equations.f90:1106.51:

G30=G30+nu_int_kernel(iq)*y(ind+3)*v**2
1
Error: Unexpected assignment statement in CONTAINS section at (1)
equations.f90:1107.15:

end if
1
Error: Expecting END MODULE statement at (1)
equations.f90:1108.12:

ind = ind+EV%lmaxnu_tau(nu_i)+1
1
Error: Unclassifiable statement at (1)
equations.f90:1109.11:

end do
1
Error: Expecting END MODULE statement at (1)
equations.f90:1111.11:

end subroutine Nu_Intvsq
1
Error: Expecting END MODULE statement at (1)
equations.f90:1585.35:

if (HighAccuracyDefault) ep = ep*0.65
1
Error: Symbol 'highaccuracydefault' at (1) has no IMPLICIT type
Fatal Error: Error count reached limit of 25.
make: *** [equations.o] Error 1



someone could guide me with this error ...

[Akhilesh Nautiyal(akhi)]

Dear Michael,

Are you planning to update modecode with new version (post Planck ) of Cosmomc. I tried to update it and I could run in successfully but unable to analyze chains generated with Getdist.
Can you please let me know when the updated version of Modcode will be released.

Regards,
Akhilesh