Help with contourf of function: 2D Plot with 3 variables

Say I have a porous geometry of connected cylindrical pores. The configuration consist of 2 layers and 3 columns. My code is a function that computes the following:

• resistance of the pore flow R(ii,jj,1,j)
• the fluid velocity which is a function of the resistance itself.

  up(ii,jj,k,j)= (p(ii,jj,j)-p(ii+1,jj,j))*(1/(8*R(ii,jj,1,j)*a(ii,jj,k,j)^2));

• ii is the index of the corresponding layer
• jj the index of the pore.
• k is the iteration over the x axis
• j is the iteration over time

I would like to plot a contourf of the attached image form. Only difference is I want the x axis (instead of apex angle) to be the Resistance of the first layer R(1,1,1,j); the y axis (instead of orientation angle) the Resistance of the second layer R(2,2,1,j). and the z axis to be the corresponding fluid velocity up(ii,jj,k,j) (instead of Potential P)

Any help or leads would be appreciated!