<code>
#数据处理的.f也不放上啦～
# caculate:Carbon_nemd_directly_200K
# 用非平衡态直接法计算200K情况下的(10，0)碳纳米管的热导率
#------------------------------------------------------------------------------------------------------------------------------------------
# 先指定单位制（units）
units metal
#--------------------
# 再设定原子类型（atom_style）
atom_style atomic
#----------------------------
# 一般在计算前将const常量先声明（赋值）了（variable 变量名 equal xxx）equal用于设定const
variable T equal 200
#--------------------------------------------------------------------------------------
# boundaray(边界条件)
# p：周期性边界条件periodic
# f：非周期性固定边界条件fixed
# s：非周期性包覆边界条件shrink-wrapped
# m：非周期性包覆最小值边界条件minimum value
# 常用p和f
boundary p p p
#------------------------------------------
# 晶体建模窝个人选择外部建模『从外部读取原子相对位置数据data建模』
read_data ./data.tube10
#---------------------------------------------------------------
# 然而感觉并不需要restart
</code>