Te acordas del programa que calcula la amplitud de las fluctuaciones para poner en el pstart.dat, bueno, algo estoy haciendo mal porque no anda. Te recuerdo que lo que haciamos era calcular el valor de sigma8 en el Z donde queriamos empesar la simulacion (d1.f) luego con ese valor corriamos el sg.f, pero resulata que desde que te fuiste el pow.f a sufrido muchisimas modificaciones de las cuales soy absolutamente responsable. La pregunta creo que es la siguiente: El factor milenario que alli aperece de donde lo sacaste/ron? y en que unidades va el espectro? que es el delta que alli aparece?.

El diego me esta matando con esta historia para colmo anda histerica con el viaje y que nadie le da bola con el mundial.

        read(32,*)aaa,ppp,cte
        pi=4.*atan(1.)

        L=128
        rLbox=600.
        h100=.65
        rL=float(L)
        delta=1.                                                                       QUE ES ESTO?
        nstep=2000
        rkmin=0.
        rkmax=pi*rL/rLbox

        drk=(rkmax-rkmin)/nstep

        ru=8./h100

        sum=0.

        do rk=drk,rkmax,drk

          p=pow(rk)

          p=p*(rlbox/2./pi)**3*delta**2                                                         (FACTOR MILENARIO)
          w=3.*(sin(rk*ru)-rk*ru*cos(rk*ru))/(rk*ru)**3.
          f=rk**2*p*w**2.
          sum=sum+f
        enddo

        write(*,*)'pi,drk,sum',pi,drk,sum
        sig=sqrt(4.*pi*sum*drk)

        write(*,*)cte/sig                                               ESTO ES LO QUE ESTOY PONIENDO EN EL PSTART.DAT

        end
 

El numero de onda lo estoy poniendo en unidades de 1/Mpc (sin h) la altura en lo que sale de las formulas
Espero que tengas un feliz reencuentro con la flia, o se quedan un tiempo mas en judeolandia?

un abraso
 
 

-- 
Manuel E. Merchan
Grupo de Investigacion en Astronomia Teorica y Experimental (IATE)
Observatorio Astronomico de la Universidad Nacional de Cordoba
Laprida 922, Cordoba, Argentina
Tel-Fax: +54-51-331066
Tel: +54-51-882402 (home)