anchor.diff:        COMMON /evolp1/ t9,hv,phie,y                   !Evolution parameters.
         COMMON /evolp2/ dt9,dhv,dphie,dydt             !Evolution parameters.
         COMMON /evolp3/ t90,hv0,phie0,y0               !Evolution parameters.
  -      COMMON /modpr/  g,tau,xnu,c(3),cosmo,xi(3)     !Model parameters.
  +      COMMON /modpr/  g,tau,xnu,c(3),cosmo,xi(3),b   !Model parameters.
         COMMON /ttime/   t,dt,dlt9dt                   !Time variables.
         COMMON /thermcb/  thm,hubcst                   !Dynamic variables.
         COMMON /endens/ rhone0,rhob0,rhob,rnb          !Energy densities.
  
positive.code: C     Contains integrands to be integrated.
  
  C----------COMMON AREAS.
        COMMON /modpr/  g,tau,xnu,c(3),cosmo,xi        !Model parameters.
        COMMON /nupar/  t9mev,tnmev,tnu,cnorm,nu,rhonu !Integration parameters.
  
  C==========================DECLARATION DIVISION=================================
  
negative.code:       COMMON /evolp2/ dt9,dhv,dphie,dydt(nnuc)       !Evolution parameters.
        COMMON /evolp3/ t90,hv0,phie0,y0               !Evolution parameters.
        COMMON /compr/  cy,ct,t9i,t9f,ytmin,inc        !Computation parameters.
        COMMON /modpr/  g,tau,xnu,c,cosmo,xi,b         !Model parameters.
        COMMON /varpr/  dt1,eta1                       !Variational parameters.
        COMMON /ttime/   t,dt,dlt9dt                    !Time variables.
        COMMON /endens/ rhone0,rhob0,rhob,rnb          !Energy densities.
  