Debugging negative pressure values using traceback in Intel oneAPI on Windows

[Niyas]

Hello everyone, 

How can I use the traceback function in my 3D numerical simulations in Fortran, while using a Windows computer and Intel oneAPI, to identify the cause of negative pressure values and what is the command for it?

[mecej4]

At the various locations in your code where new values of pressure are calculated, you can add a couple of diagnostic statements along the lines of 

pnew = ....
if (pnew < 0.0) then
   print *, i,j,k, t, vx, vy, vz, ..., ' p new is negative'
   call tracebackqq(...)
endif
p(i,j,k, m) = pnew
...

[Niyas]

@mecej4  Thanks for your reply. 

I am just getting started with traceback and debugging.

Can you please give me an idea about the 'tracebackqq' subroutine? What should it do ? 

[mecej4]

Tracebackqq() is provided by the compiler, and is documented in the compiler development reference guide  .