Welcome to RQL's Blog!
Do Better Every Day!-
putty快捷方式登录
putty利用快捷方式自动登录的两种方法。
下述操作的前提是已经在 Putty configuration 中保存了session 的各种配置（包括session名字，账户名，密码等等）。密码登录
- 桌面创建PuTTY.exe的快捷方式。
- 右键快捷方式，属性>快捷方式>目标，加入如下三个参数：
X:\yourpath\putty.exe -load "session_name" -l "username" -pw "password"session_name, username, password 替换为自己的账号信息。
- 点击应用保存。这样下次打开该快捷方式时就会自动登录属性中写的那个账户。为方便管理，可以将该快捷方式重命名为session_name
密钥登录
步骤同上，只是在修改快捷方式属性时，将参数改为
X:\yourpath\putty.exe -i "X:\mykey.ppk" username@123.45.67.89-i后的引号内为私匙路径，username 为用户名，123.45.67.89 为远程节点的IP。
如果想要调整终端的色彩和字体，更改default配置即可。
参考：https://segmentfault.com/a/1190000000639516
-
SCI论文写作与发表注意点
-
很多年轻的学者总是会质疑自己的工作是否值得发表。工作的可信度、创新性、数据、模式？此时，please appreciate your tiny achievement
-
写论文的时候，记得要站在审稿人、读者的角度，想想他们在看到这篇文章时会提出什么问题。
-
审稿人总是会比较焦躁、不耐烦的，偏见也总是存在的
-
为了发表几张图片，你可能需要画一大堆图，并可能需要花很长的时间去分析一张图
-
你的文章最好能让别人感兴趣并留下深刻印象，不要让别人读着读着睡着了
-
充分理解要呈现的结果。Always create an outline。写文章时最先准备的是图表和图注
-
文章被拒的多数是因为重复度太高，所以在写文章时避免重复性或与其他文章类似的表述
-
注意投稿期刊的格式，不同期刊的格式要求不同，例如副标题的点、作者名字、引言格式
-
写文章用新罗马字体（字号12），双倍行距，表上页码和行数；但在做PPT时用Arial字体，此外，在做PPT时注意变换图片形式，以吸引观众
- 写文章避免长度超过四行的长句以及不必要的词句和重复。
- 表述不要过于绝对或强烈，但也不要过于平淡。
- 避免一张图一句话，避免一段话只有一句话。
- 避免在没有陈述他人研究成果的前提下引用过多文献。
-
不要使用太多缩略语，不要过分夸大特征，不要随意评判别人的研究成果（他可能是你的审稿人），不要在一个句子中重复使用同一个单词
-
强调研究结果的重要影响，但也要对研究结果的caveats诚实。
- 你的图是否清晰，特征是否明显标注，你的图对读者是否友好
-
-
瞬变涡旋诊断量
折射指数
在折射指数为正的区域，波动传播。
在折射指数为负的地方，波动容易消散。
且波动倾向于往正反射指数大的区域传播。上面图a与图b分别显示的是不同情况下的纬向平均U风和折射指数的合成图。
在high index时，急流强度减弱且变得更宽，此时波动向赤道移动，离开了波动源（即eddy-driven jet区域），有利于产生eddy momentum fluxes，进而通过波流相互作用维持西风急流异常。
在low index时，急流变窄且增强，此时折射指数在波动源区达到最大，意味着波动难以离开波源，不利于eddy momentum fluxes的产生，使得波流相互作用减弱，难以维持西风急流异常。 此时若计算两种情况下的eddy momentum fluxes，会发现high index的eddy momentum fluxes高于low index。该图也表明，当急流过强时，它就像一个波导，减弱了波动的径向传播，波流相互作用减弱。
-
Linux进程管理与后台运行
常用任务管理命令
jobs //查看任务，返回任务编号n和进程号 bg %n //将编号为n的任务转后台运行 fg %n //将编号为n的任务转前台运行 ctrl+z //挂起当前任务 ctrl+c //结束当前任务 kill -9 进程号 #立即强制删除一个工作 kill -15 %n #以正常的程序方式终止一项工作不受HUP影响的后台运行
在Linux中，如果要让进程在后台运行，一般情况下，我们在命令后面加上&即可。 但这样是将命令放入到一个作业队列中了，但当网络断开或本地终端断开连接后，终端会发送hangup（HUP）信号来通知其关闭所有子进程，因此该任务也会终止。 为了使任务不受HUP（hangup）信号影响，该博文 记录了好几种方法。 本人觉得最好用的是以下两种，记录如下。
nohup
在正常命令前加入nohup，就可以让提交的命令忽略 hangup 信号。
nohup mpirun –np 8 ./wrf.exe nohup command > myout.file 2>&1 & # 指定了输出信息的文件，缺省情况下该作业的所有输出都被重定向到一个名为nohup.out的文件中用
nohup后，可以用 ` jobs ` 来查看任务， 再次登录查看任务可以用 ` top ` 或者 ` ps –a `
` ps -ef ` 可以查看主机所有运行的进程 ` ps -ef | grep ` 过滤条件
但是，请注意，以上两个命令都只能查看单个节点的运行情况。对于并行运行、集群系统、天河等电脑，则无法使用上述命令来查看任务提交情况，而应该使用qstat（四期）,yhq(tianhe)screen
通过screen来创建新的虚拟终端，在该终端中运行作业。 关于screen的进一步了解可以参考这篇博文 。https://www.ibm.com/developerworks/cn/linux/l-cn-screen/
# 最好给每个会话命名，便于区分 screen -dmS session_name # 建立一个一开始便处于断开模式下的会话（-S 表示创建screen会话时为会话指定一个名字） screen -r session_name # 重新连接指定会话 screen -d session_name # 断开其他正在运行的screen会话 screen -S session_name -X quit # 删除一个会话 screen -list # 列出所有会话名、状态 screen -wipe # 会清除会话状态为 dead 的会话 # 进入会话后，可用快捷键CTRL-a d 来暂时断开当前会话，此时对话里的任务不会终端 # 进入会话后，若选择输入 exit 来退出当前会话，则该会话就会被删除 screen -h num # 可以指定历史回滚缓冲区大小为num行
-
Linux批量修改文件名
最近发现linux有一个批量修改文件名的命令:
rename 原字符串 目标字符串 需要修改的文件根据网上资料，该命令一般有两种用法，上面的是c语言类型的用法，还有一种正则表达式用法。但个人比较喜欢前一种用法。
发现新服务器中如果用上述命令会报如下错误“”，经查找后，此时需要把语法改为
rename 's/a/b/' ./*.txt # 将a替换为b
-
关于缺测
发现再分析资料中一般是没有缺测的，即便在青藏高原，1000hPa也依然有值。
但模式输出的资料中，在有地形或地表气压低于1000hPa的地方，1000hPa等压面上就会有缺测值，因此模式陆地底层经常有大片区域的缺测。
当利用模式资料算一些比较复杂的诊断量时要注意缺测值。ncl中关于缺测的处理
ncl中的缺测值用 var@_FillValue 表示。
;缺测值赋值，气象中常用的缺测值还有-999，9.969209968386869e+36，9.96921e+36 if (any(brunt.eq.0)) then if (.not.isatt(brunt,"_FillValue")) then if (typeof(brunt).eq."double") then brunt@_FillValue = 1d20 else brunt@_FillValue = 1e20 end if end if brunt = where(brunt.eq.0, brunt@_FillValue, brunt) end if b = ismissing(a) ;当数组a中含有缺测值时，相应位置就会返回 True c = where( ismissing(a), 0, a ) ;把数组a中的缺测值用0代换 c = where( a.eq.a@_FillValue, 0, a ) ;按理来说效果同上，但实践证明这样不行 if(all(ismissing(data))) then print("Your data is all missing. Cannot create plot.") end if if(any(ismissing(data))) then print("Your data contains some missing values. Beware.") end if N = num(.not.ismissing(data)) ;统计非缺测值数量，同理也可统计缺测值数量 ;num 的作用是 Counts the number of True values in the input.使用ncl函数时一定要看清楚该函数对缺测值的处理方法。
- 自动忽略缺测值的函数：
dim_avg_n, dim_sum_n、ncl的加减乘除 - 把缺测值当正常值处理的函数：
bw_bandpass_filter
- 自动忽略缺测值的函数：