Fortran Rank mismatch in argument (Fortran 90)

  • Thread starter Thread starter MelihAltunan
  • Start date Start date
  • Tags Tags
    Argument rank
AI Thread Summary
The discussion centers on resolving a rank mismatch error in a Fortran 90 subroutine related to fluid dynamics problems using LBM methods. The user encountered a warning indicating a rank mismatch for the variable 'rho', which was attributed to discrepancies in the dimensions of the arguments between the subroutine and its call. After suggestions to check the subroutine's signature and the calling parameters, the user realized the mismatch in the number of arguments and successfully corrected the issue. The conversation highlights the importance of ensuring consistency in argument dimensions and declarations in Fortran programming. Ultimately, the user resolved the problem by aligning the subroutine call with its definition.
MelihAltunan
Messages
8
Reaction score
0
Hello everyone, i am dealing with the code which can help me to solve fluid dynamics problems with using LBM methods. Anyways, since i am beginner on Fortran i couldn't solve the rank mismatch error, i think it is easy one but i just can't fix it, i am waiting for your help. Here is the problem part;

subroutine collesion(u,v,f,feq,rho,omega,w,cx,cy,n,m,tm,tminv,stmiv)
real f(0:8,0:n,0:m)
real feq(0:8,0:n,0:m),rho(0:n,0:m)
real w(0:8), cx(0:8),cy(0:8)
real u(0:n,0:m), v(0:n,0:m)
real tm(0:8,0:8),tminv(0:8,0:8),stmiv(0:8,0:8)
real fmom(0:8,0:n,0:m),fmeq(0:8,0:n,0:m)
!calculate equilibrium moments
do i=0,n
do j=0,m
fmeq(0,i,j)=rho(i,j)
fmeq(1,i,j)=rho(i,j)*(-2.0+3.0*rho(i,j)*(u(i,j)*u(i,j)+v(i,j)*v(i,j)))
fmeq(2,i,j)=rho(i,j)*(1.0-3.0*rho(i,j)*(u(i,j)*u(i,j)+v(i,j)*v(i,j)))
fmeq(3,i,j)=rho(i,j)*u(i,j)
fmeq(4,i,j)=-rho(i,j)*u(i,j)
fmeq(5,i,j)=rho(i,j)*v(i,j)
fmeq(6,i,j)=-rho(i,j)*v(i,j)
fmeq(7,i,j)=rho(i,j)*(u(i,j)*u(i,j)-v(i,j)*v(i,j))
fmeq(8,i,j)=rho(i,j)*u(i,j)*v(i,j)
end do
end do

the error is Warning: Rank mismatch in argument 'rho' at (1) (rank-2 and scalar)

i guess it's because fmeq and rho have different ranks but how i can write in a proper way?

Thank you for now.
 
Technology news on Phys.org
I would check elsewhere in the program where the subroutine collesion is called. It is possible that the variable 'rho' has not been declared with the proper dimensions in other parts of the program. The error message suggests that there is a mismatch in the dimensions of the arguments between calling the routine collesion and the dimensions of rho in the subroutine itself.
 
i checked the code rho has same dimensions everywhere and also i am receiving 'segmentation fault-invalid memory reference' error on simply fortran then it shows the line 11 above (fmeq(0,i,j)=rho(i,j)). And when i try to run on microsoft developer studio it says 'wrong number of arguments to procedure COLLESION invoked from main: 13 found, 14 expected' without showing any lines. Do you have any idea about that? thank you for your help.
 
Asking for help like this without posting the entire program can lead to a lot of speculation and a waste of time.

Either post the entire program or a minimal one that still shows the problem...then, again, when one does the step of going to a minimal program often the problem is found.
 
By the way, I noticed that omega, n, and m do not seem to have been declared in your subroutine...this leads me to believe that you do NOT have an IMPLICIT NONE statement in your program...which opens the door for a lot of things.
 
i tried with IMPLICIT NONE but it had same problem again. Considering what you said, i added the code as an attachment. So, the best way is if you run the code with your own compiler, you will see the what is an error. I would be pleased if you can check the code that i uploaded. Thanks.
 

Attachments

No need to compile...the problem is very clear, just compare the signature of the subroutine as declared (order and quantity of arguments) and see how you are calling it:

Code:
call collesion(u,v,feq,rho,omega,w,cx,cy,n,m,tm,tminv,stmiv)
subroutine collesion(u,v,f,feq,rho,omega,w,cx,cy,n,m,tm,tminv,stmiv)

Notice anything?
 
yea, i noticed :), it is working now, thanks a lot.
 

Similar threads

Replies
8
Views
6K
Replies
17
Views
3K
Replies
2
Views
2K
Replies
4
Views
2K
Replies
1
Views
7K
Replies
14
Views
2K
Replies
2
Views
2K
Replies
20
Views
2K
Back
Top