Nextflow/Groovy处理.m8文件的循环脚本修正及检查文件生成求助

原代码主要问题 #

• 脚本块混用Groovy与Bash语法：file("$params.input_dir").list()是Groovy代码，无法被Bash解析
• Bash循环、条件判断语法不符合规范
• 未闭合循环结构与脚本引号
• 未实现处理完成后生成Kraken2_iterations.check的逻辑

修正后的完整代码 #

// 启用DSL2语法（必须放在代码最开头）
nextflow.enable.dsl=2

// 定义输入、输出目录参数
params.input_dir = "path1/path2/path3/Data"
params.output_dir = "path1/path2/path3/Output"

process RUN_KRAKEN2 {
    // 自动将输出文件复制到指定目录
    publishDir params.output_dir, mode: 'copy'

    // 声明输出文件
    output:
    path("Kraken2_iterations.check")

    // Bash执行脚本
    script:
    """
    # 创建输出目录（不存在则创建）
    mkdir -p "${params.output_dir}"

    # 初始化处理文件计数器
    count=0

    # 遍历输入目录下所有.m8文件
    for m8_file in "${params.input_dir}"/*.m8; do
        # 确保当前路径是文件（避免无匹配文件时的错误）
        if [ -f "\$m8_file" ]; then
            # 替换为实际的task1命令，示例：
            # task1 --input "\$m8_file" --output "${params.output_dir}/\$(basename \$m8_file .m8).task1.result"
            echo "Running task1 on \$m8_file"
            
            # 替换为实际的task2命令
            echo "Running task2 on \$m8_file"
            
            # 计数器自增
            ((count++))
        fi
    done

    # 所有文件处理完成后生成check文件，记录处理数量
    echo "Total processed .m8 files: \$count" > Kraken2_iterations.check
    """
}

// 启动流水线
workflow {
    RUN_KRAKEN2()
}

关键修改说明 #

• DSL2启用：nextflow.enable.dsl=2必须放在代码最顶部，否则DSL2语法不生效
• 纯Bash脚本：用Bash通配符*.m8遍历文件，避免Groovy与Bash语法混用导致的错误
• 计数器实现：用Bash原生变量count统计处理的文件数，所有文件处理完成后将结果写入Kraken2_iterations.check
• 错误防护：添加[ -f "$m8_file" ]判断，防止输入目录无.m8文件时循环执行无效命令
• 变量引用：Bash变量用$var，参数传递时用引号包裹，避免文件名含空格等特殊字符引发问题

内容的提问来源于stack exchange，提问作者Grendel