Hi,

One obvious problem is that you have defined an allocatable which is not allocated before being referenced

goodluck

sina

Thanks for your answer. Your mean this, if I clean the " allocatable" ,has the problem solved? 

Checked it..........

You need to either define real f(60,60) or real,allocatable f(:,:); allocate(f(60,60))

Thanks Bharat Kumar. 

I test it now but it did not work. 

Thanks Sina Haeri.

I test it now and it worked but the matrix does not have 3600 array. What I should do?

What do you mean it does not have 3600 array? You mean elements? If you have allocated it as f(60,60) then it "must" have 3600 elements

But it does not have 3600 elements. You can check it.

Hi,

Not sure what you mean, attached is a slightly modified version of your code and it's output as you can see (first line of the output) it "exactly" has 3600 elements, which would be utterly superrising if it didn't.

Sina

Thank you very much. Good luck.

Sina,  please accept my admiration: looking at the code of smb.  else is harder than to

write it from ground 0. 

Have you tried increasing the declared integer size - integer(kind=8), or kind=16, or kind=32 or kind=64?  You might have to use "selected_int_kind(r)" to work out the size of the integer or where it is going out of bounds. 

You did not ALLOCATE an array that is defined as ALLOCATABLE. You have to introduce a line

ALLOCATE(F(60,60))

before assign any value to the array. Now the program is running.

Hi

You must determine the size of ALLOCATABLE array by writing ALOOCATE (F(60,60)) statement.