<div dir="ltr"><div>Dear flash users,</div><div><br></div><div>I have a question regarding the data during a hydrodynamic time step.</div><div><br></div><div>I need to modify the update of the data (density , pressure , velocity ...) in some specific <br></div><div>cells inside the simulation domain.</div><div><br></div><div>For the specific cells I need to do modify the update of the solution at the end of the hydro time step Dt^n = t^n+1 - t^n <br></div><div><br></div><div>When the time step is over I will do <br></div><div><br></div><div>EINT^n+1 = EINT^n + cte*Dt^n     (1)<br></div><div><br></div><div>where EINT^n is the internal energy at time t^n  <br></div><div><br></div><div>so for those specific cells the update of the solution is not given by the hydro step <br></div><div>but by this equation (1). However I need the result of this hydro step in that cell to compute the constant cte.<br></div><div><br></div><div>The question is how do I call the data in my new physic module before and after the update.</div><div><br></div><div>For example I think that I could get the data after the time step by doing something like<br></div><div><br></div><div>use Grid_interface, ONLY : Grid_getBlkPtr, Grid_getBlkIndexLimits <br></div><div><br></div><div>implicit none<br>  <br>  integer,intent(IN)  :: blockCount<br>  integer,dimension(blockCount), intent(IN) :: blockList<br>  real,intent(IN) :: dt, time<br><br>  real,pointer,dimension(:,:,:,:) :: solnData<br>  integer, dimension(2,MDIM) :: blkLimits, blkLimitsGC</div><div><br></div><div>do lb = 1, blockCount</div><div><br>     blockID=blockList(lb)<br>     call Grid_getBlkPtr(blockID,solnData)<br>     call Grid_getBlkIndexLimits(blockID,blkLimits,blkLimitsGC)<br>     do k = blkLimits(LOW,KAXIS), blkLimits(HIGH,KAXIS)<br>        do j = blkLimits(LOW,JAXIS), blkLimits(HIGH,JAXIS)<br>           do i = blkLimits(LOW,IAXIS), blkLimits(HIGH,IAXIS)</div><div><br></div><div>              EINT^n+1 =  solnData(EINT_VAR,i,j,k)</div><div><br></div><div>....<br></div><div><br></div><div>but what about EINT^n (the internal energy before the update) ???<br></div><div><br></div><div>Thank you very much for any help on this !</div><div><br></div><div>Best regards,</div><div><br></div><div>Antoine.<br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div></div>