NUMÉRICO

Modificación de Códigos Cosmológicos Numéricos

Coordinador: Luis Ureña (lurena@fisica.ugto.mx)

Dentro de la paqueteria de LAMBDA
la liga directa a las programas que usa wmap
esta en:

http://lambda.gsfc.nasa.gov/product/map/current/m_sw.cfm

Hola a todos,

Aqui van algunas recomendaciones para el mini-curso que tendremos ya la
siguiente semana.

1) Recomendable que todos traigan laptop con alguna version de Linux (los
que tienen Mac pueden obviar este paso), y con compiladores de f77 y f90
instalados.

Para cmbfast: g77
Para CAMB: Intel Fortran 90. Pueden bajar una version gratis de la pagina
http://www.intel.com/cd/software/products/asmo-na/eng/282048.htm

Esta es mi recomendacion, ya que por alguna razon que no entiendo aun no he
podido compilar los codigos ni con gfortran ni con g95.

2) Sobre teoria de perturbaciones, recomiendo el paper 'Theoretical tools
for CMB physics' astro-ph/0412601, que es el paper mas completo que he visto
hasta ahora al respecto. Sera la guia que utilizare para el mini-curso.

3) Hay pequenios errores en la compilacion de los codigos. Voy a pedirle a
Miguel que nos de los permisos del repositorio que hizo, y actualizare los
archivos necesarios para que nadie tenga problemas para compilar.

El jueves enviare un nuevo mensaje esperando que todos tengan sus
compiladores instalados. Dare instrucciones para bajar los codigos y
compilarlos en su manera mas simple, con esto verificaremos que todos los
pueden correr en sus maquinas antes del minicurso.

4) Otros paquetes que utilizaremos seran el Healpix y el Skyview, para ver
mapas de temperatura. Aun los tengo bajo prueba, y espero el jueves poder
darles instrucciones para correrlos.

5) Para el codigo CAMB, tambien utilizaremos el codigo PICO (ajuste de
parametros) y las subrutinas para introducir modelos inflacionarios y de
quintessence.

Un programa tentativo del curso es como sigue:
Miercoles: Ecuaciones de la CMB (1 hr.)
Codigo cmbfast (1 hr.)
-Revision de estructura programatica
-Compilacion en diferentes sabores
-Ejemplos tipicos: LCDM, QCDM con w=cte
-Mapas de temperatura

Jueves: Ecuaciones de la CMB (1hr.)
Codigo CAMB (1hr.)
-Revision de estructura programatica
-Compilacion en diferentes sabores
-Ejemplos tipicos: LCDM, QCDM con w=cte
-Mapas de temperatura

Viernes: Ecuaciones de la CMB (1hr.)
Codigo CAMB (1hr.)
-Modelos de Quintessence
-Modelos inflacionarios
-PICO (ajuste de parametros)

Prometer no empobrece, pero hare todo lo posible para que se cumpla el
programa de arriba, :)

II

Hola a todos, este es otro mensaje en preparacion para el curso de la
proxima semana.

1) Comenzaremos con el codigo cmbfast, ya que es el mas sencillo de usar y
compilar (con g77 en linux).
a)Descarguen el archivo anexo cmbfast-mex.tgz a la carpeta de sus codigos de
CMB (e.g. cmb-codes).

b) Descomprimir con la orden
> tar xvzf cmbfast.tgz
Esto abrira una carpeta llamada cmbfast4.5.1 donde estaran los archivos
del codigo.

c) Cambiar de directorio
> cd cmbfast4.5-1

d) Compilar codigo
> sh configure
> make

Si nos les marca error alguno quiere decir que, en principio, ya pueden
empezar a jugar con el codigo cmbfast.

2) Si quieren ser un poco mas arriesgados, entonces pueden hacer su
primera figura de CMB anisotropies como sigue:
a) Cambiar a directorio de trabajo 'xxx'
> cd xxx

b) Correr script de prueba
> sh run_test.scr

Este script compilara el codigo y creara los archivos necesarios para
correrlo; estos ultimos los copiara al archivo de trabajo 'xxx'. Una vez
hecho esto, correra el codigo 'cmb' en el directorio 'xxx' utilizando los
parametros del archivo 'cmb_test.par'.

c) Crear figura de CMB con gnuplot
> gnuplot plot1.gp
> gv plot1.eps

3) Los de mayor experiencia en Linux pueden revisar el funcionamiento del
script e imaginarse lo que hace para que lo modifiquen a su gusto.

4) Tambien pueden revisar el archivo de parametros cmb_test.par y
modificarlo a su gusto. Para tener una nueva corrida, solo sera necesario
que hagan lo siguiente en el directorio 'xxx'
> ./cmb cmb_test.par
> gnuplot plot1.gp
> gv plot1.eps

y asi tendran una nueva figura de CMB anisotropies.

Por favor, avisenme si les fue posible compilar (y tal vez correr)
el programa.

III

Hola a todos, este es el otro mensaje en preparacion para el curso de la
proxima semana.

1) Ahora es el codigo CAMB. Para esto es necesario tener instalado un
compilador de Fortran 90.
a)Descarguen el archivo anexo cmbfast-mex.tgz a la carpeta de sus codigos de
CMB (e.g. cmb-codes).

b) Descomprimir con la orden
> tar xvzf CAMB-mex.tgz
Esto abrira una carpeta llamada CAMB donde estaran los archivos
del codigo.

c) Cambiar de directorio
> cd CAMB

d) Cambiar Makefile para ajustarlo a su compilador de preferencia (lineas
6-9 en el archivo Makefile)

d) Compilar codigo
> make

Si nos les marca error alguno quiere decir que, en principio, ya pueden
empezar a jugar con el codigo CAMB.

2) Si quieren ser un poco mas arriesgados, entonces pueden hacer su
primera figura de CMB anisotropies como sigue:
a) Cambiar a directorio de trabajo 'xxx'
> cd xxx

b) Correr script de prueba
> sh runCAMB_test.scr

Este script compilara el codigo y creara los archivos necesarios para
correrlo; estos ultimos los copiara al archivo de trabajo 'xxx'. Una vez
hecho esto, correra el codigo 'camb' en el directorio 'xxx' utilizando los
parametros del archivo 'params.ini'.

c) Crear figura de CMB con gnuplot
> gnuplot plot1.gp
> gv plot1.eps

3) Los de mayor experiencia en Linux pueden revisar el funcionamiento del
script e imaginarse lo que hace para que lo modifiquen a su gusto.

4) Tambien pueden revisar el archivo de parametros params.ini y
modificarlo a su gusto. Para tener una nueva corrida, solo sera necesario
que hagan lo siguiente en el directorio 'xxx'
> ./camb params.ini
> gnuplot plot1.gp
> gv plot1.eps

y asi tendran una nueva figura de CMB anisotropies.

5) Escribi a uno de los autores (Antony Lewis) del codigo para que me
proporcionara una rutina para usar modelos de quintessence, a lo cual
accedio. La rutina utiliza un potencial exponencial.

Si quieren probar, entonces cambiense al directorio 'xxx_quint' y corran
el script que se encuentra alli y hagan su figura con gnuplot
> sh runCAMB_test.scr
> gnuplot plot1.gp
> gv plot1.eps

Los conocedores avanzados de gnuplot pueden tal vez generar una figura
donde se comparen los resultados de una constante cosmologica y un modelo
con potencial exponencial.

Por favor, avisenme si les fue posible compilar (y tal vez correr)
el programa.

IV

Hola a todos,

He podido definir un poco mejor los horarios para vernos en el curso de
codigos de CMB, y esta es mi propuesta:

Miercoles: 11 a 12:30 hrs
Tomare un vuelo en la maniana y me dara tiempo para llegar al menos puntual
a la primera reunion; tambien tiempo para prepararme para el seminario del
IAC.

Jueves: 10 a 12 hrs.
Viernes: 10 a 12 hrs.

He podido definir un poco mas los temas a revisar, y el plan es revisar los
calculos de ciertos papers en el siguiente orden:
1) Ma & Bertshinger: astro-ph/9506072
2) Seljak & Zaldarriaga: astro-ph/9603033
3) Giovannini: astro-ph/0412061
4) Subramanian: astro-ph/0411049

Recuerden que seria recomendable que todos llevaran laptop, con los
programas que les envie anteriormente instalados y funcionando. De cualquier
manera, durante el curso podemos resolver los problemas tecnicos.

Envia un Mail a esta lista