Fortran arrays and subroutines (sub arrays)

Posted by ccook on Stack Overflow See other posts from Stack Overflow or by ccook
Published on 2010-03-22T18:37:59Z Indexed on 2010/03/23 2:21 UTC
Read the original article Hit count: 912

Filed under:

fortran

|

subroutine

|

dimensioning

I'm going through a Fortran code, and one bit has me a little puzzled.

There is a subroutine, say

SUBROUTINE SSUB(X,...)
REAL*8 X(0:N1,1:N2,0:N3-1),...
...
RETURN 
END

Which is called in another subroutine by:

CALL SSUB(W(0,1,0,1),...)

where W is a 'working array'. It appears that a specific value from W is passed to the X, however, X is dimensioned as an array. What's going on?

© Stack Overflow or respective owner

Related posts about fortran

passing a fortran int array to C++ by calling C++ function in fortran

as seen on Stack Overflow - Search for 'Stack Overflow'
hi, I am trying to call a C++ function in FORTRAN subroutine. This C+ function is supposed to update an integer array. Here is a non-working code I wrote. Can someone please let me know what the issue is: ! FORTRAN function that calls a C++ function subroutine my_function() integer(4) ar(*) … >>> More
MinGW/G++/g95 link problem - libf95 undefined reference to `MAIN_'

as seen on Stack Overflow - Search for 'Stack Overflow'
Hi folks, Summing up, my problem consists on compiling g95 objects inside a C++ application. Actually, I'm constructing an interface for an old fortran program. For this task, I'm using the wxWidgets GUI library, and calling fortran subroutines when necessary. At the beginning, I was developing the… >>> More
Fortran arrays and subroutines

as seen on Stack Overflow - Search for 'Stack Overflow'
I'm going through a Fortran code, and one bit has me a little puzzled. There is a subroutine, say SUBROUTINE SSUB(X,...) REAL*8 X(0:N1,1:N2,0:N3-1),... ... RETURN END Which is called in another subroutine by: CALL SSUB(W(0,1,0,1),...) where W is a 'working array'. It appears that a specific… >>> More
fortran 90 user defined type, passing by value ?

as seen on Stack Overflow - Search for 'Stack Overflow'
I have an issue in fortran 90. I have a user defined type and when I call one of the MPI subroutines the data looks to be passed by values (not address, as I thought it should). The output arguments aren't modified. It seems to be specific to the MPI calls I tried the same thing in a simple test… >>> More
Should External Routine Be Declared Always in Fortran?

as seen on Stack Overflow - Search for 'Stack Overflow'
In my fortran code I made the following call to dnrm2 routine d = dnrm2(n, ax, 1) Just a simple call that would return me a double precision result. The problem is whether I should declare the function at the start of my script. I found that if I don't declare it, and when I compile the code in… >>> More

Related posts about subroutine

Return data from subroutine while the subroutine is still processing

as seen on Stack Overflow - Search for 'Stack Overflow'
Is there any way to have a subroutine send data back while still processing? For instance (this example used simply to illustrate) - a subroutine reads a file. While it is reading through the file, if some condition is met, then "return" that line and keep processing. I know there are those that will… >>> More
Fortran 90 How to call a function in a subroutine in a module?

as seen on Stack Overflow - Search for 'Stack Overflow'
I have a module that "contains" a subroutine that "contains" a function. I say "use themodule" in my main program and I can "call thesubroutine", but how to I get the function that is "contained" in the subroutine? The code looks like this: module useful integer, parameter :: N=2 double precision… >>> More
vim opens my folds on paste. how to prevent that ?

as seen on Stack Overflow - Search for 'Stack Overflow'
There's something weird going on in my fortran folding. This is the example file module foo contains subroutine hello() end subroutine hello subroutine hello() end subroutine subroutine hello() end subroutine end module foo subroutine hello() end subroutine subroutine hello() end subroutine subroutine… >>> More
GLSL subroutine not being used

as seen on Game Development - Search for 'Game Development'
I'm using a gaussian blur fragment shader. In it, I thought it would be concise to include 2 subroutines: one for selecting the horizontal texture coordinate offsets, and another for the vertical texture coordinate offsets. This way, I just have one gaussian blur shader to manage. Here is the code… >>> More
Perl :Getting the name of the current subroutine

as seen on Stack Overflow - Search for 'Stack Overflow'
Hi all, In perl we can get the name of the current package and current line number Using the predefined variables like _PACKAGE and __LINE . Like this I want to get the name of the current subroutine. Example use strict; use warnings; print __PACKAGE__; sub test() { print __LINE__; } &test(); In… >>> More