printSTDERR"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){
printSTDERR"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");
printSTDERR"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){
printSTDERR"BZ is mandatory\n";
$bz=get_answer("Vertical comp. of the Earth's mag. field (MAGNET) [downwards,muT]",0,"BZ");
printSTDERR"BZ can't be zero\n";
$bz=get_answer("Vertical comp. of the Earth's mag. field (MAGNET) [downwards,muT]",0.1,"BZ");
$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){
printSTDERR"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");
printSTDERR"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){
printSTDERR"BZ is mandatory\n";
$bz=get_answer("Vertical comp. of the Earth's mag. field (MAGNET) [downwards,muT]",0,"BZ");
printSTDERR"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