| 查看: 363 | 回复: 0 | ||
[求助]
高手帮帮忙啊
|
|
这个程序应该是计算密度分布的 但是看不懂，请教高手 而且编译的时候还有警告！ c program calculate density profile c density profiles : (1) polymer; (2) nano; (3) solvent; (4) total parameter (NdpdMax=1400000) parameter (Ndenpro=500) parameter (nGn=15) implicit double precision (a-h,o-z) character*25 file_in,file_out real*8 rdpd(3,NdpdMax),vdpd(3,NdpdMax),rdpdcont(3,NdpdMax) real*8 denpro(4,Ndenpro) real*8 dz real*8 s_size(3),cell(3),sd,rho,D0 Integer*4 species(NdpdMax) Integer*4 n_dpd_polymer,n_dpd_ps,n_dpd_nano integer*4 mdenpro Real*8 Ene(-1:3,-1:3),elsc,d_els,rm,T_solv,r_eq open(2, file='density.in') read(2,*)file_in, file_out, dz close(2) open(40, file=file_in) read(40,*)itt,ndpd,npoly,lchain,lchainA,n_dpd_ps & ,n_dpd_nano read(40,*)T_solv read(40,*)ENE(-1,-1),ENE(1,1),ENE(-1,1),ENE(0,-1) & , ENE(0,1),ENE(0,0),ENE(2,2),ENE(1,2) read(40,*)ENE(-1,2), ENE(0,2),ENE(3,3),ENE(2,3) & , ENE(1,3),ENE(-1,3),ENE(0,3) read(40,*)s_size,D0 D0=D0+2.75 dv=s_size(1)*s_size(2)*dz n_dpd_polymer=npoly*lchain mdenpro=D0/dz+0.5 do jn=1,mdenpro denpro(1,jn)=0.0 denpro(2,jn)=0.0 denpro(3,jn)=0.0 denpro(4,jn)=0.0 enddo do i=1,n_dpd_ps read(40,*)rdpd(1,i),rdpd(2,i),rdpd(3,i) > ,vdpd(1,i),vdpd(2,i),vdpd(3,i) & ,rdpdcont(1,i),rdpdcont(2,i),rdpdcont(3,i) & ,species(i) m=rdpd(3,i)/dz+1 denpro(4,m)=denpro(4,m)+1 if(i.le.n_dpd_polymer) denpro(1,m)=denpro(1,m)+1 if(i.gt.n_dpd_polymer.and.i.le.(n_dpd_polymer+n_dpd_nano)) then denpro(2,m)=denpro(2,m)+1 endif if(i.gt.(n_dpd_polymer+n_dpd_nano).and.i.le.n_dpd_ps) then denpro(3,m)=denpro(3,m)+1 endif enddo open(13, file=file_out) do k=1,mdenpro write(13, '(f10.3,4f18.7)')(k-0.5)*dz, & denpro(1,k)/dv, & denpro(2,k)/dv, & denpro(3,k)/dv, & denpro(4,k)/dv enddo stop end |
回复此楼» 猜你喜欢
今天务委会开完了，明天出结果吗
已经有17人回复
-
明天应该可查了！？
已经有6人回复
-
有没有大神帮我看看基金代码
已经有10人回复
-
放榜前的不淡定
已经有19人回复
-
某些机构，以效率低为荣，以效率低作为存在感
已经有10人回复
-
2026国自然函评费到账
已经有23人回复
-
在坚冰还盖着北海的时候，我看到了怒放的梅花。
已经有8人回复
-
今日不放榜？网传国自然预计 8 月 27 日可查结果
已经有22人回复
-
什么时候开奖？
已经有13人回复
-
如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》
已经有8人回复
找到一些相关的精华帖子，希望有用哦~
急…高手帮帮忙，谢谢了
已经有5人回复
- 已有核磁图，但有一个峰不知是什么，求高手帮忙，谢谢！
已经有17人回复
- PCR电泳图分析（有杂带，有拖尾）求高手们帮忙分析一下
已经有4人回复
- 请高手帮忙看是什么意思？
已经有5人回复
- SDS-PAGE的问题请教，大家帮帮忙啊！求高手啊！
已经有14人回复
- COD达十多万的含油废水怎么测COD？？头疼啊，求各路高手帮帮忙。。
已经有18人回复
- 论文latex突然出了一个问题, 高手帮忙解决啊~
已经有23人回复
- 关于联合培养法国CEA和美国LLNL实验室问题的咨询，希望高手帮帮忙!!!
已经有3人回复
- 请高手帮帮忙...对乙酰氨基酚氧化问题
已经有10人回复
- 急！！！！叔丁基肼都是以盐酸盐形式存在吗？高手请帮忙
已经有8人回复
- 不好溶解的配体怎么培养晶体啊？各位高手帮帮忙~~~
已经有12人回复
- 柱子堵了怎么办？？高手们帮帮忙啊~~~
已经有7人回复
- 【求助】请虫友们帮帮解一下这个化合物啊，有氢谱，碳谱和质谱数据
已经有9人回复
- 【求助/交流】请求高手指点，帮忙解释下TMHMM预测结果
已经有14人回复
- 【求助】求助各位高手油脂碘值测定韦氏试剂怎么配啊？
已经有19人回复
- 【求助】急求高手帮忙！只有cif如何晶体画图（绝对没造假！）
已经有10人回复
- 【求助】请高手帮忙看看，派送大红包。dmol3 模块
已经有18人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
28