在NGS基础:测序原始数据下载一文中提到可以使用SRA-toolkit中的命令fastq-dump从NCBI下载原始测序数据,命令如下。

nohup fastq-dump -v --split-3 --gzip SRR5908360 &
nohup fastq-dump -v --split-3 --gzip SRR5908361 &

这个代码,给我们4个提示:

  1. fastq-dump不只可以转换下载好的sra文件为fastq文件,还可以顺带下载sra文件。

    只需提供SRR号,就可以获得FASTQ序列。

    不需要先调用prefetch下载,然后再转换。

    其它参数解释见引用文章。

  2. 每一行命令后面&号表示把命令放入后台运行,当前终端可以继续输入其它命令;

    此处也相当于实现了一个手动并行下载多样本,配合for可以自动并行下载。

  3. nohup表示让程序在终端因人为原因或网络原因断开后不挂断,适用于运行时间比较长的命令,一般与&连用,形式如nohup 你的命令 & (注意空格的存在)。

    如果程序运行输出错误信息,则会写入当前目录下nohup.out文件里面,供后续查看和调试。

  4. 经常会有一些培训班“拿来主义”比较严重,以上推文和生信宝典的其它推文都被发现过直接用于某些培训班的教材,但从未申请过授权,也未引用过出处。

    更有甚者,盗版易生信早期培训教案和视频,用于自己的课程或在全网发布,希望大家多多举报。

言归正传,通常我们运行程序前,会有个预判,如前面那个例子,运行时间比较长,会使用nohup 我的命令 &的形式进行运行,从而保证程序不受网络或终端异常退出的影响。

但有时也会有误判,如没想到某个程序运行了半个小时还没结束,或数据传输时网太慢,需要传输很久,这时怎么办?中止程序,然后加上nohup再从头运行?还是有更好的办法?

下面看这个例子:马上要去吃午饭了,把文件同步到另一个服务器,饭后回来继续操作:

ysx@ehbio:~/test/Bigwig$ rsync -av * ysx@46.93.19.14:/tmp
ysx@46.93.19.14's password:
sending incremental file list
test1Y_DK10.bw

输入密码后,发现同步速度太慢了,1分钟只同步了1个文件,后面还有99个文件,待会离开后,如果网断了,终端退出,程序终止怎么办?同步不能完成,饭后怎么愉快的工作?

还好我们有下面的方案,一步步跟着操作,补救一下。

第一步,按ctrl+z把程序挂起,操作后屏幕会出现如下提示([1]中的1表示命令的作业号,后面会用到):

^Z
[1]+  已停止               rsync -av * ysx@46.93.19.14:/tmp

第二步(可选),用jobs命令查看下任务状态,跟刚才的屏幕提示一致,程序被暂时终止,作业号还是1

ysx@ehbio:~/test/Bigwig$ jobs
[1]+  已停止               rsync -av * ysx@46.93.19.14:/tmp

第三步,使用bg %1命令把作业号为1的任务放入后台,并从停止状态变为运行状态,相当于加了&后接着运行。再用jobs查看,任务状态变成了运行中,这一步很关键。如果没有运行bg %1则程序处于停止状态,一直不会运行,吃几顿饭都不会运行。

ysx@ehbio:~/test/Bigwig$ bg %1
[1]+ rsync -av * ysx@46.93.19.14:/tmp &
ysx@ehbio:~/test/Bigwig$ jobs
[1]+  运行中               rsync -av * ysx@46.93.19.14:/tmp &

第四步,运行disown -h %1,表示在终端关闭时不对作业号为1的程序发送终止信号,外部因素将不影响程序的运行。通过ps命令查看下任务进程 (可选)。

ysx@ehbio:~/test/Bigwig$ disown -h %1
ysx@ehbio:~/test/Bigwig$ ps -auwx | grep 'rsync'
ysx       18214  0.0  0.0 117844  1720 ?        S    09:43   0:01 rsync -av *.bw ysx@46.93.19.14:/tmp
ysx       18215  0.1  0.0 182376  8360 ?        S    09:43   0:04 ssh -l ysx 46.93.19.14 rsync --server -vlogDtpre.iLsfxC . /tmp
ysx       18340  0.0  0.0 112724   984 pts/1    S+   10:17   0:00 grep --color=auto rsync

通过以上4步就完成了对这次操作的事后补救。以后遇到同类问题,试一试这个新方案吧!

同时还有5点提示:

  1. 例子中使用的是rsync同步,从节省时间来看,不是一个很好的例子。

    因为把命令停掉再运行一次时,已经同步完整的数据不会再同步,时间损失不会太大。

    这也是使用同步命令rsync相比于scp的一个好处。

    更多同步方式见(Linux服务器数据定期同步和备份方式。

  2. 例子中的rsync或其它涉及两个服务器交互的命令,都需要我们人为输入登录密码,因此直接加nohup &运行是行不通的,无法接受密码的输入。

    因此通过上面这个操作先在前台启动运行、输入密码,再放入后台不挂断运行。

    从这个角度看,是一个不错的例子。

    当然解决这个问题也有其它方式,具体见ssh免密码登录远程服务器。

  3. 如果程序运行时,已加了&号,放入后台了,则只需运行jobs获得作业号,再运行disown不挂断即可。

  4. 程序作业号不一定都是1,如果之前就有程序在后台运行,作业号相应的会自加。

    后面用到作业号时也需要相应修改,不要刻板总用1

  5. nohupdisown都可以使程序不挂断,可以获得一样的效果,但原理不太一致。

    nohup可以使程序忽略挂断信号(SIGHUP)或者使程序脱离终端的控制,从而终端不能再对其发送挂断信号(SIGHUP);

    disown则是内生于shell,告诉shell在终止时不对对应程序发送挂断信号(SIGHUP)。

猜你喜欢

10000+:菌群分析 宝宝与猫狗 梅毒狂想曲 提DNA发Nature Cell专刊 肠道指挥大脑

系列教程:微生物组入门 Biostar 微生物组  宏基因组

专业技能:学术图表 高分文章 生信宝典 不可或缺的人

一文读懂:宏基因组 寄生虫益处 进化树

必备技能:提问 搜索  Endnote

文献阅读 热心肠 SemanticScholar Geenmedical

扩增子分析:图表解读 分析流程 统计绘图

16S功能预测   PICRUSt  FAPROTAX  Bugbase Tax4Fun

在线工具:16S预测培养基 生信绘图

科研经验:云笔记  云协作 公众号

编程模板: Shell  R Perl

生物科普:  肠道细菌 人体上的生命 生命大跃进  细胞暗战 人体奥秘

写在后面

为鼓励读者交流、快速解决科研困难,我们建立了“宏基因组”专业讨论群,目前己有国内外5000+ 一线科研人员加入。参与讨论,获得专业解答,欢迎分享此文至朋友圈,并扫码加主编好友带你入群,务必备注“姓名-单位-研究方向-职称/年级”。PI请明示身份,另有海内外微生物相关PI群供大佬合作交流。技术问题寻求帮助,首先阅读《如何优雅的提问》学习解决问题思路,仍未解决群内讨论,问题不私聊,帮助同行。

学习16S扩增子、宏基因组科研思路和分析实战,关注“宏基因组”

点击阅读原文,跳转最新文章目录阅读

耗时很长的程序忘加nohup就运行了怎么办?相关推荐

  1. 送书 | 耗时很长的程序忘加nohup就运行了怎么办?

    在NGS基础:测序原始数据下载一文中提到可以使用SRA-toolkit中的命令fastq-dump从NCBI下载原始测序数据,命令如下. nohup fastq-dump -v --split-3 - ...

  2. 开发那些事儿:如何解决RK芯片视频处理编解码耗时很长的问题?

    流媒体视频直播包括以下几个步骤:采集->处理->编码和封装->推流到服务器->服务器流分发->播放器流播放. 在流媒体处理编码的过程中,会有硬解码和软解码两种播放方式.两 ...

  3. JAVA入门级教学之(JAVA程序的加载和运行)

    JAVA程序的加载和运行 多思考多动脑(边参考文章最后的示意图,边按步骤理解) 1.JAVA程序的加载和运行包括两个非常重要的阶段: 编译阶段 运行阶段 2.我们先来了解一下什么是编译阶段: 首先,我 ...

  4. c语言中组件出现错误,错误:R6034 运行时错误! 应用程序尝试加载 C 运行库时出现错误。 有关详细信息,请联系应用程序的支持团队...

    错误:R6034 运行时错误! 应用程序尝试加载 C 运行库时出现错误. 有关详细信息,请联系应用程序的支持团队 错误消息 ArcGIS Engine 应用程序 10.4 和 10.4.1 版本会返回 ...

  5. 耗时很长的服务器端事件中让客户端得到中间过程信息的合理解决方案

    需求:B/S结构的系统里,用户点一个按钮系统开始发送上千封邮件,要求把发送信息(发送成功数,失败数,剩余数量...)动态实时的反馈给客户. 分析和实施过程当中遇到的问题: 一:最低级的问题 由于客户催 ...

  6. 解决magento保存产品时耗时很长的问题

    以前我在更新产品属性值(拿price为例)的时候,通常会这样做: foreach($product_ids as $id){$product = Mage::getModel('catalog/pro ...

  7. VC6.0中ObjectARX应用程序的加载与运行

    使用  ObjectARX ,首先要确定你的目标平台,获得适当的开发环境.如果你是在 AutoCAD2002  平台上开发,你就需要具备下面的工具和软件: � AutoCAD 2002 中文版或英文版 ...

  8. java运行原理_Java程序的加载与运行原理详解

    Java程序的运行包括两个非常重要的阶段: 一.编译阶段 第一步:程序员需要在计算机硬盘中任意位置创建一个.java扩展名的文件,该文件被称为 java源文件,源文件当中编写的是java源代码/源程序 ...

  9. linux udt 源码,UDT linux下关闭链接耗时过长

    不知道使用UDT的人多不多,很大可能会石沉大海啊 ========================================================= 我在使用UDT库的时候,linux ...

最新文章

  1. Xcode9新特性介绍-中文篇
  2. 软件工程个人项目——买书的最低价格
  3. autojsui界面关闭_autojs 第九次 ui界面交互获取
  4. autoconf.mk详析
  5. HEC-RAS二维溃坝洪水模拟(尾矿库)
  6. 有效的单元测试--总结
  7. java loadjs_Javarscript中模块(module)、加载(load)与捆绑(bundle)详解
  8. 细数音频放大器的分类
  9. LINUX基本指令、NGINX负载均衡基本操作与解释
  10. 计算机二级用的ms什么版本,计算机二级ms office用的哪个版本
  11. 服务器排队系统怎么做,多服务器排队系统的平均排队时间的近似分析
  12. OpenStack计算节点nove启动报错
  13. Python条件语句和循环
  14. Ubuntu16.04在4K显示器中,字体太小,进行调整
  15. 【转】ATF中SMC深入理解
  16. python pexpect使用介绍
  17. PATA 题解 1002
  18. 7-65 字符串替换 (15 分) 本题要求编写程序,将给定字符串中的大写英文字母按以下对应规则替换: 原字母 对应字母 A Z B Y C X D W … … X C Y B Z A
  19. Open3D键盘切换上下帧显示点云
  20. Excel解密——okfone解密大师

热门文章

  1. 互联网寒冬裁员潮,这些人可能有点慌!
  2. 100%的BAT招聘岗位都考的知识,你精通了吗?
  3. springboot 加载配置信息(静态)
  4. 多级反馈队列调度算法具体原理
  5. 卫星覆盖分析品质 FOM
  6. ICON: 从单张图片重建穿衣服人体模型
  7. ICRA 2021 | UPSLAM:联合全景SLAM
  8. 大场景三维点云的语义分割综述
  9. 谷歌最新论文:手机端的实时3-D目标检测
  10. 怎么两边同时取ln_脏辫发型怎么编编发教程图解简单易学!