COMPUTE INTERFACE VELOCITIES AND LOAD UP GHOST CELLS FOR FLUID STATE VARIABLES
Note that it is done on a per species basis 5/23/2015 - may need to be changed for 2D/1D sims which have only one element in the x2 direction...
v1i(lx1+1,:,:)=v1i(lx1,:,:) !avoids issues with top boundary velocity spikes which may arise vadvx1(1,:,6)=mean(vadvx1(1,:,6)); !for eq sims
FOR INTERNAL ENERGY NOW DEAL WITH ADVECTION ALONG X3; FIRST IDENTIFY MY NEIGHBORS
PASS X3 BOUNDARY CONDITIONS WITH GENERIC HALOING ROUTINES we only need one ghost cell to compute interface velocities AFTER HALOING CAN COMPUTE THE X3 INTERFACE VELOCITIES NORMALLY
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| integer, | intent(in) | :: | isp | |||
| logical, | intent(in) | :: | isperiodic | |||
| real(kind=wp), | intent(inout), | dimension(-1:,-1:,-1:,:) | :: | ns | ||
| real(kind=wp), | intent(inout), | dimension(-1:,-1:,-1:,:) | :: | rhovs1 | ||
| real(kind=wp), | intent(inout), | dimension(-1:,-1:,-1:,:) | :: | vs1 | ||
| real(kind=wp), | intent(inout), | dimension(-1:,-1:,-1:,:) | :: | vs2 | ||
| real(kind=wp), | intent(inout), | dimension(-1:,-1:,-1:,:) | :: | vs3 | ||
| real(kind=wp), | intent(inout), | dimension(-1:,-1:,-1:,:) | :: | rhoes | ||
| real(kind=wp), | intent(out), | dimension(1:size(vs1,1)-3,1:size(vs1,2)-4,1:size(vs1,3)-4) | :: | v1i | ||
| real(kind=wp), | intent(out), | dimension(1:size(vs1,1)-4,1:size(vs1,2)-3,1:size(vs1,3)-4) | :: | v2i | ||
| real(kind=wp), | intent(out), | dimension(1:size(vs1,1)-4,1:size(vs1,2)-4,1:size(vs1,3)-3) | :: | v3i |