24小时热门版块排行榜

>论坛更新日志 (3661)
>导师招生 (543)
>考博 (431)
>虫友互识 (230)
>招聘信息布告栏 (199)
>休闲灌水 (186)
>物理 (140)
>文献求助 (132)
>硕博家园 (35)
>海外博后 (33)
>学术会议 (28)
>博后之家 (27)
>教师之家 (23)
>基金申请 (22)
>论文投稿 (20)
>留学生活 (12)

返回列表

本帖产生 1 个模拟EPI ，点击这里进行查看

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

flowerstudy

木虫 (正式写手)

应助: 4 (幼儿园)
金币: 4402.6
帖子: 854
在线: 137.4小时
虫号: 162807

[交流] 【求助】求高人给我发一份GEMC的模拟程序

小弟，最近在学习MC，觉得GEMC对我比较有用，可是在网上，我怎么都找不着。
如果有人有的话，能否给我发一份呢？万分感谢！
我的邮箱地址：zw198204@163.com

回复此楼

» 猜你喜欢

越来越担心硕士毕业后没有竞争力了已经有17人回复
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急已经有6人回复
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急已经有6人回复
【长期招聘】复旦大学 AI 制药实验室博士后招聘已经有11人回复
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急已经有7人回复
2027博士申请-临床检验诊断学、病原生物学、微生物方向已经有5人回复
各位基础老师们，看完中央《从严治党》的决定后，对学术界更加有信心了吧已经有4人回复
申博自荐 2篇一区TOP+1篇EI+2公开发明专利 2027博士申请已经有3人回复

» 本主题相关价值贴推荐，对您同样有帮助:

求助sybyl关于添加形式电荷的问题已经有8人回复
正在做大学生创新实践项目，车身方面，研究方向一直反复不定，跪求高人指点啊已经有4人回复
苯环上的烷基取代。求高人指点。谢谢！已经有23人回复
除掉反应生成的乙醇，拉动反应平衡。求高人指点已经有15人回复
求各位高人给小弟支个招：我是否需要跨专业考研？已经有24人回复
【求助】求将c语言改写为C++，求教高人啊！已经有6人回复
【求助】求高人给分析一下三方相的xrd，帮我算算晶胞参数啊,自己不在行这个啊，赞助50 已经有8人回复
【求助】100金币求助高人一个程序解方程已经有15人回复

» 抢金币啦！回帖就可以得到:

查看全部散金贴

1楼 2010-11-28 18:52:27

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

qphll

金虫 (正式写手)

模拟EPI: 10
应助: 18 (小学生)
金币: 2708.7
帖子: 554
在线: 227.1小时
虫号: 89654

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
zh1987hs(金币+10):谢谢 2010-12-04 23:13:56
ghcacj(金币+3, 模拟EPI+1):精彩答疑 2010-12-06 12:16:56

答应过的事情, 虽然觉得其实对楼主帮助也不大. 因为我不可能给你写一个完整的GibbsMC, 如果你自己不努力去理解, 那么什么都是白搭.

这是大体结构. 按照需要填充就是了. 应该此贴不会再更新了.

文件名: GibbsMC.f90

文件内容如下:

PROGRAM MC

  USE global
  IMPLICIT NONE

  ! Local
  INTEGER :: ibox, iblock, istep
  INTEGER :: ts,tf,rate
  DOUBLE PRECISION :: engtst, virtst, engerr

  ! Initialize the random number generator
  CALL seed_random(idum)

  ! Read the input file
  CALL read_input

  ! Get the start time
  CALL system_clock(count=ts, count_rate=rate)

  ! Perform minimization
  CALL relax

  ! Calculate energy and virial
  DO ibox = 1,n_box
CALL eng_total(ibox,energy(ibox),virial(ibox))
CALL write_config(ibox,0)
  ENDDO

  ! Loop over the blocks
  DO iblock = 1,n_blocks

! Reset basic statistics
CALL stats(1,0,0)

! Reset optional statistics
IF(opt .EQ. 1) CALL sample(1,0,0)

! Loop over the steps per block
DO istep = 1,n_steps

   ! Perform a trial move
   CALL trial_move(iblock,istep)

   ! Accumulate basic statistics
   CALL stats(2,iblock,istep)

   ! Accumulate optional statistics
   IF(opt .EQ. 1) call sample(2,iblock,istep)

! End loop over the steps per block
ENDDO

! Average basic statistics for iblock
CALL stats(3,iblock,0)

! Average optional statistics for iblock
IF(opt .EQ. 1) call sample(3,iblock,0)

! Update the max values (disp, volchng, etc)
CALL update_max

! Print status
IF(iblock .LE. n_equil) THEN
   IF(iblock .EQ. 1) WRITE(*,'(A)')'Equilibration stage... '
   IF(MOD(NINT(100.0d0*DBLE(iblock)),10*n_equil) .EQ. 0)&
      &WRITE(*,'(I0,A3)')NINT(100.0d0*DBLE(iblock)/DBLE(n_equil)),'%  '
   IF(iblock .EQ. n_equil) WRITE(*,*)''
ELSE
   IF(iblock .EQ. n_equil + 1)  WRITE(*,'(A)')'Production stage... '
   IF(MOD(NINT(100.0d0*DBLE(iblock-n_equil)),10*(n_blocks-n_equil)).EQ. 0)&
      &WRITE(*,'(I0,A3)')NINT(100.0d0*DBLE(iblock-n_equil)/DBLE(n_blocks-n_equil)),'%  '
   IF(iblock .EQ. n_blocks) WRITE(*,*)''
ENDIF

! Check the total and updated energy consistency for each box (every 10 blocks)
IF(MOD(iblock,1) .EQ. 0) THEN
   DO ibox = 1,n_box
      CALL eng_total(ibox,engtst,virtst)
      engerr = dABS((engtst - energy(ibox)))
      IF (engerr .GE. 1.0d-5) THEN
      WRITE(*,*) 'WARNING: LARGE ERROR BETWEEN THE TOTAL AND UPDATED ENERGY'
      WRITE(*,'(A7,I0,A6,I0,A10,F20.8,A10,F20.8)')&
      &'BLOCK: ',iblock,' BOX: ',ibox,' UPD ENG: ',energy(ibox),' TOT ENG: ',engtst
      ENDIF
   ENDDO
ENDIF

  ! End loop over the blocks
  ENDDO

  ! Get the finish time
  CALL system_clock(count=tf, count_rate=rate)

  ! Write the elapsed time to screen
  CALL time_elapsed(tf,ts,rate)

  ! Write basic statistics
  CALL stats(4,0,0)

  ! Write optional statistics
  IF(opt .EQ. 1) CALL sample(4,n_blocks+1,0)

  ! Write the final configurations
  DO ibox = 1,n_box
CALL write_config(ibox,n_blocks+1)
  ENDDO

END PROGRAM MC