如果你是跟着“Ubuntu WSL2环境下从下载CORSIKA到跑通第一次模拟”这个需求点进来的大概率已经知道CORSIKA是什么了它全称是COsmic Ray SImulations for KAscade是大气簇射蒙特卡洛模拟里绕不开的标准工具KASCADE、KCDC以及不少宇宙线实验都在用它做参考基准。这篇文章我想用跑通一次真实模拟的完整过程把从环境搭建、源码下载、编译到写输入卡、看结果这条路捋一遍。目标读者就是两类人一是刚接触宇宙线物理、想在WSL2里跑一个最小例子验证流程的学生二是做探测器模拟、急需拿到CORSIKA输出数据继续往下走的工程向选手。下面每一步都是我在Ubuntu 22.04 WSL2环境里实际用过的方案尽量把容易卡住的地方一并讲清楚。1. 为什么要先选WSL2CORSIKA 对运行环境到底有多挑1.1 在 Windows 上直接跑 CORSIKA 的难点在哪里CORSIKA 的经典版本是 Fortran 程序官方发布的是 Linux 源码你拿到 tar 包后要自己配置、编译、链接粒子物理模型库。虽然有一些第三方打包好的 Windows 可执行版但版本老、模型不全跑起来一旦出问题很难排查因为你不知道编译时用了哪个模型、哪个Fortran编译选项。所以多数人最终都会回到 Linux 环境。那有人会问用 VMware 虚拟机或者直接装双系统不行吗行但 WSL2 在“复制粘贴命令—看到输出—改参数—再跑”这个循环里体验最好。虚拟机需要单独分配 CPU 和内存图形界面要装增强工具文件共享有时候还很别扭双系统更麻烦你只是为了跑个模拟不想每次重启切换系统。WSL2 本质上是一个轻量级虚拟机但集成度做得高Windows 和 Linux 文件可以互相访问网络共享也做了处理对跑 CORSIKA 这种命令行程序来说非常顺手。1.2 WSL2 相比虚拟机和双系统的取舍说句公道话WSL2 不是完全没有代价。它是基于 Hyper-V 虚拟化平台的第一次启用后可能需要重启另外如果你做的是需要用到 CUDA 的 GPU 加速任务WSL2 虽然官方支持 NVIDIA CUDA但配置复杂度明显高于纯 Linux 主机。CORSIKA 本身是 CPU 密集型的蒙特卡洛模拟主要吃主频和内存带宽GPU 不是必需所以 WSL2 对它来说反而是很合适的环境。如果你是 Windows 10 用户请先确认系统版本在 2004 以上并且 BIOS 里开启了虚拟化。装好之后我先跑一个简单的性能感受我用一台 i5-12400、16GB 内存的机器在 WSL2 里跑 10 个 10^4 到 10^5 GeV 的质子簇射用时在几十秒量级跟物理机上的体验没什么差别。但如果把数据放在 /mnt/c 下运行速度会明显变慢这个坑后面单独说。1.3 安装 WSL2 与 Ubuntu 22.04 的具体操作安装步骤不复杂管理员权限打开 PowerShell 或 Windows Terminal依次执行wsl --install wsl --install -d Ubuntu-22.04第一条命令会帮你启用需要的 Windows 功能并安装 WSL2 内核第二条直接装 Ubuntu 22.04。装完重启首次启动会让你设 Linux 用户名和密码。随后进入 Ubuntu 终端先把包装源和基础工具刷新一遍sudo apt update sudo apt upgrade -y我习惯顺手确认一下 WSL 版本避免后面出现奇怪问题wsl -l -v看到 Ubuntu 那行的 VERSION 是 2 就没问题。如果显示 1用wsl --set-version Ubuntu-22.04 2转一下。2. CORSIKA 下载与版本选择2.1 CORSIKA 到底是干什么的在写下载步骤之前我稍微解释一下它解决的问题。高能宇宙线进入大气层后会和原子核碰撞产生一堆次级粒子这些粒子继续碰撞或者衰变形成一簇覆盖面很大的“空气簇射”。CORSIKA 就是用一个非常详细的蒙特卡洛模型逐粒子跟踪这个过程从入射核子开始一路模拟到地面输出每个粒子的种类、能量、位置、方向、到达时间等信息。它的核心价值在于“逐粒子、分模型、可微调”你可以选择不同的高能强子模型比如 QGSJET-II-04、SIBYLL 2.3d、EPOS-LHC也可以控制能量阈值、观测高度、大气模型、切伦科夫光输出等。这正好是 WSL2 这种命令行友好环境最擅长的场景写一个输入卡跑一遍改参数再跑一遍。2.2 到哪里拿源码CORSIKA 的官方发布站点是卡尔斯鲁厄理工学院KIT维护的页面直接搜“CORSIKA KIT”就能找到。下载不是随意点击就能拿到需要填写一个申请表格写清楚单位、用途、所属研究组。审核通过后官方会邮件回复一个带密码的下载链接。整个过程通常一两个工作日。经历过的人都知道这个流程对国内用户稍微有点考验有时收件延迟有时下载链接的服务器速度一般。我的建议是尽早提交申请同时看下实验室有没有人保留过 7.74x 的包。需要注意CORSIKA 的许可协议要求仅用于科研和教育不能随便外传所以我不建议在这里贴任何第三方下载链接你自己从官方渠道申请最稳妥。拿到 tar 包之后文件名一般长这样corsika-77400.tar.gz7.74 对应 774007.75 对应 77500。解压mkdir -p ~/CORSIKA tar xvzf corsika-77400.tar.gz -C ~/CORSIKA cd ~/CORSIKA/corsika-774002.3 版本怎么选7.x 与 8.x 的区别经典路线选 7.x也就是 CORSIKA 7.74 或 7.75。这类版本使用 configure makefile 的构建方式输入卡格式非常成熟文档和网上资料最多绝大多数后处理工具也是围绕 7.x 输出格式写的。第一次跑模拟我强烈建议用 7.x。CORSIKA 8 是官方新的模块化重写版本用 C17 和 CMake 构建架构更清爽也在持续增加新物理模型。但它迭代很快接口还在变很多人光是把 CORSIKA 8 编译出来就折腾了半天更别提不同版本的输入参数经常对不上。如果你不是非要用某种新特性现阶段别拿 8.x 做第一个项目。先把 7.x 跑通理解了粒子数据流再切换到 8 会轻松很多。3. Ubuntu 侧依赖安装与目录规划3.1 需要装哪些依赖包CORSIKA 的编译器和库依赖并不多Ubuntu 22.04 默认源里就能装齐。我一般先装这一套sudo apt install -y gfortran gcc g make sudo apt install -y libpng-dev libgsl-dev libx11-dev sudo apt install -y zlib1g-dev perl逐个说下用途gfortranCORSIKA 7.x 的粒子传输核心是 Fortran 写的必须靠它编译。gcc / g部分辅助程序、接口工具和绘图相关组件需要 C/C 编译器。make驱动整个编译流程。libpng-dev生成绘图结果时需要写 PNG 文件。libgsl-devGNU 科学计算库某些版本的可选功能会用到。libx11-dev如果你的编译选项里包含了 X11 显示相关功能需要这个头文件库。perlconfigure 脚本和一些辅助工具会用到。这里有个小提醒Ubuntu 22.04 默认的 gfortran 是 11 版对 CORSIKA 7.74 完全没问题。如果以后碰到老版本源码在较新 gfortran 下编译报错优先想到-fno-range-check这类兼容性选项不要一上来就想卸载重装。3.2 目录规划与文件系统注意事项WSL2 的 Linux 文件系统和 Windows 文件系统是分开放置的。你的 Linux 家目录在\\wsl$\Ubuntu-22.04\home\用户名\这样的隐藏位置而 Windows 盘符挂载为/mnt/c、/mnt/d等。CORSIKA 最好放在 Linux 文件系统里也就是/home/用户名/xxx不要图方便放到/mnt/c/Users/xxx下面。原因有两个第一drvfs 文件系统对符号链接、文件锁和权限的支持没有本地 ext4 那么自然configure 过程会生成很多链接和脚本在 /mnt/c 下容易出问题第二跨文件系统 I/O 性能差蒙特卡洛程序要频繁写中间状态和输出文件放在 /mnt/c 上可能比本地 Linux 目录慢一个数量级。我把源码放在~/CORSIKA/下运行目录单独留一个run/这样所有产物都在同一个地方后续清理也方便。mkdir -p ~/CORSIKA/corsika-77400/run/outrun/out目录是给模拟输出数据用的后面写 OUTFILE 参数时会直接引用这个路径。4. configure 到 make真正把 CORSIKA 编译出来4.1 configure 的工作原理CORSIKA 7.x 的构建流程和很多经典 Fortran 程序类似源码根目录下有一个configure脚本它负责检测你系统里的编译器、是否安装了 X11、GSL、libpng根据检测结果生成 makefile然后把一份干净的运行目录复制到你指定的位置。进入源码目录后先看看帮助信息./configure --help不同小版本的选项有差异但核心逻辑一致。我通常直接运行./configure运行过程中脚本会问一些问题比如安装目录、要不要启用某些耦合输出、使用哪种强子模型。如果你不确定全部回车用默认值即可。也可以用非交互方式指定前缀./configure --prefix$HOME/CORSIKA/corsika-77400--prefix是最终可执行文件所在的位置。configure 结束后源码目录里会多出一个makefile同时会生成一个类似run的目录里面放着一份默认输入卡default.input和一套编译好的可执行文件模板。这里稍微注意不同版本复制出来的运行目录名字可能不同有的叫run有的直接在源码树里生成corsika目录以实际输出为准。4.2 编译几个模型库的产物接下来就是经典的 makemake -f makefile这条命令会把选定的模型库和 CORSIKA 主程序一起编译。具体编译哪些模型取决于 configure 阶段的选项。比如选了 QGSJET-II-04 和 SIBYLL 2.3dmake 之后在 run 目录下通常能看到类似这些名字的文件corsika77400Linux_QGSII_gfortran corsika77400Linux_SIBYLL23_gfortran corsika77400Linux_QGSJETII_proton_gfortran命名规律一般是corsika版本_平台_模型_编译器。不同小版本和后缀命名不完全一致但按这个规律去找不会错。第一次 make 需要几分钟因为要把强子模型的大段 Fortran 代码和主程序一起编译。看到类似Generation of CORSIKA program version 7.7400 successful的提示就说明成功了然后去 run 目录里确认可执行文件存在。ls -lh ~/CORSIKA/corsika-77400/run/4.3 编译报错怎么排查如果 make 过程报错最常见的是三类第一类找不到编译器。报错信息里直接出现gfortran: command not found这种就是前面依赖没装全回 3.1 把包装上。第二类Fortran 语法兼容性问题。Ubuntu 22.04 的 gfortran 11 对语言标准检查比较严格有些老代码会触发Error: Unclassifiable statement之类的错误。这时候可以给 configure 或 make 加上宽松选项。比如export FCFLAGS-ffixed-line-length-none -fno-range-check ./configure make -f makefile这是我对付老 Fortran 程序最常用的组合拳。-ffixed-line-length-none允许固定格式源码行长度超过 72 字符-fno-range-check允许比如sqrt(-0.0)这类不影响实际结果的数值检查通过。第三类链接阶段缺少数学库、png 库。报错末尾通常会有undefined reference to字样顺着末尾去补对应的 dev 包即可。5. 编写第一个输入卡并跑通模拟5.1 输入卡参数逐行说明CORSIKA 的输入文件是纯文本扩展名随意我习惯叫proton10.inp。每一行一个关键词加若干参数。下面这个例子是我在一个 7.74 版本上实际跑通过的最小配置模拟 10 个能量在 10^4 到 10^5 GeV 范围内的垂直入射质子RUNNR 1 EVTNR 1 NSHOW 10 PRMPAR 14 ESLOPE -2.0 ERANGE 1.E04 1.E05 THETAP 0. 0. PHIP 0. 0. ECUTS 0.3 0.3 0.3 0.3 ATMOD 1 SEED 1 0 0 SEED 2 0 0 SEED 3 0 0 SEED 4 0 0 SEED 5 0 0 SEED 6 0 0 OUTFILE /home/你的用户名/CORSIKA/corsika-77400/run/out/proton10.dat逐行解释RUNNR这次运行的编号输出文件里会用它做标识默认 1。EVTNR起始事件号跟 RUNNR 配合使用一般也是 1。NSHOW模拟多少个簇射。我建议第一次只跑 10目的就是验证流程别一上来就 5000 个。PRMPAR入射粒子类型。14 是质子如果你是氦核或铁核要换对应代码最常见的几个是 1gamma、14proton、56iron。ESLOPE能谱指数。CORSIKA 会按你指定的谱型抽取能量这里 -2.0 表示平坦一点如果你只要离散能量点可以研究ERANGE配合能量抽样模式但第一次不用管那么多。ERANGE能量范围下限和上限单位 GeV。1.E04 到 1.E05 对 WSL2 来说是安全选择既能体现簇射过程又不会慢到怀疑人生。THETAP天顶角范围和抽样单位度。0. 0.表示固定垂直入射。PHIP方位角范围和抽样。0. 0.表示固定。ECUTS能量截断阈值四个数分别是 e、e-、gamma、mu 的截断能量单位 GeV。0.3 是常用经验值低于这个能量的粒子不再往下跟踪。ATMOD大气模型1 表示美国标准大气也是最常用的默认选项。SEED随机数种子。CORSIKA 允许很多个SEED行分别对应强子相互作用、衰变、底强子截面等不同的随机流。第二三个数一般写成0 0第一个数每个 SEED 必须不同否则粒子流高度相关。OUTFILE二进制输出文件路径。记得把路径里的“你的用户名”换成实际值而且确保run/out目录已经存在。5.2 实际运行与判断成功的方法把输入文件放到 run 目录然后执行可执行文件把输入重定向进去cd ~/CORSIKA/corsika-77400/run ./corsika77400Linux_QGSII_gfortran proton10.inp程序启动后会打印一连串初始化信息包括大气模型参数、能量范围、粒子种类、随机数种子等。接着开始逐事件模拟。10 个质子在这个能量区间通常很快可能半分钟左右就结束了。如果想看耗时加timetime ./corsika77400Linux_QGSII_gfortran proton10.inp跑完怎么判断成功了我的方法是看三处一是终端有没有出现 Fortran runtime error、Segmentation fault、STOP这类中断信息。正常的结尾会输出类似“end of run”和 CPU 时间统计的段。二是检查 OUTFILE 指定的文件是否生成ls -lh ~/CORSIKA/corsika-77400/run/out/proton10.dat文件体积在几十 KB 到几百 KB 之间通常就说明有数据。三是看跑完的事件数是不是你输入的 NSHOW。如果只跑了 3 个就退出多半是随机数种子或者能量抽样设置有问题可以回头检查 SEED 行。5.3 输出文件和后续读取CORSIKA 的二进制输出格式不是简单的 CSV里面有事件头、粒子记录、运行结尾标记等结构。用记事本打开是乱码这正常。后续处理一般有以下几种方式用 CORSIKA 自带的读数据例程去解析二进制流官方源码里能找到读文件模板转换成 ROOT 文件后用 ROOT 分析先用corsika前处理工具把二进制输出转成你需要的文本格式。第一次跑通后我建议先确认“有文件生成”这件事本身暂时别急着解析全部粒子。你可以用hexdump -C proton10.dat | head看一眼文件头看到非全零的内容就行。等你需要做纵向发展曲线、地面粒子分布、Cherenkov 光子输出时再来针对具体 OUTFILE 格式写解析器。6. WSL2 图形显示给后处理程序一个 X11 环境6.1 CORSIKA 本体不需要图形但后处理可能要CORSIKA 的主模拟阶段是纯命令行不需要窗口环境。但跑完之后很多人想画个簇射纵向发展图、粒子横向分布图如果用到经典 X11 风格的可视化工具或者你想在 WSL2 里跑一些别的 Linux 图形软件就需要一个 X Server。最新版的 Windows 11 自带 WSLgLinux GUI 程序可以直接弹出窗口一般不用额外配置。但不少 Windows 10 用户或特殊构建的 WSL2 环境没有 WSLg这时推荐用 VcXsrv 来补一个 X Server。我不建议以“让 CORSIKA 自带一个图形界面”为目标因为它本来就没有这种界面。图形显示主要用于数据可视化环节比如临时用 Python 画图时需要一个 Qt 后端或者用 xool、ROOT 的 TBrowser 看数据。配置好 X11 后这些都能用。6.2 实际配置步骤Windows 侧安装 VcXsrv一路默认。启动 XLaunch 时注意几个选项Display settings 选 Multiple windowsSelect how to start clients 选 Start no client勾选 Disable access control否则连接时经常被拒绝。在 Ubuntu 侧设置显示地址。WSL2 里最简单的方式是直接用 :0export DISPLAY:0为了测试是否生效先安一个小工具sudo apt install -y x11-apps xeyes如果屏幕上出现一个跟着鼠标转的眼睛说明 X11 通道已经通了。要是:0连不上可以试试把宿主机 IP 填进去。在 Ubuntu 里查看ip route | grep default cat /etc/resolv.conf | grep nameserver有些版本需要用 nameserver 那个地址作为 DISPLAY类似export DISPLAY192.168.x.x:0.0每次新开终端都要重新 export 很烦可以在~/.bashrc末尾加一行export DISPLAY:0然后source ~/.bashrc。配置结束后再启动你的可视化后处理工具就不会报cannot connect to X server了。7. 实操中比较常见的坑编译期和运行期7.1 文件系统带来的坑我在前面反复强调源码和运行目录不要放 /mnt/c这不是玄学是我踩过的坑。有一次我把整个 CORSIKA 解压到 Windows 桌面configure 一切正常但 make 时总报permission denied或者符号链接创建失败。原因就是 drvfs 对 Linux 符号链接的支持不一致Fortran 编译过程要生成的 Mod 文件也会受文件锁机制影响时好时坏。如果你已经报错了不要犹豫直接拷贝回 Linux 目录重新 configure 一遍正常情况下能省下很多排查时间。接收文件可以用cp -r /mnt/c/Users/xxx/Desktop/corsika-77400 ~/CORSIKA/删掉源码目录里之前残留的 makefile 和编译中间产物再重新走流程。7.2 Fortran 与编译宏的坑Ubuntu 22.04 的 gfortran 11 对老代码比较“严格”编译时会报一些标准不匹配的问题。除了前面提的FCFLAGS还有一个小技巧是在 configure 之前清空环境避免你之前 export 过乱七八糟的变量unset FC F77 FFLAGS FCFLAGS如果你要换不同版本的 gfortran还可以装sudo apt install -y gfortran-10然后用./configure --with-fortran/usr/bin/gfortran-10但注意 configure 是否支持这个参数要看版本不支持的话就用FC/usr/bin/gfortran-10 ./configure的方式临时指定。7.3 模拟运行效率的坑第一次跑模拟最容易犯的错误是参数给太大。有人直接把 ERANGE 写到 10^8 GeV还跑 5000 个事件在 WSL2 里跑了一个多小时没结束还以为是程序卡死了。实际上高能量簇射粒子数会指数级膨胀计算量暴增非常明显。我的建议是从小到大测试先跑 10 个事件能量 10^4 到 10^5 GeV验证流程再跑 100 个事件看文件输出速度和体积最后再按物理需求逐步放开能量范围、天顶角范围和事件数。另外WSL2 默认会占用你 Windows 一半内存但实际可用内存不够时它会启动回收机制。你可以自己在 Windows 用户目录下建一个.wslconfig文件限制 WSL2 的资源比如[wsl2] memory8GB processors4 swap2GB改完在 PowerShell 里执行wsl --shutdown再重新进 Ubuntu 生效。这样能避免 WSL2 把宿主机内存占满也可以在内存充足时给模拟预留更多空间。7.4 随机数种子引发的结果重复CORSIKA 对随机数种子非常敏感多个 SEED 行之间第一个参数不能重复。我以前偷懒把所有 SEED 都写成1 0 0结果跑出来的簇射高度相关数据完全不能用。正确做法是让每个 SEED 互不相同比如1 0 0、2 0 0、3 0 0这样依次递增。想做可重复实验时记录下这次输入卡下次跑同样的 SEED 就能复现同一条随机流。最后再分享一个我的运行习惯每次模拟的输入卡命名不要用input.txt这种无意义名字改成p_1e4_1e5_th0_10evt.inp这种能自解释的文件名输出文件也按事件和能量归档。CORSIKA 的输入卡本身不长但等你要跑几十组参数扫谱时一个清晰的文件命名体系能帮你省掉大量找文件的时间。我现在跑 CORSIKA 基本就在这个 WSL2 环境里进行日常也不需要额外开虚拟机。先把这条“下载-编译-跑通-出数据”的最小链路打通后面无论切换到 ROOT 分析、做切伦科夫光模拟还是升级到 CORSIKA 8都是在它基础上继续加东西。如果你在跑通的路上卡住了回看一下第 4 节和第 7 节的报错场景大部分问题都能在那里找到答案。