Script to test uniform grid memory access time

cartesian_grid.py

from fipy import *
from scipy import interpolate
import numpy as np

##Create mesh in the Cartesian Coordinate system (X,Y)
nx = 100
ny = nx
dx = 1.0
dy = dx
CL = nx*dx
Cmesh = Grid2D(dx=dx, dy=dy, nx=nx, ny=ny)

C_V_xy = CellVariable(name = "C_V in XY plane",mesh=Cmesh,value = 0.0)

cvxy_res = C_V_xy((20,20))