ERROR CHECKING TO MAKE SURE DIFFRENCING IS DONE OVER A CONSISTENTLY-SIZED GRID just bail on it and let the user figure it out NOW EXECUTE THE FINITE DIFFERENCES - NOTE THAT LOOP INDICES ARE MEANT TO INDEX ARRAY BEING DIFFERENCED AND NOT THE MESH STRUCTURE, WHICH USES INPUT BOUNDS. TO KEEP THE CODE CLEAN I'VE ALIASED THE GRID VARS SO THAT THEY MAY BE ACCESSED BY LOOP INDEX. fwd diff. at beginning, note that h1 is cell-centered centered diff. in the middleq backward diff. at end
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
real(kind=wp), | intent(in), | dimension(:,:,:) | :: | f | ||
type(curvmesh), | intent(in) | :: | x | |||
integer, | intent(in) | :: | lbnd1 | |||
integer, | intent(in) | :: | ubnd1 | |||
integer, | intent(in) | :: | lbnd2 | |||
integer, | intent(in) | :: | ubnd2 | |||
integer, | intent(in) | :: | lbnd3 | |||
integer, | intent(in) | :: | ubnd3 |