Merge branch 'master' into dev-asoreyh

README.md

@@ -66,7 +66,7 @@ Even more, by using ARTI, it is possible to calculate in a very precise way the
 
 When using ARTI, please cite the reference reported below:
 
-Christian Sarmiento-Cano, Mauricio Suárez-Durán, Rolando Calderón-Ardila, Adriana Vásquez-Ramírez, Andrei Jaimes-Motta, Luis A. Núñez, Sergio Dasso, Iván Sidelnik and Hernán Asorey (for the LAGO Collaboration), _"The ARTI Framework: Cosmic Rays Atmospheric Background Simulations"_, _European Journal of Physics_ **C**, submitted, [arXiv:2010.14591[astro-ph.IM]](https://arxiv.org/abs/2010.14591), [doi:10.48550/arXiv.2010.14591](https://doi.org/10.48550/arXiv.2010.14591) (2022).
+Christian Sarmiento-Cano, Mauricio Suárez-Durán, Rolando Calderón-Ardila, Adriana Vásquez-Ramírez, Andrei Jaimes-Motta, Luis A. Núñez, Sergio Dasso, Iván Sidelnik and Hernán Asorey (for the LAGO Collaboration), _"The ARTI Framework: Cosmic Rays Atmospheric Background Simulations"_, _European Journal of Physics_ **C**, ***accepted***, [arXiv:2010.14591[astro-ph.IM]](https://arxiv.org/abs/2010.14591), [doi:10.48550/arXiv.2010.14591](https://doi.org/10.48550/arXiv.2010.14591) (2022).
 [(bibtex)](docs/references/arti.bib)
 
 <p align="right">(<a href="#top">back to top</a>)</p>
@@ -88,6 +88,8 @@ ARTI requires the installation of the following packages:
 
 * bash
 * perl 5.0 
+* perl LWD::simple module
+* perl JSON module
 * gcc
 * make
 * screen
@@ -97,7 +99,7 @@ ARTI requires the installation of the following packages:
 As a one-liner for Ubuntu/Debian:
 
 ```bash 
-sudo apt install build-essential screen rsync git
+sudo apt install build-essential perl screen rsync git
 ```
 (git is optional).
 

docs/references/arti.bib

@@ -2,7 +2,7 @@
 	title={The ARTI Framework: Cosmic Rays Atmospheric Background Simulations},
 	author={Sarmiento-Cano, Christian and Su{\'a}rez-Dur{\'a}n, Mauricio and Calder{\'o}n-Ardila, Rolando and V{\'a}squez-Ram{\'\i}rez, Adriana and Jaimes-Motta, Andrei and N{\'u}{\~n}ez, Luis Alberto and Dasso, Sergio and Sidelnik, Iv{\'a}n and Asorey, Hern{\'a}n},
 	journal={European Physics Journal C},
-	volume={submitted},
+	volume={accepted},
 	preprint={arXiv:2010.14591},
 	doi={https://doi.org/10.48550/arXiv.2010.14591},
 	year={2022}

eas/rain.pl

@@ -397,39 +397,39 @@ for (my $i=0; $i < $nofruns; $i++) {
       }
     }
     else {
-      $altitude = get_answer("Observation level above sea level [cm]", 0, "OBSLEV");
+      $altitude = get_answer("Observation level above sea level [cm]", 1, "OBSLEV");
       while ($altitude == 0) {
-        print STDERR "Observation level is mandatory\n";
+        print STDERR "Observation level can't be zero due to CORSIKA requirements. For sea level use 1 (1 cm ~ sea level!)\n";
       }
       $modatm = get_answer("Atmospheric model selection. Start number with 'E' to use external atmospheres module, or 'G' for GDAS module", 19, "$atmcrd");
-      $bx = get_answer("Horizontal comp. of the Earth's mag. field (MAGNET) [North,muT],\nsee values at http://www.ngdc.noaa.gov/geomagmodels/struts/calcIGRFWMM", 0, "BX");
+      $bx = get_answer("Horizontal comp. of the Earth's mag. field (MAGNET) [North,muT],\nget the values from https://geomag.bgs.ac.uk/data_service/models_compass/igrf_calc.html", 0.1, "BX");
       while ($bx == 0) {
-        print STDERR "BX is mandatory\n";
-        $bx = get_answer("Horizontal comp. of the Earth's mag. field (MAGNET) [North,muT],\nsee values at http://www.ngdc.noaa.gov/geomagmodels/struts/calcIGRFWMM", 0, "BX");
+        print STDERR "BX can't be zero\n";
+        $bx = get_answer("Horizontal comp. of the Earth's mag. field (MAGNET) [North,muT],\nget the values from https://geomag.bgs.ac.uk/data_service/models_compass/igrf_calc.html", 0.1, "BX");
       }
-      $bz = get_answer("Vertical comp. of the Earth's mag. field (MAGNET) [downwards,muT]", 0, "BZ");
+      $bz = get_answer("Vertical comp. of the Earth's mag. field (MAGNET) [downwards,muT]", 0.1, "BZ");
       while ($bz == 0) {
-        print STDERR "BZ is mandatory\n";
-        $bz = get_answer("Vertical comp. of the Earth's mag. field (MAGNET) [downwards,muT]", 0, "BZ");
+        print STDERR "BZ can't be zero\n";
+        $bz = get_answer("Vertical comp. of the Earth's mag. field (MAGNET) [downwards,muT]", 0.1, "BZ");
       }
     }
   }
   else {
-    $altitude = get_answer("Observation level above sea level [cm]", 0, "OBSLEV");
+    $altitude = get_answer("Observation level above sea level [cm]", 1, "OBSLEV");
     while ($altitude == 0) {
-      print STDERR "Observation level is mandatory\n";
-      $altitude = get_answer("Observation level above sea level [cm]", 0, "OBSLEV");
+      print STDERR "Observation level can't be zero due to CORSIKA requirements. For sea level use 1 (1 cm ~ sea level!)\n";
+      $altitude = get_answer("Observation level above sea level [cm]", 1, "OBSLEV");
     }
     $modatm = get_answer("Atmospheric model selection. Start number with 'E' to use external atmospheres module, or 'G' for GDAS module", 19, "$atmcrd");
-    $bx = get_answer("Horizontal comp. of the Earth's mag. field (MAGNET) [North,muT],\nsee values at http://www.ngdc.noaa.gov/geomagmodels/struts/calcIGRFWMM", 0, "BX");
+    $bx = get_answer("Horizontal comp. of the Earth's mag. field (MAGNET) [North,muT],\nget the values from https://geomag.bgs.ac.uk/data_service/models_compass/igrf_calc.html", 0.1, "BX");
     while ($bx == 0) {
-      print STDERR "BX is mandatory\n";
-      $bx = get_answer("Horizontal comp. of the Earth's mag. field (MAGNET) [North,muT],\nsee values at http://www.ngdc.noaa.gov/geomagmodels/struts/calcIGRFWMM", 0, "BX");
+      print STDERR "BX can't be zero\n";
+      $bx = get_answer("Horizontal comp. of the Earth's mag. field (MAGNET) [North,muT],\nget the values from https://geomag.bgs.ac.uk/data_service/models_compass/igrf_calc.html", 0.1, "BX");
     }
-    $bz = get_answer("Vertical comp. of the Earth's mag. field (MAGNET) [downwards,muT]", 0, "BZ");
+    $bz = get_answer("Vertical comp. of the Earth's mag. field (MAGNET) [downwards,muT]", 0.1, "BZ");
     while ($bz == 0) {
-      print STDERR "BZ is mandatory\n";
-      $bz = get_answer("Vertical comp. of the Earth's mag. field (MAGNET) [downwards,muT]", 0, "BZ");
+      print STDERR "BZ can't be zero\n";
+      $bz = get_answer("Vertical comp. of the Earth's mag. field (MAGNET) [downwards,muT]", 0.1, "BZ");
     }
   }
   if (uc(substr($modatm,0,1)) eq "E") { # using external atmospheres bernlhor