Question

I’m trying to check how many times certain lines are executed in few timesteps, here is part of my code:
1)write “countd” every time:

real(kind=8) function dergfm(jp,ip,lp)
integer :: jp,ip,lp,countd
real(kind=8) press
    .
    .
    .
   countd=countd+1
   !if (countd < 5) then 
   print*, "countd= ", countd
   !endif
    .
    .
    .
end function dergfm

result:

countd=            1
.
.
.
countd=            21504

2)write “countd” first 4 times:

    real(kind=8) function dergfm(jp,ip,lp)
integer :: jp,ip,lp,countd
real(kind=8) press
    .
    .
    .
   countd=countd+1
   if (countd < 5) then 
   print*, "countd= ", countd
   endif
    .
    .
    .
end function dergfm

results:

countd=            1
countd=            2
countd=            3
countd=            4

3)write “countd” larger than 5

    real(kind=8) function dergfm(jp,ip,lp)
integer :: jp,ip,lp,countd
real(kind=8) press
    .
    .
    .
   countd=countd+1
   if (countd > 5) then 
   print*, "countd= ", countd
   endif
    .
    .
    .
end function dergfm

results:

[none]

It looks like comparing these integers works for .lt. but doesn’t for .gt. (nor .eq.)

Answer 1

countd is a local variable to a function. As such, its value it not necessarily retained between calls, unless you add the “save” attribute to the declaration. Try this and see if the code starts behaving. And how do you initialize countd? I would use the following:

integer, save :: countd = 0

In this case the “save” is optional because it is implied by the initialization in the declaration.