<meta http-equiv="Content-Type" content="text/html; charset=utf-8"><div dir="auto"><div>Hi Victorien,<div dir="auto"><br></div><div dir="auto">You just need to put the allocate and </div>getCellCoords call outside of the j loop. And in between the 'end do's, put a deallocate(yCoord).</div><div dir="auto"><br><br><div class="gmail_quote" dir="auto"><div dir="ltr" class="gmail_attr">On Tue, Sep 8, 2020, 6:45 AM Victorien Bouffetier <<a href="mailto:victorien.bouffetier@u-bordeaux.fr">victorien.bouffetier@u-bordeaux.fr</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div style="font-family:arial,helvetica,sans-serif;font-size:12pt;color:#000000"><div><div style="font-family:arial,helvetica,sans-serif;font-size:12pt;color:#000000"><div>Hi everyone!<br></div><br><div>I would need to have acces to the physical coordinates of the mesh in the file Simulation_adjustEvolution.F90, but I don't  know how to do so. I have been trying to get the physical coodinates using Grid_getCellCoords(...), but it keeps crashing.<br></div><br><div>here is the piece of codes which fails. It shows no problem at the compilation but  crashes while executing the code when entering the do procedure.<br></div><br><br><div>do lb = 1, blkcnt<br>        call Grid_getBlkIndexLimits(blklst(lb),blkLimits,blkLimitsGC)<br>        call Grid_getBlkPtr(blklst(lb), blkPtr)<br><br>           do j = blkLimits(LOW,JAXIS), blkLimits(HIGH,JAXIS)<br></div><br><div>               allocate(yCoord(blkLimitsGC(HIGH,JAXIS)))</div><div>               call Grid_getCellCoords(JAXIS, blockID,CENTER,gcell,yCoord, blkLimitsGC(HIGH,JAXIS))<br></div><br><div>                     if (yCoord(j)<=y_threshold) then<br></div><br><div>                              ...<br></div><br><div>                    end if<br></div><div>          end do <br></div><div>end do<br></div><br><br><br><div>When the code crashes it returns this error pointing at the line where I have the Grid_getCellCoords(...).<br></div><br><br><div>Program received signal SIGSEGV: Segmentation fault - invalid memory reference.<br></div><br><br><div>So my questions are: <br></div><div>_ Am I using the right subroutine to get the physical coordinates?<br></div><div>_ Do you have any idea of why it says me 'invalid memory reference' while my variables and subroutine are declared?<br></div><br><br><div>Thank you so much for your patience and attention,</div><br><div>Best regards,<br></div><br><br>----------------------------------------------------<br><div>Victorien Bouffetier<br></div><div>phd student<br></div><div>Centre Lasers Intenses et Applications<br></div><div>Bordeaux University</div></div><br></div></div></div></blockquote></div></div></div>