Corsika install - ciemat-tic/codec GitHub Wiki

Corsika install and configuration for ACME

Required libs:

yum install perl-Switch

Corsika install

tar -xzvf corsika-73500-lago-single.tar.gz
cd corsika-73500-lago-single
./coconut -b

Add corsika to user path

#simplify app name
cd /home/localsoft/corsika/corsika-73500-lago-single/run
ln -s corsika73500Linux_QGSII_gheisha corsika

#add it to profile
#more /etc/profile
#CORSIKA
export PATH=$PATH:/home/localsoft/corsika/corsika-73500-lago-single/run

Install aux libs

git clone https://github.com/supermanue/arti.git
cd arti
make
#create simpler name && add it to profile
ln -s do_grid.sh corsika_create_input.sh
#more /etc/profile
export PATH=$PATH:/home/localsoft/corsika/arti/sims

Ejecución

./do_sims.sh -w ../../corsika-73500-lago-single/run/ -p prueba -t 3600 -v 75600 -u manuel -y -k 55000

(te pregunta cosas. Darle a OK a todo excepto...

• Horizontal comp. of the Earth's mag. field (MAGNET) [North,muT], see values at http://www.ngdc.noaa.gov/geomagmodels/struts/calcIGRFWMM?
  • <BX 0>: 21.39
• Vertical comp. of the Earth's mag. field (MAGNET) [downwards,muT]?
  • <BZ 0>: -11.91

La salida en la terminal nos muestra las opciones que hemos elegido y a partir de ahí tenemos que darle algunos otros parámetros. Les recomiendo que usen un rango de energía baja, recuerden que van a correr sobre sus maquinas.

Al final esto va a generar varios archivos llamados por ejemplo:

go-prueba-all_01.sh alojados en la corsikaxxx/run/

Ejecutarlos todos, por ejemplo con el siguiente script.

fo

uego de esto en la carpeta corsikaxxx/run/prueba/ se encontraran todos los archivos binarios de las simulaciones corridas.

Ahora para analizar los archivos se usará el código de corsikareader para extraer los datos importantes y el analysis para hacer esto usamos la siguiente linea de bash que ejecutamos directamente sobre la terminal:

for i in DAT??????.bz2; do j=$(echo $i | sed -e 's/.bz2//'); u=$(echo $j | sed -e 's/DAT//'); bzip2 -d -k $i; echo $j | ../../arti/analysis/lagocrkread | ../../arti/analysis/analysis -p -v $u; rm $j; done

y luego,

bzcat *sec.bz2 | ../../arti/analysis/showers -a 20 -d 20 -c 5200. -v salida

Al final tendremos tres archivos llamados:

salida.dst salida.hst salida.shw.bz2