Lamento joderte apenas llegando de victoria but i need you.
Tengo el siguiente problema:
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)