Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

regarding dimensions #1

Open
omkantnirala opened this issue May 20, 2019 · 0 comments
Open

regarding dimensions #1

omkantnirala opened this issue May 20, 2019 · 0 comments

Comments

@omkantnirala
Copy link

I have a fortran code for finding the number of hydrogen bond from md simulation trajectory. But still i did not understand about the dimension.

    "dimension it(6000,5000)
    dimension cx(6000,5000),cy(6000,5000),cz(6000,5000)
    dimension vx(6000,5000),vy(6000,5000),vz(6000,5000)
    dimension cmx(6000,5000),cmy(6000,5000),cmz(6000,5000)
    dimension ax(6000),ay(6000),az(6000),at(6000)
    dimension snx(6000,5000),sny(6000,5000),snz(6000,5000)
    dimension shx(6000,5000,5),shy(6000,5000,5),shz(6000,5000,5)

c
dimension histt(800),histth(800),sumh(800),ph(10),phnew(10)
dimension histt1(800),histth1(800),sumh1(800)
dimension h(5000,108,108)
dimension h1(5000,108,108)
dimension avau1(10000),avau2(10000),an(10000) "
c
real mas,masinv,mn,md
do i=1,800
histt(i)=0.
histth(i)=0.
sumh(i)=0.
enddo
c
do i=1,5000
do j=1,76
do k=1,76
h(i,j,k)=0.
h1(i,j,k)=0.
enddo
delt=1./5000
do i=1,ntotal
read(51,)at(i),ax(i),ay(i),az(i)
c write (1,
) ax(i), ay(i), az(i)
do j=1,natom
c read(51,)it(i,j),cx(i,j),cy(i,j),cz(i,j),vx(i,j),vy(i,j),vz(i,j)
read(51,
)it(i,j),cx(i,j),cy(i,j),cz(i,j)
enddo
enddo

c Calculate centre-of-mass coordinates
do 12 i=1,ntotal
ax(i)=ax(i)
ay(i)=ay(i)
az(i)=az(i)
Just i want to know that about ".....". Can anyone help me to define the dimension terms.

Thanks

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant