 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 1 | 1/1 | 返回列表 |
| 查看: 207 | 回復: 0 | |||
[交流]
【討論】麻煩各位高手給看一下我這程序有問題嗎？多謝
|
|
麻煩高手幫忙給看一下我這段程序有問題嗎？為什么開始幾十步的時候，能量基本保持不變，但是后來，動能急劇增加呢，勢能也在增加，有勞各位了，多謝了 SUBROUTINE L_J (n,a,ahalf,rcut,position_x,position_y,position_z,u,f_x,f_y,f_z,Ep) IMPLICIT NONE ! the number of particles INTEGER,INTENT(IN)::n ! the potential cutoff REAL,INTENT(IN)::rcut ! the length and half length of box REAL,INTENT(IN)::a,ahalf !potential energy REAL,INTENT(OUT)::Ep ! the position of particles REAL(8),DIMENSION(n),INTENT(INOUT)::position_x,position_y,position_z ! the potential and force of particles REAL,DIMENSION(n),INTENT(OUT)::u,f_x,f_y,f_z ! the contral variable of loop INTEGER::i,j,k ! the label of particles INTEGER::ijk ! the distance of every two particles along the three axis REAL::distance_x,distance_y,distance_z ! the distance's square of every two particles REAL::distance_sq REAL::distance6,distance12,rcut_sq,rcut6,rcut12 REAL::Ec,ff ! initialize the potential and force u=0. Ep=0 f_x=0. f_y=0. f_z=0. ijk=0 !calculate the rcut_sq rcut_sq=rcut*rcut rcut6=rcut_sq**3 rcut12=rcut6*rcut6 Ec=(1./rcut12-1./rcut6) position_i: DO i=1,n-1 distance_ij: DO j=i+1,n distance_x=position_x(i)-position_x(j) distance_y=position_y(i)-position_y(j) distance_z=position_z(i)-position_z(j) !x_distance: distance_x=distance_x-a*REAL(ANINT(distance_x/a)) !y_distance: distance_y=distance_y-a*REAL(ANINT(distance_y/a)) !z_distance: distance_z=distance_z-a*REAL(ANINT(distance_z/a)) ! calculate the distance_sq distance_sq=distance_x**2+distance_y**2+distance_z**2 ! calculate potential and force force:IF (distance_sq distance6=distance_sq**3 distance12=distance6*distance6 ff=(distance_x/distance_sq)*((12./distance12)-(6./distance6)) ! u(ijk)=u(ijk)+(1./distance12-1./distance6)-Ec f_x(i)=f_x(i)+4.*ff f_y(i)=f_y(i)+4.*ff f_z(i)=f_z(i)+4.*ff f_x(j)=f_x(j)-4.*ff f_y(j)=f_y(j)-4.*ff f_z(j)=f_z(j)-4.*ff Ep=Ep+(1./distance12-1./distance6)-Ec END IF force END DO distance_ij ! u(ijk)=4.*u(ijk) END DO position_i Ep=4.*Ep Ep=Ep/real(n) END SUBROUTINE L_J [ Last edited by ghcacj on 2010-12-16 at 13:08 ] |
回復此樓» 猜你喜歡
07化學280分求調(diào)劑
已經(jīng)有3人回復
-
考研化學308分求調(diào)劑
已經(jīng)有8人回復
-
070300化學求調(diào)劑
已經(jīng)有7人回復
-
生物學一志愿985，分數(shù)349求調(diào)劑
已經(jīng)有9人回復
-
341求調(diào)劑(一志愿湖南大學070300)
已經(jīng)有4人回復
-
336化工調(diào)劑
已經(jīng)有5人回復
-
279分求調(diào)劑 一志愿211
已經(jīng)有17人回復
-
291求調(diào)劑
已經(jīng)有4人回復
-
328求調(diào)劑，英語六級551，有科研經(jīng)歷
已經(jīng)有12人回復
-
一志愿上海交大生物與醫(yī)藥專碩324分，求調(diào)劑
已經(jīng)有5人回復
» 搶金幣啦！回帖就可以得到:
-
紹興大學博士教師招聘1位（化學化工相關）
+5/255
-
2026年國基求好運！
+5/190
-
西北大學化學與材料科學學院博士招生（還有兩個名額，3月30日截至）
+1/63
-
散金幣，求基金中標��！
+1/60
-
天津理工大學國家杰青王鐵團隊招聘
+1/39
-
2026設施園藝作物環(huán)境調(diào)控與綠色生產(chǎn)團隊招收相關專業(yè)調(diào)劑研究生
+1/36
-
哈爾濱工業(yè)大學（深圳）-何自開教授團隊誠招化學工程與技術專業(yè)博士研究生
+1/35
-
催化方向推薦一個極好的博導
+1/32
-
招收熒光探針/有機合成/生物學 博士后--薪金: 30萬-50萬
+1/32
-
從專科一路到211碩士，到TOP5博士，也想進高校體制
+1/30
-
哈工大鄭州高等研究院劉紹琴與查正寶教授課題組誠招合成生物學背景相關的優(yōu)秀人才加盟
+1/29
-
湖南大學微生物結(jié)構(gòu)與功能實驗室2026年計劃招收博士研究生
+1/17
-
北京石油化工學院智能防腐涂層課題組招收2026級調(diào)劑研究生
+1/16
-
招聘科研助理——儲能方向（儲能、流體傳熱、CFD模擬）
+1/15
-
寧波東方理工大學余鵬課題組2026年上海交大聯(lián)培博士生、博士后、科研助理招聘
+1/8
-
雙一流天津工業(yè)大學電信學院李鴻強教授招收2026年申請審核制博士
+1/7
-
大連海事大學船舶潔凈能源研究中心2026年博士研究生招生啟事
+1/2
-
“機械、材料”招生
+1/1
-
招收生物與醫(yī)藥（專碩）調(diào)劑生（廣西某高校、B區(qū)）
+1/1
-
武漢紡織大學紡織纖維及制品教育部重點實驗室招收博士研究生
+1/1
相關版塊跳轉(zhuǎn)
我要訂閱樓主 leigp 的主題更新
| 1 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 070300化學求調(diào)劑 +7 | 苑豆豆 2026-03-20 | 7/350 |
|
|---|---|---|---|---|
|
|
[考研] 生物學一志愿985，分數(shù)349求調(diào)劑 +6 | zxts12 2026-03-21 | 9/450 |
|
|
[考研] 316求調(diào)劑 +7 | 梁茜雯 2026-03-19 | 7/350 |
|
|
|
[考研]
 材料與化工考研調(diào)劑
+4
|
孅華 2026-03-22 | 4/200 |
|
|
|
[考研] 276求調(diào)劑 +3 | YNRYG 2026-03-21 | 4/200 |
|
|
|
[考研] 一志愿西安交通大學材料工程專業(yè) 282分求調(diào)劑 +11 | 楓橋ZL 2026-03-18 | 13/650 |
|
|
|
[考研] 287求調(diào)劑 +8 | 晨昏線與星海 2026-03-19 | 9/450 |
|
|
|
[考研] 260求調(diào)劑 +3 | 朱芷琳 2026-03-20 | 4/200 |
|
|
|
[考研] 材料學碩301分求調(diào)劑 +7 | Liyouyumairs 2026-03-21 | 7/350 |
|
|
|
[考研] 307求調(diào)劑 +3 | 余意卿 2026-03-18 | 3/150 |
|
|
|
[考研] 一志愿武理材料305分求調(diào)劑 +6 | 想上岸的鯉魚 2026-03-18 | 7/350 |
|
|
|
[考研] 一志愿西南交大，求調(diào)劑 +5 | 材化逐夢人 2026-03-18 | 5/250 |
|
|
|
[考研] 296求調(diào)劑 +6 | www_q 2026-03-18 | 10/500 |
|
|
|
[考研] 308求調(diào)劑 +3 | 阿姐阿姐家啊 2026-03-18 | 3/150 |
|
|
|
[考研] 一志愿中海洋材料工程專碩330分求調(diào)劑 +8 | 小材化本科 2026-03-18 | 8/400 |
|
|
|
[考研] 一志愿蘇州大學材料求調(diào)劑，總分315（英一） +5 | sbdksD 2026-03-19 | 5/250 |
|
|
|
[考研] A區(qū)線材料學調(diào)劑 +5 | 周周無極 2026-03-20 | 5/250 |
|
|
|
[考研] 086500 325 求調(diào)劑 +3 | 領帶小熊 2026-03-19 | 3/150 |
|
|
|
[考研] 085600材料與化工調(diào)劑 324分 +10 | llllkkkhh 2026-03-18 | 12/600 |
|
|
|
[考研] 085600材料與化工求調(diào)劑 +6 | 緒幸與子 2026-03-17 | 6/300 |
|